Changeset 22811

Timestamp:

05/29/18 14:33:48 (7 years ago)

Author:

adhikari

Message:

minor

File:

• issm/trunk-jpl/examples/EsaWahr/runme.m (modified) (10 diffs)

issm/trunk-jpl/examples/EsaWahr/runme.m

@@ -1 +1 @@
 
 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 +16 @@
         %export_fig('Fig0.pdf'); 
 
-end
-
-if any(steps==1)
+end % }}} 
+
+if any(steps==1) % {{{ Anisotropic mesh creation  
         disp('   Step 1: Anisotropic mesh creation');
 
@@ -36 +38 @@
         %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 +63 @@
         %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 +78 @@
         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 +91 @@
         % 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 +113 @@
         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 +150 @@
         %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 +191 @@
                 % 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 +198 @@
         %export_fig('Fig6.pdf'); 
 
-end
-
+end % }}} 
+