Index: /issm/trunk-jpl/examples/Greenland/CheatSheet.m
===================================================================
--- /issm/trunk-jpl/examples/Greenland/CheatSheet.m	(revision 18263)
+++ /issm/trunk-jpl/examples/Greenland/CheatSheet.m	(revision 18264)
@@ -9,5 +9,5 @@
 
 	%Save surface mass balance, by looping through 200 years, or 1000 steps
-	surfmb=[]; for i=1:100; surfmb=[surfmb ...
+	surfmb=[]; for i=1:1000; surfmb=[surfmb ...
 		md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
 
@@ -19,5 +19,5 @@
 
 	%Save velocity by looping through 1000 steps
-	vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end
+	vel=[]; for i=1:1000; vel=[vel md.results.TransientSolution(i).Vel]; end
 
 	%Plot velocity time series in second subplot
@@ -28,5 +28,5 @@
 
 	%Save Ice Volume by looping through 1000 steps
-	volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end
+	volume=[]; for i=1:1000; volume=[volume md.results.TransientSolution(i).IceVolume]; end
 
 	%Plot volume time series in third subplot
@@ -35,4 +35,3 @@
 	%Title this plot Mean Velocity and add an x label of years
 	title('Ice Volume'); xlabel('years');
-
 end
Index: /issm/trunk-jpl/examples/Greenland/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Greenland/runme.m	(revision 18263)
+++ /issm/trunk-jpl/examples/Greenland/runme.m	(revision 18264)
@@ -146,4 +146,5 @@
 if any(steps==6)
 	disp('   Step 6: Extract Box SMB');
+	md = loadmodel('./Models/Greenland.Transient');
 
 	ncbox='../Data/Box_Greenland_SMB_monthly_1840-2012_5km_ver20140421.nc';
@@ -164,5 +165,5 @@
 	%Area of grid for 5km box
 	area_of_grid=5000*5000;
-	totalsmb=reshape(sum(sum(smbbox,1),2),length(t),1);
+	totalsmb=reshape(sum(sum(smbbox/1000,1),2),length(t),1)*area_of_grid;
 
 	%save surface mass balance mat dataset
@@ -171,7 +172,7 @@
 
 	%plot a time series of total SMB
-	plot(t,totalsmb); title('Total Surface mass balance, m^3');
-
-	clear smbbox smbmean totalsmb x1 y1 t
+	plot(t,totalsmb/1e9); title('Total Surface mass balance, Gt'); xlabel('year'); ylabel('Gt/yr');
+
+	clear smbbox
 
 end
@@ -211,36 +212,34 @@
 
 	%load historic transient model
-	md = loadmodel('./Models/Greenland.HistoricTransient');
 
 	%Create Line Plots of relaxation run.  Create a figure.
-	figure
+
 
 	%Save surface mass balance, by looping through 200 years, or 1000 steps
-	surfmb=[]; for i=1:100; surfmb=[surfmb ...
-		md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
+
 
 	%Plot surface mass balance time series in first subplot
-	subplot(3,1,1); plot([0.2:0.2:200],mean(surfmb)); 
+
 	
 	%Title this plot Mean surface mass balance
-	title('Mean Surface mass balance');
+
 
 	%Save velocity by looping through 1000 steps
-	vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end
+
 
 	%Plot velocity time series in second subplot
-	subplot(3,1,2); plot([0.2:0.2:200],mean(vel)); 
+
 	
 	%Title this plot Mean Velocity
-	title('Mean Velocity');
+
 
 	%Save Ice Volume by looping through 1000 steps
-	volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end
+
 
 	%Plot volume time series in third subplot
-	subplot(3,1,3); plot([0.2:0.2:200],volume); 
+
 	
 	%Title this plot Mean Velocity and add an x label of years
-	title('Ice Volume'); xlabel('years');
+
 end