Index: /issm/trunk-jpl/examples/SlrGRACE/runme.m
===================================================================
--- /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 27028)
+++ /issm/trunk-jpl/examples/SlrGRACE/runme.m	(revision 27029)
@@ -1,11 +1,11 @@
-clear all;
+%clear all;
 addpath('../Data','../Functions');
 
-steps=[1:7];
+%steps=[1:7];
 
 if any(steps==1) % {{{ 
 	disp('   Step 1: Global mesh creation');
 
-	numrefine=5;
+	numrefine=1;
 	resolution=150*1e3;			% inital resolution [m]
 	radius = 6.371012*10^6;		% mean radius of Earth, m
@@ -50,5 +50,6 @@
 	save ./Models/SlrGRACE_Mesh md;
 
-	plotmodel (md,'data',md.mask.ocean_levelset,'edgecolor','k');
+	plotmodel (md,'data',md.mask.ocean_levelset,'edgecolor','k','view',[45 45]);
+
 end % }}} 
 if any(steps==2) % {{{ 
@@ -70,5 +71,7 @@
 	md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
-	md.masstransport.spcthickness = [md.geometry.thickness + ice_load; 0];
+	md.masstransport.spcthickness = repmat(md.geometry.thickness,1,2); 
+	md.masstransport.spcthickness(:,end) = md.masstransport.spcthickness(:,end) + ice_load; 
+	md.masstransport.spcthickness(end+1,:) = [0 1]; % dummy time. 
 
 	md.smb.mass_balance=zeros(md.mesh.numberofvertices,1); 
@@ -140,5 +143,5 @@
 	md = loadmodel('./Models/SlrGRACE_Solution');
 
-	sol1 = (md.masstransport.spcthickness(1:end-1)-md.geometry.thickness)*100; % [cm]
+	sol1 = (md.masstransport.spcthickness(1:end-1,2)-md.masstransport.spcthickness(1:end-1,1))*100; % [cm]
 	sol2 = (md.results.TransientSolution.Sealevel-md.results.TransientSolution.Bed)*1000;	% [mm]
 
@@ -208,6 +211,7 @@
 	% masstransport evalulates diff between the successive times, so we should cumsum. 
 	num_time = size(ice_load,2); 
-	md.masstransport.spcthickness = [repmat(md.geometry.thickness,1,size(ice_load,2)) + ice_load];
-	md.masstransport.spcthickness =[ md.masstransport.spcthickness; 0:num_time-1]; 
+	md.masstransport.spcthickness = repmat(md.geometry.thickness,1,num_time+1); 
+	md.masstransport.spcthickness(:,2:end) = md.masstransport.spcthickness(:,2:end) + ice_load;
+	md.masstransport.spcthickness(end+1,:) = 0:num_time; % dummy time. 
 
 	%Physics 
@@ -241,13 +245,13 @@
 	md = loadmodel('./Models/SlrGRACE_Transient');
 
-	time = md.masstransport.spcthickness(end,:);
+	time = md.masstransport.spcthickness(end,2:end);
 
 	for tt=1:length(time)
 		gmsl(tt) = md.results.TransientSolution(tt).CumBslc*1000; % [mm]
-		sol1(:,tt) = (md.masstransport.spcthickness(1:end-1,tt)-md.geometry.thickness)*100; % [cm]
+		sol1(:,tt) = (md.masstransport.spcthickness(1:end-1,tt+1)-md.geometry.thickness)*100; % [cm]
 		sol2(:,tt) = (md.results.TransientSolution(tt).Sealevel-md.results.TransientSolution(tt).Bed)*1000;	% [mm]
 	end
 	sol_name = {'Change in water equivalent height [cm]', 'Relative sea level [mm]'};
-	movie_name = {'Movie_dH.mp4','Movie_slr.mp4'};
+	movie_name = {'Movie_dH','Movie_slr'};
 	
 	res = 1.0;
@@ -273,5 +277,4 @@
 		vidObj = VideoWriter(movie_name{kk});
 		vidObj.FrameRate=2; % frames per second
-		vifObj.FileFormat='mp4';
 		open(vidObj);
 
@@ -311,5 +314,4 @@
 		close(vidObj);
 	end
-	!open *.mp4;
 
 	% plot GMSL time series