Changeset 13686

Timestamp:

10/16/12 08:42:06 (12 years ago)

Author:

jschierm

Message:

CHG: Update surfaceforcings.py corresponding to changeset [13521].

File:

• issm/trunk-jpl/src/m/classes/surfaceforcings.py (modified) (4 diffs)

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

@@ -22 +22 @@
                 self.isdelta18o = 0
                 self.hc = float('NaN')
+                self.href = float('NaN')
+                self.smbref = float('NaN')
                 self.smb_pos_max = float('NaN')
                 self.smb_pos_min = float('NaN')
@@ -56 +58 @@
                 string="%s\n%s"%(string,fielddisplay(self,'issmbgradients','is smb gradients method activated (0 or 1, default is 0)'))
                 string="%s\n%s"%(string,fielddisplay(self,'hc',' elevation of intersection between accumulation and ablation regime required if smb gradients is activated'))
+                string="%s\n%s"%(string,fielddisplay(self,'href',' reference elevation from which deviation is used to calculate SMB adjustment in smb gradients method'))
+                string="%s\n%s"%(string,fielddisplay(self,'smbref',' reference smb from which deviation is calculated in smb gradients method'))
                 string="%s\n%s"%(string,fielddisplay(self,'smb_pos_max',' maximum value of positive smb required if smb gradients is activated'))
                 string="%s\n%s"%(string,fielddisplay(self,'smb_pos_min',' minimum value of positive smb required if smb gradients is activated'))
@@ -94 +98 @@
                         elif self.issmbgradients:
                                 md = checkfield(md,'surfaceforcings.hc','forcing',1,'NaN',1)
+                                md = checkfield(md,'surfaceforcings.href','forcing',1,'NaN',1)
+                                md = checkfield(md,'surfaceforcings.smbref','forcing',1,'NaN',1)
                                 md = checkfield(md,'surfaceforcings.smb_pos_max','forcing',1,'NaN',1)
                                 md = checkfield(md,'surfaceforcings.smb_pos_min','forcing',1,'NaN',1)
@@ -130 +136 @@
                 if self.issmbgradients:
                         WriteData(fid,'object',self,'fieldname','hc','format','DoubleMat','mattype',1)
+                        WriteData(fid,'object',self,'fieldname','href','format','DoubleMat','mattype',1)
+                        WriteData(fid,'object',self,'fieldname','smbref','format','DoubleMat','mattype',1)
                         WriteData(fid,'object',self,'fieldname','smb_pos_max','format','DoubleMat','mattype',1)
                         WriteData(fid,'object',self,'fieldname','smb_pos_min','format','DoubleMat','mattype',1)