Changeset 12393

Timestamp:

06/08/12 13:44:01 (13 years ago)

Author:

seroussi

Message:

better check for spctemperature (enthalpy)

File:

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

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

@@ -48 +48 @@
                         checkfield(md,'thermal.spctemperature','forcing',1);
                         if (ismember(EnthalpyAnalysisEnum,analyses) & md.thermal.isenthalpy),
-                        checkfield(md,'thermal.spctemperature','<',md.materials.meltingpoint-md.materials.beta*md.materials.rho_ice*md.constants.g*md.geometry.thickness,'message','spctemperature should be below the adjusted melting point');
+                        checkfield(md,'thermal.spctemperature','<',md.materials.meltingpoint-md.materials.beta*md.materials.rho_ice*md.constants.g*(md.geometry.surface-md.mesh.z),'message','spctemperature should be below the adjusted melting point');
                         checkfield(md,'thermal.isenthalpy','numel',1,'values',[0 1]);
                         end