Changeset 22860

Timestamp:

06/21/18 06:51:41 (7 years ago)

Author:

sjohnsen

Message:

CHG: Testing Transmisivity based on basal meltingrates

File:

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

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

@@ -597 +597 @@
         IssmDouble sediment_transmitivity;
         IssmDouble FullLayer_transmitivity;
-        IssmDouble T;
+        IssmDouble meltingrate;
+        IssmDouble groundedice;
         IssmDouble base_elev,prestep_head,water_sheet;
         IssmDouble sediment_thickness       = element->GetMaterialParameter(HydrologydcSedimentThicknessEnum);
@@ -628 +629 @@
         /*Let's deal with the frozen parts for which transmitivity is zero*/
         if (isthermal){
-                element->GetInputValue(&T,gauss,TemperatureEnum);
-                if (T<273.1) sediment_transmitivity=0.;
+                element->GetInputValue(&meltingrate,gauss,BasalforcingsGroundediceMeltingRateEnum);
+                element->GetInputValue(&groundedice,gauss,MaskGroundediceLevelsetEnum); 
+                if ((meltingrate<=0.0) && (groundedice>0)) sediment_transmitivity=0.;
         }