- Timestamp:
- 08/20/12 17:39:30 (13 years ago)
- Location:
- issm/branches/trunk-jpl-damage
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/branches/trunk-jpl-damage
- Property svn:ignore
-
old new 1 projects 1 2 autom4te.cache 2 3 aclocal.m4
-
- Property svn:mergeinfo changed
/issm/trunk-jpl merged: 12948-13099
- Property svn:ignore
-
issm/branches/trunk-jpl-damage/src/m/classes/balancethickness.py
r12168 r13101 1 1 #module imports 2 2 from fielddisplay import fielddisplay 3 class balancethickness: 3 from EnumDefinitions import * 4 from checkfield import * 5 from WriteData import * 6 7 class balancethickness(object): 8 """ 9 BALANCETHICKNESS class definition 10 11 Usage: 12 balancethickness=balancethickness(); 13 """ 14 4 15 #properties 5 16 def __init__(self): 6 17 # {{{ Properties 7 self.spcthickness = float('NaN')8 self.thickening_rate 9 self.stabilization 18 self.spcthickness = float('NaN') 19 self.thickening_rate = float('NaN') 20 self.stabilization = 0 10 21 11 22 #set defaults … … 13 24 14 25 #}}} 15 def __repr__( obj):26 def __repr__(self): 16 27 # {{{ Display 17 28 18 29 string=' balance thickness solution parameters:' 19 30 20 string="%s\n\n%s"%(string,fielddisplay( obj,'spcthickness','thickness constraints (NaN means no constraint)'))21 string="%s\n%s"%(string,fielddisplay( obj,'thickening_rate','ice thickening rate used in the mass conservation (dh/dt)'))22 string="%s\n%s"%(string,fielddisplay( obj,'stabilization','0: None, 1: SU, 2: MacAyeal''s artificial diffusivity, 3:DG'))31 string="%s\n\n%s"%(string,fielddisplay(self,'spcthickness','thickness constraints (NaN means no constraint)')) 32 string="%s\n%s"%(string,fielddisplay(self,'thickening_rate','ice thickening rate used in the mass conservation (dh/dt)')) 33 string="%s\n%s"%(string,fielddisplay(self,'stabilization',"0: None, 1: SU, 2: MacAyeal's artificial diffusivity, 3:DG")) 23 34 return string 24 35 #}}} 25 def setdefaultparameters( obj):36 def setdefaultparameters(self): 26 37 # {{{setdefaultparameters 27 38 28 39 #Type of stabilization used 29 obj.stabilization=140 self.stabilization=1 30 41 31 return obj42 return self 32 43 #}}} 33 44 45 def checkconsistency(self,md,solution,analyses): # {{{ 46 #Early return 47 if not solution==BalancethicknessSolutionEnum(): 48 return md 49 50 md = checkfield(md,'balancethickness.spcthickness','forcing',1) 51 md = checkfield(md,'balancethickness.thickening_rate','size',[md.mesh.numberofvertices],'NaN',1) 52 md = checkfield(md,'balancethickness.stabilization','size',[1],'values',[0,1,2,3]) 53 54 return md 55 # }}} 56 57 def marshall(self,fid): # {{{ 58 WriteData(fid,'object',self,'fieldname','spcthickness','format','DoubleMat','mattype',1) 59 WriteData(fid,'object',self,'fieldname','thickening_rate','format','DoubleMat','mattype',1) 60 WriteData(fid,'object',self,'fieldname','stabilization','format','Integer') 61 # }}} 62
Note:
See TracChangeset
for help on using the changeset viewer.