Changeset 14640 for issm/trunk-jpl/src/m/classes/initialization.py
- Timestamp:
- 04/18/13 10:20:16 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/initialization.py
r14577 r14640 1 #module imports2 1 import numpy 3 2 from fielddisplay import fielddisplay … … 14 13 """ 15 14 16 #properties 17 def __init__(self): 18 # {{{ Properties 15 def __init__(self): # {{{ 19 16 20 17 self.vx = float('NaN') … … 29 26 self.watercolumn = float('NaN') 30 27 self.sediment_head = float('NaN') 31 self.water_pressure = float('NaN')32 28 33 29 #set defaults … … 35 31 36 32 #}}} 37 def __repr__(self): 38 # {{{ Display 33 def __repr__(self): # {{{ 39 34 string=' initial field values:' 40 35 41 string="%s\n%s"%(string,fielddisplay(self,'vx','x component of velocity '))42 string="%s\n%s"%(string,fielddisplay(self,'vy','y component of velocity '))43 string="%s\n%s"%(string,fielddisplay(self,'vz','z component of velocity '))44 string="%s\n%s"%(string,fielddisplay(self,'vel','velocity norm '))45 string="%s\n%s"%(string,fielddisplay(self,'pressure','pressure field'))46 string="%s\n%s"%(string,fielddisplay(self,'temperature','temperature in Kelvins'))47 string="%s\n%s"%(string,fielddisplay(self,'surfacetemp','surface temperature in Kelvins'))48 string="%s\n%s"%(string,fielddisplay(self,'basaltemp','basal temperature in Kelvins'))36 string="%s\n%s"%(string,fielddisplay(self,'vx','x component of velocity [m/yr]')) 37 string="%s\n%s"%(string,fielddisplay(self,'vy','y component of velocity [m/yr]')) 38 string="%s\n%s"%(string,fielddisplay(self,'vz','z component of velocity [m/yr]')) 39 string="%s\n%s"%(string,fielddisplay(self,'vel','velocity norm [m/yr]')) 40 string="%s\n%s"%(string,fielddisplay(self,'pressure','pressure [Pa]')) 41 string="%s\n%s"%(string,fielddisplay(self,'temperature','temperature [K]')) 42 string="%s\n%s"%(string,fielddisplay(self,'surfacetemp','surface temperature [K]')) 43 string="%s\n%s"%(string,fielddisplay(self,'basaltemp','basal temperature [K]')) 49 44 string="%s\n%s"%(string,fielddisplay(self,'waterfraction','fraction of water in the ice')) 50 string="%s\n%s"%(string,fielddisplay(self,'watercolumn','thickness of subglacial water')) 51 string="%s\n%s"%(string,fielddisplay(self,'sediment_head','sediment water head of subglacial system (meters)')) 52 string="%s\n%s"%(string,fielddisplay(self,'water_pressure','water pressure due to the sediment head (Pa)')) 45 string="%s\n%s"%(string,fielddisplay(self,'watercolumn','thickness of subglacial water [m]')) 46 string="%s\n%s"%(string,fielddisplay(self,'sediment_head','sediment water head of subglacial system [m]')) 53 47 54 48 return string 55 49 #}}} 56 57 def setdefaultparameters(self): 58 # {{{setdefaultparameters 50 def setdefaultparameters(self): # {{{ 59 51 return self 60 52 #}}} 61 62 53 def checkconsistency(self,md,solution,analyses): # {{{ 63 54 if DiagnosticHorizAnalysisEnum() in analyses: … … 86 77 return md 87 78 # }}} 88 89 79 def marshall(self,fid): # {{{ 90 80 WriteData(fid,'data',self.vx,'format','DoubleMat','mattype',1,'enum',VxEnum()) … … 97 87 WriteData(fid,'data',self.sediment_head,'format','DoubleMat','mattype',1,'enum',SedimentHeadEnum()) 98 88 # }}} 99
Note:
See TracChangeset
for help on using the changeset viewer.