Changeset 19527 for issm/trunk-jpl/src/m/classes/SMBgradients.py

Timestamp:

08/29/15 11:55:39 (10 years ago)

Author:

Eric.Larour

Message:

CHG: moved md.surfaceforcings to md.smb.
By doing so, had to rename the SMB class to SMBforcing class (it's just that, a mass_balance forcing inside
a SMB class, hence the name).
We also now have an smb_core solution, taken out of the mass transport core. Makes more sense long term.
Synced all enums according to the new changes, and operated the adjustments in all the test decks.

In addition, progressing in terms of GEMB integration into ISSM, specifically at the SMBgemb level (which
is spurring all the changes described above). Brought the class up to the level of the GEMB.m call in Alex's
code. Starting the C integration now.

File:

• issm/trunk-jpl/src/m/classes/SMBgradients.py (modified) (5 diffs)

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

@@ -18 +18 @@
                 self.b_pos   = float('NaN')
                 self.b_neg   = float('NaN')
+                self.requested_outputs      = []
                 #}}}
         def __repr__(self): # {{{
@@ -27 +28 @@
                 string="%s\n%s"%(string,fielddisplay(self,'b_pos',' slope of hs - smb regression line for accumulation regime required if smb gradients is activated'))
                 string="%s\n%s"%(string,fielddisplay(self,'b_neg',' slope of hs - smb regression line for ablation regime required if smb gradients is activated'))
+                string="%s\n%s"%(string,fielddisplay(self,'requested_outputs','additional outputs requested'))
 
                 return string
@@ -34 +36 @@
                 #Nothing for now
                 return self
+        #}}}
+        def defaultoutputs(self,md): # {{{
+                return []
         #}}}
         def initialize(self,md): # {{{
@@ -44 +49 @@
 
                 if MasstransportAnalysisEnum() in analyses:
-                        md = checkfield(md,'fieldname','surfaceforcings.href','timeseries',1,'NaN',1)
-                        md = checkfield(md,'fieldname','surfaceforcings.smbref','timeseries',1,'NaN',1)
-                        md = checkfield(md,'fieldname','surfaceforcings.b_pos','timeseries',1,'NaN',1)
-                        md = checkfield(md,'fieldname','surfaceforcings.b_neg','timeseries',1,'NaN',1)
+                        md = checkfield(md,'fieldname','smb.href','timeseries',1,'NaN',1)
+                        md = checkfield(md,'fieldname','smb.smbref','timeseries',1,'NaN',1)
+                        md = checkfield(md,'fieldname','smb.b_pos','timeseries',1,'NaN',1)
+                        md = checkfield(md,'fieldname','smb.b_neg','timeseries',1,'NaN',1)
 
+                md = checkfield(md,'fieldname','masstransport.requested_outputs','stringrow',1)
                 return md
         # }}}
@@ -55 +61 @@
                 yts=365.0*24.0*3600.0
 
-                WriteData(fid,'enum',SurfaceforcingsEnum(),'data',SMBgradientsEnum(),'format','Integer');
-                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','href','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1)
-                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
-                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
-                WriteData(fid,'object',self,'class','surfaceforcings','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
+                WriteData(fid,'enum',SmbEnum(),'data',SMBgradientsEnum(),'format','Integer');
+                WriteData(fid,'object',self,'class','smb','fieldname','href','format','DoubleMat','mattype',1,'timeserieslength',md.mesh.numberofvertices+1)
+                WriteData(fid,'object',self,'class','smb','fieldname','smbref','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
+                WriteData(fid,'object',self,'class','smb','fieldname','b_pos','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
+                WriteData(fid,'object',self,'class','smb','fieldname','b_neg','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1)
+                
+                #process requested outputs
+                outputs = self.requested_outputs
+                indices = [i for i, x in enumerate(outputs) if x == 'default']
+                if len(indices) > 0:
+                        outputscopy=outputs[0:max(0,indices[0]-1)]+self.defaultoutputs(md)+outputs[indices[0]+1:]
+                        outputs    =outputscopy
+                WriteData(fid,'data',outputs,'enum',SmbRequestedOutputsEnum(),'format','StringArray')
+
         # }}}