Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 22844)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 22845)
@@ -597,9 +597,13 @@
 	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,4 +625,11 @@
 		_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;
 }/*}}}*/