Changeset 19897


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

CHG: added check on Inf values in checkconsistency

Location:
issm/trunk-jpl/src/m
Files:
74 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk-jpl/src/m/classes/SMBcomponents.m

    r19527 r19897  
    5050
    5151                        if ismember(MasstransportAnalysisEnum(),analyses),
    52                                 md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1);
     52                                md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1,'Inf',1);
    5353                        end
    5454                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    55                                 md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices 1],'NaN',1);
     55                                md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    5656                        end
    5757                        if ismember(MasstransportAnalysisEnum(),analyses),
    58                                 md = checkfield(md,'fieldname','smb.runoff','timeseries',1,'NaN',1);
     58                                md = checkfield(md,'fieldname','smb.runoff','timeseries',1,'NaN',1,'Inf',1);
    5959                        end
    6060                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    61                                 md = checkfield(md,'fieldname','smb.runoff','size',[md.mesh.numberofvertices 1],'NaN',1);
     61                                md = checkfield(md,'fieldname','smb.runoff','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    6262                        end
    6363                        if ismember(MasstransportAnalysisEnum(),analyses),
    64                                 md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1);
     64                                md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1,'Inf',1);
    6565                        end
    6666                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    67                                 md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices 1],'NaN',1);
     67                                md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    6868                        end
    6969                        md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);

  • issm/trunk-jpl/src/m/classes/SMBcomponents.py

    r19527 r19897  
    5656
    5757                if MasstransportAnalysisEnum() in analyses:
    58                         md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1)
     58                        md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1,'Inf',1)
    5959
    6060                if BalancethicknessAnalysisEnum() in analyses:
    61                         md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices],'NaN',1)
     61                        md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    6262
    6363                if MasstransportAnalysisEnum() in analyses:
    64                         md = checkfield(md,'fieldname','smb.runoff','timeseries',1,'NaN',1)
     64                        md = checkfield(md,'fieldname','smb.runoff','timeseries',1,'NaN',1,'Inf',1)
    6565
    6666                if BalancethicknessAnalysisEnum() in analyses:
    67                         md = checkfield(md,'fieldname','smb.runoff','size',[md.mesh.numberofvertices],'NaN',1)
     67                        md = checkfield(md,'fieldname','smb.runoff','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    6868
    6969                if MasstransportAnalysisEnum() in analyses:
    70                         md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1)
     70                        md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1,'Inf',1)
    7171
    7272                if BalancethicknessAnalysisEnum() in analyses:
    73                         md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices],'NaN',1)
     73                        md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    7474               
    7575                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)

  • issm/trunk-jpl/src/m/classes/SMBd18opdd.m

    r19527 r19897  
    6969                        if ismember(MasstransportAnalysisEnum(),analyses),
    7070                                md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',1);
    71                                 md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'size',[md.mesh.numberofvertices 1]);
    72                                 md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'size',[md.mesh.numberofvertices 1]);
     71                                md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     72                                md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    7373                                md = checkfield(md,'fieldname','smb.rlaps','>=',0,'numel',1);
    7474                                md = checkfield(md,'fieldname','smb.rlapslgm','>=',0,'numel',1);
    7575                                if(self.isd18opd==1)
    76                                         md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    77                                         md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    78                                         md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'size',[2,NaN],'singletimeseries',1);
     76                                        md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     77                                        md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     78                                        md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
    7979                                        md = checkfield(md,'fieldname','smb.dpermil','>=',0,'numel',1);
    8080                                end

  • issm/trunk-jpl/src/m/classes/SMBd18opdd.py

    r19527 r19897  
    9191                if MasstransportAnalysisEnum() in analyses:
    9292                        md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',[1])
    93                         md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'size',[md.mesh.numberofvertices,1])
    94                         md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'size',[md.mesh.numberofvertices,1])
     93                        md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
     94                        md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
    9595                        md = checkfield(md,'fieldname','smb.rlaps','>=',0,'numel',[1])
    9696                        md = checkfield(md,'fieldname','smb.rlapslgm','>=',0,'numel',[1])
    9797
    9898                        if self.isd18opd:
    99                                 md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    100                                 md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    101                                 md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
     99                                md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     100                                md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     101                                md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
    102102                                md = checkfield(md,'fieldname','smb.dpermil','>=',0,'numel',[1])
    103103               

  • issm/trunk-jpl/src/m/classes/SMBforcing.m

    r19879 r19897  
    5555
    5656                        if ismember(MasstransportAnalysisEnum(),analyses),
    57                                 md = checkfield(md,'fieldname','smb.mass_balance','timeseries',1,'NaN',1);
     57                                md = checkfield(md,'fieldname','smb.mass_balance','timeseries',1,'NaN',1,'Inf',1);
    5858                        end
    5959                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    60                                 md = checkfield(md,'fieldname','smb.mass_balance','size',[md.mesh.numberofvertices 1],'NaN',1);
     60                                md = checkfield(md,'fieldname','smb.mass_balance','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    6161                        end
    6262                        md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);

  • issm/trunk-jpl/src/m/classes/SMBforcing.py

    r19527 r19897  
    4343
    4444                if MasstransportAnalysisEnum() in analyses:
    45                         md = checkfield(md,'fieldname','smb.mass_balance','timeseries',1,'NaN',1)
     45                        md = checkfield(md,'fieldname','smb.mass_balance','timeseries',1,'NaN',1,'Inf',1)
    4646
    4747                if BalancethicknessAnalysisEnum() in analyses:
    48                         md = checkfield(md,'fieldname','smb.mass_balance','size',[md.mesh.numberofvertices],'NaN',1)
     48                        md = checkfield(md,'fieldname','smb.mass_balance','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    4949
    5050                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)

  • issm/trunk-jpl/src/m/classes/SMBgemb.m

    r19587 r19897  
    153153                md = checkfield(md,'fieldname','smb.isturbulentflux','values',[0 1]);
    154154               
    155                 md = checkfield(md,'fieldname','smb.Ta','timeseries',1,'NaN',1,'>',273-60,'<',273+60); %60 celsius max value
    156                 md = checkfield(md,'fieldname','smb.V','timeseries',1,'NaN',1,'>=',0,'<',45); %max 500 km/h
    157                 md = checkfield(md,'fieldname','smb.dswrf','timeseries',1,'NaN',1,'>=',0,'<=',1400);
    158                 md = checkfield(md,'fieldname','smb.dlwrf','timeseries',1,'NaN',1,'>=',0);
    159                 md = checkfield(md,'fieldname','smb.P','timeseries',1,'NaN',1,'>=',0,'<=',100);
    160                 md = checkfield(md,'fieldname','smb.eAir','timeseries',1,'NaN',1);
    161                
    162                 md = checkfield(md,'fieldname','smb.Tmean','size',[md.mesh.numberofelements 1],'NaN',1,'>',273-60,'<',273+60); %60 celsius max value
    163                 md = checkfield(md,'fieldname','smb.C','size',[md.mesh.numberofelements 1],'NaN',1,'>=',0);
    164                 md = checkfield(md,'fieldname','smb.Tz','size',[md.mesh.numberofelements 1],'NaN',1,'>=',0,'<=',5000);
    165                 md = checkfield(md,'fieldname','smb.Vz','size',[md.mesh.numberofelements 1],'NaN',1,'>=',0,'<=',5000);
    166                
    167                 md = checkfield(md,'fieldname','smb.aIdx','NaN',1,'values',[1,2,3,4]);
    168                 md = checkfield(md,'fieldname','smb.swIdx','NaN',1,'values',[0,1]);
    169                 md = checkfield(md,'fieldname','smb.denIdx','NaN',1,'values',[1,2,3,4,5]);
    170                
    171                 md = checkfield(md,'fieldname','smb.zTop','NaN',1,'>=',0);
    172                 md = checkfield(md,'fieldname','smb.dzTop','NaN',1,'>',0);
    173                 md = checkfield(md,'fieldname','smb.dzMin','NaN',1,'>',0);
    174                 md = checkfield(md,'fieldname','smb.zY','NaN',1,'>=',1);
    175                 md = checkfield(md,'fieldname','smb.outputFreq','NaN',1,'>',0,'<',10*365); %10 years max
    176                 md = checkfield(md,'fieldname','smb.InitDensityScaling','NaN',1,'>=',0,'<=',1);
     155                md = checkfield(md,'fieldname','smb.Ta','timeseries',1,'NaN',1,'Inf',1,'>',273-60,'<',273+60); %60 celsius max value
     156                md = checkfield(md,'fieldname','smb.V','timeseries',1,'NaN',1,'Inf',1,'>=',0,'<',45); %max 500 km/h
     157                md = checkfield(md,'fieldname','smb.dswrf','timeseries',1,'NaN',1,'Inf',1,'>=',0,'<=',1400);
     158                md = checkfield(md,'fieldname','smb.dlwrf','timeseries',1,'NaN',1,'Inf',1,'>=',0);
     159                md = checkfield(md,'fieldname','smb.P','timeseries',1,'NaN',1,'Inf',1,'>=',0,'<=',100);
     160                md = checkfield(md,'fieldname','smb.eAir','timeseries',1,'NaN',1,'Inf',1);
     161               
     162                md = checkfield(md,'fieldname','smb.Tmean','size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1,'>',273-60,'<',273+60); %60 celsius max value
     163                md = checkfield(md,'fieldname','smb.C','size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1,'>=',0);
     164                md = checkfield(md,'fieldname','smb.Tz','size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1,'>=',0,'<=',5000);
     165                md = checkfield(md,'fieldname','smb.Vz','size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1,'>=',0,'<=',5000);
     166               
     167                md = checkfield(md,'fieldname','smb.aIdx','NaN',1,'Inf',1,'values',[1,2,3,4]);
     168                md = checkfield(md,'fieldname','smb.swIdx','NaN',1,'Inf',1,'values',[0,1]);
     169                md = checkfield(md,'fieldname','smb.denIdx','NaN',1,'Inf',1,'values',[1,2,3,4,5]);
     170               
     171                md = checkfield(md,'fieldname','smb.zTop','NaN',1,'Inf',1,'>=',0);
     172                md = checkfield(md,'fieldname','smb.dzTop','NaN',1,'Inf',1,'>',0);
     173                md = checkfield(md,'fieldname','smb.dzMin','NaN',1,'Inf',1,'>',0);
     174                md = checkfield(md,'fieldname','smb.zY','NaN',1,'Inf',1,'>=',1);
     175                md = checkfield(md,'fieldname','smb.outputFreq','NaN',1,'Inf',1,'>',0,'<',10*365); %10 years max
     176                md = checkfield(md,'fieldname','smb.InitDensityScaling','NaN',1,'Inf',1,'>=',0,'<=',1);
    177177               
    178178                switch self.aIdx,
    179179                case {1 2}
    180                         md = checkfield(md,'fieldname','smb.aSnow','NaN',1,'>=',.64,'<=',.89);
    181                         md = checkfield(md,'fieldname','smb.aIce','NaN',1,'>=',.27,'<=',.58);
     180                        md = checkfield(md,'fieldname','smb.aSnow','NaN',1,'Inf',1,'>=',.64,'<=',.89);
     181                        md = checkfield(md,'fieldname','smb.aIce','NaN',1,'Inf',1,'>=',.27,'<=',.58);
    182182                case 3
    183                         md = checkfield(md,'fieldname','smb.cldFrac','NaN',1,'>=',0,'<=',1);
     183                        md = checkfield(md,'fieldname','smb.cldFrac','NaN',1,'Inf',1,'>=',0,'<=',1);
    184184                case 4
    185                         md = checkfield(md,'fieldname','smb.t0wet','NaN',1,'>=',15,'<=',21.9);
    186                         md = checkfield(md,'fieldname','smb.t0dry','NaN',1,'>=',30,'<=',30);
    187                         md = checkfield(md,'fieldname','smb.K','NaN',1,'>=',7,'<=',7);
     185                        md = checkfield(md,'fieldname','smb.t0wet','NaN',1,'Inf',1,'>=',15,'<=',21.9);
     186                        md = checkfield(md,'fieldname','smb.t0dry','NaN',1,'Inf',1,'>=',30,'<=',30);
     187                        md = checkfield(md,'fieldname','smb.K','NaN',1,'Inf',1,'>=',7,'<=',7);
    188188                end
    189189               

  • issm/trunk-jpl/src/m/classes/SMBgradients.m

    r19527 r19897  
    4242
    4343                        if ismember(MasstransportAnalysisEnum(),analyses),
    44                                 md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1);
    45                                 md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1);
    46                                 md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1);
    47                                 md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1);
     44                                md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1,'Inf',1);
     45                                md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1,'Inf',1);
     46                                md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1,'Inf',1);
     47                                md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1,'Inf',1);
    4848                        end
    4949                        md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);

  • issm/trunk-jpl/src/m/classes/SMBgradients.py

    r19527 r19897  
    4949
    5050                if MasstransportAnalysisEnum() in analyses:
    51                         md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1)
    52                         md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1)
    53                         md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1)
    54                         md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1)
     51                        md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1,'Inf',1)
     52                        md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1,'Inf',1)
     53                        md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1,'Inf',1)
     54                        md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1,'Inf',1)
    5555
    5656                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)

  • issm/trunk-jpl/src/m/classes/SMBhenning.m

    r19527 r19897  
    3636
    3737                        if ismember(MasstransportAnalysisEnum(),analyses),
    38                                 md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1);
     38                                md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1,'Inf',1);
    3939                        end
    4040                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    41                                 md = checkfield(md,'fieldname','smb.smbref','size',[md.mesh.numberofvertices 1],'NaN',1);
     41                                md = checkfield(md,'fieldname','smb.smbref','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    4242                        end
    4343                        md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);

  • issm/trunk-jpl/src/m/classes/SMBmeltcomponents.m

    r19527 r19897  
    5454
    5555                        if ismember(MasstransportAnalysisEnum(),analyses),
    56                                 md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1);
     56                                md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1,'Inf',1);
    5757                        end
    5858                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    59                                 md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices 1],'NaN',1);
     59                                md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    6060                        end
    6161                        if ismember(MasstransportAnalysisEnum(),analyses),
    62                                 md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1);
     62                                md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1,'Inf',1);
    6363                        end
    6464                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    65                                 md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices 1],'NaN',1);
     65                                md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    6666                        end
    6767                        if ismember(MasstransportAnalysisEnum(),analyses),
    68                                 md = checkfield(md,'fieldname','smb.refreeze','timeseries',1,'NaN',1);
     68                                md = checkfield(md,'fieldname','smb.refreeze','timeseries',1,'NaN',1,'Inf',1);
    6969                        end
    7070                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    71                                 md = checkfield(md,'fieldname','smb.refreeze','size',[md.mesh.numberofvertices 1],'NaN',1);
     71                                md = checkfield(md,'fieldname','smb.refreeze','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    7272                        end
    7373                        if ismember(MasstransportAnalysisEnum(),analyses),
    74                                 md = checkfield(md,'fieldname','smb.melt','timeseries',1,'NaN',1);
     74                                md = checkfield(md,'fieldname','smb.melt','timeseries',1,'NaN',1,'Inf',1);
    7575                        end
    7676                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    77                                 md = checkfield(md,'fieldname','smb.melt','size',[md.mesh.numberofvertices 1],'NaN',1);
     77                                md = checkfield(md,'fieldname','smb.melt','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    7878                        end
    7979                        md = checkfield(md,'fieldname','smb.requested_outputs','stringrow',1);

  • issm/trunk-jpl/src/m/classes/SMBmeltcomponents.py

    r19527 r19897  
    6262
    6363                if MasstransportAnalysisEnum() in analyses:
    64                         md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1)
     64                        md = checkfield(md,'fieldname','smb.accumulation','timeseries',1,'NaN',1,'Inf',1)
    6565
    6666                if BalancethicknessAnalysisEnum() in analyses:
    67                         md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices],'NaN',1)
     67                        md = checkfield(md,'fieldname','smb.accumulation','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    6868
    6969                if MasstransportAnalysisEnum() in analyses:
    70                         md = checkfield(md,'fieldname','smb.melt','timeseries',1,'NaN',1)
     70                        md = checkfield(md,'fieldname','smb.melt','timeseries',1,'NaN',1,'Inf',1)
    7171
    7272                if BalancethicknessAnalysisEnum() in analyses:
    73                         md = checkfield(md,'fieldname','smb.melt','size',[md.mesh.numberofvertices],'NaN',1)
     73                        md = checkfield(md,'fieldname','smb.melt','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    7474
    7575                if MasstransportAnalysisEnum() in analyses:
    76                         md = checkfield(md,'fieldname','smb.refreeze','timeseries',1,'NaN',1)
     76                        md = checkfield(md,'fieldname','smb.refreeze','timeseries',1,'NaN',1,'Inf',1)
    7777
    7878                if BalancethicknessAnalysisEnum() in analyses:
    79                         md = checkfield(md,'fieldname','smb.refreeze','size',[md.mesh.numberofvertices],'NaN',1)
     79                        md = checkfield(md,'fieldname','smb.refreeze','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    8080
    8181                if MasstransportAnalysisEnum() in analyses:
    82                         md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1)
     82                        md = checkfield(md,'fieldname','smb.evaporation','timeseries',1,'NaN',1,'Inf',1)
    8383
    8484                if BalancethicknessAnalysisEnum() in analyses:
    85                         md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices],'NaN',1)
     85                        md = checkfield(md,'fieldname','smb.evaporation','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    8686
    8787                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)

  • issm/trunk-jpl/src/m/classes/SMBpdd.m

    r19527 r19897  
    7878                        if ismember(MasstransportAnalysisEnum(),analyses),
    7979                                md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',1);
    80                                 md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'size',[md.mesh.numberofvertices 1]);
    81                                 md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'size',[md.mesh.numberofvertices 1]);
     80                                md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     81                                md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    8282                                md = checkfield(md,'fieldname','smb.rlaps','>=',0,'numel',1);
    8383                                md = checkfield(md,'fieldname','smb.rlapslgm','>=',0,'numel',1);
    8484                                if(self.isdelta18o==0 & self.ismungsm==0)
    85                                         md = checkfield(md,'fieldname','smb.monthlytemperatures','timeseries',1,'NaN',1);
    86                                         md = checkfield(md,'fieldname','smb.precipitation','timeseries',1,'NaN',1);
     85                                        md = checkfield(md,'fieldname','smb.monthlytemperatures','timeseries',1,'NaN',1,'Inf',1);
     86                                        md = checkfield(md,'fieldname','smb.precipitation','timeseries',1,'NaN',1,'Inf',1);
    8787                                elseif(self.isdelta18o==1)
    88                                         md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'size',[2,NaN],'singletimeseries',1);
    89                                         md = checkfield(md,'fieldname','smb.delta18o_surface','NaN',1,'size',[2,NaN],'singletimeseries',1);
    90                                         md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    91                                         md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    92                                         md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    93                                         md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);                                       
    94                                         md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'size',[2,NaN],'singletimeseries',1);
    95                                         md = checkfield(md,'fieldname','smb.sealev','NaN',1,'size',[2,NaN],'singletimeseries',1);
     88                                        md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
     89                                        md = checkfield(md,'fieldname','smb.delta18o_surface','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
     90                                        md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     91                                        md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     92                                        md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     93                                        md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);                                       
     94                                        md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
     95                                        md = checkfield(md,'fieldname','smb.sealev','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
    9696                                elseif(self.ismungsm==1)
    97                                         md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    98                                         md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    99                                         md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);
    100                                         md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'timeseries',1);                                       
    101                                         md = checkfield(md,'fieldname','smb.Pfac','NaN',1,'size',[2,NaN],'singletimeseries',1);
    102                                         md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'size',[2,NaN],'singletimeseries',1);
    103                                         md = checkfield(md,'fieldname','smb.sealev','NaN',1,'size',[2,NaN],'singletimeseries',1);
     97                                        md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     98                                        md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     99                                        md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);
     100                                        md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1 12],'NaN',1,'Inf',1,'timeseries',1);                                       
     101                                        md = checkfield(md,'fieldname','smb.Pfac','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
     102                                        md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
     103                                        md = checkfield(md,'fieldname','smb.sealev','NaN',1,'Inf',1,'size',[2,NaN],'singletimeseries',1);
    104104                                end
    105105                        end

  • issm/trunk-jpl/src/m/classes/SMBpdd.py

    r19527 r19897  
    120120                if MasstransportAnalysisEnum() in analyses:
    121121                        md = checkfield(md,'fieldname','smb.desfac','<=',1,'numel',[1])
    122                         md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'size',[md.mesh.numberofvertices,1])
    123                         md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'size',[md.mesh.numberofvertices,1])
     122                        md = checkfield(md,'fieldname','smb.s0p','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
     123                        md = checkfield(md,'fieldname','smb.s0t','>=',0,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
    124124                        md = checkfield(md,'fieldname','smb.rlaps','>=',0,'numel',[1])
    125125                        md = checkfield(md,'fieldname','smb.rlapslgm','>=',0,'numel',[1])
    126126
    127127                        if (self.isdelta18o==0 and self.ismungsm==0):
    128                                 md = checkfield(md,'fieldname','smb.monthlytemperatures','NaN',1,'timeseries',1)
    129                                 md = checkfield(md,'fieldname','smb.precipitation','NaN',1,'timeseries',1)
     128                                md = checkfield(md,'fieldname','smb.monthlytemperatures','NaN',1,'Inf',1,'timeseries',1)
     129                                md = checkfield(md,'fieldname','smb.precipitation','NaN',1,'Inf',1,'timeseries',1)
    130130                        elif self.isdelta18o:
    131                                 md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
    132                                 md = checkfield(md,'fieldname','smb.delta18o_surface','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
    133                                 md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    134                                 md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    135                                 md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    136                                 md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)                                       
    137                                 md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
    138                                 md = checkfield(md,'fieldname','smb.sealev','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
     131                                md = checkfield(md,'fieldname','smb.delta18o','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
     132                                md = checkfield(md,'fieldname','smb.delta18o_surface','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
     133                                md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     134                                md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     135                                md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     136                                md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)                                       
     137                                md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
     138                                md = checkfield(md,'fieldname','smb.sealev','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
    139139                        elif self.ismungsm:
    140                                 md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    141                                 md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    142                                 md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)
    143                                 md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'timeseries',1)                                       
    144                                 md = checkfield(md,'fieldname','smb.Pfac','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
    145                                 md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
    146                                 md = checkfield(md,'fieldname','smb.sealev','NaN',1,'size',[2,numpy.nan],'singletimeseries',1)
     140                                md = checkfield(md,'fieldname','smb.temperatures_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     141                                md = checkfield(md,'fieldname','smb.temperatures_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     142                                md = checkfield(md,'fieldname','smb.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)
     143                                md = checkfield(md,'fieldname','smb.precipitations_lgm','size',[md.mesh.numberofvertices+1,12],'NaN',1,'Inf',1,'timeseries',1)                                       
     144                                md = checkfield(md,'fieldname','smb.Pfac','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
     145                                md = checkfield(md,'fieldname','smb.Tdiff','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
     146                                md = checkfield(md,'fieldname','smb.sealev','NaN',1,'Inf',1,'size',[2,numpy.nan],'singletimeseries',1)
    147147
    148148                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)

  • issm/trunk-jpl/src/m/classes/adinversion.m

    r19042 r19897  
    7373                                {'BalancethicknessThickeningRate' 'FrictionCoefficient' 'MaterialsRheologyBbar' 'DamageDbar',...
    7474                                'Vx' 'Vy' 'Thickness' 'BalancethicknessOmega' 'BalancethicknessApparentMassbalance'});
    75                         md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
     75                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1,'Inf',1);
    7676                        md = checkfield(md,'fieldname','inversion.maxsteps','numel',1,'>=',0);
    7777                        md = checkfield(md,'fieldname','inversion.maxiter','numel',1,'>=',0);
    … …  
    8484
    8585                        if solution==BalancethicknessSolutionEnum()
    86                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
    87                                 md = checkfield(md,'fieldname','inversion.surface_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     86                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
     87                                md = checkfield(md,'fieldname','inversion.surface_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    8888                        elseif solution==BalancethicknessSoftSolutionEnum()
    89                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     89                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    9090                        else
    91                                 md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     91                                md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    9292                                if ~strcmp(domaintype(md.mesh),'2Dvertical'),
    93                                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     93                                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    9494                                end
    9595                        end

  • issm/trunk-jpl/src/m/classes/balancethickness.m

    r19879 r19897  
    3232
    3333                        md = checkfield(md,'fieldname','balancethickness.spcthickness');
    34                         md = checkfield(md,'fieldname','balancethickness.thickening_rate','size',[md.mesh.numberofvertices 1],'NaN',1);
     34                        md = checkfield(md,'fieldname','balancethickness.thickening_rate','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    3535                        md = checkfield(md,'fieldname','balancethickness.stabilization','size',[1 1],'values',[0 1 2 3]);
    3636
    37                         %md = checkfield(md,'fieldname','balancethickness.omega','size',[md.mesh.numberofvertices 1],'NaN',1,'>=',0);
     37                        %md = checkfield(md,'fieldname','balancethickness.omega','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1,'>=',0);
    3838                end % }}}
    3939                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/balancethickness.py

    r17497 r19897  
    4343
    4444                md = checkfield(md,'fieldname','balancethickness.spcthickness')
    45                 md = checkfield(md,'fieldname','balancethickness.thickening_rate','size',[md.mesh.numberofvertices],'NaN',1)
     45                md = checkfield(md,'fieldname','balancethickness.thickening_rate','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    4646                md = checkfield(md,'fieldname','balancethickness.stabilization','size',[1],'values',[0,1,2,3])
    4747

  • issm/trunk-jpl/src/m/classes/basalforcings.m

    r19879 r19897  
    5050
    5151                        if ismember(MasstransportAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.ismasstransport==0),
    52                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
    53                                 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1);
     52                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
     53                                md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
    5454                        end
    5555                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    56                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices 1]);
    57                                 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices 1]);
     56                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     57                                md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    5858                        end
    5959                        if ismember(ThermalAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.isthermal==0),
    60                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
    61                                 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1);
    62                                 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0);
     60                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
     61                                md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
     62                                md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0);
    6363                        end
    6464                end % }}}

  • issm/trunk-jpl/src/m/classes/basalforcings.py

    r19048 r19897  
    5555
    5656                if MasstransportAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.ismasstransport):
    57                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
    58                         md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1)
     57                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
     58                        md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
    5959
    6060                if BalancethicknessAnalysisEnum() in analyses:
    61                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices])
    62                         md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices])
     61                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     62                        md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    6363
    6464                if ThermalAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.isthermal):
    65                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
    66                         md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1)
    67                         md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0)
     65                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
     66                        md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
     67                        md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0)
    6868
    6969                return md

  • issm/trunk-jpl/src/m/classes/calving.m

    r19879 r19897  
    4747                        md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
    4848                        md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
    49                         md = checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1);
    50                         md = checkfield(md,'fieldname','calving.meltingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1);
     49                        md = checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1,'Inf',1);
     50                        md = checkfield(md,'fieldname','calving.meltingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1,'Inf',1);
    5151                end % }}}
    5252                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/calving.py

    r19158 r19897  
    5555                md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1)
    5656                md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
    57                 md = checkfield(md,'fieldname','calving.calvingrate','>=',0,'timeseries',1,'NaN',1);
    58                 md = checkfield(md,'fieldname','calving.meltingrate','>=',0,'timeseries',1,'NaN',1);
     57                md = checkfield(md,'fieldname','calving.calvingrate','>=',0,'timeseries',1,'NaN',1,'Inf',1);
     58                md = checkfield(md,'fieldname','calving.meltingrate','>=',0,'timeseries',1,'NaN',1,'Inf',1);
    5959
    6060                return md

  • issm/trunk-jpl/src/m/classes/calvingdev.m

    r19158 r19897  
    5050                        md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
    5151                        md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
    52                         md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'size',[md.mesh.numberofvertices 1],'>=',0);
     52                        md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>=',0);
    5353                end % }}}
    5454                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/calvinglevermann.m

    r19167 r19897  
    5050                        md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
    5151                        md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
    52                         md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'size',[md.mesh.numberofvertices 1],'>=',0);
     52                        md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>=',0);
    5353                end % }}}
    5454                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/calvinglevermann.py

    r19167 r19897  
    5656                md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
    5757                md = checkfield(md,'fieldname','calving.coeff','size',[md.mesh.numberofvertices],'>',0)
    58                 md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'size',[md.mesh.numberofvertices],'>=',0)
     58                md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices],'>=',0)
    5959                return md
    6060        # }}}

  • issm/trunk-jpl/src/m/classes/calvingpi.m

    r19158 r19897  
    5050                        md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
    5151                        md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
    52                         md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'size',[md.mesh.numberofvertices 1],'>=',0);
     52                        md = checkfield(md,'fieldname','calving.meltingrate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>=',0);
    5353                end % }}}
    5454                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/friction.m

    r19879 r19897  
    4040                        %Early return
    4141                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    42                         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1);
    43                         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
    44                         md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]);
     42                        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1);
     43                        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     44                        md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    4545                end % }}}
    4646                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/friction.py

    r19707 r19897  
    4545                        return md
    4646
    47                 md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1)
    48                 md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements])
    49                 md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements])
     47                md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1)
     48                md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements])
     49                md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements])
    5050
    5151                return md

  • issm/trunk-jpl/src/m/classes/frictioncoulomb.m

    r19707 r19897  
    4242                        %Early return
    4343                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    44                         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1);
    45                         md = checkfield(md,'fieldname','friction.coefficientcoulomb','timeseries',1,'NaN',1);
    46                         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
    47                         md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]);
     44                        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1);
     45                        md = checkfield(md,'fieldname','friction.coefficientcoulomb','timeseries',1,'NaN',1,'Inf',1);
     46                        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     47                        md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    4848                end % }}}
    4949                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/frictioncoulomb.py

    r19851 r19897  
    4848            return md
    4949
    50         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1)
    51         md = checkfield(md,'fieldname','friction.coefficientcoulomb','timeseries',1,'NaN',1)
    52         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements])
    53         md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements])
     50        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1)
     51        md = checkfield(md,'fieldname','friction.coefficientcoulomb','timeseries',1,'NaN',1,'Inf',1)
     52        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements])
     53        md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements])
    5454
    5555        return md

  • issm/trunk-jpl/src/m/classes/frictionhydro.m

    r19195 r19897  
    2929                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    3030                        md = checkfield(md,'fieldname','friction.Coupling','numel',[1],'values',[0 1]);
    31                         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
    32                         md = checkfield(md,'fieldname','friction.C','NaN',1,'size',[md.mesh.numberofelements 1]);
    33                         md = checkfield(md,'fieldname','friction.As','NaN',1,'size',[md.mesh.numberofelements 1]);
     31                        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     32                        md = checkfield(md,'fieldname','friction.C','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     33                        md = checkfield(md,'fieldname','friction.As','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    3434                        if self.Coupling==0,
    35                                 md = checkfield(md,'fieldname','friction.effective_pressure','NaN',1,'timeseries',1);
     35                                md = checkfield(md,'fieldname','friction.effective_pressure','NaN',1,'Inf',1,'timeseries',1);
    3636            end
    3737                end % }}}

  • issm/trunk-jpl/src/m/classes/frictionsommers.m

    r19738 r19897  
    2929                        %Early return
    3030                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    31                         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1);
     31                        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1);
    3232                end % }}}
    3333                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/frictiontemp.m

    r19042 r19897  
    3434                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    3535
    36                         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1);
    37                         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
    38                         md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]);
    39                         md = checkfield(md,'fieldname','friction.gamma','NaN',1,'numel',1,'>',0.);
     36                        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1);
     37                        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     38                        md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     39                        md = checkfield(md,'fieldname','friction.gamma','NaN',1,'Inf',1,'numel',1,'>',0.);
    4040
    4141                        %Check that temperature is provided
    42                         md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'size',[md.mesh.numberofvertices 1]);
     42                        md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    4343                end % }}}
    4444                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/frictionwaterlayer.m

    r19373 r19897  
    3131                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    3232
    33                         md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1);
    34                         md = checkfield(md,'fieldname','friction.f','size',[1 1],'NaN',1);
    35                         md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]);
    36                         md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]);
     33                        md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1,'Inf',1);
     34                        md = checkfield(md,'fieldname','friction.f','size',[1 1],'NaN',1,'Inf',1);
     35                        md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
     36                        md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    3737                        md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1,'>=',0.);
    3838

  • issm/trunk-jpl/src/m/classes/frictionweertman.m

    r19048 r19897  
    3838                        %Early return
    3939                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    40                         md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1);
    41                         md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements 1]);
     40                        md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1,'Inf',1);
     41                        md = checkfield(md,'fieldname','friction.m','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    4242                end % }}}
    4343                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/frictionweertman.py

    r19048 r19897  
    3737                        return md
    3838
    39                 md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1)
    40                 md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements])
     39                md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1,'Inf',1)
     40                md = checkfield(md,'fieldname','friction.m','NaN',1,'Inf',1,'size',[md.mesh.numberofelements])
    4141
    4242                return md

  • issm/trunk-jpl/src/m/classes/frictionweertmantemp.m

    r19040 r19897  
    2626                        %Early return
    2727                        if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end
    28                         md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1);
    29                         md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements 1]);
     28                        md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1,'Inf',1);
     29                        md = checkfield(md,'fieldname','friction.m','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
    3030                end % }}}
    3131                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/geometry.m

    r19879 r19897  
    6565
    6666                        if (solution==TransientSolutionEnum() & md.transient.isgia) | (solution==GiaSolutionEnum()),
    67                                 md = checkfield(md,'fieldname','geometry.thickness','timeseries',1,'NaN',1,'>=',0);
     67                                md = checkfield(md,'fieldname','geometry.thickness','timeseries',1,'NaN',1,'Inf',1,'>=',0);
    6868                        else
    69                                 md = checkfield(md,'fieldname','geometry.surface'  ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
    70                                 md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'size',[md.mesh.numberofvertices 1]);
    71                                 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
     69                                md = checkfield(md,'fieldname','geometry.surface'  ,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     70                                md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     71                                md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>',0);
    7272                                if any(abs(self.thickness-self.surface+self.base)>10^-9),
    7373                                        md = checkmessage(md,['equality thickness=surface-base violated']);
    7474                                end
    7575                                if solution==TransientSolutionEnum() & md.transient.isgroundingline,
    76                                         md = checkfield(md,'fieldname','geometry.bed','NaN',1,'size',[md.mesh.numberofvertices 1]);
     76                                        md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    7777                                end
    7878                        end

  • issm/trunk-jpl/src/m/classes/geometry.py

    r19573 r19897  
    4646        def checkconsistency(self,md,solution,analyses):    # {{{
    4747
    48                 md = checkfield(md,'fieldname','geometry.surface'  ,'NaN',1,'size',[md.mesh.numberofvertices])
    49                 md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'size',[md.mesh.numberofvertices])
    50                 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices],'>',0,'timeseries',1)
     48                md = checkfield(md,'fieldname','geometry.surface'  ,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     49                md = checkfield(md,'fieldname','geometry.base'      ,'NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     50                md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices],'>',0,'timeseries',1)
    5151                if any(abs(self.thickness-self.surface+self.base)>10**-9):
    5252                        md.checkmessage("equality thickness=surface-base violated")
    5353                if solution==TransientSolutionEnum() and md.transient.isgroundingline:
    54                         md = checkfield(md,'fieldname','geometry.bed','NaN',1,'size',[md.mesh.numberofvertices])
     54                        md = checkfield(md,'fieldname','geometry.bed','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    5555
    5656                return md

  • issm/trunk-jpl/src/m/classes/gia.m

    r19879 r19897  
    4444
    4545                        if ~ismember(GiaAnalysisEnum(),analyses), return; end
    46                         md = checkfield(md,'fieldname','gia.mantle_viscosity','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
    47                         md = checkfield(md,'fieldname','gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0);
     46                        md = checkfield(md,'fieldname','gia.mantle_viscosity','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>',0);
     47                        md = checkfield(md,'fieldname','gia.lithosphere_thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1],'>',0);
    4848                        md = checkfield(md,'fieldname','gia.cross_section_shape','numel',[1],'values',[1,2]);
    4949

  • issm/trunk-jpl/src/m/classes/gia.py

    r19048 r19897  
    4848                        return md
    4949               
    50                 md = checkfield(md,'fieldname','gia.mantle_viscosity','NaN',1,'size',[md.mesh.numberofvertices,1],'>',0)
    51                 md = checkfield(md,'fieldname','gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices,1],'>',0)
     50                md = checkfield(md,'fieldname','gia.mantle_viscosity','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1],'>',0)
     51                md = checkfield(md,'fieldname','gia.lithosphere_thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1],'>',0)
    5252                md = checkfield(md,'fieldname','gia.cross_section_shape','numel',[1],'values',[1,2])
    5353

  • issm/trunk-jpl/src/m/classes/hydrologydc.m

    r19048 r19897  
    114114                                md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0,'numel',1);
    115115                        end
    116                         md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'timeseries',1);
     116                        md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'Inf',1,'timeseries',1);
    117117
    118118                        md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1);

  • issm/trunk-jpl/src/m/classes/hydrologydc.py

    r19048 r19897  
    163163                        md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0.,'numel',[1])
    164164
    165                 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'timeseries',1)
     165                md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'Inf',1,'timeseries',1)
    166166                md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1)
    167167                md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0.,'numel',[1])

  • issm/trunk-jpl/src/m/classes/hydrologysommers.m

    r19757 r19897  
    3838                        end
    3939
    40                         md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices 1],'NaN',1);
    41                         md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1);
    42                         md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements 1],'NaN',1);
    43                         md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1);
    44                         md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'timeseries',1);
    45                         md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'timeseries',1);
    46                         md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements 1],'NaN',1);
    47                         md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1);
     40                        md = checkfield(md,'fieldname','hydrology.head','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
     41                        md = checkfield(md,'fieldname','hydrology.gap_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1);
     42                        md = checkfield(md,'fieldname','hydrology.bump_spacing','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1);
     43                        md = checkfield(md,'fieldname','hydrology.bump_height','>=',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1);
     44                        md = checkfield(md,'fieldname','hydrology.englacial_input','>=',0,'NaN',1,'Inf',1,'timeseries',1);
     45                        md = checkfield(md,'fieldname','hydrology.moulin_input','>=',0,'NaN',1,'Inf',1,'timeseries',1);
     46                        md = checkfield(md,'fieldname','hydrology.reynolds','>',0,'size',[md.mesh.numberofelements 1],'NaN',1,'Inf',1);
     47                        md = checkfield(md,'fieldname','hydrology.neumannflux','timeseries',1,'NaN',1,'Inf',1);
    4848                        md = checkfield(md,'fieldname','hydrology.spchead','size',[md.mesh.numberofvertices 1]);
    4949                end % }}}

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

    r19879 r19897  
    6565                        if ismember(StressbalanceAnalysisEnum(),analyses)
    6666                                if ~(isnan(md.initialization.vx) | isnan(md.initialization.vy)),
    67                                         md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices 1]);
    68                                         md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices 1]);
     67                                        md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     68                                        md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    6969                                end
    7070                        end
    7171                        if ismember(MasstransportAnalysisEnum(),analyses),
    72                                 md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices 1]);
    73                                 md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices 1]);
     72                                md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     73                                md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    7474                        end
    7575                        if ismember(BalancethicknessAnalysisEnum(),analyses) & solution==BalancethicknessSolutionEnum,
    76                                 md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices 1]);
    77                                 md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices 1]);
     76                                md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     77                                md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    7878                                %Triangle with zero velocity
    7979                                if any(sum(abs(md.initialization.vx(md.mesh.elements)),2)==0 & sum(abs(md.initialization.vy(md.mesh.elements)),2)==0)
    … …  
    8282                        end
    8383                        if ismember(ThermalAnalysisEnum(),analyses),
    84                                 md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices 1]);
    85                                 md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices 1]);
     84                                md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     85                                md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    8686                                if dimension(md.mesh)==3
    87                                         md = checkfield(md,'fieldname','initialization.vz','NaN',1,'size',[md.mesh.numberofvertices 1]);
     87                                        md = checkfield(md,'fieldname','initialization.vz','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    8888                                end
    89                                 md = checkfield(md,'fieldname','initialization.pressure','NaN',1,'size',[md.mesh.numberofvertices 1]);
    90                                 md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'size',[md.mesh.numberofvertices 1]);
     89                                md = checkfield(md,'fieldname','initialization.pressure','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     90                                md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    9191                        end
    9292                        if (ismember(EnthalpyAnalysisEnum(),analyses) & md.thermal.isenthalpy)
    … …  
    9696                        if ismember(HydrologyShreveAnalysisEnum(),analyses),
    9797                                if isa(md.hydrology,'hydrologyshreve'),
    98                                         md = checkfield(md,'fieldname','initialization.watercolumn','NaN',1,'size',[md.mesh.numberofvertices 1]);
     98                                        md = checkfield(md,'fieldname','initialization.watercolumn','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    9999                                end
    100100                        end
    101101                        if ismember(HydrologyDCInefficientAnalysisEnum(),analyses),
    102102                                if isa(md.hydrology,'hydrologydc'),
    103                                         md = checkfield(md,'fieldname','initialization.sediment_head','NaN',1,'size',[md.mesh.numberofvertices 1]);
     103                                        md = checkfield(md,'fieldname','initialization.sediment_head','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    104104                                end
    105105                        end
    … …  
    107107                                if isa(md.hydrology,'hydrologydc'),
    108108                                        if md.hydrology.isefficientlayer==1,
    109                                                 md = checkfield(md,'fieldname','initialization.epl_head','NaN',1,'size',[md.mesh.numberofvertices 1]);
    110                                                 md = checkfield(md,'fieldname','initialization.epl_thickness','NaN',1,'size',[md.mesh.numberofvertices 1]);
     109                                                md = checkfield(md,'fieldname','initialization.epl_head','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     110                                                md = checkfield(md,'fieldname','initialization.epl_thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    111111                                        end
    112112                                end

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

    r19048 r19897  
    7171                if StressbalanceAnalysisEnum() in analyses:
    7272                        if not numpy.any(numpy.logical_or(numpy.isnan(md.initialization.vx),numpy.isnan(md.initialization.vy))):
    73                                 md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices])
    74                                 md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices])
     73                                md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     74                                md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    7575                if MasstransportAnalysisEnum() in analyses:
    76                         md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices])
    77                         md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices])
     76                        md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     77                        md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    7878                if BalancethicknessAnalysisEnum() in analyses:
    79                         md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices])
    80                         md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices])
     79                        md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     80                        md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    8181                        #Triangle with zero velocity
    8282                        if numpy.any(numpy.logical_and(numpy.sum(numpy.abs(md.initialization.vx[md.mesh.elements-1]),axis=1)==0,\
    … …  
    8484                                md.checkmessage("at least one triangle has all its vertices with a zero velocity")
    8585                if ThermalAnalysisEnum() in analyses:
    86                         md = checkfield(md,'fieldname','initialization.vx','NaN',1,'size',[md.mesh.numberofvertices])
    87                         md = checkfield(md,'fieldname','initialization.vy','NaN',1,'size',[md.mesh.numberofvertices])
    88                         md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'size',[md.mesh.numberofvertices])
     86                        md = checkfield(md,'fieldname','initialization.vx','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     87                        md = checkfield(md,'fieldname','initialization.vy','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     88                        md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    8989                        if md.mesh.dimension()==3:
    90                                 md = checkfield(md,'fieldname','initialization.vz','NaN',1,'size',[md.mesh.numberofvertices])
    91                         md = checkfield(md,'fieldname','initialization.pressure','NaN',1,'size',[md.mesh.numberofvertices])
     90                                md = checkfield(md,'fieldname','initialization.vz','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     91                        md = checkfield(md,'fieldname','initialization.pressure','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    9292                        if (EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy):
    9393                                md = checkfield(md,'fieldname','initialization.waterfraction','>=',0,'size',[md.mesh.numberofvertices])
    … …  
    9595                if HydrologyShreveAnalysisEnum() in analyses:
    9696                        if hasattr(md.hydrology,'hydrologyshreve'):
    97                                 md = checkfield(md,'fieldname','initialization.watercolumn','NaN',1,'size',[md.mesh.numberofvertices])
     97                                md = checkfield(md,'fieldname','initialization.watercolumn','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    9898                if HydrologyDCInefficientAnalysisEnum() in analyses:
    9999                        if hasattr(md.hydrology,'hydrologydc'):
    100                                 md = checkfield(md,'fieldname','initialization.sediment_head','NaN',1,'size',[md.mesh.numberofvertices,1])
     100                                md = checkfield(md,'fieldname','initialization.sediment_head','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
    101101                if HydrologyDCEfficientAnalysisEnum() in analyses:
    102102                        if hasattr(md.hydrology,'hydrologydc'):
    103103                                if md.hydrology.isefficientlayer==1:
    104                                         md = checkfield(md,'fieldname','initialization.epl_head','NaN',1,'size',[md.mesh.numberofvertices,1])
    105                                         md = checkfield(md,'fieldname','initialization.epl_thickness','NaN',1,'size',[md.mesh.numberofvertices,1])
     104                                        md = checkfield(md,'fieldname','initialization.epl_head','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
     105                                        md = checkfield(md,'fieldname','initialization.epl_thickness','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices,1])
    106106
    107107                return md

  • issm/trunk-jpl/src/m/classes/inversion.m

    r19048 r19897  
    159159
    160160                        if solution==BalancethicknessSolutionEnum()
    161                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     161                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    162162                        elseif solution==BalancethicknessSoftSolutionEnum()
    163                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     163                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    164164                        else
    165                                 md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
    166                                 md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     165                                md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
     166                                md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    167167                        end
    168168                end % }}}

  • issm/trunk-jpl/src/m/classes/inversion.py

    r19048 r19897  
    151151
    152152                if solution==BalancethicknessSolutionEnum():
    153                         md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices],'NaN',1)
     153                        md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    154154                else:
    155                         md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices],'NaN',1)
    156                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices],'NaN',1)
     155                        md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
     156                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    157157
    158158                return md

  • issm/trunk-jpl/src/m/classes/inversionvalidation.m

    r19048 r19897  
    6767                        md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0 1]);
    6868                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',supportedcontrols());
    69                         md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
     69                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1,'Inf',1);
    7070                        md = checkfield(md,'fieldname','inversion.cost_functions','size',[1 num_costfunc],'values',supportedcostfunctions());
    7171                        md = checkfield(md,'fieldname','inversion.cost_functions_coefficients','size',[md.mesh.numberofvertices num_costfunc],'>=',0);
    … …  
    7474
    7575                        if solution==BalancethicknessSolutionEnum()
    76                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     76                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    7777                        elseif solution==BalancethicknessSoftSolutionEnum()
    78                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     78                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    7979                        else
    80                                 md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     80                                md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    8181                                if ~strcmp(domaintype(md.mesh),'2Dvertical'),
    82                                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     82                                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    8383                                end
    8484                        end

  • issm/trunk-jpl/src/m/classes/linearbasalforcings.m

    r19042 r19897  
    4949
    5050                        if ismember(MasstransportAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.ismasstransport==0),
    51                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
     51                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
    5252                                md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0,'numel',1);
    5353                                md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<','basalforcings.upperwater_elevation','numel',1);
    … …  
    5555                        end
    5656                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    57                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices 1]);
     57                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    5858                                md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0,'numel',1);
    5959                                md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<','basalforcings.upperwater_elevation','numel',1);
    … …  
    6161                        end
    6262                        if ismember(ThermalAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.isthermal==0),
    63                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
     63                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
    6464                                md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0,'numel',1);
    6565                                md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<','basalforcings.upperwater_elevation','numel',1);
    6666                                md = checkfield(md,'fieldname','basalforcings.upperwater_elevation','<',0,'numel',1);
    67                                 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0);
     67                                md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0);
    6868                        end
    6969                end % }}}

  • issm/trunk-jpl/src/m/classes/linearbasalforcings.py

    r19027 r19897  
    6969
    7070                if MasstransportAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.ismasstransport):
    71                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
     71                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
    7272                        md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0);
    7373                        md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<',md.basalforcings.upperwater_elevation);
    … …  
    7575
    7676                if BalancethicknessAnalysisEnum() in analyses:
    77                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices])
     77                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
    7878                        md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0);
    7979                        md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<',md.basalforcings.upperwater_elevation);
    … …  
    8181
    8282                if ThermalAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.isthermal):
    83                         md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1)
     83                        md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1)
    8484                        md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0);
    8585                        md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<',md.basalforcings.upperwater_elevation);
    8686                        md = checkfield(md,'fieldname','basalforcings.upperwater_elevation','<',0);
    87                         md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0)
     87                        md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0)
    8888
    8989                return md

  • issm/trunk-jpl/src/m/classes/m1qn3inversion.m

    r19879 r19897  
    8484                        md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0 1]);
    8585                        md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',supportedcontrols());
    86                         md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1);
     86                        md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[1 num_controls],'>',0,'NaN',1,'Inf',1);
    8787                        md = checkfield(md,'fieldname','inversion.maxsteps','numel',1,'>=',0);
    8888                        md = checkfield(md,'fieldname','inversion.maxiter','numel',1,'>=',0);
    … …  
    9595
    9696                        if solution==BalancethicknessSolutionEnum()
    97                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
    98                                 md = checkfield(md,'fieldname','inversion.surface_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     97                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
     98                                md = checkfield(md,'fieldname','inversion.surface_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    9999                        elseif solution==BalancethicknessSoftSolutionEnum()
    100                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     100                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    101101                        else
    102                                 md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     102                                md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    103103                                if ~strcmp(domaintype(md.mesh),'2Dvertical'),
    104                                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     104                                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    105105                                end
    106106                        end

  • issm/trunk-jpl/src/m/classes/m1qn3inversion.py

    r19048 r19897  
    145145                md = checkfield(md,'fieldname','inversion.incomplete_adjoint','values',[0,1])
    146146                md = checkfield(md,'fieldname','inversion.control_parameters','cell',1,'values',supportedcontrols())
    147                 md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[num_controls],'>',0,'NaN',1)
     147                md = checkfield(md,'fieldname','inversion.control_scaling_factors','size',[num_controls],'>',0,'NaN',1,'Inf',1)
    148148                md = checkfield(md,'fieldname','inversion.maxsteps','numel',[1],'>=',0)
    149149                md = checkfield(md,'fieldname','inversion.maxiter','numel',[1],'>=',0)
    … …  
    156156
    157157                if solution==BalancethicknessSolutionEnum():
    158                         md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices],'NaN',1)
     158                        md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    159159                else:
    160                         md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices],'NaN',1)
    161                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices],'NaN',1)
     160                        md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
     161                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1)
    162162
    163163                return md

  • issm/trunk-jpl/src/m/classes/masscon.m

    r19071 r19897  
    4343                        md = checkfield(md,'fieldname','self.definitionenum','field',self.definitionenum,'values',[Outputdefinition1Enum:Outputdefinition100Enum]);
    4444
    45                         md = checkfield(md,'fieldname','self.levelset','field',self.levelset,'timeseries',1,'NaN',1);
     45                        md = checkfield(md,'fieldname','self.levelset','field',self.levelset,'timeseries',1,'NaN',1,'Inf',1);
    4646
    4747                end % }}}

  • issm/trunk-jpl/src/m/classes/massconaxpby.m

    r19071 r19897  
    5353                        md = checkfield(md,'fieldname','self.definitionenum','field',self.definitionenum,'values',[Outputdefinition1Enum:Outputdefinition100Enum]);
    5454
    55                         md = checkfield(md,'fieldname','self.alpha','field',self.alpha,'NaN',1,'size',[1 1]);
    56                         md = checkfield(md,'fieldname','self.betaa','field',self.beta,'NaN',1,'size',[1 1]);
     55                        md = checkfield(md,'fieldname','self.alpha','field',self.alpha,'NaN',1,'Inf',1,'size',[1 1]);
     56                        md = checkfield(md,'fieldname','self.betaa','field',self.beta,'NaN',1,'Inf',1,'size',[1 1]);
    5757
    5858                end % }}}

  • issm/trunk-jpl/src/m/classes/matice.m

    r19879 r19897  
    147147                        md = checkfield(md,'fieldname','materials.rho_freshwater','>',0);
    148148                        md = checkfield(md,'fieldname','materials.mu_water','>',0);
    149                         md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1);
     149                        md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1,'Inf',1);
    150150                        md = checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements 1]);
    151151                        md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'Cuffey' 'Paterson' 'Arrhenius' 'LliboutryDuval'});

  • issm/trunk-jpl/src/m/classes/matice.py

    r19048 r19897  
    123123                md = checkfield(md,'fieldname','materials.rho_freshwater','>',0)
    124124                md = checkfield(md,'fieldname','materials.mu_water','>',0)
    125                 md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1)
     125                md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1,'Inf',1)
    126126                md = checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements])
    127127                md = checkfield(md,'fieldname','materials.rheology_law','values',['None','Cuffey','Paterson','Arrhenius','LliboutryDuval'])

  • issm/trunk-jpl/src/m/classes/mesh2d.m

    r19879 r19897  
    8181                function md = checkconsistency(self,md,solution,analyses) % {{{
    8282
    83                         md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices 1]);
    84                         md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices 1]);
    85                         md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',1:md.mesh.numberofvertices);
     83                        md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     84                        md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     85                        md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.mesh.numberofvertices);
    8686                        md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements 3]);
    8787                        if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))));
    … …  
    9191                        md = checkfield(md,'fieldname','mesh.numberofvertices','>',0);
    9292                        md = checkfield(md,'fieldname','mesh.average_vertex_connectivity','>=',9,'message','''mesh.average_vertex_connectivity'' should be at least 9 in 2d');
    93                         md = checkfield(md,'fieldname','mesh.segments','NaN',1,'>',0,'size',[NaN 3]);
     93                        md = checkfield(md,'fieldname','mesh.segments','NaN',1,'Inf',1,'>',0,'size',[NaN 3]);
    9494
    9595                        switch(solution),

  • issm/trunk-jpl/src/m/classes/mesh2d.py

    r18558 r19897  
    8484        def checkconsistency(self,md,solution,analyses):    # {{{
    8585
    86                 md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices])
    87                 md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices])
    88                 md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',numpy.arange(1,md.mesh.numberofvertices+1))
     86                md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     87                md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     88                md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',numpy.arange(1,md.mesh.numberofvertices+1))
    8989                md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements,3])
    9090                if numpy.any(numpy.logical_not(m.ismember(numpy.arange(1,md.mesh.numberofvertices+1),md.mesh.elements))):

  • issm/trunk-jpl/src/m/classes/mesh2dvertical.m

    r19040 r19897  
    113113                function md = checkconsistency(self,md,solution,analyses) % {{{
    114114
    115                         md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices 1]);
    116                         md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices 1]);
    117                         md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',1:md.mesh.numberofvertices);
     115                        md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     116                        md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     117                        md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.mesh.numberofvertices);
    118118                        md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements 3]);
    119119                        if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))));

  • issm/trunk-jpl/src/m/classes/mesh3dprisms.m

    r19040 r19897  
    144144                function md = checkconsistency(self,md,solution,analyses) % {{{
    145145
    146                         md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices 1]);
    147                         md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices 1]);
    148                         md = checkfield(md,'fieldname','mesh.z','NaN',1,'size',[md.mesh.numberofvertices 1]);
    149                         md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',1:md.mesh.numberofvertices);
     146                        md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     147                        md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     148                        md = checkfield(md,'fieldname','mesh.z','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     149                        md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.mesh.numberofvertices);
    150150                        md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements 6]);
    151151                        if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))));

  • issm/trunk-jpl/src/m/classes/mesh3dprisms.py

    r18558 r19897  
    106106        def checkconsistency(self,md,solution,analyses):    # {{{
    107107
    108                 md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices])
    109                 md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices])
    110                 md = checkfield(md,'fieldname','mesh.z','NaN',1,'size',[md.mesh.numberofvertices])
    111                 md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',numpy.arange(1,md.mesh.numberofvertices+1))
     108                md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     109                md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     110                md = checkfield(md,'fieldname','mesh.z','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices])
     111                md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',numpy.arange(1,md.mesh.numberofvertices+1))
    112112                md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements,6])
    113113                if numpy.any(numpy.logical_not(m.ismember(numpy.arange(1,md.mesh.numberofvertices+1),md.mesh.elements))):

  • issm/trunk-jpl/src/m/classes/mesh3dsurface.m

    r19076 r19897  
    8181                function md = checkconsistency(obj,md,solution,analyses) % {{{
    8282
    83                         md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices 1]);
    84                         md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices 1]);
    85                         md = checkfield(md,'fieldname','mesh.z','NaN',1,'size',[md.mesh.numberofvertices 1]);
    86                         md = checkfield(md,'fieldname','mesh.lat','NaN',1,'size',[md.mesh.numberofvertices 1]);
    87                         md = checkfield(md,'fieldname','mesh.long','NaN',1,'size',[md.mesh.numberofvertices 1]);
    88                         md = checkfield(md,'fieldname','mesh.r','NaN',1,'size',[md.mesh.numberofvertices 1]);
    89                         md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',1:md.mesh.numberofvertices);
     83                        md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     84                        md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     85                        md = checkfield(md,'fieldname','mesh.z','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     86                        md = checkfield(md,'fieldname','mesh.lat','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     87                        md = checkfield(md,'fieldname','mesh.long','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     88                        md = checkfield(md,'fieldname','mesh.r','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     89                        md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.mesh.numberofvertices);
    9090                        md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements 3]);
    9191                        if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))));

  • issm/trunk-jpl/src/m/classes/mesh3dtetras.m

    r19040 r19897  
    144144                function md = checkconsistency(self,md,solution,analyses) % {{{
    145145
    146                         md = checkfield(md,'fieldname','mesh.x','NaN',1,'size',[md.mesh.numberofvertices 1]);
    147                         md = checkfield(md,'fieldname','mesh.y','NaN',1,'size',[md.mesh.numberofvertices 1]);
    148                         md = checkfield(md,'fieldname','mesh.z','NaN',1,'size',[md.mesh.numberofvertices 1]);
    149                         md = checkfield(md,'fieldname','mesh.elements','NaN',1,'>',0,'values',1:md.mesh.numberofvertices);
     146                        md = checkfield(md,'fieldname','mesh.x','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     147                        md = checkfield(md,'fieldname','mesh.y','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     148                        md = checkfield(md,'fieldname','mesh.z','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
     149                        md = checkfield(md,'fieldname','mesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.mesh.numberofvertices);
    150150                        md = checkfield(md,'fieldname','mesh.elements','size',[md.mesh.numberofelements 4]);
    151151                        if any(~ismember(1:md.mesh.numberofvertices,sort(unique(md.mesh.elements(:)))));

  • issm/trunk-jpl/src/m/classes/misfit.m

    r19067 r19897  
    6464                                error('misfit error message: ''timeinterpolation'' field should be a string!');
    6565                        end
    66                         md = checkfield(md,'fieldname','self.observation','field',self.observation,'timeseries',1,'NaN',1);
     66                        md = checkfield(md,'fieldname','self.observation','field',self.observation,'timeseries',1,'NaN',1,'Inf',1);
    6767                        md = checkfield(md,'fieldname','self.timeinterpolation','field',self.timeinterpolation,'values',{'nearestneighbor'});
    68                         md = checkfield(md,'fieldname','self.weights','field',self.weights,'timeseries',1,'NaN',1);
     68                        md = checkfield(md,'fieldname','self.weights','field',self.weights,'timeseries',1,'NaN',1,'Inf',1);
    6969
    7070                end % }}}

  • issm/trunk-jpl/src/m/classes/mismipbasalforcings.m

    r19713 r19897  
    4949
    5050                        if ismember(MasstransportAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.ismasstransport==0),
    51                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
     51                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
    5252                                md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',1);
    5353                                md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',1);
    … …  
    5555                        end
    5656                        if ismember(BalancethicknessAnalysisEnum(),analyses),
    57                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'size',[md.mesh.numberofvertices 1]);
     57                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
    5858                                md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',1);
    5959                                md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',1);
    … …  
    6161                        end
    6262                        if ismember(ThermalAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.isthermal==0),
    63                                 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1);
     63                                md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'Inf',1,'timeseries',1);
    6464                                md = checkfield(md,'fieldname','basalforcings.meltrate_factor','>=',0,'numel',1);
    6565                                md = checkfield(md,'fieldname','basalforcings.threshold_thickness','>=',0,'numel',1);
    6666                                md = checkfield(md,'fieldname','basalforcings.upperdepth_melt','<=',0,'numel',1);
    67                                 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0);
     67                                md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'Inf',1,'timeseries',1,'>=',0);
    6868                        end
    6969                        if isnan(md.geometry.bed),

  • issm/trunk-jpl/src/m/classes/mismipbasalforcings.py

    r19783 r19897  
    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    # }}}

  • issm/trunk-jpl/src/m/classes/spheremesh.m

    r19040 r19897  
    4141                function md = checkconsistency(self,md,solution,analyses) % {{{
    4242
    43                         md = checkfield(md,'fieldname','spheremesh.x','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    44                         md = checkfield(md,'fieldname','spheremesh.y','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    45                         md = checkfield(md,'fieldname','spheremesh.z','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    46                         md = checkfield(md,'fieldname','spheremesh.r','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    47                         md = checkfield(md,'fieldname','spheremesh.theta','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    48                         md = checkfield(md,'fieldname','spheremesh.phi','NaN',1,'size',[md.spheremesh.numberofvertices 1]);
    49                         md = checkfield(md,'fieldname','spheremesh.elements','NaN',1,'>',0,'values',1:md.spheremesh.numberofvertices);
     43                        md = checkfield(md,'fieldname','spheremesh.x','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     44                        md = checkfield(md,'fieldname','spheremesh.y','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     45                        md = checkfield(md,'fieldname','spheremesh.z','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     46                        md = checkfield(md,'fieldname','spheremesh.r','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     47                        md = checkfield(md,'fieldname','spheremesh.theta','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     48                        md = checkfield(md,'fieldname','spheremesh.phi','NaN',1,'Inf',1,'size',[md.spheremesh.numberofvertices 1]);
     49                        md = checkfield(md,'fieldname','spheremesh.elements','NaN',1,'Inf',1,'>',0,'values',1:md.spheremesh.numberofvertices);
    5050                        md = checkfield(md,'fieldname','spheremesh.elements','size',[md.spheremesh.numberofelements 3]);
    5151                        if any(~ismember(1:md.spheremesh.numberofvertices,sort(unique(md.spheremesh.elements(:)))));
    … …  
    5555                        md = checkfield(md,'fieldname','spheremesh.numberofelements','>',0);
    5656                        md = checkfield(md,'fieldname','spheremesh.numberofvertices','>',0);
    57                         md = checkfield(md,'fieldname','spheremesh.elementconnectivity','size',[md.spheremesh.numberofelements 3],'NaN',1);
     57                        md = checkfield(md,'fieldname','spheremesh.elementconnectivity','size',[md.spheremesh.numberofelements 3],'NaN',1,'Inf',1);
    5858                end % }}}
    5959                function disp(self) % {{{

  • issm/trunk-jpl/src/m/classes/stressbalance.m

    r19892 r19897  
    137137                        md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1);
    138138                        md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1);
    139                         md = checkfield(md,'fieldname','stressbalance.restol','size',[1 1],'>',0,'NaN',1);
     139                        md = checkfield(md,'fieldname','stressbalance.restol','size',[1 1],'>',0,'NaN',1,'Inf',1);
    140140                        md = checkfield(md,'fieldname','stressbalance.reltol','size',[1 1]);
    141141                        md = checkfield(md,'fieldname','stressbalance.abstol','size',[1 1]);
    142142                        md = checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0 1 2]);
    143                         md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1 1],'NaN',1);
    144                         md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1 1],'NaN',1);
     143                        md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1 1],'NaN',1,'Inf',1);
     144                        md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1 1],'NaN',1,'Inf',1);
    145145                        md = checkfield(md,'fieldname','stressbalance.maxiter','size',[1 1],'>=',1);
    146146                        md = checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices 6]);

  • issm/trunk-jpl/src/m/classes/stressbalance.py

    r19294 r19897  
    141141                md = checkfield(md,'fieldname','stressbalance.abstol','size',[1])
    142142                md = checkfield(md,'fieldname','stressbalance.isnewton','numel',[1],'values',[0,1,2])
    143                 md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1],'NaN',1)
    144                 md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1],'NaN',1)
     143                md = checkfield(md,'fieldname','stressbalance.FSreconditioning','size',[1],'NaN',1,'Inf',1)
     144                md = checkfield(md,'fieldname','stressbalance.viscosity_overshoot','size',[1],'NaN',1,'Inf',1)
    145145                md = checkfield(md,'fieldname','stressbalance.maxiter','size',[1],'>=',1)
    146146                md = checkfield(md,'fieldname','stressbalance.referential','size',[md.mesh.numberofvertices,6])

  • issm/trunk-jpl/src/m/classes/taoinversion.m

    r19048 r19897  
    118118
    119119                        if solution==BalancethicknessSolutionEnum()
    120                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     120                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    121121                        elseif solution==BalancethicknessSoftSolutionEnum()
    122                                 md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     122                                md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    123123                        else
    124                                 md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
    125                                 md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1);
     124                                md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
     125                                md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices 1],'NaN',1,'Inf',1);
    126126                        end
    127127                end % }}}

  • issm/trunk-jpl/src/m/classes/taoinversion.py

    r19379 r19897  
    156156
    157157                if solution==BalancethicknessSolutionEnum():
    158                         md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices, 1],'NaN',1)
     158                        md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices, 1],'NaN',1,'Inf',1)
    159159                elif solution==BalancethicknessSoftSolutionEnum():
    160                         md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices, 1],'NaN',1)
     160                        md = checkfield(md,'fieldname','inversion.thickness_obs','size',[md.mesh.numberofvertices, 1],'NaN',1,'Inf',1)
    161161                else:
    162                         md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices, 1],'NaN',1)
    163                         md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices, 1],'NaN',1)
     162                        md = checkfield(md,'fieldname','inversion.vx_obs','size',[md.mesh.numberofvertices, 1],'NaN',1,'Inf',1)
     163                        md = checkfield(md,'fieldname','inversion.vy_obs','size',[md.mesh.numberofvertices, 1],'NaN',1,'Inf',1)
    164164
    165165                def marshall(self, md, fid):

  • issm/trunk-jpl/src/m/classes/timestepping.m

    r19894 r19897  
    5353                function md = checkconsistency(self,md,solution,analyses) % {{{
    5454
    55                         md = checkfield(md,'fieldname','timestepping.start_time','numel',[1],'NaN',1);
    56                         md = checkfield(md,'fieldname','timestepping.final_time','numel',[1],'NaN',1);
    57                         md = checkfield(md,'fieldname','timestepping.time_step','numel',[1],'>=',0,'NaN',1);
     55                        md = checkfield(md,'fieldname','timestepping.start_time','numel',[1],'NaN',1,'Inf',1);
     56                        md = checkfield(md,'fieldname','timestepping.final_time','numel',[1],'NaN',1,'Inf',1);
     57                        md = checkfield(md,'fieldname','timestepping.time_step','numel',[1],'>=',0,'NaN',1,'Inf',1);
    5858                        md = checkfield(md,'fieldname','timestepping.time_adapt','numel',[1],'values',[0 1]);
    5959                        md = checkfield(md,'fieldname','timestepping.cfl_coefficient','numel',[1],'>',0,'<=',1);

  • issm/trunk-jpl/src/m/classes/timestepping.py

    r19707 r19897  
    5353        def checkconsistency(self,md,solution,analyses):    # {{{
    5454
    55                 md = checkfield(md,'fieldname','timestepping.start_time','numel',[1],'NaN',1)
    56                 md = checkfield(md,'fieldname','timestepping.final_time','numel',[1],'NaN',1)
    57                 md = checkfield(md,'fieldname','timestepping.time_step','numel',[1],'>=',0,'NaN',1)
     55                md = checkfield(md,'fieldname','timestepping.start_time','numel',[1],'NaN',1,'Inf',1)
     56                md = checkfield(md,'fieldname','timestepping.final_time','numel',[1],'NaN',1,'Inf',1)
     57                md = checkfield(md,'fieldname','timestepping.time_step','numel',[1],'>=',0,'NaN',1,'Inf',1)
    5858                md = checkfield(md,'fieldname','timestepping.time_adapt','numel',[1],'values',[0,1])
    5959                md = checkfield(md,'fieldname','timestepping.cfl_coefficient','numel',[1],'>',0,'<=',1)

  • issm/trunk-jpl/src/m/consistency/checkfield.m

    r19560 r19897  
    88%   Available options:
    99%      - NaN: 1 if check that there is no NaN
     10%      - Inf: 1 if check that there is no Inf
    1011%      - size: [lines cols], NaN for non checked dimensions
    1112%      - >:  greater than provided value
    … …  
    9192        end
    9293end
     94
     95%check Inf
     96if getfieldvalue(options,'Inf',0);
     97        field2=reshape(field,prod(size(field)),1);
     98        if any(isinf(field2)),
     99                md = checkmessage(md,getfieldvalue(options,'message',...
     100                        ['Inf values found in field ''' fieldname '''']));
     101        end
     102end
     103
    93104
    94105%check cell
Note: See TracChangeset for help on using the changeset viewer.