Changeset 26403

Timestamp:

08/20/21 14:44:54 (4 years ago)

Author:

schlegel

Message:

CHG: update to newest sea level class

File:

• issm/trunk-jpl/examples/SlrGRACE/runme.m (modified) (4 diffs)

issm/trunk-jpl/examples/SlrGRACE/runme.m

@@ -89 +89 @@
         md.solidearth.settings.reltol=NaN;
         md.solidearth.settings.abstol=1e-3;
-        md.solidearth.settings.computesealevelchange=1;
+        md.solidearth.settings.sealevelloading=1;
         md.solidearth.settings.isgrd=1;
         md.solidearth.settings.grdmodel=1;
@@ -115 +115 @@
         md = loadmodel('./Models/SlrGRACE_Parameterization');
 
-        md.solidearth.settings.rigid=1;
+        md.solidearth.settings.viscous=0;
         md.solidearth.settings.elastic=1;
         md.solidearth.settings.rotation=1;
@@ -207 +207 @@
         % masstransport evalulates diff between the successive times, so we should cumsum. 
         num_time = size(ice_load,2); 
-        md.masstransport.spcthickness = [md.geometry.thickness + ice_load];
-        md.masstransport.spcthickness(md.mesh.numberofvertices+1,:)=[0:num_time-1]; 
+        md.masstransport.spcthickness = [repmat(md.geometry.thickness,1,size(ice_load,2)) + ice_load];
+        md.masstransport.spcthickness =[ md.masstransport.spcthickness; 0:num_time-1]; 
 
         %Physics 
@@ -217 +217 @@
         md.transient.isslc=1;
 
-        md.solidearth.settings.rigid=1;
+        md.solidearth.settings.viscous=0;
         md.solidearth.settings.elastic=1;
         md.solidearth.settings.rotation=1;