Ignore:
Timestamp:
12/16/15 14:34:56 (9 years ago)
Author:
Mathieu Morlighem
Message:

CHG: added check on Inf values in checkconsistency

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/trunk-jpl/src/py3/classes/mismipbasalforcings.py

    r19895 r19898  
    5858        if MasstransportAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and md.transient.ismasstransport==0):
    5959
    60             md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
     60            md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
    6161            md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',[1])
    6262            md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',[1])
     
    6565        if BalancethicknessAnalysisEnum() in analyses:
    6666
    67             md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices])
     67            md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    6868            md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',[1])
    6969            md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',[1])
     
    7272        if ThermalAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and md.transient.isthermal==0):
    7373
    74             md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
     74            md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
    7575            md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',[1])
    7676            md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',[1])
    7777            md = checkfield(md,'fieldname','basalforcings.upperdepth_melt','<=',0,'numel',[1])
    78             md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0)
     78            md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0)
    7979        return md
    8080    # }}}
Note: See TracChangeset for help on using the changeset viewer.