Changeset 15863

Timestamp:

08/22/13 01:15:22 (12 years ago)

Author:

jbondzio

Message:

CHG changed value of diffusive constant for enthalpy to value from Aschwanden 2012

File:

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

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

@@ -398 +398 @@
 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
         }
 }