Index: /issm/trunk-jpl/examples/SlrGRACE/runme.m
===================================================================
--- /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 25937)
+++ /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 25938)
@@ -2,5 +2,5 @@
 addpath('../Data','../Functions');
 
-steps=[6];
+steps=[7];
 
 if any(steps==1) %{{{
@@ -199,4 +199,5 @@
 	md.transient.isstressbalance=0;
 	md.transient.isthermal=0;
+	md.transient.isgia=1; 
 	md.transient.isslr=1;
 
@@ -210,4 +211,6 @@
 	md.verbose=verbose('111111111');
 
+	md.solidearth.requested_outputs = {'Sealevel','SealevelRSL'};
+	
 	md=solve(md,'tr');
 
@@ -219,10 +222,10 @@
 	md = loadmodel('./Models/SlrGRACE_Transient');
 
-	time = md.slr.deltathickness(end,:);
+	time = md.solidearth.surfaceload.icethicknesschange(end,:);
 
 	for tt=1:length(time)
-		gmsl(tt) = md.results.TransientSolution(tt).SealevelEustatic*1000;	% GMSL rate mm/yr
-		sol1(:,tt) = md.slr.deltathickness(1:end-1,tt)*100;					% ice equivalent height [cm/yr]
-		sol2(:,tt) = md.results.TransientSolution(tt+1).Sealevel*1000;		% mm/yr
+		gmsl(tt) = md.results.TransientSolution(tt).Bslr*1000; % GMSL rate mm/yr
+		sol1(:,tt) = md.solidearth.surfaceload.icethicknesschange(1:end-1,tt)*100;	% ice equivalent height [cm/yr]
+		sol2(:,tt) = md.results.TransientSolution(tt+1).SealevelRSL*1000;	% mm/yr
 	end
 	sol_name = {'Change in water equivalent height [cm]', 'Relative sea-level [mm]'};