Index: /issm/trunk-jpl/examples/SlrGRACE/runme.m
===================================================================
--- /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 26402)
+++ /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 26403)
@@ -89,5 +89,5 @@
 	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,5 +115,5 @@
 	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,6 +207,6 @@
 	% 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,5 +217,5 @@
 	md.transient.isslc=1;
 
-	md.solidearth.settings.rigid=1;
+	md.solidearth.settings.viscous=0;
 	md.solidearth.settings.elastic=1;
 	md.solidearth.settings.rotation=1;