Changeset 18623

Timestamp:

10/14/14 02:51:17 (10 years ago)

Author:

jbondzio

Message:

BUG: H=c*(T-Tref)+… for nodes at pressure melting point

File:

• issm/trunk-jpl/src/m/classes/initialization.m (modified) (1 diff)

issm/trunk-jpl/src/m/classes/initialization.m

@@ -130 +130 @@
                         if md.thermal.isenthalpy,
                                 tpmp = md.materials.meltingpoint - md.materials.beta*md.initialization.pressure;
-                                pos  = find(md.initialization.temperature>tpmp);
+                                pos  = find(md.initialization.temperature>=tpmp);
                                 enthalpy      = md.materials.heatcapacity*(md.initialization.temperature-md.constants.referencetemperature);
-                                enthalpy(pos) = md.materials.heatcapacity*tpmp(pos) - md.constants.referencetemperature + md.materials.latentheat*md.initialization.waterfraction(pos);
+                                enthalpy(pos) = md.materials.heatcapacity*(tpmp(pos) - md.constants.referencetemperature) + md.materials.latentheat*md.initialization.waterfraction(pos);
                                 WriteData(fid,'data',enthalpy,'format','DoubleMat','mattype',1,'enum',EnthalpyEnum());
                         end