Index: /issm/trunk-jpl/examples/EsaGRACE/runme.m
===================================================================
--- /issm/trunk-jpl/examples/EsaGRACE/runme.m	(revision 22814)
+++ /issm/trunk-jpl/examples/EsaGRACE/runme.m	(revision 22815)
@@ -1,5 +1,7 @@
 
 clear all;
-steps=[0:5]; % [1:5]; 
+addpath('../Data','../Functions');
+
+steps=[0]; % [0:5]; 
 
 if any(steps==0) % Download GRACE land_mass data {{{
@@ -12,4 +14,6 @@
    mget(f,'GRCTellus.JPL.200204_201701.LND.RL05_1.DSTvSCS1411.nc');
 
+	!mv *.nc '../Data/'
+	
 	% display the content of the netcdf file. 
 	%ncdisp('GRCTellus.JPL.200204_201701.LND.RL05_1.DSTvSCS1411.nc')
@@ -122,6 +126,8 @@
 	sol3 = md.results.EsaSolution.EsaNmotion*1000; % [mm] 
 	sol4 = md.results.EsaSolution.EsaEmotion*1000; % [mm] 
+
 	sol_name={'Change in water equivalent height [cm]', 'Vertical displacement [mm]',...
 		'Horizontal (NS) displacement [mm]', 'Horizontal (EW) displacement [mm]'}; 
+	fig_name={'Fig_dH.pdf','Fig_vert.pdf','Fig_horzNS.pdf','Fig_horzEW.pdf'}; 
 
 	res = 1.0; % degree 
@@ -160,18 +166,23 @@
 		set(0,'DefaultAxesFontSize',18,'DefaultAxesLineWidth',1,'DefaultTextFontSize',18,'DefaultLineMarkerSize',8)
 		figure1=figure('Position', [100, 100, 1000, 500]); 
-		gcf; 
-		load coast; 
-		cla; 
+		gcf; load coast; cla; 
 		pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on; 
+		% mask out oceans {{{
+		if (kk==1)
+			geoshow(flipud(lat),flipud(long),'DisplayType','polygon','FaceColor','white'); 
+		end % }}}
 		plot(long,lat,'k'); hold off; 
+		% colormap, xlim etc {{{
 		c1=colorbar;
-		colormap(jet); 
+		colormap('haxby');
+		caxis([-min(abs(min(sol)),abs(max(sol))) min(abs(min(sol)),abs(max(sol)))]); 
 		xlim([-180 180]); 
 		ylim([-90 90]); 
+		% }}}
 		grid on; 
 		title(sol_name(kk)); 
 		set(gcf,'color','w');
-
-		%export_fig('Fig5.pdf'); 
+		
+		%export_fig(fig_name{kk}); 
 	end
 
Index: /issm/trunk-jpl/examples/EsaWahr/runme.m
===================================================================
--- /issm/trunk-jpl/examples/EsaWahr/runme.m	(revision 22814)
+++ /issm/trunk-jpl/examples/EsaWahr/runme.m	(revision 22815)
@@ -3,5 +3,5 @@
 addpath('../Functions');
 
-steps=[0:6]; % [0:6]; 
+steps=[0]; % [0:6]; 
 
 if any(steps==0) % Simple mesh creation {{{
Index: /issm/trunk-jpl/examples/SlrFarrell/runme.m
===================================================================
--- /issm/trunk-jpl/examples/SlrFarrell/runme.m	(revision 22814)
+++ /issm/trunk-jpl/examples/SlrFarrell/runme.m	(revision 22815)
@@ -1,5 +1,6 @@
 
 clear all;
-steps=[5]; % 
+
+steps=[1]; % [1:5]; 
 
 if any(steps==1) % Global mesh creation {{{ 
@@ -128,4 +129,7 @@
 	md.verbose=verbose('111111111');
 
+	% Choose different convergence threshold. [10% 1% 0.1%] to match Farrell 3 panels in Fig. 1
+	md.slr.reltol = 0.1/100; % per cent change in solution 
+
 	% Solve 
 	md=solve(md,'Slr');
@@ -141,6 +145,6 @@
 	% solutions. 
 	sol = md.results.SealevelriseSolution.Sealevel*100; % per cent normalized by GMSL (which 1 m)  
-	
-	res = 1.0; % degree 
+
+	res = 1; % degree 
 
 	% Make a grid of lats and lons, based on the min and max of the original vectors
@@ -150,6 +154,6 @@
 	% Make a interpolation object
 	F = scatteredInterpolant(md.mesh.lat,md.mesh.long,sol); 
-	F.Method = 'linear';
-	F.ExtrapolationMethod = 'linear'; 
+	F.Method = 'natural'; % for smooth contour 
+	F.ExtrapolationMethod = 'none'; 
 	
 	% Do the interpolation to get gridded solutions... 
@@ -161,17 +165,18 @@
 	figure1=figure('Position', [100, 100, 1000, 500]); 
 	gcf; load coast; cla; 
-	%pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on;
-	contourf(lon_grid,lat_grid,sol_grid,[96 98 100 102 104 105]); shading flat; hold on;
-	geoshow(lat,long,'DisplayType','polygon','FaceColor','white'); 
+	pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on;
+	[C,h]=contour(lon_grid,lat_grid,sol_grid,[96 98 100 102 104 105],'-k','linewidth',2);
+	clabel(C,h,'FontSize',18,'Color','red','LabelSpacing',500); 
+	geoshow(lat,long,'DisplayType','polygon','FaceColor',[.82 .82 .82]); 
 	plot(long,lat,'k'); hold off; 
 	% define colormap, caxis, xlim etc {{{
 	c1=colorbar;
-	%colormap('haxby'); 
-	%caxis([96 104]); 
-	xlim([-180 180]); 
-	ylim([-90 90]); 
+	colormap(flipud(haxby)); 
+	caxis([96 105]); 
+	xlim([-170 170]); 
+	ylim([-85 85]); 
 	% }}} 
 	grid on; 
-	title('Relative sea-level [mm]'); 
+	title('Relative sea-level [% of GMSL]'); 
 	set(gcf,'color','w');