Changeset 19027
- Timestamp:
- 01/20/15 15:09:59 (10 years ago)
- Location:
- issm/trunk-jpl/src
- Files:
-
- 52 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified issm/trunk-jpl/src/c/classes/Inputs/TransientInput.h ¶
r18450 r19027 42 42 Input* SpawnTriaInput(int index1,int index2,int index3); 43 43 Input* SpawnSegInput(int index1,int index2); 44 Input* PointwiseDivide(Input* forcingB){_error_("not implemented yet");};45 Input* PointwiseMin(Input* forcingB){_error_("not implemented yet");};46 Input* PointwiseMax(Input* forcingB){_error_("not implemented yet");};44 Input* PointwiseDivide(Input* input_in){_error_("not implemented yet");}; 45 Input* PointwiseMin(Input* input_in){_error_("not implemented yet");}; 46 Input* PointwiseMax(Input* input_in){_error_("not implemented yet");}; 47 47 int GetResultInterpolation(void); 48 48 int GetResultNumberOfNodes(void); -
TabularUnified issm/trunk-jpl/src/m/classes/SMB.m ¶
r17720 r19027 38 38 39 39 if ismember(MasstransportAnalysisEnum(),analyses), 40 md = checkfield(md,'fieldname','surfaceforcings.mass_balance',' forcing',1,'NaN',1);40 md = checkfield(md,'fieldname','surfaceforcings.mass_balance','timeseries',1,'NaN',1); 41 41 end 42 42 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 53 53 54 54 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBEnum(),'format','Integer'); 55 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);55 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 56 56 end % }}} 57 57 end -
TabularUnified issm/trunk-jpl/src/m/classes/SMB.py ¶
r17497 r19027 38 38 39 39 if MasstransportAnalysisEnum() in analyses: 40 md = checkfield(md,'fieldname','surfaceforcings.mass_balance',' forcing',1,'NaN',1)40 md = checkfield(md,'fieldname','surfaceforcings.mass_balance','timeseries',1,'NaN',1) 41 41 42 42 if BalancethicknessAnalysisEnum() in analyses: … … 50 50 51 51 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBEnum(),'format','Integer'); 52 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)52 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 53 53 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/SMBcomponents.m ¶
r18000 r19027 44 44 45 45 if ismember(MasstransportAnalysisEnum(),analyses), 46 md = checkfield(md,'fieldname','surfaceforcings.accumulation',' forcing',1,'NaN',1);46 md = checkfield(md,'fieldname','surfaceforcings.accumulation','timeseries',1,'NaN',1); 47 47 end 48 48 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 50 50 end 51 51 if ismember(MasstransportAnalysisEnum(),analyses), 52 md = checkfield(md,'fieldname','surfaceforcings.runoff',' forcing',1,'NaN',1);52 md = checkfield(md,'fieldname','surfaceforcings.runoff','timeseries',1,'NaN',1); 53 53 end 54 54 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 56 56 end 57 57 if ismember(MasstransportAnalysisEnum(),analyses), 58 md = checkfield(md,'fieldname','surfaceforcings.evaporation',' forcing',1,'NaN',1);58 md = checkfield(md,'fieldname','surfaceforcings.evaporation','timeseries',1,'NaN',1); 59 59 end 60 60 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 73 73 74 74 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBcomponentsEnum(),'format','Integer'); 75 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);76 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','runoff','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);77 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);75 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 76 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','runoff','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 77 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 78 78 end % }}} 79 79 end -
TabularUnified issm/trunk-jpl/src/m/classes/SMBcomponents.py ¶
r18007 r19027 51 51 52 52 if MasstransportAnalysisEnum() in analyses: 53 md = checkfield(md,'fieldname','surfaceforcings.accumulation',' forcing',1,'NaN',1)53 md = checkfield(md,'fieldname','surfaceforcings.accumulation','timeseries',1,'NaN',1) 54 54 55 55 if BalancethicknessAnalysisEnum() in analyses: … … 57 57 58 58 if MasstransportAnalysisEnum() in analyses: 59 md = checkfield(md,'fieldname','surfaceforcings.runoff',' forcing',1,'NaN',1)59 md = checkfield(md,'fieldname','surfaceforcings.runoff','timeseries',1,'NaN',1) 60 60 61 61 if BalancethicknessAnalysisEnum() in analyses: … … 63 63 64 64 if MasstransportAnalysisEnum() in analyses: 65 md = checkfield(md,'fieldname','surfaceforcings.evaporation',' forcing',1,'NaN',1)65 md = checkfield(md,'fieldname','surfaceforcings.evaporation','timeseries',1,'NaN',1) 66 66 67 67 if BalancethicknessAnalysisEnum() in analyses: … … 75 75 76 76 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBcomponentsEnum(),'format','Integer'); 77 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)78 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','runoff','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)79 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)77 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 78 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','runoff','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 79 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 80 80 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/SMBgradients.m ¶
r18265 r19027 38 38 39 39 if ismember(MasstransportAnalysisEnum(),analyses), 40 md = checkfield(md,'fieldname','surfaceforcings.href',' forcing',1,'NaN',1);41 md = checkfield(md,'fieldname','surfaceforcings.smbref',' forcing',1,'NaN',1);42 md = checkfield(md,'fieldname','surfaceforcings.b_pos',' forcing',1,'NaN',1);43 md = checkfield(md,'fieldname','surfaceforcings.b_neg',' forcing',1,'NaN',1);40 md = checkfield(md,'fieldname','surfaceforcings.href','timeseries',1,'NaN',1); 41 md = checkfield(md,'fieldname','surfaceforcings.smbref','timeseries',1,'NaN',1); 42 md = checkfield(md,'fieldname','surfaceforcings.b_pos','timeseries',1,'NaN',1); 43 md = checkfield(md,'fieldname','surfaceforcings.b_neg','timeseries',1,'NaN',1); 44 44 end 45 45 end % }}} … … 59 59 60 60 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBgradientsEnum(),'format','Integer'); 61 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','href','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);62 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);63 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);64 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);61 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','href','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 62 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 63 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 64 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 65 65 66 66 end % }}} -
TabularUnified issm/trunk-jpl/src/m/classes/SMBgradients.py ¶
r17497 r19027 44 44 45 45 if MasstransportAnalysisEnum() in analyses: 46 md = checkfield(md,'fieldname','surfaceforcings.href',' forcing',1,'NaN',1)47 md = checkfield(md,'fieldname','surfaceforcings.smbref',' forcing',1,'NaN',1)48 md = checkfield(md,'fieldname','surfaceforcings.b_pos',' forcing',1,'NaN',1)49 md = checkfield(md,'fieldname','surfaceforcings.b_neg',' forcing',1,'NaN',1)46 md = checkfield(md,'fieldname','surfaceforcings.href','timeseries',1,'NaN',1) 47 md = checkfield(md,'fieldname','surfaceforcings.smbref','timeseries',1,'NaN',1) 48 md = checkfield(md,'fieldname','surfaceforcings.b_pos','timeseries',1,'NaN',1) 49 md = checkfield(md,'fieldname','surfaceforcings.b_neg','timeseries',1,'NaN',1) 50 50 51 51 return md … … 56 56 57 57 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBgradientsEnum(),'format','Integer'); 58 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','href','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)59 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)60 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)61 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)58 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','href','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 59 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 60 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 61 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 62 62 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/SMBhenning.m ¶
r17087 r19027 32 32 33 33 if ismember(MasstransportAnalysisEnum(),analyses), 34 md = checkfield(md,'fieldname','surfaceforcings.smbref',' forcing',1,'NaN',1);34 md = checkfield(md,'fieldname','surfaceforcings.smbref','timeseries',1,'NaN',1); 35 35 end 36 36 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 47 47 48 48 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBhenningEnum(),'format','Integer'); 49 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);49 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 50 50 end % }}} 51 51 end -
TabularUnified issm/trunk-jpl/src/m/classes/SMBmeltcomponents.m ¶
r18000 r19027 50 50 51 51 if ismember(MasstransportAnalysisEnum(),analyses), 52 md = checkfield(md,'fieldname','surfaceforcings.accumulation',' forcing',1,'NaN',1);52 md = checkfield(md,'fieldname','surfaceforcings.accumulation','timeseries',1,'NaN',1); 53 53 end 54 54 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 56 56 end 57 57 if ismember(MasstransportAnalysisEnum(),analyses), 58 md = checkfield(md,'fieldname','surfaceforcings.evaporation',' forcing',1,'NaN',1);58 md = checkfield(md,'fieldname','surfaceforcings.evaporation','timeseries',1,'NaN',1); 59 59 end 60 60 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 62 62 end 63 63 if ismember(MasstransportAnalysisEnum(),analyses), 64 md = checkfield(md,'fieldname','surfaceforcings.refreeze',' forcing',1,'NaN',1);64 md = checkfield(md,'fieldname','surfaceforcings.refreeze','timeseries',1,'NaN',1); 65 65 end 66 66 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 68 68 end 69 69 if ismember(MasstransportAnalysisEnum(),analyses), 70 md = checkfield(md,'fieldname','surfaceforcings.melt',' forcing',1,'NaN',1);70 md = checkfield(md,'fieldname','surfaceforcings.melt','timeseries',1,'NaN',1); 71 71 end 72 72 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 86 86 87 87 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBmeltcomponentsEnum(),'format','Integer'); 88 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);89 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);90 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','melt','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);91 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','refreeze','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);88 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 89 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 90 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','melt','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 91 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','refreeze','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 92 92 end % }}} 93 93 end -
TabularUnified issm/trunk-jpl/src/m/classes/SMBmeltcomponents.py ¶
r18007 r19027 57 57 58 58 if MasstransportAnalysisEnum() in analyses: 59 md = checkfield(md,'fieldname','surfaceforcings.accumulation',' forcing',1,'NaN',1)59 md = checkfield(md,'fieldname','surfaceforcings.accumulation','timeseries',1,'NaN',1) 60 60 61 61 if BalancethicknessAnalysisEnum() in analyses: … … 63 63 64 64 if MasstransportAnalysisEnum() in analyses: 65 md = checkfield(md,'fieldname','surfaceforcings.melt',' forcing',1,'NaN',1)65 md = checkfield(md,'fieldname','surfaceforcings.melt','timeseries',1,'NaN',1) 66 66 67 67 if BalancethicknessAnalysisEnum() in analyses: … … 69 69 70 70 if MasstransportAnalysisEnum() in analyses: 71 md = checkfield(md,'fieldname','surfaceforcings.refreeze',' forcing',1,'NaN',1)71 md = checkfield(md,'fieldname','surfaceforcings.refreeze','timeseries',1,'NaN',1) 72 72 73 73 if BalancethicknessAnalysisEnum() in analyses: … … 75 75 76 76 if MasstransportAnalysisEnum() in analyses: 77 md = checkfield(md,'fieldname','surfaceforcings.evaporation',' forcing',1,'NaN',1)77 md = checkfield(md,'fieldname','surfaceforcings.evaporation','timeseries',1,'NaN',1) 78 78 79 79 if BalancethicknessAnalysisEnum() in analyses: … … 87 87 88 88 WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBmeltcomponentsEnum(),'format','Integer'); 89 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)90 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)91 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','melt','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)92 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','refreeze','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)89 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','accumulation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 90 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','evaporation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 91 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','melt','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 92 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','refreeze','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 93 93 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/SMBpdd.m ¶
r18984 r19027 75 75 md = checkfield(md,'fieldname','surfaceforcings.rlapslgm','>=',0,'numel',1); 76 76 if(obj.isdelta18o==0 & obj.ismungsm==0) 77 md = checkfield(md,'fieldname','surfaceforcings.monthlytemperatures',' forcing',1,'NaN',1);78 md = checkfield(md,'fieldname','surfaceforcings.precipitation',' forcing',1,'NaN',1);77 md = checkfield(md,'fieldname','surfaceforcings.monthlytemperatures','timeseries',1,'NaN',1); 78 md = checkfield(md,'fieldname','surfaceforcings.precipitation','timeseries',1,'NaN',1); 79 79 elseif(obj.isdelta18o==1) 80 80 md = checkfield(md,'fieldname','surfaceforcings.delta18o','NaN',1); … … 146 146 if(obj.isdelta18o==0 & obj.ismungsm==0) 147 147 %WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','monthlytemperatures','format','DoubleMat','mattype',1); 148 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','monthlytemperatures','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);149 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','precipitation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);148 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','monthlytemperatures','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 149 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','precipitation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 150 150 elseif obj.isdelta18o 151 151 WriteData(fid,'object',obj,'class','surfaceforcings','fieldname','temperatures_presentday','format','DoubleMat','mattype',1); -
TabularUnified issm/trunk-jpl/src/m/classes/SMBpdd.py ¶
r18984 r19027 153 153 if (self.isdelta18o==0 and self.ismungsm==0): 154 154 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','monthlytemperatures','format','DoubleMat','mattype',1) 155 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','precipitation','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)155 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','precipitation','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 156 156 elif self.isdelta18o: 157 157 WriteData(fid,'object',self,'class','surfaceforcings','fieldname','temperatures_presentday','format','DoubleMat','mattype',1) -
TabularUnified issm/trunk-jpl/src/m/classes/basalforcings.m ¶
r18093 r19027 45 45 46 46 if ismember(MasstransportAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.ismasstransport==0), 47 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1);48 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,' forcing',1);47 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1); 48 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1); 49 49 end 50 50 if ismember(BalancethicknessAnalysisEnum(),analyses), … … 53 53 end 54 54 if ismember(ThermalAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.isthermal==0), 55 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1);56 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,' forcing',1);57 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,' forcing',1,'>=',0);55 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1); 56 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1); 57 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0); 58 58 end 59 59 end % }}} … … 71 71 72 72 WriteData(fid,'enum',BasalforcingsEnum(),'data',FloatingMeltRateEnum(),'format','Integer'); 73 WriteData(fid,'object',obj,'fieldname','groundedice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)74 WriteData(fid,'object',obj,'fieldname','floatingice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)75 WriteData(fid,'object',obj,'fieldname','geothermalflux','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);73 WriteData(fid,'object',obj,'fieldname','groundedice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 74 WriteData(fid,'object',obj,'fieldname','floatingice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 75 WriteData(fid,'object',obj,'fieldname','geothermalflux','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 76 76 end % }}} 77 77 end -
TabularUnified issm/trunk-jpl/src/m/classes/basalforcings.py ¶
r18093 r19027 48 48 49 49 if MasstransportAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.ismasstransport): 50 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1)51 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,' forcing',1)50 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1) 51 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1) 52 52 53 53 if BalancethicknessAnalysisEnum() in analyses: … … 56 56 57 57 if ThermalAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.isthermal): 58 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1)59 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,' forcing',1)60 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,' forcing',1,'>=',0)58 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1) 59 md = checkfield(md,'fieldname','basalforcings.floatingice_melting_rate','NaN',1,'timeseries',1) 60 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0) 61 61 62 62 return md … … 67 67 68 68 WriteData(fid,'enum',BasalforcingsEnum(),'data',FloatingMeltRateEnum(),'format','Integer'); 69 WriteData(fid,'object',self,'fieldname','groundedice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)70 WriteData(fid,'object',self,'fieldname','floatingice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)71 WriteData(fid,'object',self,'fieldname','geothermalflux','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)69 WriteData(fid,'object',self,'fieldname','groundedice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 70 WriteData(fid,'object',self,'fieldname','floatingice_melting_rate','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 71 WriteData(fid,'object',self,'fieldname','geothermalflux','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 72 72 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/calving.m ¶
r19024 r19027 40 40 41 41 md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]); 42 md = checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,' forcing',1,'NaN',1);43 md = checkfield(md,'fieldname','calving.meltingrate(1:md.mesh.numberofvertices,:)','>=',0,' forcing',1,'NaN',1);42 md = checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1); 43 md = checkfield(md,'fieldname','calving.meltingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1); 44 44 end % }}} 45 45 function disp(obj) % {{{ … … 53 53 WriteData(fid,'enum',CalvingLawEnum(),'data',DefaultCalvingEnum(),'format','Integer'); 54 54 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',obj.stabilization,'format','Integer'); 55 WriteData(fid,'object',obj,'fieldname','calvingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts);56 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts);55 WriteData(fid,'object',obj,'fieldname','calvingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts); 56 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts); 57 57 end % }}} 58 58 end -
TabularUnified issm/trunk-jpl/src/m/classes/calving.py ¶
r19024 r19027 45 45 46 46 md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]); 47 md = checkfield(md,'fieldname','calving.calvingrate','>=',0,' forcing',1,'NaN',1);48 md = checkfield(md,'fieldname','calving.meltingrate','>=',0,' forcing',1,'NaN',1);47 md = checkfield(md,'fieldname','calving.calvingrate','>=',0,'timeseries',1,'NaN',1); 48 md = checkfield(md,'fieldname','calving.meltingrate','>=',0,'timeseries',1,'NaN',1); 49 49 50 50 return md … … 56 56 WriteData(fid,'enum',CalvingLawEnum(),'data',DefaultCalvingEnum(),'format','Integer'); 57 57 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',self.stabilization,'format','Integer'); 58 WriteData(fid,'object',self,'fieldname','calvingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts)59 WriteData(fid,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts)58 WriteData(fid,'object',self,'fieldname','calvingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts) 59 WriteData(fid,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts) 60 60 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/calvinglevermann.m ¶
r19024 r19027 57 57 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',obj.stabilization,'format','Integer'); 58 58 WriteData(fid,'enum',CalvinglevermannCoeffEnum(),'data',obj.coeff,'format','DoubleMat','mattype',1); 59 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts);59 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts); 60 60 end % }}} 61 61 end -
TabularUnified issm/trunk-jpl/src/m/classes/calvinglevermann.py ¶
r19024 r19027 54 54 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',self.stabilization,'format','Integer'); 55 55 WriteData(fid,'enum',CalvinglevermannCoeffEnum(),'data',self.coeff,'format','DoubleMat','mattype',1) 56 WriteData(fid,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts)56 WriteData(fid,'object',self,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts) 57 57 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/calvingpi.m ¶
r19024 r19027 57 57 WriteData(fid,'enum',LevelsetStabilizationEnum(),'data',obj.stabilization,'format','Integer'); 58 58 WriteData(fid,'enum',CalvingpiCoeffEnum(),'data',obj.coeff,'format','DoubleMat','mattype',1); 59 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1,'scale',1./yts);59 WriteData(fid,'object',obj,'fieldname','meltingrate','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1,'scale',1./yts); 60 60 end % }}} 61 61 end -
TabularUnified issm/trunk-jpl/src/m/classes/damage.m ¶
r18279 r19027 150 150 md = checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2]); 151 151 md = checkfield(md,'fieldname','damage.D','>=',0,'<=',obj.max_damage,'size',[md.mesh.numberofvertices 1]); 152 md = checkfield(md,'fieldname','damage.spcdamage',' forcing',1);152 md = checkfield(md,'fieldname','damage.spcdamage','timeseries',1); 153 153 md = checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0); 154 154 md = checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0 1 2]); … … 213 213 WriteData(fid,'object',obj,'fieldname','law','format','Integer'); 214 214 WriteData(fid,'object',obj,'fieldname','D','format','DoubleMat','mattype',1); 215 WriteData(fid,'object',obj,'fieldname','spcdamage','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);215 WriteData(fid,'object',obj,'fieldname','spcdamage','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 216 216 WriteData(fid,'object',obj,'fieldname','max_damage','format','Double'); 217 217 -
TabularUnified issm/trunk-jpl/src/m/classes/damage.py ¶
r18278 r19027 131 131 md = checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0) 132 132 md = checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2,3]) 133 md = checkfield(md,'fieldname','damage.spcdamage',' forcing',1)133 md = checkfield(md,'fieldname','damage.spcdamage','timeseries',1) 134 134 md = checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0,1,2]) 135 135 md = checkfield(md,'fieldname','damage.maxiter','>=0',0) … … 159 159 WriteData(fid,'object',self,'fieldname','D','format','DoubleMat','mattype',1) 160 160 WriteData(fid,'object',self,'fieldname','law','format','Integer') 161 WriteData(fid,'object',self,'fieldname','spcdamage','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)161 WriteData(fid,'object',self,'fieldname','spcdamage','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 162 162 WriteData(fid,'object',self,'fieldname','max_damage','format','Double') 163 163 WriteData(fid,'object',self,'fieldname','stabilization','format','Integer') -
TabularUnified issm/trunk-jpl/src/m/classes/friction.m ¶
r17952 r19027 35 35 %Early return 36 36 if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end 37 md = checkfield(md,'fieldname','friction.coefficient',' forcing',1,'NaN',1);37 md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1); 38 38 md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]); 39 39 md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]); … … 49 49 50 50 WriteData(fid,'enum',FrictionLawEnum,'data',1,'format','Integer'); 51 WriteData(fid,'object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);51 WriteData(fid,'object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 52 52 %WriteData(fid,'object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1); 53 53 WriteData(fid,'object',obj,'fieldname','p','format','DoubleMat','mattype',2); -
TabularUnified issm/trunk-jpl/src/m/classes/friction.py ¶
r17952 r19027 38 38 return md 39 39 40 md = checkfield(md,'fieldname','friction.coefficient',' forcing',1,'NaN',1)40 md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1) 41 41 md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements]) 42 42 md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements]) -
TabularUnified issm/trunk-jpl/src/m/classes/frictionhydro.m ¶
r19001 r19027 43 43 WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',2); 44 44 WriteData(fid,'class','friction','object',obj,'fieldname','As','format','DoubleMat','mattype',2); 45 WriteData(fid,'class','friction','object',obj,'fieldname','effective_pressure','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);45 WriteData(fid,'class','friction','object',obj,'fieldname','effective_pressure','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 46 46 end % }}} 47 47 end -
TabularUnified issm/trunk-jpl/src/m/classes/frictiontemp.m ¶
r18732 r19027 34 34 if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end 35 35 36 md = checkfield(md,'fieldname','friction.coefficient',' forcing',1,'NaN',1);36 md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1); 37 37 md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]); 38 38 md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]); … … 53 53 WriteData(fid,'enum',FrictionLawEnum,'data',4,'format','Integer'); 54 54 WriteData(fid,'class','friction','object',obj,'fieldname','gamma','format','Double'); 55 WriteData(fid,'class','friction','object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);55 WriteData(fid,'class','friction','object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 56 56 WriteData(fid,'class','friction','object',obj,'fieldname','p','format','DoubleMat','mattype',2); 57 57 WriteData(fid,'class','friction','object',obj,'fieldname','q','format','DoubleMat','mattype',2); -
TabularUnified issm/trunk-jpl/src/m/classes/frictionwaterlayer.m ¶
r18772 r19027 30 30 if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end 31 31 32 md = checkfield(md,'fieldname','friction.coefficient',' forcing',1,'NaN',1);32 md = checkfield(md,'fieldname','friction.coefficient','timeseries',1,'NaN',1); 33 33 md = checkfield(md,'fieldname','friction.q','NaN',1,'size',[md.mesh.numberofelements 1]); 34 34 md = checkfield(md,'fieldname','friction.p','NaN',1,'size',[md.mesh.numberofelements 1]); 35 md = checkfield(md,'fieldname','thermal.spctemperature',' forcing',1,'>=',0.);35 md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1,'>=',0.); 36 36 37 37 end % }}} … … 46 46 47 47 WriteData(fid,'enum',FrictionLawEnum,'data',5,'format','Integer'); 48 WriteData(fid,'class','friction','object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);48 WriteData(fid,'class','friction','object',obj,'fieldname','coefficient','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 49 49 WriteData(fid,'class','friction','object',obj,'fieldname','p','format','DoubleMat','mattype',2); 50 50 WriteData(fid,'class','friction','object',obj,'fieldname','q','format','DoubleMat','mattype',2); 51 WriteData(fid,'class','friction','object',obj,'fieldname','water_layer','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);51 WriteData(fid,'class','friction','object',obj,'fieldname','water_layer','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 52 52 end % }}} 53 53 end -
TabularUnified issm/trunk-jpl/src/m/classes/frictionweertman.m ¶
r18800 r19027 34 34 %Early return 35 35 if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end 36 md = checkfield(md,'fieldname','friction.C',' forcing',1,'NaN',1);36 md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1); 37 37 md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements 1]); 38 38 end % }}} … … 51 51 52 52 WriteData(fid,'enum',FrictionLawEnum,'data',2,'format','Integer'); 53 WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);53 WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 54 54 WriteData(fid,'class','friction','object',obj,'fieldname','m','format','DoubleMat','mattype',2); 55 55 -
TabularUnified issm/trunk-jpl/src/m/classes/frictionweertman.py ¶
r18800 r19027 36 36 return md 37 37 38 md = checkfield(md,'fieldname','friction.C',' forcing',1,'NaN',1)38 md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1) 39 39 md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements]) 40 40 -
TabularUnified issm/trunk-jpl/src/m/classes/frictionweertmantemp.m ¶
r18804 r19027 26 26 %Early return 27 27 if ~ismember(StressbalanceAnalysisEnum(),analyses) & ~ismember(ThermalAnalysisEnum(),analyses), return; end 28 md = checkfield(md,'fieldname','friction.C',' forcing',1,'NaN',1);28 md = checkfield(md,'fieldname','friction.C','timeseries',1,'NaN',1); 29 29 md = checkfield(md,'fieldname','friction.m','NaN',1,'size',[md.mesh.numberofelements 1]); 30 30 end % }}} … … 42 42 WriteData(fid,'enum',FrictionLawEnum,'data',6,'format','Integer'); 43 43 WriteData(fid,'class','friction','object',obj,'fieldname','gamma','format','Double'); 44 WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);44 WriteData(fid,'class','friction','object',obj,'fieldname','C','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 45 45 WriteData(fid,'class','friction','object',obj,'fieldname','m','format','DoubleMat','mattype',2); 46 46 -
TabularUnified issm/trunk-jpl/src/m/classes/geometry.m ¶
r18492 r19027 63 63 64 64 if (solution==TransientSolutionEnum() & md.transient.isgia) | (solution==GiaSolutionEnum()), 65 md = checkfield(md,'fieldname','geometry.thickness',' forcing',1,'NaN',1,'>=',0);65 md = checkfield(md,'fieldname','geometry.thickness','timeseries',1,'NaN',1,'>=',0); 66 66 else 67 67 md = checkfield(md,'fieldname','geometry.surface' ,'NaN',1,'size',[md.mesh.numberofvertices 1]); … … 87 87 function marshall(obj,md,fid) % {{{ 88 88 WriteData(fid,'data',obj.surface,'format','DoubleMat','mattype',1,'enum',SurfaceEnum()); 89 WriteData(fid,'data',obj.thickness,'format','DoubleMat','mattype',1,'enum',ThicknessEnum(),' forcinglength',md.mesh.numberofvertices+1);89 WriteData(fid,'data',obj.thickness,'format','DoubleMat','mattype',1,'enum',ThicknessEnum(),'timeserieslength',md.mesh.numberofvertices+1); 90 90 WriteData(fid,'data',obj.base,'format','DoubleMat','mattype',1,'enum',BaseEnum()); 91 91 WriteData(fid,'data',obj.bed,'format','DoubleMat','mattype',1,'enum',BedEnum()); -
TabularUnified issm/trunk-jpl/src/m/classes/geometry.py ¶
r17590 r19027 39 39 md = checkfield(md,'fieldname','geometry.surface' ,'NaN',1,'size',[md.mesh.numberofvertices]) 40 40 md = checkfield(md,'fieldname','geometry.base' ,'NaN',1,'size',[md.mesh.numberofvertices]) 41 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices],'>',0,' forcing',1)41 md = checkfield(md,'fieldname','geometry.thickness','NaN',1,'size',[md.mesh.numberofvertices],'>',0,'timeseries',1) 42 42 if any((self.thickness-self.surface+self.base)>10**-9): 43 43 md.checkmessage("equality thickness=surface-base violated") … … 49 49 def marshall(self,md,fid): # {{{ 50 50 WriteData(fid,'data',self.surface,'format','DoubleMat','mattype',1,'enum',SurfaceEnum()) 51 WriteData(fid,'data',self.thickness,'format','DoubleMat','mattype',1,'enum',ThicknessEnum(),' forcinglength',md.mesh.numberofvertices+1)51 WriteData(fid,'data',self.thickness,'format','DoubleMat','mattype',1,'enum',ThicknessEnum(),'timeserieslength',md.mesh.numberofvertices+1) 52 52 WriteData(fid,'data',self.base,'format','DoubleMat','mattype',1,'enum',BaseEnum()) 53 53 WriteData(fid,'data',self.bed,'format','DoubleMat','mattype',1,'enum',BedEnum()) -
TabularUnified issm/trunk-jpl/src/m/classes/hydrologydc.m ¶
r18982 r19027 105 105 md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0,'numel',1); 106 106 end 107 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,' forcing',1);108 109 md = checkfield(md,'fieldname','hydrology.spcsediment_head',' forcing',1);107 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'timeseries',1); 108 109 md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1); 110 110 md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0,'numel',1); 111 111 md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0,'numel',1); … … 114 114 115 115 if obj.isefficientlayer==1, 116 md = checkfield(md,'fieldname','hydrology.spcepl_head',' forcing',1);116 md = checkfield(md,'fieldname','hydrology.spcepl_head','timeseries',1); 117 117 md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices 1],'values',[0 1]); 118 118 md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0,'numel',1); … … 195 195 WriteData(fid,'object',obj,'fieldname','leakage_factor','format','Double'); 196 196 end 197 WriteData(fid,'object',obj,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)198 199 WriteData(fid,'object',obj,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);197 WriteData(fid,'object',obj,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 198 199 WriteData(fid,'object',obj,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 200 200 WriteData(fid,'object',obj,'fieldname','sediment_compressibility','format','Double'); 201 201 WriteData(fid,'object',obj,'fieldname','sediment_porosity','format','Double'); … … 204 204 205 205 if obj.isefficientlayer==1, 206 WriteData(fid,'object',obj,'fieldname','spcepl_head','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);206 WriteData(fid,'object',obj,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 207 207 WriteData(fid,'object',obj,'fieldname','mask_eplactive_node','format','DoubleMat','mattype',1); 208 208 WriteData(fid,'object',obj,'fieldname','epl_compressibility','format','Double'); -
TabularUnified issm/trunk-jpl/src/m/classes/hydrologydc.py ¶
r18986 r19027 154 154 md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0.,'numel',[1]) 155 155 156 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,' forcing',1)157 md = checkfield(md,'fieldname','hydrology.spcsediment_head',' forcing',1)156 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'timeseries',1) 157 md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1) 158 158 md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0.,'numel',[1]) 159 159 md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0.,'numel',[1]) … … 161 161 md = checkfield(md,'fieldname','hydrology.sediment_transmitivity','>=',0,'size',[md.mesh.numberofvertices,1]) 162 162 if self.isefficientlayer==1: 163 md = checkfield(md,'fieldname','hydrology.spcepl_head',' forcing',1)163 md = checkfield(md,'fieldname','hydrology.spcepl_head','timeseries',1) 164 164 md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices,1],'values',[0,1]) 165 165 md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0.,'numel',[1]) … … 190 190 WriteData(fid,'object',self,'fieldname','leakage_factor','format','Double') 191 191 192 WriteData(fid,'object',self,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)193 WriteData(fid,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)192 WriteData(fid,'object',self,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 193 WriteData(fid,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 194 194 WriteData(fid,'object',self,'fieldname','sediment_compressibility','format','Double') 195 195 WriteData(fid,'object',self,'fieldname','sediment_porosity','format','Double') … … 198 198 199 199 if self.isefficientlayer==1: 200 WriteData(fid,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)200 WriteData(fid,'object',self,'fieldname','spcepl_head','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 201 201 WriteData(fid,'object',self,'fieldname','mask_eplactive_node','format','DoubleMat','mattype',1) 202 202 WriteData(fid,'object',self,'fieldname','epl_compressibility','format','Double') -
TabularUnified issm/trunk-jpl/src/m/classes/hydrologyshreve.m ¶
r17881 r19027 43 43 end 44 44 45 md = checkfield(md,'fieldname','hydrology.spcwatercolumn',' forcing',1);45 md = checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1); 46 46 md = checkfield(md,'fieldname','hydrology.stabilization','>=',0); 47 47 end % }}} … … 54 54 function marshall(obj,md,fid) % {{{ 55 55 WriteData(fid,'enum',HydrologyModelEnum(),'data',HydrologyshreveEnum(),'format','Integer'); 56 WriteData(fid,'object',obj,'fieldname','spcwatercolumn','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);56 WriteData(fid,'object',obj,'fieldname','spcwatercolumn','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 57 57 WriteData(fid,'object',obj,'fieldname','stabilization','format','Double'); 58 58 end % }}} -
TabularUnified issm/trunk-jpl/src/m/classes/hydrologyshreve.py ¶
r17881 r19027 40 40 return md 41 41 42 md = checkfield(md,'fieldname','hydrology.spcwatercolumn',' forcing',1)42 md = checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1) 43 43 md = checkfield(md,'fieldname','hydrology.stabilization','>=',0) 44 44 … … 47 47 def marshall(self,md,fid): # {{{ 48 48 WriteData(fid,'enum',HydrologyModelEnum(),'data',HydrologyshreveEnum(),'format','Integer'); 49 WriteData(fid,'object',self,'fieldname','spcwatercolumn','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)49 WriteData(fid,'object',self,'fieldname','spcwatercolumn','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 50 50 WriteData(fid,'object',self,'fieldname','stabilization','format','Double') 51 51 # }}} -
TabularUnified issm/trunk-jpl/src/m/classes/linearbasalforcings.m ¶
r18096 r19027 49 49 50 50 if ismember(MasstransportAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.ismasstransport==0), 51 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1);51 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1); 52 52 md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0,'numel',1); 53 53 md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<','basalforcings.upperwater_elevation','numel',1); … … 61 61 end 62 62 if ismember(ThermalAnalysisEnum(),analyses) & ~(solution==TransientSolutionEnum() & md.transient.isthermal==0), 63 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1);63 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1); 64 64 md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0,'numel',1); 65 65 md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<','basalforcings.upperwater_elevation','numel',1); 66 66 md = checkfield(md,'fieldname','basalforcings.upperwater_elevation','<',0,'numel',1); 67 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,' forcing',1,'>=',0);67 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0); 68 68 end 69 69 end % }}} … … 87 87 floatingice_melting_rate(pos)=md.basalforcings.deepwater_melting_rate*(md.geometry.base(pos)-md.basalforcings.upperwater_elevation)/(md.basalforcings.deepwater_elevation-md.basalforcings.upperwater_elevation); 88 88 WriteData(fid,'enum',BasalforcingsEnum(),'data',LinearFloatingMeltRateEnum(),'format','Integer'); 89 WriteData(fid,'data',floatingice_melting_rate,'format','DoubleMat','enum',BasalforcingsFloatingiceMeltingRateEnum(),'mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)90 WriteData(fid,'object',obj,'fieldname','groundedice_melting_rate','format','DoubleMat','enum',BasalforcingsGroundediceMeltingRateEnum(),'mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)91 WriteData(fid,'object',obj,'fieldname','geothermalflux','enum',BasalforcingsGeothermalfluxEnum(),'format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);89 WriteData(fid,'data',floatingice_melting_rate,'format','DoubleMat','enum',BasalforcingsFloatingiceMeltingRateEnum(),'mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 90 WriteData(fid,'object',obj,'fieldname','groundedice_melting_rate','format','DoubleMat','enum',BasalforcingsGroundediceMeltingRateEnum(),'mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 91 WriteData(fid,'object',obj,'fieldname','geothermalflux','enum',BasalforcingsGeothermalfluxEnum(),'format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 92 92 WriteData(fid,'object',obj,'fieldname','deepwater_melting_rate','format','Double','enum',BasalforcingsDeepwaterMeltingRateEnum(),'scale',1./yts) 93 93 WriteData(fid,'object',obj,'fieldname','deepwater_elevation','format','Double','enum',BasalforcingsDeepwaterElevationEnum()) -
TabularUnified issm/trunk-jpl/src/m/classes/linearbasalforcings.py ¶
r18116 r19027 69 69 70 70 if MasstransportAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.ismasstransport): 71 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1)71 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1) 72 72 md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0); 73 73 md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<',md.basalforcings.upperwater_elevation); … … 81 81 82 82 if ThermalAnalysisEnum() in analyses and not (solution==TransientSolutionEnum() and not md.transient.isthermal): 83 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,' forcing',1)83 md = checkfield(md,'fieldname','basalforcings.groundedice_melting_rate','NaN',1,'timeseries',1) 84 84 md = checkfield(md,'fieldname','basalforcings.deepwater_melting_rate','>=',0); 85 85 md = checkfield(md,'fieldname','basalforcings.deepwater_elevation','<',md.basalforcings.upperwater_elevation); 86 86 md = checkfield(md,'fieldname','basalforcings.upperwater_elevation','<',0); 87 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,' forcing',1,'>=',0)87 md = checkfield(md,'fieldname','basalforcings.geothermalflux','NaN',1,'timeseries',1,'>=',0) 88 88 89 89 return md … … 100 100 101 101 WriteData(fid,'enum',BasalforcingsEnum(),'data',LinearFloatingMeltRateEnum(),'format','Integer'); 102 WriteData(fid,'object',self,'fieldname','groundedice_melting_rate','enum',BasalforcingsGroundediceMeltingRateEnum(),'format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)103 WriteData(fid,'data',floatingice_melting_rate,'enum',BasalforcingsFloatingiceMeltingRateEnum(),'format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)104 WriteData(fid,'object',self,'fieldname','geothermalflux','enum',BasalforcingsGeothermalfluxEnum(),'format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)102 WriteData(fid,'object',self,'fieldname','groundedice_melting_rate','enum',BasalforcingsGroundediceMeltingRateEnum(),'format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 103 WriteData(fid,'data',floatingice_melting_rate,'enum',BasalforcingsFloatingiceMeltingRateEnum(),'format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 104 WriteData(fid,'object',self,'fieldname','geothermalflux','enum',BasalforcingsGeothermalfluxEnum(),'format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 105 105 WriteData(fid,'object',self,'fieldname','deepwater_melting_rate','enum',BasalforcingsDeepwaterMeltingRateEnum(),'format','Double','scale',1./yts) 106 106 WriteData(fid,'object',self,'fieldname','deepwater_elevation','enum',BasalforcingsDeepwaterElevationEnum(),'format','Double') -
TabularUnified issm/trunk-jpl/src/m/classes/masscon.m ¶
r18881 r19027 43 43 md = checkfield(md,'fieldname','obj.definitionenum','field',obj.definitionenum,'values',[Outputdefinition1Enum,Outputdefinition2Enum,Outputdefinition3Enum,Outputdefinition4Enum,Outputdefinition5Enum,Outputdefinition6Enum,Outputdefinition7Enum,Outputdefinition8Enum,Outputdefinition9Enum,Outputdefinition10Enum]); 44 44 45 md = checkfield(md,'fieldname','obj.levelset','field',obj.levelset,' forcing',1,'NaN',1);45 md = checkfield(md,'fieldname','obj.levelset','field',obj.levelset,'timeseries',1,'NaN',1); 46 46 47 47 end % }}} -
TabularUnified issm/trunk-jpl/src/m/classes/masstransport.m ¶
r18757 r19027 87 87 if ~ismember(MasstransportAnalysisEnum(),analyses) | (solution==TransientSolutionEnum() & md.transient.ismasstransport==0), return; end 88 88 89 md = checkfield(md,'fieldname','masstransport.spcthickness',' forcing',1);89 md = checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1); 90 90 md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0 1]); 91 91 md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',{'Absolute' 'Incremental'}); … … 113 113 yts=365.*24.*3600.; 114 114 115 WriteData(fid,'object',obj,'fieldname','spcthickness','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);115 WriteData(fid,'object',obj,'fieldname','spcthickness','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 116 116 WriteData(fid,'object',obj,'fieldname','isfreesurface','format','Boolean'); 117 117 WriteData(fid,'object',obj,'fieldname','min_thickness','format','Double'); -
TabularUnified issm/trunk-jpl/src/m/classes/masstransport.py ¶
r18756 r19027 67 67 return md 68 68 69 md = checkfield(md,'fieldname','masstransport.spcthickness',' forcing',1)69 md = checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1) 70 70 md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0,1]) 71 71 md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',['Absolute','Incremental']) … … 80 80 yts=365.*24.*3600. 81 81 82 WriteData(fid,'object',self,'fieldname','spcthickness','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)82 WriteData(fid,'object',self,'fieldname','spcthickness','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 83 83 WriteData(fid,'object',self,'fieldname','isfreesurface','format','Boolean') 84 84 WriteData(fid,'object',self,'fieldname','min_thickness','format','Double') -
TabularUnified issm/trunk-jpl/src/m/classes/matice.m ¶
r19014 r19027 143 143 md = checkfield(md,'fieldname','materials.rho_freshwater','>',0); 144 144 md = checkfield(md,'fieldname','materials.mu_water','>',0); 145 md = checkfield(md,'fieldname','materials.rheology_B','>',0,' forcing',1,'NaN',1);145 md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1); 146 146 md = checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements 1]); 147 147 md = checkfield(md,'fieldname','materials.rheology_law','values',{'None' 'Cuffey' 'Paterson' 'Arrhenius' 'LliboutryDuval'}); … … 192 192 WriteData(fid,'object',obj,'class','materials','fieldname','mixed_layer_capacity','format','Double'); 193 193 WriteData(fid,'object',obj,'class','materials','fieldname','thermal_exchange_velocity','format','Double'); 194 WriteData(fid,'object',obj,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);194 WriteData(fid,'object',obj,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 195 195 WriteData(fid,'object',obj,'class','materials','fieldname','rheology_n','format','DoubleMat','mattype',2); 196 196 WriteData(fid,'data',StringToEnum(obj.rheology_law),'enum',MaterialsRheologyLawEnum(),'format','Integer'); -
TabularUnified issm/trunk-jpl/src/m/classes/matice.py ¶
r19014 r19027 117 117 md = checkfield(md,'fieldname','materials.rho_freshwater','>',0) 118 118 md = checkfield(md,'fieldname','materials.mu_water','>',0) 119 md = checkfield(md,'fieldname','materials.rheology_B','>',0,' forcing',1,'NaN',1)119 md = checkfield(md,'fieldname','materials.rheology_B','>',0,'timeseries',1,'NaN',1) 120 120 md = checkfield(md,'fieldname','materials.rheology_n','>',0,'size',[md.mesh.numberofelements]) 121 121 md = checkfield(md,'fieldname','materials.rheology_law','values',['None','Cuffey','Paterson','Arrhenius','LliboutryDuval']) … … 140 140 WriteData(fid,'object',self,'class','materials','fieldname','mixed_layer_capacity','format','Double') 141 141 WriteData(fid,'object',self,'class','materials','fieldname','thermal_exchange_velocity','format','Double') 142 WriteData(fid,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)142 WriteData(fid,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 143 143 WriteData(fid,'object',self,'class','materials','fieldname','rheology_B','format','DoubleMat','mattype',1) 144 144 WriteData(fid,'object',self,'class','materials','fieldname','rheology_n','format','DoubleMat','mattype',2) -
TabularUnified issm/trunk-jpl/src/m/classes/misfit.m ¶
r18881 r19027 60 60 error('misfit error message: ''timeinterpolation'' field should be a string!'); 61 61 end 62 md = checkfield(md,'fieldname','obj.observation','field',obj.observation,' forcing',1,'NaN',1);62 md = checkfield(md,'fieldname','obj.observation','field',obj.observation,'timeseries',1,'NaN',1); 63 63 md = checkfield(md,'fieldname','obj.timeinterpolation','field',obj.timeinterpolation,'values',{'nearestneighbor'}); 64 md = checkfield(md,'fieldname','obj.weights','field',obj.weights,' forcing',1,'NaN',1);64 md = checkfield(md,'fieldname','obj.weights','field',obj.weights,'timeseries',1,'NaN',1); 65 65 66 66 end % }}} … … 84 84 WriteData(fid,'object',obj,'fieldname','definitionenum','format','Integer'); 85 85 WriteData(fid,'object',obj,'fieldname','model_enum','format','Integer'); 86 WriteData(fid,'object',obj,'fieldname','observation','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);86 WriteData(fid,'object',obj,'fieldname','observation','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 87 87 WriteData(fid,'object',obj,'fieldname','observation_enum','format','Integer'); 88 88 WriteData(fid,'object',obj,'fieldname','timeinterpolation','format','String'); 89 WriteData(fid,'object',obj,'fieldname','weights','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);89 WriteData(fid,'object',obj,'fieldname','weights','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 90 90 WriteData(fid,'object',obj,'fieldname','weights_enum','format','Integer'); 91 91 -
TabularUnified issm/trunk-jpl/src/m/classes/snowpack.m ¶
r16764 r19027 455 455 yts=365.0*24.0*3600.0; 456 456 457 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);458 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);459 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);457 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 458 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 459 WriteData(fid,'object',obj,'class','snowpack','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 460 460 WriteData(fid,'object',obj,'class','snowpack','fieldname','restol','format','Double'); 461 461 WriteData(fid,'object',obj,'class','snowpack','fieldname','reltol','format','Double'); -
TabularUnified issm/trunk-jpl/src/m/classes/stressbalance.m ¶
r18652 r19027 126 126 if ~ismember(StressbalanceAnalysisEnum(),analyses), return; end 127 127 128 md = checkfield(md,'fieldname','stressbalance.spcvx',' forcing',1);129 md = checkfield(md,'fieldname','stressbalance.spcvy',' forcing',1);130 md = checkfield(md,'fieldname','stressbalance.spcvz',' forcing',1);128 md = checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1); 129 md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1); 130 md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1); 131 131 md = checkfield(md,'fieldname','stressbalance.restol','size',[1 1],'>',0,'NaN',1); 132 132 md = checkfield(md,'fieldname','stressbalance.reltol','size',[1 1]); … … 218 218 yts=365.0*24.0*3600.0; 219 219 220 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);221 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);222 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1);220 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 221 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 222 WriteData(fid,'object',obj,'class','stressbalance','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1); 223 223 WriteData(fid,'object',obj,'class','stressbalance','fieldname','restol','format','Double'); 224 224 WriteData(fid,'object',obj,'class','stressbalance','fieldname','reltol','format','Double'); -
TabularUnified issm/trunk-jpl/src/m/classes/stressbalance.py ¶
r18619 r19027 123 123 return md 124 124 125 md = checkfield(md,'fieldname','stressbalance.spcvx',' forcing',1)126 md = checkfield(md,'fieldname','stressbalance.spcvy',' forcing',1)125 md = checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1) 126 md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1) 127 127 if m.strcmp(md.mesh.domaintype(),'3D'): 128 md = checkfield(md,'fieldname','stressbalance.spcvz',' forcing',1)128 md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1) 129 129 md = checkfield(md,'fieldname','stressbalance.restol','size',[1],'>',0) 130 130 md = checkfield(md,'fieldname','stressbalance.reltol','size',[1]) … … 168 168 yts=365.0*24.0*3600.0 169 169 170 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)171 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)172 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,' forcinglength',md.mesh.numberofvertices+1)170 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvx','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 171 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvy','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 172 WriteData(fid,'object',self,'class','stressbalance','fieldname','spcvz','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 173 173 WriteData(fid,'object',self,'class','stressbalance','fieldname','restol','format','Double') 174 174 WriteData(fid,'object',self,'class','stressbalance','fieldname','reltol','format','Double') -
TabularUnified issm/trunk-jpl/src/m/classes/thermal.m ¶
r18754 r19027 91 91 92 92 md = checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0 1 2]); 93 md = checkfield(md,'fieldname','thermal.spctemperature',' forcing',1);93 md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1); 94 94 if (ismember(EnthalpyAnalysisEnum(),analyses) & md.thermal.isenthalpy & dimension(md.mesh)==3), 95 95 pos=find(md.thermal.spctemperature(1:md.mesh.numberofvertices,:)~=NaN); … … 124 124 end % }}} 125 125 function marshall(obj,md,fid) % {{{ 126 WriteData(fid,'object',obj,'fieldname','spctemperature','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1);126 WriteData(fid,'object',obj,'fieldname','spctemperature','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1); 127 127 WriteData(fid,'object',obj,'fieldname','penalty_threshold','format','Integer'); 128 128 WriteData(fid,'object',obj,'fieldname','stabilization','format','Integer'); -
TabularUnified issm/trunk-jpl/src/m/classes/thermal.py ¶
r18620 r19027 86 86 87 87 md = checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0,1,2]) 88 md = checkfield(md,'fieldname','thermal.spctemperature',' forcing',1)88 md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1) 89 89 if EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy and md.mesh.dimension()==3: 90 90 pos=numpy.nonzero(numpy.logical_not(numpy.isnan(md.thermal.spctemperature[0:md.mesh.numberofvertices]))) … … 102 102 # }}} 103 103 def marshall(self,md,fid): # {{{ 104 WriteData(fid,'object',self,'fieldname','spctemperature','format','DoubleMat','mattype',1,' forcinglength',md.mesh.numberofvertices+1)104 WriteData(fid,'object',self,'fieldname','spctemperature','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1) 105 105 WriteData(fid,'object',self,'fieldname','penalty_threshold','format','Integer') 106 106 WriteData(fid,'object',self,'fieldname','stabilization','format','Integer') -
TabularUnified issm/trunk-jpl/src/m/consistency/checkfield.m ¶
r18077 r19027 14 14 % - <=: smaller or equal to provided value 15 15 % - < vec: smallerthan provided values on each vertex 16 % - forcing: 1 if check forcingconsistency (size and time)16 % - timeseries: 1 if check time series consistency (size and time) 17 17 % - values: cell of strings or vector of acceptable values 18 18 % - numel: list of acceptable number of elements … … 198 198 199 199 %Check forcings (size and times) 200 if getfieldvalue(options,' forcing',0),200 if getfieldvalue(options,'timeseries',0), 201 201 if size(field,1)==md.mesh.numberofvertices, 202 202 if ~size(field,2)==1, -
TabularUnified issm/trunk-jpl/src/m/consistency/checkfield.py ¶
r18077 r19027 21 21 - <=: smaller or equal to provided value 22 22 - < vec: smallerthan provided values on each vertex 23 - forcing: 1 if check forcingconsistency (size and time)23 - timeseries: 1 if check time series consistency (size and time) 24 24 - values: cell of strings or vector of acceptable values 25 25 - numel: list of acceptable number of elements … … 151 151 152 152 #Check forcings (size and times) 153 if options.getfieldvalue(' forcing',0):153 if options.getfieldvalue('timeseries',0): 154 154 if numpy.size(field,0)==md.mesh.numberofvertices: 155 155 if numpy.ndim(field)>1 and not numpy.size(field,1)==1: -
TabularUnified issm/trunk-jpl/src/m/solve/WriteData.m ¶
r17945 r19027 27 27 format = getfieldvalue(options,'format'); 28 28 mattype = getfieldvalue(options,'mattype',0); %only required for matrices 29 forcinglength = getfieldvalue(options,'forcinglength',-1);29 timeserieslength = getfieldvalue(options,'timeserieslength',-1); 30 30 31 31 %Process sparse matrices … … 37 37 if exist(options,'scale'), 38 38 scale = getfieldvalue(options,'scale'); 39 if size(data,1)== forcinglength,39 if size(data,1)==timeserieslength, 40 40 data(1:end-1,:) = scale.*data(1:end-1,:); 41 41 else … … 43 43 end 44 44 end 45 if(size(data,1)== forcinglength),45 if(size(data,1)==timeserieslength), 46 46 yts=365.0*24.0*3600.0; 47 47 data(end,:) = data(end,:)*yts; -
TabularUnified issm/trunk-jpl/src/m/solve/WriteData.py ¶
r17945 r19027 35 35 format = options.getfieldvalue('format') 36 36 mattype = options.getfieldvalue('mattype',0) #only required for matrices 37 forcinglength = options.getfieldvalue('forcinglength',-1)37 timeserieslength = options.getfieldvalue('timeserieslength',-1) 38 38 39 39 #Process sparse matrices … … 46 46 scale = options.getfieldvalue('scale') 47 47 if numpy.size(data) > 1 : 48 if numpy.size(data,0)== forcinglength:48 if numpy.size(data,0)==timeserieslength: 49 49 data=numpy.array(data) 50 50 data[0:-1,:] = scale*data[0:-1,:] … … 54 54 data = scale*data 55 55 if numpy.size(data) > 1 : 56 if numpy.size(data,0)== forcinglength:56 if numpy.size(data,0)==timeserieslength: 57 57 yts=365.0*24.0*3600.0 58 58 data[-1,:] = yts*data[-1,:]
Note:
See TracChangeset
for help on using the changeset viewer.