Changeset 11290

Timestamp:

02/01/12 13:41:36 (13 years ago)

Author:

seroussi

Message:

different enthalpy diffusion parameter depending on enthalpy value

File:

• issm/trunk-jpl/src/c/objects/Materials/Matpar.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/objects/Materials/Matpar.cpp

@@ -416 +416 @@
 /*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*/