Changeset 18754

Timestamp:

11/07/14 10:14:00 (10 years ago)

Author:

seroussi

Message:

CHG: minor

File:

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

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

@@ -95 +95 @@
                                 pos=find(md.thermal.spctemperature(1:md.mesh.numberofvertices,:)~=NaN);
                                 replicate=repmat(md.geometry.surface-md.mesh.z,1,size(md.thermal.spctemperature,2));
-                                md = checkfield(md,'fieldname','thermal.spctemperature(find(md.thermal.spctemperature(1:md.mesh.numberofvertices,:)~=NaN))','<=',md.materials.meltingpoint-md.materials.beta*md.materials.rho_ice*md.constants.g*replicate(pos),'message','spctemperature should be less or equal than the adjusted melting point');
+                                md = checkfield(md,'fieldname','thermal.spctemperature(find(md.thermal.spctemperature(1:md.mesh.numberofvertices,:)~=NaN))','<=',md.materials.meltingpoint-md.materials.beta*md.materials.rho_ice*md.constants.g*replicate(pos)+10^-5,'message','spctemperature should be less or equal than the adjusted melting point');
                                 md = checkfield(md,'fieldname','thermal.isenthalpy','numel',[1],'values',[0 1]);
                                 md = checkfield(md,'fieldname','thermal.isdynamicbasalspc','numel', [1],'values',[0 1]);