Changeset 15143
- Timestamp:
- 05/29/13 14:50:26 (12 years ago)
- Location:
- issm/trunk-jpl/src/m/classes
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/gia.m
r15131 r15143 8 8 mantle_viscosity = NaN; 9 9 lithosphere_thickness = NaN; 10 output_rates = 0;11 10 cross_section_shape = 0; 12 11 end … … 21 20 end % }}} 22 21 function obj = setdefaultparameters(obj) % {{{ 23 obj.output_rates=1; %don't output rates by default (see irate in GiaDeflectionCorex)24 22 obj.cross_section_shape=1; %square as default (see iedge in GiaDeflectionCorex) 25 23 end % }}} … … 29 27 md = checkfield(md,'gia.mantle_viscosity','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0); 30 28 md = checkfield(md,'gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices 1],'>',0); 31 md = checkfield(md,'gia.output_rates','numel',[1],'values',[1,2]);32 29 md = checkfield(md,'gia.cross_section_shape','numel',[1],'values',[1,2]); 33 30 … … 53 50 fielddisplay(obj,'mantle_viscosity','mantle viscosity[Pa s]'); 54 51 fielddisplay(obj,'lithosphere_thickness','lithosphere thickness (km)'); 55 fielddisplay(obj,'output_rates','1: fetch w solution (m) (default). 2: fetch dw/dt (mm/yr). See irate in GiaDeflectionCore');56 52 fielddisplay(obj,'cross_section_shape','1: square-edged (default). 2: elliptical. See iedge in GiaDeflectionCore'); 57 53 … … 60 56 WriteData(fid,'object',obj,'fieldname','mantle_viscosity','format','DoubleMat','mattype',1); 61 57 WriteData(fid,'object',obj,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1,'scale',10^3); %from km to m 62 WriteData(fid,'object',obj,'fieldname','output_rates','format','Integer');63 58 WriteData(fid,'object',obj,'fieldname','cross_section_shape','format','Integer'); 64 59 end % }}} -
issm/trunk-jpl/src/m/classes/gia.py
r15131 r15143 15 15 self.mantle_viscosity = float('NaN'); 16 16 self.lithosphere_thickness = float('NaN'); 17 self.output_rates = 0;18 17 self.cross_section_shape = 0; 19 18 … … 28 27 string="%s\n%s"%(string,fielddisplay(self,'mantle_viscosity','mantle viscosity constraints (NaN means no constraint) (Pa s)')) 29 28 string="%s\n%s"%(string,fielddisplay(self,'lithosphere_thickness','lithosphere thickness constraints (NaN means no constraint) (m)')) 30 string="%s\n%s"%(string,fielddisplay(self,'output_rates','1: deflection (w) (m), 2: rate of deflection (dw/dt) (mm/yr)'))31 29 string="%s\n%s"%(string,fielddisplay(self,'cross_section_shape',"1: square-edged, 2: elliptical-edged surface")) 32 30 return string … … 34 32 def setdefaultparameters(self): # {{{ 35 33 36 self.output_rates=1;37 34 self.cross_section_shape=1; 38 35 … … 47 44 md = checkfield(md,'gia.mantle_viscosity','NaN',1,'size',[md.mesh.numberofvertices,1],'>',0) 48 45 md = checkfield(md,'gia.lithosphere_thickness','NaN',1,'size',[md.mesh.numberofvertices,1],'>',0) 49 md = checkfield(md,'gia.output_rates','numel',[1],'values',[1,2])50 46 md = checkfield(md,'gia.cross_section_shape','numel',[1],'values',[1,2]) 51 47 … … 59 55 WriteData(fid,'object',self,'fieldname','mantle_viscosity','format','DoubleMat','mattype',1); 60 56 WriteData(fid,'object',self,'fieldname','lithosphere_thickness','format','DoubleMat','mattype',1,'scale',10.**3.); 61 WriteData(fid,'object',self,'fieldname','output_rates','format','Integer');62 57 WriteData(fid,'object',self,'fieldname','cross_section_shape','format','Integer'); 63 58 # }}}
Note:
See TracChangeset
for help on using the changeset viewer.