Changeset 25981

Timestamp:

02/04/21 15:49:51 (4 years ago)

Author:

jdquinn

Message:

BUG: Last plot was not working

File:

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

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

@@ -1 @@
-%PIG Uncertainty Quantification Application
+%PIG Uncertainty Quantification Tutorial
 steps=[1];
 
@@ -115 +115 @@
         md.qmu.mass_flux_profile_directory='./MassFluxes/';
 
-        %%  sampling analysis
+        %sampling analysis
         md.qmu.method           =dakota_method('nond_samp');
         md.qmu.method(end)      =dmeth_params_set(...
@@ -124 +124 @@
         ); %random or lhs
 
-        %%  a variety of parameters
+        %a variety of parameters
         md.qmu.params.direct=true;
         md.qmu.params.analysis_driver='';
@@ -136 +136 @@
         md.verbose=verbose(0);
 
-        %Here, we choose to run with 4 processors, 3 for DAKOTA
+        %Here, we choose to run with 4 processors, 3 for Dakota
         % while one serves as the master
         md.cluster=generic('name',oshostname,'np',4);
@@ -226 +226 @@
         md.stressbalance.restol=10^-5; %tighten for qmu analyses
 
-        %Here, we choose to run with 2 processors, 1 for DAKOTA
+        %Here, we choose to run with 2 processors, 1 for Dakota
         % while one serves as the master
         md.cluster=generic('name',oshostname,'np',2);
@@ -250 +250 @@
         md = loadmodel('./Models/PIG.Sampling');
 
-        plotmodel(md,'data','mesh','partitionedges','on',...
+        plotmodel(md,'data','mesh','partition',partition,...
         'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
         'title','Partition Edges on ISSM mesh','grid','on');
@@ -301 +301 @@
 
         %To plot importance factors
-        ifa=importancefactors(md,'scaled_FrictionCoefficient',['indexed_MassFlux_' num2str(index)]);
-        ifb=importancefactors(md,'scaled_MaterialsRheologyB',['indexed_MassFlux_' num2str(index)]);
-        ifh=importancefactors(md,'scaled_Thickness',['indexed_MassFlux_' num2str(index)]);
+        ifa=importancefactors(md,'scaled_FrictionCoefficient',['indexed_MassFlux_' num2str(index)],partition);
+        ifb=importancefactors(md,'scaled_MaterialsRheologyB',['indexed_MassFlux_' num2str(index)],partition);
+        ifh=importancefactors(md,'scaled_Thickness',['indexed_MassFlux_' num2str(index)],partition);
 
         plotmodel(md,'data',ifh,'data',ifa,'data',ifb,'expdisp#all',...