Changeset 25987 for issm/trunk-jpl/examples/UncertaintyQuantification/runme.m

Timestamp:

02/08/21 14:40:10 (4 years ago)

Author:

jdquinn

Message:

CHG: Reordering of steps; added caxis bounds to sensitivity plots

File:

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

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

@@ -1 +1 @@
 %PIG Uncertainty Quantification Tutorial
-steps=[1];
-
-if any(steps==1) 
+steps=[1:7];
+
+if any(steps==1)
         disp('   Step 1: plot flux gates');
 
@@ -28 +28 @@
                 '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 +58 @@
 
         save Models/PIG.CrossOvers DeltaHH
-end 
-
-if any(steps==3) 
+end
+
+if any(steps==3)
         disp('   Step 3: sampling analysis');
 
@@ -141 +141 @@
 
         %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 +150 @@
 
         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 +242 @@
 
         %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 +253 @@
 
         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 +272 @@
         '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 +294 @@
                 '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 +313 @@
                 'colorbartitle#3','If_{B}','unit#all','km','figure',2,...
                 'title','Importance Factors: H, \alpha, B');
-end 
+end