Index: /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17082)
+++ /issm/trunk-jpl/src/c/analyses/FreeSurfaceTopAnalysis.cpp	(revision 17083)
@@ -16,5 +16,9 @@
 	if (iomodel->meshtype==Mesh2DhorizontalEnum)return;
 
+	int smb_model;
 	int finiteelement = P1Enum;
+
+	/*Fetch data needed: */
+	iomodel->Constant(&smb_model,SurfaceforcingsEnum);
 
 	/*Update elements: */
@@ -30,5 +34,4 @@
 	iomodel->FetchDataToInput(elements,SurfaceEnum);
 	iomodel->FetchDataToInput(elements,MaskIceLevelsetEnum);
-	iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum,0.);
 	iomodel->FetchDataToInput(elements,VxEnum);
 	iomodel->FetchDataToInput(elements,MeshVertexonsurfaceEnum);
@@ -37,4 +40,12 @@
 		iomodel->FetchDataToInput(elements,MeshElementonbedEnum);
 		iomodel->FetchDataToInput(elements,MeshElementonsurfaceEnum);
+	}
+	switch(smb_model){
+		case SMBEnum:
+			iomodel->FetchDataToInput(elements,SurfaceforcingsMassBalanceEnum,0.);
+			break;
+		default:
+			/*Nothing for now*/
+			;
 	}
 }/*}}}*/