Changeset 18585
- Timestamp:
- 10/07/14 09:07:21 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/hydrologydc.py
r17497 r18585 23 23 self.transfer_flag = 0 24 24 self.leakage_factor = 0 25 self.basal_moulin_input = float('NaN') 25 26 26 27 self.spcsediment_head = float('NaN') … … 36 37 self.epl_porosity = 0 37 38 self.epl_initial_thickness = 0 39 self.epl_max_thickness = 0 38 40 self.epl_conductivity = 0 39 41 40 42 #set defaults 41 43 self.setdefaultparameters() 44 42 45 #}}} 43 46 … … 59 62 if self.sedimentlimit_flag==1: 60 63 string="%s\n%s"%(string,fielddisplay(self,'sedimentlimit','user defined upper limit for the inefficient layer [m]')) 61 string="%s\n%s"%(string,fielddisplay(self,'transfer_flag',['what kind of transfer method is applied between the layers'])) 64 string="%s\n%s"%(string,fielddisplay(self,'basal_moulin_input','water flux at a given point [m3 s-1]')) 65 string="%s\n%s"%(string,fielddisplay(self,'transfer_flag','what kind of transfer method is applied between the layers')) 62 66 string='%55s 0: no transfer',' ' 63 67 string='%55s 1: constant leakage factor: %s',' ','leakage_factor' … … 103 107 self.epl_porosity = 0.4 104 108 self.epl_initial_thickness = 1.0 109 self.epl_max_thickness = 5.0 105 110 self.epl_conductivity = 8.0e-02 106 111 … … 108 113 # }}} 109 114 115 def initialize(self): # {{{ 116 if numpy.all(numpy.isnan(self.basal_moulin_input): 117 self.basal_moulin_input=numpy.zeros(md.mesh.numberofvertices,1) 118 print" no hydrology.basal_moulin_input specified: values set as zero" 119 120 return self 121 122 # }}} 110 123 def checkconsistency(self,md,solution,analyses): #{{{ 111 124 … … 127 140 if self.transfer_flag==1: 128 141 md = checkfield(md,'fieldname','hydrology.leakage_factor','>',0,'numel',1) 129 142 143 md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'forcing',1) 130 144 md = checkfield(md,'fieldname','hydrology.spcsediment_head','forcing',1) 131 145 md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0,'numel',1) … … 159 173 WriteData(fid,'object',self,'fieldname','leakage_factor','format','Double') 160 174 175 WriteData(fid,'object',self,'fieldname','basal_moulin_input','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1) 161 176 WriteData(fid,'object',self,'fieldname','spcsediment_head','format','DoubleMat','mattype',1,'forcinglength',md.mesh.numberofvertices+1) 162 177 WriteData(fid,'object',self,'fieldname','sediment_compressibility','format','Double')
Note:
See TracChangeset
for help on using the changeset viewer.