Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 22859)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 22860)
@@ -597,5 +597,6 @@
 	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,6 +629,7 @@
 	/*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.;
 	}