Changeset 14118
- Timestamp:
- 12/09/12 21:25:20 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/workshop/2012/Talks/10_PIG_UQ/runme.m
r14029 r14118 69 69 70 70 %load model and cross over errors 71 md = loadmodel('. /Models/Pig.Control_drag');71 md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag'); 72 72 load -mat Models/Pig.CrossOvers 73 73 74 74 %partition the mesh 75 75 md.qmu.numberofpartitions=50; 76 md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on'); 76 md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,... 77 'weighting','on'); 77 78 md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing 78 79 … … 144 145 145 146 %load model 146 md = loadmodel('. /Models/Pig.Control_drag');147 md = loadmodel('../09_PIG_Application/Models/Pig.Control_drag'); 147 148 148 149 %partition the mesh … … 155 156 %variables 156 157 md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05); 157 md.qmu.variables.DragCoefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.05); 158 md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,0.05); 158 md.qmu.variables.DragCoefficient=normal_uncertain(... 159 'scaled_FrictionCoefficient',1,0.05); 160 md.qmu.variables.rheology_B=normal_uncertain(... 161 'scaled_MaterialsRheologyB',1,0.05); 159 162 160 163 %responses … … 218 221 md = loadmodel('./Models/Pig.Sampling'); 219 222 220 plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,'linewidth',2,'axis#all','image','unit','km','colorbar','off','title','','meshcolor','b','grid','on'); 223 plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,... 224 'linewidth',2, 'axis#all','image','unit','km','colorbar','off',... 225 'title','','meshcolor','b','grid','on'); 221 226 222 227 end … … 232 237 %retrieve results for this specific profile 233 238 result=md.results.dakota.dresp_dat(md.qmu.numberofpartitions+index); 234 error;235 239 236 240 %plot histogram 237 plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off','xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor','none','EdgeColor','red'); 241 plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off',... 242 'xlabelplt','M (Gt/yr)','ylabelplt','F','FontSize',8,'FaceColor',... 243 'none','EdgeColor','red'); 238 244 239 245 end
Note:
See TracChangeset
for help on using the changeset viewer.