Index: /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 15862)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matpar.cpp	(revision 15863)
@@ -398,8 +398,8 @@
 IssmDouble Matpar::GetEnthalpyDiffusionParameter(IssmDouble enthalpy,IssmDouble pressure){
 	if(enthalpy<PureIceEnthalpy(pressure)){
-		return thermalconductivity/(rho_ice*heatcapacity);
+	  return thermalconductivity/(rho_ice*heatcapacity);
 	}
 	else{
-		return 0.1*thermalconductivity/(rho_ice*heatcapacity);
+	  return 1.045*pow(10,-4.)/rho_ice; // K0=1.045*1e-4 from Aschwanden 2012. TODO: fetch K0 from model
 	}
 }