Changeset 16040

Timestamp:

08/30/13 11:14:55 (12 years ago)

Author:

seroussi

Message:

CHG: check that you have watercolumn in model

Location:

issm/trunk-jpl/src/m/classes

Files:

• initialization.m (modified) (1 diff)
• initialization.py (modified) (1 diff)

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

@@ -56 +56 @@
                         if (ismember(EnthalpyAnalysisEnum(),analyses) & md.thermal.isenthalpy) | solution==EnthalpySolutionEnum(),
                                 md = checkfield(md,'initialization.waterfraction','>=',0,'size',[md.mesh.numberofvertices 1]);
+                                md = checkfield(md,'initialization.watercolumn'  ,'>=',0,'size',[md.mesh.numberofvertices 1]);
                         end
                         if ismember(HydrologyShreveAnalysisEnum(),analyses),

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

@@ -69 +69 @@
                 if (EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy) or solution==EnthalpySolutionEnum():
                         md = checkfield(md,'initialization.waterfraction','>=',0,'size',[md.mesh.numberofvertices])
+                        md = checkfield(md,'initialization.watercolumn'  ,'>=',0,'size',[md.mesh.numberofvertices])
                 if HydrologyShreveAnalysisEnum() in analyses:
                         md = checkfield(md,'initialization.watercolumn','NaN',1,'size',[md.mesh.numberofvertices])