Index: /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp =================================================================== --- /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (revision 25957) +++ /issm/trunk-jpl/src/c/analyses/MasstransportAnalysis.cpp (revision 25958) @@ -795,7 +795,8 @@ void MasstransportAnalysis::InputUpdateFromSolution(IssmDouble* solution,Element* element){/*{{{*/ + /*Only update if on base*/ if(!element->IsOnBase()) return; - + /*Fetch dof list and allocate solution vector*/ int *doflist = NULL; @@ -805,4 +806,7 @@ IssmDouble* newthickness = xNew(numnodes); IssmDouble* thicknessresidual = xNew(numnodes); + + /*recover time step:*/ + IssmDouble dt = element->FindParam(TimesteppingTimeStepEnum); /*Use the dof list to index into the solution vector: */ @@ -840,5 +844,5 @@ IssmDouble* oldthickness = xNew(numvertices); IssmDouble* cumdeltathickness = xNew(numvertices); - IssmDouble* deltathickness = xNew(numvertices); + IssmDouble* icethicknessrate = xNew(numvertices); IssmDouble* newbase = xNew(numvertices); IssmDouble* bed = xNew(numvertices); @@ -866,5 +870,5 @@ for(int i=0;iAddBasalInput(BaseEnum,newbase,P1Enum); element->AddBasalInput(SealevelchangeCumDeltathicknessEnum,cumdeltathickness,P1Enum); - element->AddBasalInput(SurfaceloadIceThicknessChangeEnum,deltathickness,P1Enum); + element->AddBasalInput(SurfaceloadIceThicknessRateEnum,icethicknessrate,P1Enum); /*Free ressources:*/ @@ -911,5 +915,5 @@ xDelete(newsurface); xDelete(oldthickness); - xDelete(deltathickness); + xDelete(icethicknessrate); xDelete(oldbase); xDelete(oldsurface); Index: /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp =================================================================== --- /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp (revision 25957) +++ /issm/trunk-jpl/src/c/analyses/SealevelchangeAnalysis.cpp (revision 25958) @@ -40,8 +40,9 @@ iomodel->FetchDataToInput(inputs,elements,"md.mask.ice_levelset",MaskIceLevelsetEnum); iomodel->FetchData(&geodetic,"md.solidearth.settings.computesealevelchange"); - iomodel->FetchDataToInput(inputs,elements,"md.solidearth.surfaceload.icethicknesschange",SurfaceloadIceThicknessChangeEnum); iomodel->FetchDataToInput(inputs,elements,"md.solidearth.initialsealevel",SealevelEnum,0); iomodel->FetchDataToInput(inputs,elements,"md.geometry.bed",BedEnum); - iomodel->FetchDataToInput(inputs,elements,"md.solidearth.surfaceload.waterheightchange",SurfaceloadWaterHeightChangeEnum); + iomodel->FetchDataToInput(inputs,elements,"md.solidearth.surfaceload.icethicknesschange",SurfaceloadIceThicknessRateEnum); + iomodel->FetchDataToInput(inputs,elements,"md.solidearth.surfaceload.waterheightchange",SurfaceloadWaterHeightRateEnum); + iomodel->FetchDataToInput(inputs,elements,"md.solidearth.surfaceload.other",SurfaceloadOtherRateEnum); /*dynamic sea level: */ Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp =================================================================== --- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 25957) +++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp (revision 25958) @@ -5523,4 +5523,6 @@ /*}}}*/ void Tria::SealevelchangeMomentOfInertia(IssmDouble* dI_list,IssmDouble* Sg_old, SealevelMasks* masks){/*{{{*/ + + /*early return if we are not on an ice cap OR ocean:*/ if(!masks->isiceonly[this->lid] && !masks->isoceanin[this->lid]){ @@ -5599,17 +5601,20 @@ } else if(masks->isiceonly[this->lid]){ - IssmDouble rho_ice, I; - - /*recover material parameters: */ + IssmDouble rho_ice, dIdt, dt; + + + /*recover parameters: */ rho_ice=FindParam(MaterialsRhoIceEnum); + dt=FindParam(TimesteppingTimeStepEnum); /*Compute ice thickness change: */ - Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessChangeEnum); + Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessRateEnum); if (!deltathickness_input)_error_("delta thickness input needed to compute sea level change!"); - deltathickness_input->GetInputAverage(&I); - - dI_list[0] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/planetarea; - dI_list[1] = -4*PI*(rho_ice*I*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/planetarea; - dI_list[2] = +4*PI*(rho_ice*I*area)*pow(re,4)*(1-pow(sin(late),2))/planetarea; + deltathickness_input->GetInputAverage(&dIdt); + + + dI_list[0] = -4*PI*(rho_ice*dIdt*dt*area)*pow(re,4)*(sin(late)*cos(late)*cos(longe))/planetarea; + dI_list[1] = -4*PI*(rho_ice*dIdt*dt*area)*pow(re,4)*(sin(late)*cos(late)*sin(longe))/planetarea; + dI_list[2] = +4*PI*(rho_ice*dIdt*dt*area)*pow(re,4)*(1-pow(sin(late),2))/planetarea; } @@ -5819,9 +5824,10 @@ IssmDouble area; IssmDouble phi=1.0; //WARNING: do not touch this, default is entire elemnt contributes barystatic - IssmDouble I; //change in ice thickness or water level(Farrel and Clarke, Equ. 4) + IssmDouble Idot; //ice thickness rate (Farrel and Clarke, Equ. 4) bool notfullygrounded=false; bool scaleoceanarea= false; bool computerigid= false; int glfraction=1; + IssmDouble dt=0; /*output: */ @@ -5868,4 +5874,5 @@ rho_ice=FindParam(MaterialsRhoIceEnum); rho_water=FindParam(MaterialsRhoSeawaterEnum); + dt=FindParam(TimesteppingTimeStepEnum); this->parameters->FindParam(&computerigid,SolidearthSettingsRigidEnum); this->parameters->FindParam(&scaleoceanarea,SolidearthSettingsOceanAreaScalingEnum); @@ -5892,9 +5899,9 @@ /*Retrieve ice thickness at vertices: */ - Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessChangeEnum); + Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessRateEnum); if (!deltathickness_input)_error_("delta thickness input needed to compute sea level change!"); /*/Average ice thickness over grounded area of the element only: {{{*/ - if(!notfullygrounded)deltathickness_input->GetInputAverage(&I); + if(!notfullygrounded)deltathickness_input->GetInputAverage(&Idot); else{ IssmDouble total_weight=0; @@ -5909,12 +5916,12 @@ /* Start looping on the number of gaussian points and average over these gaussian points: */ total_weight=0; - I=0; + Idot=0; while(gauss->next()){ - IssmDouble Ig=0; - deltathickness_input->GetInputValue(&Ig,gauss); - I+=Ig*gauss->weight; + IssmDouble Idotg=0; + deltathickness_input->GetInputValue(&Idotg,gauss); + Idot+=Idotg*gauss->weight; total_weight+=gauss->weight; } - I=I/total_weight; + Idot=Idot/total_weight; delete gauss; } @@ -5924,13 +5931,13 @@ _assert_(oceanarea>0.); if(scaleoceanarea) oceanarea=3.619e+14; // use true ocean area, m^2 - bslcice = rho_ice*area*phi*I/(oceanarea*rho_water); + bslcice = rho_ice*area*phi*Idot*dt/(oceanarea*rho_water); _assert_(!xIsNan(bslcice)); if(computerigid){ /*convert from m to kg/m^2:*/ - I=I*rho_ice*phi; + Idot=Idot*rho_ice*phi; /*convolve:*/ - for(int i=0;iparameters->FindParam(&computeelastic,SolidearthSettingsElasticEnum); this->parameters->FindParam(&scaleoceanarea,SolidearthSettingsOceanAreaScalingEnum); @@ -5993,20 +6002,20 @@ /*Retrieve water height at vertices: */ - Input* deltathickness_input=this->GetInput(SurfaceloadWaterHeightChangeEnum); - if (!deltathickness_input)_error_("SurfaceloadWaterHeightChangeEnum input needed to compute sea level change!"); - deltathickness_input->GetInputAverage(&W); + Input* deltathickness_input=this->GetInput(SurfaceloadWaterHeightRateEnum); + if (!deltathickness_input)_error_("SurfaceloadWaterHeightRateEnum input needed to compute sea level change!"); + deltathickness_input->GetInputAverage(&Wdot); /*Compute barystatic component:*/ _assert_(oceanarea>0.); if(scaleoceanarea) oceanarea=3.619e+14; // use true ocean area, m^2 - bslchydro = rho_freshwater*area*phi*W/(oceanarea*rho_water); + bslchydro = rho_freshwater*area*phi*Wdot*dt/(oceanarea*rho_water); _assert_(!xIsNan(bslchydro)); if(computeelastic){ /*convert from m to kg/m^2:*/ - W=W*rho_freshwater*phi; + Wdot=Wdot*rho_freshwater*phi; /*convolve:*/ - for(int i=0;iisiceonly[this->lid]){ + IssmDouble dt=0; + dt=FindParam(TimesteppingTimeStepEnum); /*Compute ice thickness change: */ - Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessChangeEnum); + Input* deltathickness_input=this->GetInput(SurfaceloadIceThicknessRateEnum); if (!deltathickness_input)_error_("delta thickness input needed to compute sea level change!"); - deltathickness_input->GetInputAverage(&I); + deltathickness_input->GetInputAverage(&Idot); /*convert to kg/m^2*/ - I=I*rho_ice; + Idot=Idot*rho_ice; for(int i=0;i((strlen(field)+1)); xMemCpy(fieldname,field,(strlen(field)+1)); } - else if(strcmp(string_in,"SurfaceloadIceThicknessChange")==0){ + else if(strcmp(string_in,"SurfaceloadIceThicknessRate")==0){ const char* field = "md.solidearth.surfaceload.icethicknesschange"; - input_enum = SurfaceloadIceThicknessChangeEnum; + input_enum = SurfaceloadIceThicknessRateEnum; fieldname=xNew((strlen(field)+1)); xMemCpy(fieldname,field,(strlen(field)+1)); } @@ -178,7 +178,7 @@ fieldname=xNew((strlen(field)+1)); xMemCpy(fieldname,field,(strlen(field)+1)); } - else if(strcmp(string_in,"SurfaceloadWaterHeightChange")==0){ + else if(strcmp(string_in,"SurfaceloadWaterHeightRate")==0){ const char* field = "md.solidearth.surfaceload.waterheightchange"; - input_enum = SurfaceloadWaterHeightChangeEnum; + input_enum = SurfaceloadWaterHeightRateEnum; fieldname=xNew((strlen(field)+1)); xMemCpy(fieldname,field,(strlen(field)+1)); } Index: /issm/trunk-jpl/src/m/classes/mmeadditionalsolidearthsolution.m =================================================================== --- /issm/trunk-jpl/src/m/classes/mmeadditionalsolidearthsolution.m (revision 25957) +++ /issm/trunk-jpl/src/m/classes/mmeadditionalsolidearthsolution.m (revision 25958) @@ -69,13 +69,40 @@ function marshall(self,prefix,md,fid) % {{{ WriteData(fid,prefix,'object',self,'data',3,'name','md.solidearth.external.nature','format','Integer'); %code 3 for mmeadditionalsolidearthsolution class - WriteData(fid,prefix,'name','md.solidearth.external.nummodels','data',length(self.displacementeast),'format','Integer'); + nummodels=length(self.displacementeast); + WriteData(fid,prefix,'name','md.solidearth.external.nummodels','data',nummodels,'format','Integer'); WriteData(fid,prefix,'object',self,'fieldname','modelid','format','Double'); - WriteData(fid,prefix,'object',self,'fieldname','displacementeast','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','displacementup','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','displacementnorth','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','geoid','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','barystaticsealevel','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); + + %transform our cell array of time series into cell array of time series of rates + for i=1:nummodels, + displacementeast=self.displacementeast{i}; + displacementnorth=self.displacementnorth{i}; + displacementup=self.displacementup{i}; + geoid=self.geoid{i}; + barystaticsealevel=self.barystaticsealevel{i}; + + time=displacementeast(end,:); + dt=diff(time,1,2); + + displacementeast_rate=diff(displacementeast(1:end-1,:),1,2)./dt; + displacementnorth_rate=diff(displacementnorth(1:end-1,:),1,2)./dt; + displacementup_rate=diff(displacementup(1:end-1,:),1,2)./dt; + geoid_rate=diff(geoid(1:end-1,:),1,2)./dt; + barystaticsealevel_rate=diff(barystaticsealevel(1:end-1,:),1,2)./dt; + + self.displacementeast{i}=displacementeast_rate; + self.displacementnorth{i}=displacementnorth_rate; + self.displacementup{i}=displacementup_rate; + self.geoid{i}=geoid_rate; + self.barystaticsealevel{i}=barystaticsealevel_rate; + end + + WriteData(fid,prefix,'object',self,'fieldname','displacementeast','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','displacementup','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','displacementnorth','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','geoid','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','barystaticsealevel','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); end % }}} + function savemodeljs(self,fid,modelname) % {{{ error('mmeadditionalsolidearthsolution error message: not implemented yet'); Index: /issm/trunk-jpl/src/m/classes/mmeofflinesolidearthsolution.m =================================================================== --- /issm/trunk-jpl/src/m/classes/mmeofflinesolidearthsolution.m (revision 25957) +++ /issm/trunk-jpl/src/m/classes/mmeofflinesolidearthsolution.m (revision 25958) @@ -70,10 +70,36 @@ WriteData(fid,prefix,'object',self,'data',4,'name','md.solidearth.external.nature','format','Integer'); %code 4 for mmeofflinesolidearthsolution class WriteData(fid,prefix,'object',self,'fieldname','modelid','format','Double'); - WriteData(fid,prefix,'name','md.solidearth.external.nummodels','data',length(self.displacementeast),'format','Integer'); - WriteData(fid,prefix,'object',self,'fieldname','displacementeast','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','displacementup','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','displacementnorth','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','geoid','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); - WriteData(fid,prefix,'object',self,'fieldname','barystaticsealevel','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1e-3); + nummodels=length(self.displacementeast); + WriteData(fid,prefix,'name','md.solidearth.external.nummodels','data',nummodels,'format','Integer'); + + %transform our cell array of time series into cell array of time series of rates + for i=1:nummodels, + displacementeast=self.displacementeast{i}; + displacementnorth=self.displacementnorth{i}; + displacementup=self.displacementup{i}; + geoid=self.geoid{i}; + barystaticsealevel=self.barystaticsealevel{i}; + + time=displacementeast(end,:); + dt=diff(time,1,2); + + displacementeast_rate=diff(displacementeast(1:end-1,:),1,2)./dt; + displacementnorth_rate=diff(displacementnorth(1:end-1,:),1,2)./dt; + displacementup_rate=diff(displacementup(1:end-1,:),1,2)./dt; + geoid_rate=diff(geoid(1:end-1,:),1,2)./dt; + barystaticsealevel_rate=diff(barystaticsealevel(1:end-1,:),1,2)./dt; + + self.displacementeast{i}=displacementeast_rate; + self.displacementnorth{i}=displacementnorth_rate; + self.displacementup{i}=displacementup_rate; + self.geoid{i}=geoid_rate; + self.barystaticsealevel{i}=barystaticsealevel_rate; + end + + WriteData(fid,prefix,'object',self,'fieldname','displacementeast','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','displacementup','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','displacementnorth','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','geoid','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); + WriteData(fid,prefix,'object',self,'fieldname','barystaticsealevel','format','MatArray','timeserieslength',md.mesh.numberofvertices+1,'yts',md.constants.yts,'scale',1/yts); end % }}} Index: /issm/trunk-jpl/src/m/classes/surfaceload.m =================================================================== --- /issm/trunk-jpl/src/m/classes/surfaceload.m (revision 25957) +++ /issm/trunk-jpl/src/m/classes/surfaceload.m (revision 25958) @@ -65,26 +65,89 @@ end % }}} function marshall(self,prefix,md,fid) % {{{ - + + %deal with ice thickness change: {{{ if isempty(self.icethicknesschange), self.icethicknesschange=zeros(md.mesh.numberofelements+1,1); end + + if isa(self.icethicknesschange,'cell'), + %transform our cell array of time series into cell array of time series of rates + nummodels=length(self.icethicknesschange); + for i=1:nummodels, + icethicknesschange=self.icethicknesschange{i}; + time=icethicknesschange(end,:); + dt=diff(time,1,2); + icethicknesschange_rate=diff(icethicknesschange(1:end-1,:),1,2)./dt; + self.icethicknesschange{i}=icethicknesschange_rate; + end + WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',... + 'format','MatArray','timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); + else + icethicknesschange=self.icethicknesschange; + time=icethicknesschange(end,:); + dt=diff(time,1,2); + icethicknesschange_rate=diff(icethicknesschange(1:end-1,:),1,2)./dt; + self.icethicknesschange=icethicknesschange_rate; + + WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',... + 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); + end + %}}} + %deal with water height change: {{{ if isempty(self.waterheightchange), self.waterheightchange=zeros(md.mesh.numberofelements+1,1); end + + if isa(self.waterheightchange,'cell'), + %transform our cell array of time series into cell array of time series of rates + nummodels=length(self.waterheightchange); + for i=1:nummodels, + waterheightchange=self.waterheightchange{i}; + time=waterheightchange(end,:); + dt=diff(time,1,2); + waterheightchange_rate=diff(waterheightchange(1:end-1,:),1,2)./dt; + self.waterheightchange{i}=waterheightchange_rate; + end + WriteData(fid,prefix,'object',self,'fieldname','waterheightchange','name','md.solidearth.surfaceload.waterheightchange',... + 'format','MatArray','timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); + else + waterheightchange=self.waterheightchange; + time=waterheightchange(end,:); + dt=diff(time,1,2); + waterheightchange_rate=diff(waterheightchange(1:end-1,:),1,2)./dt; + self.waterheightchange=waterheightchange_rate; + + WriteData(fid,prefix,'object',self,'fieldname','waterheightchange','name','md.solidearth.surfaceload.waterheightchange',... + 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); + end + %}}} + %deal with other: {{{ if isempty(self.other), self.other=zeros(md.mesh.numberofelements+1,1); end - if isa(self.icethicknesschange,'cell'), - WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',... - 'format','MatArray','timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); + + if isa(self.otherchange,'cell'), + %transform our cell array of time series into cell array of time series of rates + nummodels=length(self.otherchange); + for i=1:nummodels, + otherchange=self.otherchange{i}; + time=otherchange(end,:); + dt=diff(time,1,2); + otherchange_rate=diff(otherchange(1:end-1,:),1,2)./dt; + self.otherchange{i}=otherchange_rate; + end + WriteData(fid,prefix,'object',self,'fieldname','otherchange','name','md.solidearth.surfaceload.otherchange',... + 'format','MatArray','timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); else - WriteData(fid,prefix,'object',self,'fieldname','icethicknesschange','name','md.solidearth.surfaceload.icethicknesschange',... - 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); + otherchange=self.otherchange; + time=otherchange(end,:); + dt=diff(time,1,2); + otherchange_rate=diff(otherchange(1:end-1,:),1,2)./dt; + self.otherchange=otherchange_rate; + + WriteData(fid,prefix,'object',self,'fieldname','otherchange','name','md.solidearth.surfaceload.otherchange',... + 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts,'scale',1/yts); end - - WriteData(fid,prefix,'object',self,'fieldname','waterheightchange','name','md.solidearth.surfaceload.waterheightchange',... - 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); - WriteData(fid,prefix,'object',self,'fieldname','other','name','md.solidearth.surfaceload.other',... - 'format','DoubleMat','mattype',2,'timeserieslength',md.mesh.numberofelements+1,'yts',md.constants.yts); + %}}} end % }}}