Index: /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 22967)
+++ /issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp	(revision 22968)
@@ -52,22 +52,13 @@
 	iomodel->FetchDataToInput(elements,"md.slr.Ngia",SealevelNGiaRateEnum);
 	iomodel->FetchDataToInput(elements,"md.slr.Ugia",SealevelUGiaRateEnum);
-
-	/*Initialize cumdeltalthickness and sealevel rise rate input: unfortunately, we don't have femmodel, so we 
-	 * need to iterate on elements, otherwise we would have used InputUpdateFromConstantx: */
-	counter=0;
-	for(int i=0;i<iomodel->numberofelements;i++){
-		if(iomodel->my_elements[i]){
-			Element* element=(Element*)elements->GetObjectByOffset(counter);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelriseCumDeltathicknessEnum);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelNEsaRateEnum);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelUEsaRateEnum);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelRSLRateEnum);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelEustaticMaskEnum);
-			element->InputUpdateFromConstant((IssmDouble)0.,SealevelEustaticOceanMaskEnum);
-			counter++;
-		}
-	}
-	
 	iomodel->FetchDataToInput(elements,"md.slr.steric_rate",SealevelriseStericRateEnum);
+
+	/*Initialize cumdeltalthickness and sealevel rise rate input*/
+	InputUpdateFromConstantx(elements,0.,SealevelriseCumDeltathicknessEnum);
+	InputUpdateFromConstantx(elements,0.,SealevelNEsaRateEnum);
+	InputUpdateFromConstantx(elements,0.,SealevelUEsaRateEnum);
+	InputUpdateFromConstantx(elements,0.,SealevelRSLRateEnum);
+	InputUpdateFromConstantx(elements,0.,SealevelEustaticMaskEnum);
+	InputUpdateFromConstantx(elements,0.,SealevelEustaticOceanMaskEnum);
 
 }/*}}}*/