Changeset 23670 for issm/trunk-jpl/src/py3/classes/SMBforcing.py
- Timestamp:
- 01/31/19 07:34:11 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/py3/classes/SMBforcing.py
r19898 r23670 1 import numpy 1 import numpy as np 2 2 from fielddisplay import fielddisplay 3 from EnumDefinitions import *4 3 from checkfield import checkfield 5 4 from WriteData import WriteData … … 34 33 def initialize(self,md): # {{{ 35 34 36 if n umpy.all(numpy.isnan(self.mass_balance)):37 self.mass_balance=n umpy.zeros((md.mesh.numberofvertices,1))35 if np.all(np.isnan(self.mass_balance)): 36 self.mass_balance=np.zeros((md.mesh.numberofvertices)) 38 37 print(" no SMBforcing.mass_balance specified: values set as zero") 39 38 … … 42 41 def checkconsistency(self,md,solution,analyses): # {{{ 43 42 44 if MasstransportAnalysisEnum()in analyses:43 if 'MasstransportAnalysis' in analyses: 45 44 md = checkfield(md,'fieldname','smb.mass_balance','timeseries',1,'NaN',1,'Inf',1) 46 45 47 if BalancethicknessAnalysisEnum()in analyses:46 if 'BalancethicknessAnalysis' in analyses: 48 47 md = checkfield(md,'fieldname','smb.mass_balance','size',[md.mesh.numberofvertices],'NaN',1,'Inf',1) 49 48 … … 51 50 return md 52 51 # }}} 53 def marshall(self, md,fid): # {{{52 def marshall(self,prefix,md,fid): # {{{ 54 53 55 yts= 365.0*24.0*3600.054 yts=md.constants.yts 56 55 57 WriteData(fid,'enum',SmbEnum(),'data',SMBforcingEnum(),'format','Integer'); 58 WriteData(fid,'object',self,'class','smb','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1) 56 WriteData(fid,prefix,'name','md.smb.model','data',1,'format','Integer'); 57 WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','DoubleMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts) 58 #WriteData(fid,prefix,'object',self,'class','smb','fieldname','mass_balance','format','CompressedMat','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); 59 59 60 60 #process requested outputs … … 64 64 outputscopy=outputs[0:max(0,indices[0]-1)]+self.defaultoutputs(md)+outputs[indices[0]+1:] 65 65 outputs =outputscopy 66 WriteData(fid, 'data',outputs,'enum',SmbRequestedOutputsEnum(),'format','StringArray')66 WriteData(fid,prefix,'data',outputs,'name','md.smb.requested_outputs','format','StringArray') 67 67 68 68 # }}}
Note:
See TracChangeset
for help on using the changeset viewer.