Changeset 14198
- Timestamp:
- 12/19/12 10:19:22 (12 years ago)
- Location:
- issm/trunk-jpl
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/surfaceforcings.py
r14141 r14198 21 21 self.issmbgradients = 0 22 22 self.isdelta18o = 0 23 self.hc = float('NaN')24 23 self.href = float('NaN') 25 24 self.smbref = float('NaN') 26 self.smb_pos_max = float('NaN')27 self.smb_pos_min = float('NaN')28 self.a_pos = float('NaN')29 25 self.b_pos = float('NaN') 30 self.a_neg = float('NaN')31 26 self.b_neg = float('NaN') 32 27 self.monthlytemperatures = float('NaN') … … 57 52 string="%s\n%s"%(string,fielddisplay(self,'delta18o_surface','surface elevation of the delta18o site, required if pdd is activated and delta18o activated')) 58 53 string="%s\n%s"%(string,fielddisplay(self,'issmbgradients','is smb gradients method activated (0 or 1, default is 0)')) 59 string="%s\n%s"%(string,fielddisplay(self,'hc',' elevation of intersection between accumulation and ablation regime required if smb gradients is activated'))60 54 string="%s\n%s"%(string,fielddisplay(self,'href',' reference elevation from which deviation is used to calculate SMB adjustment in smb gradients method')) 61 55 string="%s\n%s"%(string,fielddisplay(self,'smbref',' reference smb from which deviation is calculated in smb gradients method')) 62 string="%s\n%s"%(string,fielddisplay(self,'smb_pos_max',' maximum value of positive smb required if smb gradients is activated'))63 string="%s\n%s"%(string,fielddisplay(self,'smb_pos_min',' minimum value of positive smb required if smb gradients is activated'))64 string="%s\n%s"%(string,fielddisplay(self,'a_pos',' intercept of hs - smb regression line for accumulation regime required if smb gradients is activated'))65 56 string="%s\n%s"%(string,fielddisplay(self,'b_pos',' slope of hs - smb regression line for accumulation regime required if smb gradients is activated')) 66 string="%s\n%s"%(string,fielddisplay(self,'a_neg',' intercept of hs - smb regression line for ablation regime required if smb gradients is activated'))67 57 string="%s\n%s"%(string,fielddisplay(self,'b_neg',' slope of hs - smb regression line for ablation regime required if smb gradients is activated')) 68 58 … … 96 86 md = checkfield(md,'surfaceforcings.precipitations_presentday','size',[md.mesh.numberofvertices+1,12],'NaN',1) 97 87 elif self.issmbgradients: 98 md = checkfield(md,'surfaceforcings.hc','forcing',1,'NaN',1)99 88 md = checkfield(md,'surfaceforcings.href','forcing',1,'NaN',1) 100 89 md = checkfield(md,'surfaceforcings.smbref','forcing',1,'NaN',1) 101 md = checkfield(md,'surfaceforcings.smb_pos_max','forcing',1,'NaN',1)102 md = checkfield(md,'surfaceforcings.smb_pos_min','forcing',1,'NaN',1)103 md = checkfield(md,'surfaceforcings.a_pos','forcing',1,'NaN',1)104 90 md = checkfield(md,'surfaceforcings.b_pos','forcing',1,'NaN',1) 105 md = checkfield(md,'surfaceforcings.a_neg','forcing',1,'NaN',1)106 91 md = checkfield(md,'surfaceforcings.b_neg','forcing',1,'NaN',1) 107 92 else: … … 134 119 135 120 if self.issmbgradients: 136 WriteData(fid,'object',self,'fieldname','hc','format','DoubleMat','mattype',1)137 121 WriteData(fid,'object',self,'fieldname','href','format','DoubleMat','mattype',1) 138 122 WriteData(fid,'object',self,'fieldname','smbref','format','DoubleMat','mattype',1) 139 WriteData(fid,'object',self,'fieldname','smb_pos_max','format','DoubleMat','mattype',1)140 WriteData(fid,'object',self,'fieldname','smb_pos_min','format','DoubleMat','mattype',1)141 WriteData(fid,'object',self,'fieldname','a_pos','format','DoubleMat','mattype',1)142 123 WriteData(fid,'object',self,'fieldname','b_pos','format','DoubleMat','mattype',1) 143 WriteData(fid,'object',self,'fieldname','a_neg','format','DoubleMat','mattype',1)144 124 WriteData(fid,'object',self,'fieldname','b_neg','format','DoubleMat','mattype',1) 145 125 # }}} -
issm/trunk-jpl/test/NightlyRun/test328.m
r14179 r14198 7 7 md.surfaceforcings.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y; 8 8 md.surfaceforcings.href=md.geometry.surface; 9 md.surfaceforcings.smbref= 1000 - 0.001*md.mesh.x - 0.005*md.mesh.y;9 md.surfaceforcings.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y; 10 10 md.transient.requested_outputs=TotalSmbEnum(); 11 11 md.cluster=generic('name',oshostname(),'np',3); -
issm/trunk-jpl/test/NightlyRun/test328.py
r14190 r14198 19 19 md.transient.requested_outputs=TotalSmbEnum() 20 20 md.surfaceforcings.href=copy.deepcopy(md.geometry.surface).reshape(-1) 21 md.surfaceforcings.smbref= 1000 - 0.001*md.mesh.x - 0.005*md.mesh.y;21 md.surfaceforcings.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y; 22 22 md.cluster=generic('name',oshostname(),'np',3) 23 23 md=solve(md,TransientSolutionEnum()) -
issm/trunk-jpl/test/NightlyRun/test329.m
r14179 r14198 8 8 md.surfaceforcings.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y; 9 9 md.surfaceforcings.href=md.geometry.surface; 10 md.surfaceforcings.smbref= 1000 - 0.001*md.mesh.x - 0.005*md.mesh.y;10 md.surfaceforcings.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y; 11 11 md.transient.requested_outputs=TotalSmbEnum(); 12 12 md.cluster=generic('name',oshostname(),'np',3); -
issm/trunk-jpl/test/NightlyRun/test329.py
r14190 r14198 19 19 md.surfaceforcings.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y 20 20 md.surfaceforcings.href=copy.deepcopy(md.geometry.surface).reshape(-1) 21 md.surfaceforcings.smbref= 1000 - 0.001*md.mesh.x - 0.005*md.mesh.y;21 md.surfaceforcings.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y; 22 22 md.transient.requested_outputs=TotalSmbEnum() 23 23 md.cluster=generic('name',oshostname(),'np',3)
Note:
See TracChangeset
for help on using the changeset viewer.