Index: /issm/trunk-jpl/examples/EsaWahr/runme.m
===================================================================
--- /issm/trunk-jpl/examples/EsaWahr/runme.m	(revision 22810)
+++ /issm/trunk-jpl/examples/EsaWahr/runme.m	(revision 22811)
@@ -1,7 +1,9 @@
 
 clear all;
-steps=[1:6]; % [1:6]; 
-
-if any(steps==0)
+addpath('../Functions');
+
+steps=[0:6]; % [0:6]; 
+
+if any(steps==0) % Simple mesh creation {{{
 	disp('   Step 0: Mesh creation');
 
@@ -14,7 +16,7 @@
 	%export_fig('Fig0.pdf'); 
 
-end
-
-if any(steps==1)
+end % }}} 
+
+if any(steps==1) % {{{ Anisotropic mesh creation  
 	disp('   Step 1: Anisotropic mesh creation');
 
@@ -36,7 +38,7 @@
 	%export_fig('Fig1.pdf'); 
 
-end
-
-if any(steps==2)
+end % }}} 
+
+if any(steps==2) % Define loads {{{ 
 	disp('   Step 2: Define loads');
 	md = loadmodel('./Models/EsaWahr.Mesh');
@@ -61,7 +63,7 @@
 	%export_fig('Fig2.pdf'); 
 
-end
-
-if any(steps==3)
+end % }}} 
+
+if any(steps==3) % Parameterization {{{ 
 	disp('   Step 3: Parameterization');
 	md = loadmodel('./Models/EsaWahr.Loads');
@@ -76,5 +78,5 @@
 	md.mask.groundedice_levelset = ones(md.mesh.numberofvertices,1); % 1 = ice is grounnded 
 
-	%% IGNORE BUT DO NOT DELETE %% 
+	%% IGNORE BUT DO NOT DELETE %% {{{ 
 	% Geometry: Important only when you want to couple with Ice Flow Model 
 	di=md.materials.rho_ice/md.materials.rho_water;
@@ -89,11 +91,11 @@
 	% Miscellaneous: 
 	md.miscellaneous.name='EsaWahr';
-	%% IGNORE BUT DO NOT DELETE %% 
+	%% IGNORE BUT DO NOT DELETE %% }}}  
 	
 	save ./Models/EsaWahr.Parameterization md; 
 
-end
-
-if any(steps==4)
+end % }}} 
+
+if any(steps==4) % Solve {{{ 
 	disp('   Step 4: Solve Esa solver');
 	md = loadmodel('./Models/EsaWahr.Parameterization');
@@ -111,7 +113,7 @@
 	save ./Models/EsaWahr.Solution md; 
 
-end
-
-if any(steps==5)
+end % }}} 
+
+if any(steps==5) % Plot solutions {{{ 
 	disp('   Step 5: Plot solutions');
 	md = loadmodel('./Models/EsaWahr.Solution');
@@ -148,7 +150,7 @@
 	%export_fig('Fig5.pdf'); 
 
-end
-
-if any(steps==6)
+end % }}} 
+
+if any(steps==6) % Compare results against semi-analytic solutions {{{ 
 	disp('   Step 6: Compare results against Wahr semi-analytic solutions');
 	md = loadmodel('./Models/EsaWahr.Solution');
@@ -189,5 +191,5 @@
 		% first box 
 		ag1 = gca;
-		leg1a = legend(ag1,[h3,h1,h4,h2],'Vertical (Wahr)','Vertical (ISSM)','Horizontal (Wahr)','Horizontal (ISSM)',1);
+		leg1a = legend(ag1,[h3,h1,h4,h2],'Vertical (Wahr)','Vertical (ISSM)','Horizontal (Wahr)','Horizontal (ISSM)');
 		set(leg1a,'location','east','Orientation','Vertical','Box','Off','FontSize',16); 
 		% 
@@ -196,4 +198,4 @@
 	%export_fig('Fig6.pdf'); 
 
-end
-
+end % }}} 
+