Changeset 14118

Timestamp:

12/09/12 21:25:20 (12 years ago)

Author:

Eric.Larour

Message:

CHG: minor

File:

• issm/workshop/2012/Talks/10_PIG_UQ/runme.m (modified) (5 diffs)

issm/workshop/2012/Talks/10_PIG_UQ/runme.m

@@ -69 +69 @@
         
         %load model and cross over errors
-        md = loadmodel('./Models/Pig.Control_drag');
+        md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag');
         load -mat Models/Pig.CrossOvers
 
         %partition the mesh
         md.qmu.numberofpartitions=50;
-        md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+        md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,...
+        'weighting','on');
         md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing
 
@@ -144 +145 @@
         
         %load model
-        md = loadmodel('./Models/Pig.Control_drag');
+        md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag');
 
         %partition the mesh
@@ -155 +156 @@
         %variables
         md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);
-        md.qmu.variables.DragCoefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.05);
-        md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,0.05);
+        md.qmu.variables.DragCoefficient=normal_uncertain(...
+        'scaled_FrictionCoefficient',1,0.05);
+        md.qmu.variables.rheology_B=normal_uncertain(...
+        'scaled_MaterialsRheologyB',1,0.05);
 
         %responses
@@ -218 +221 @@
         md = loadmodel('./Models/Pig.Sampling');
         
-        plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,'linewidth',2,'axis#all','image','unit','km','colorbar','off','title','','meshcolor','b','grid','on');
+        plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,...
+        'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
+        'title','','meshcolor','b','grid','on');
 
 end 
@@ -232 +237 @@
         %retrieve results for this specific profile
         result=md.results.dakota.dresp_dat(md.qmu.numberofpartitions+index); 
-        error;
                 
         %plot histogram
-        plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off','xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor','none','EdgeColor','red');
+        plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off',...
+        'xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor',...
+        'none','EdgeColor','red');
 
 end