Changeset 21405

Timestamp:

11/21/16 12:00:04 (8 years ago)

Author:

jbondzio

Message:

BUG: ensure that temperate base gets proper enthalpy value. Rounding errors would cause test (T>=Tpmp) to fail.

File:

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

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

@@ -129 +129 @@
                         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-1e-12);
                                 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);