Changeset 16972 for issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp

Timestamp:

11/29/13 10:08:46 (11 years ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing hydro tests

File:

• issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp (modified) (2 diffs)

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

@@ -187 +187 @@
                 D[0][0]=D_scalar;
                 D[1][1]=D_scalar;
-                GetB(B,element,xyz_list,gauss); 
+                GetB(B,basalelement,xyz_list,gauss); 
                 TripleMultiply(B,2,numnodes,1,
                                         &D[0][0],2,2,0,
@@ -249 +249 @@
         IssmDouble sediment_storing = SedimentStoring(basalelement);
         basalelement->FindParam(&dt,TimesteppingTimeStepEnum);
-        Input* water_input       = element->GetInput(BasalforcingsMeltingRateEnum); _assert_(water_input);
-        Input* transfer_input    = element->GetInput(WaterTransferEnum);            _assert_(transfer_input);
-        if(dt!= 0.){old_wh_input = element->GetInput(SedimentHeadOldEnum);          _assert_(old_wh_input);}
+        Input* water_input       = basalelement->GetInput(BasalforcingsMeltingRateEnum); _assert_(water_input);
+        Input* transfer_input    = basalelement->GetInput(WaterTransferEnum);            _assert_(transfer_input);
+        if(dt!= 0.){old_wh_input = basalelement->GetInput(SedimentHeadOldEnum);          _assert_(old_wh_input);}
 
         /* Start  looping on the number of gaussian points: */