Changeset 22845

Timestamp:

06/15/18 08:52:56 (7 years ago)

Author:

bdef

Message:

NEW:Testing transmitivity limiting based on frozen ground

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 base_elev,prestep_head,water_sheet;
         IssmDouble sediment_thickness       = element->GetMaterialParameter(HydrologydcSedimentThicknessEnum);
-
+        bool isthermal;
+
+        element->FindParam(&isthermal,TransientIsthermalEnum);
         element->FindParam(&unconf_scheme,HydrologydcUnconfinedFlagEnum);
         SedTrans_input->GetInputValue(&FullLayer_transmitivity,gauss);
+
         switch(unconf_scheme){
         case 0:
@@ -621 +625 @@
                 _error_("UnconfinedFlag is 0 or 1");
         }
+
+        /*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.;
+        }
+
         return sediment_transmitivity;
 }/*}}}*/