Changeset 18272

Timestamp:

07/21/14 14:56:29 (11 years ago)

Author:

schlegel

Message:

CHG: update cheatsheet

Location:

issm/trunk-jpl/examples/Inversion

Files:

• CheatSheet.m (moved) (moved from issm/trunk-jpl/examples/Inversion/runme2.m ) (5 diffs)
• runme.m (modified) (4 diffs)

issm/trunk-jpl/examples/Inversion/CheatSheet.m ¶

@@ -4 +4 @@
         md = model;
         md = triangle(md,'DomainOutline.exp',100000);
+
+        %Hands on: STEP 1
+        md.geometry.base    = md.geometry.base+100;
+        md.geometry.surface = md.geometry.surface+100;
+
+        %Hands on: STEP 2
         md = setmask(md,'','');
+
         md = parameterize(md,'Square.par');
         md = setflowequation(md,'SSA','all');
         md.cluster = generic('np',2);
 
-        %START
-        md.geometry.base    = md.geometry.base+100;
-        md.geometry.surface = md.geometry.surface+100;
+        %Hands on: STEP 3
+        md.materials.rheology_B(:) = 1.8*10^8;
+
+        %Hands on: STEP 4
         md.friction.coefficient(:)=50;
         md.friction.coefficient(find(md.mesh.x<600000 & md.mesh.x>400000))=10;
-        %END
 
         md = solve(md,StressbalanceSolutionEnum);
@@ -23 +30 @@
         loadmodel('model1.mat');
 
+        %Hands on: STEP 1
         md.friction.coefficient(:)=50;
 
@@ -32 +40 @@
 
         md = solve(md,StressbalanceSolutionEnum);
+
         save model2 md
 end
@@ -41 +50 @@
         maxsteps = 20;
         md.inversion.iscontrol = 1;
+
+        %Hands on: STEP 1
         md.inversion.control_parameters = {'FrictionCoefficient'};
         md.inversion.maxsteps = maxsteps;
+
+        %Hands on: STEP 2
         md.inversion.cost_functions = 101;
         md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,1);
+
+        %Hands on: STEP 3
         md.inversion.min_parameters    = 1*ones(md.mesh.numberofvertices,1);
         md.inversion.max_parameters    = 100*ones(md.mesh.numberofvertices,1);
@@ -51 +66 @@
         md.verbose=verbose(0);
         md=solve(md,StressbalanceSolutionEnum);
+
         plotmodel(md,'data',md.results.StressbalanceSolution.FrictionCoefficient,'figure',1);
         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  

issm/trunk-jpl/examples/Inversion/runme.m ¶

@@ -9 +9 @@
         md.cluster = generic('np',2);
         md = solve(md,StressbalanceSolutionEnum);
+
         plotmodel(md,'data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'figure',1);
         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  
@@ -25 +26 @@
 
         md = solve(md,StressbalanceSolutionEnum);
+
         plotmodel(md,'data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'figure',1);
-        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  
+        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2, 'gridded#all',1);  
+
         save model2 md
 end
@@ -46 +49 @@
         md.verbose=verbose(0);
         md=solve(md,StressbalanceSolutionEnum);
+
         plotmodel(md,'data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'figure',1);
-        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  
+        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2,'gridded#all',1);  
 end
 if step==4
@@ -67 +71 @@
         md.verbose=verbose(0);
         md=solve(md,StressbalanceSolutionEnum);
+
         plotmodel(md,'data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'figure',1);
         plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);