Changeset 20712

Timestamp:

06/07/16 16:08:16 (9 years ago)

Author:

Mathieu Morlighem

Message:

CHG: minor adjustments

File:

• issm/trunk-jpl/examples/Inversion/CheatSheet.m (modified) (5 diffs)

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

@@ -1 +1 @@
 step=1;
 if step==1
-        %Generate observation
+        %Generate observations
         md = model;
         md = triangle(md,'DomainOutline.exp',100000);
@@ -23 +23 @@
         md.friction.coefficient(find(md.mesh.x<600000 & md.mesh.x>400000))=10;
 
-        md = solve(md,StressbalanceSolutionEnum);
+        md = solve(md,StressbalanceSolutionEnum());
 
-        plotmodel(md,'data',md.friction.coefficient,'figure',1);
-        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);
-
+        plotmodel(md,'axis#all','tight','data',md.friction.coefficient,'caxis',[0 100],'title','"True" \alpha',...
+                'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"')
         save model1 md
 end
@@ -42 +41 @@
         md.inversion.vel_obs = md.results.StressbalanceSolution.Vel;
 
-        md = solve(md,StressbalanceSolutionEnum);
+        md = solve(md,StressbalanceSolutionEnum());
 
-        plotmodel(md,'data',md.friction.coefficient,'figure',1);
-        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);
+        plotmodel(md,'axis#all','tight','data',md.friction.coefficient,'caxis',[0 100],'title','\alpha first guess',...
+                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 
         save model2 md
@@ -62 +61 @@
 
         %Hands on: STEP 2
-        md.inversion.cost_functions = 101;
-        md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,1);
+        md.inversion.cost_functions = [101 103];
+        md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,2);
 
         %Hands on: STEP 3
@@ -71 +70 @@
         %Go solve!
         md.verbose=verbose(0);
-        md=solve(md,StressbalanceSolutionEnum);
+        md=solve(md,StressbalanceSolutionEnum());
 
-        plotmodel(md,'data',md.results.StressbalanceSolution.FrictionCoefficient,'figure',1,'caxis',[10 50]);
-        plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2);  
+        plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.FrictionCoefficient,'caxis',[0 100],'title','inferred \alpha',...
+                'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 end