Index: /issm/trunk-jpl/examples/Inversion/runme.m
===================================================================
--- /issm/trunk-jpl/examples/Inversion/runme.m	(revision 20707)
+++ /issm/trunk-jpl/examples/Inversion/runme.m	(revision 20708)
@@ -1,5 +1,5 @@
 step=1;
 if step==1
-	%Generate observation
+	%Generate observations
 	md = model;
 	md = triangle(md,'DomainOutline.exp',100000);
@@ -8,10 +8,7 @@
 	md = setflowequation(md,'SSA','all');
 	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,'gridded#all',1);
-	plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'figure',2,'gridded#all',1);  
-
+	md = solve(md,StressbalanceSolutionEnum());
+	plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','"True" B',...
+		'data',md.results.StressbalanceSolution.Vel,'title','"observed velocities"')
 	save model1 md
 end
@@ -27,9 +24,7 @@
 	md.inversion.vel_obs = md.results.StressbalanceSolution.Vel;
 
-	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);  
-
+	md = solve(md,StressbalanceSolutionEnum());
+	plotmodel(md,'axis#all','tight','data',md.materials.rheology_B,'caxis',[ 1.3 1.9]*10^8,'title','B first guess',...
+		'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 	save model2 md
 end
@@ -45,13 +40,12 @@
 	md.inversion.cost_functions = 101;
 	md.inversion.cost_functions_coefficients = ones(md.mesh.numberofvertices,1);
-	md.inversion.min_parameters    = paterson(273)*ones(md.mesh.numberofvertices,1);
-	md.inversion.max_parameters    = paterson(200)*ones(md.mesh.numberofvertices,1);
+	md.inversion.min_parameters    = cuffey(273)*ones(md.mesh.numberofvertices,1);
+	md.inversion.max_parameters    = cuffey(200)*ones(md.mesh.numberofvertices,1);
 
 	%Go solve!
 	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);  
+	md=solve(md,StressbalanceSolutionEnum());
+	plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
+		'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 end
 if step==4
@@ -67,12 +61,11 @@
 	md.inversion.cost_functions_coefficients      = ones(md.mesh.numberofvertices,1);
 	md.inversion.cost_functions_coefficients(:,2) = 10^-16*ones(md.mesh.numberofvertices,1);
-	md.inversion.min_parameters    = paterson(273)*ones(md.mesh.numberofvertices,1);
-	md.inversion.max_parameters    = paterson(200)*ones(md.mesh.numberofvertices,1);
+	md.inversion.min_parameters    = cuffey(273)*ones(md.mesh.numberofvertices,1);
+	md.inversion.max_parameters    = cuffey(200)*ones(md.mesh.numberofvertices,1);
 
 	%Go solve!
 	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);  
+	md=solve(md,StressbalanceSolutionEnum());
+	plotmodel(md,'axis#all','tight','data',md.results.StressbalanceSolution.MaterialsRheologyBbar,'caxis',[ 1.3 1.9]*10^8,'title','inferred B',...
+		'data',md.results.StressbalanceSolution.Vel,'title','modeled velocities')
 end