Changeset 26274 for issm/trunk-jpl/src/m/classes/solidearthsolution.m
- Timestamp:
- 05/17/21 14:08:16 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/classes/solidearthsolution.m
r26229 r26274 10 10 displacementup=[]; 11 11 geoid=[]; 12 barystaticsealevel=[];13 12 end 14 13 methods … … 27 26 self.displacementup=[]; 28 27 self.geoid=[]; 29 self.barystaticsealevel=[];30 28 31 29 end % }}} … … 36 34 md = checkfield(md,'fieldname','solidearth.external.displacementup','Inf',1,'timeseries',1); 37 35 md = checkfield(md,'fieldname','solidearth.external.geoid','Inf',1,'timeseries',1); 38 md = checkfield(md,'fieldname','solidearth.external.barystaticsealevel','Inf',1,'timeseries',1);39 36 40 37 end % }}} … … 45 42 fielddisplay(self,'displacementup','solid-Earth bedrock uplift time series (m)'); 46 43 fielddisplay(self,'geoid','solid-Earth geoid time series (m)'); 47 fielddisplay(self,'barystaticsealevel','barystatic absolute sea-level time series (m)');48 44 49 45 end % }}} … … 58 54 displacementup_rate=[self.displacementup;0]; 59 55 geoid_rate=[self.geoid;0]; 60 barystaticsealevel_rate=[self.barystaticsealevel;0];61 56 else 62 57 time=self.displacementeast(end,:); … … 70 65 geoid_rate=diff(self.geoid(1:end-1,:),1,2)./dt; 71 66 geoid_rate(end+1,:)=time(1:end-1); 72 barystaticsealevel_rate=diff(self.barystaticsealevel(1:end-1,:),1,2)./dt;73 barystaticsealevel_rate(end+1,:)=time(1:end-1);74 67 end 75 68 WriteData(fid,prefix,'object',self,'fieldname','displacementeast','data',displacementeast_rate,'format','DoubleMat','name', 'md.solidearth.external.displacementeast','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); … … 77 70 WriteData(fid,prefix,'object',self,'fieldname','displacementnorth','data',displacementnorth_rate,'format','DoubleMat','name', 'md.solidearth.external.displacementnorth','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); 78 71 WriteData(fid,prefix,'object',self,'fieldname','geoid','data',geoid_rate,'format','DoubleMat','name', 'md.solidearth.external.geoid','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts); 79 WriteData(fid,prefix,'object',self,'fieldname','barystaticsealevel','data',barystaticsealevel_rate,'format','DoubleMat','name', 'md.solidearth.external.barystaticsealevel','mattype',1,'scale',1./yts,'timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts);80 72 81 73 end % }}} … … 85 77 writejs1Darray(fid,[modelname '.solidearth.external.displacementup'],self.displacementup); 86 78 writejs1Darray(fid,[modelname '.solidearth.external.geoid'],self.geoid); 87 writejs1Darray(fid,[modelname '.solidearth.external.barystaticsealevel'],self.barystaticsealevel);88 79 end % }}} 89 80 function self = extrude(self,md) % {{{
Note:
See TracChangeset
for help on using the changeset viewer.