source: issm/oecreview/Archive/11250-11268/ISSM-11264-11265.diff

proj/ice/larour/issm-uci-clean/trunk-jpl/src/m/classes/thermal.m

@@ -45 +45 @@
                 function checkconsistency(obj,md,solution,analyses) % {{{
 
                         %Early return
-                        if ~ismember(ThermalAnalysisEnum,analyses) | (solution==TransientSolutionEnum & md.transient.isthermal==0), return; end
+                        if (~ismember(ThermalAnalysisEnum,analyses) & ~ismember(EnthalpyAnalysisEnum,analyses)) | (solution==TransientSolutionEnum & md.transient.isthermal==0), return; end
 
                         checkfield(md,'thermal.stabilization','numel',1,'values',[0 1 2]);
                         checkfield(md,'thermal.spctemperature','forcing',1);
+                        if ismember(EnthalpyAnalysisEnum,analyses),
+                        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');
+                        end
                 end % }}}
                 function disp(obj) % {{{
                         disp(sprintf('   Thermal solution parameters:'));