Changeset 22968

Timestamp:

07/17/18 13:22:53 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: fixing AD

File:

• issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/SealevelriseAnalysis.cpp

@@ -52 +52 @@
         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);
 
 }/*}}}*/