Index: /issm/trunk-jpl/src/c/objects/Materials/Matpar.cpp
===================================================================
--- /issm/trunk-jpl/src/c/objects/Materials/Matpar.cpp	(revision 11289)
+++ /issm/trunk-jpl/src/c/objects/Materials/Matpar.cpp	(revision 11290)
@@ -416,5 +416,10 @@
 /*FUNCTION Matpar::GetEnthalpyDiffusionParameter{{{1*/
 double Matpar::GetEnthalpyDiffusionParameter(double enthalpy,double pressure){
-	return thermalconductivity/(rho_ice*heatcapacity);
+	if(enthalpy<PureIceEnthalpy(pressure)){
+		return thermalconductivity/(rho_ice*heatcapacity);
+	}
+	else{
+		return 0.1*thermalconductivity/(rho_ice*heatcapacity);
+	}
 }
 /*}}}1*/