Index: /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m
===================================================================
--- /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 25980)
+++ /issm/trunk-jpl/examples/UncertaintyQuantification/runme.m	(revision 25981)
@@ -1,3 +1,3 @@
-%PIG Uncertainty Quantification Application
+%PIG Uncertainty Quantification Tutorial
 steps=[1];
 
@@ -115,5 +115,5 @@
 	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,5 +124,5 @@
 	); %random or lhs
 
-	%%  a variety of parameters
+	%a variety of parameters
 	md.qmu.params.direct=true;
 	md.qmu.params.analysis_driver='';
@@ -136,5 +136,5 @@
 	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,5 +226,5 @@
 	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,5 +250,5 @@
 	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,7 +301,7 @@
 
 	%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',...