Index: /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m
===================================================================
--- /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 25986)
+++ /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 25987)
@@ -1,6 +1,6 @@
 %PIG Uncertainty Quantification Tutorial
-steps=[1];
-
-if any(steps==1) 
+steps=[1:7];
+
+if any(steps==1)
 	disp('   Step 1: plot flux gates');
 
@@ -28,7 +28,7 @@
 		'8','9','10','11','12','13'},...
 		'textposition',textpositions);
-end 
-
-if any(steps==2) 
+end
+
+if any(steps==2)
 	disp('   Step 2: compute cross overs from CRESIS');
 
@@ -58,7 +58,7 @@
 
 	save Models/PIG.CrossOvers DeltaHH
-end 
-
-if any(steps==3) 
+end
+
+if any(steps==3)
 	disp('   Step 3: sampling analysis');
 
@@ -141,5 +141,5 @@
 
 	%Dakota runs in parallel with a master/slave configuration.
-	% At least 2 cpu's are needed to run the UQ
+	% At least 2 CPUs are needed to run the UQ
 	md.qmu.params.evaluation_scheduling='master';
 	md.qmu.params.processors_per_evaluation=md.cluster.np-1;
@@ -150,8 +150,19 @@
 
 	save ./Models/PIG.Sampling md;
-end 
-
-if any(steps==4) 
-	disp('   Step 4: sensitivity analysis');
+end
+
+if any(steps==4)
+	disp('   Step 4: plot partition');
+
+	%load model
+	md = loadmodel('./Models/PIG.Sampling');
+
+	plotmodel(md,'data','mesh','partition',partition,...
+	'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
+	'title','Partition Edges on ISSM mesh','grid','on');
+end
+
+if any(steps==5)
+	disp('   Step 5: sensitivity analysis');
 
 	%load model
@@ -231,5 +242,5 @@
 
 	%Dakota runs in parallel with a master/slave configuration.
-	% At least 2 cpu's are needed to run the UQ
+	% At least 2 CPUs are needed to run the UQ
 	md.qmu.params.evaluation_scheduling='master';
 	md.qmu.params.processors_per_evaluation=md.cluster.np-1;
@@ -242,18 +253,7 @@
 
 	save ./Models/PIG.Sensitivity md;
-end 
-
-if any(steps==5) 
-	disp('   Step 5: plot partition');
-
-	%load model
-	md = loadmodel('./Models/PIG.Sampling');
-
-	plotmodel(md,'data','mesh','partition',partition,...
-	'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
-	'title','Partition Edges on ISSM mesh','grid','on');
-end 
-
-if any(steps==6) 
+end
+
+if any(steps==6)
 	disp('   Step 6: plot histogram');
 
@@ -272,8 +272,7 @@
 	'xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor',...
 	'none','EdgeColor','red');
-
-end 
-
-if any(steps==7) 
+end
+
+if any(steps==7)
 	disp('   Step 7: plot sensitivity');
 
@@ -295,4 +294,5 @@
 		'expstyle#all','b-','linewidth#all',2,...
 		'nlines',3,'ncols',1, 'axis#all','image',...
+		'caxis#1',[0 150],'caxis#2',[-20 0],'caxis#3',[-25 0],...
 		'colorbar#all','on','colorbarfontsize#all',10,...
 		'colorbartitle#1','S_{H}', 'colorbartitle#2','S_{\alpha}',...
@@ -313,3 +313,3 @@
 		'colorbartitle#3','If_{B}','unit#all','km','figure',2,...
 		'title','Importance Factors: H, \alpha, B');
-end 
+end