Changeset 20739
- Timestamp:
- 06/15/16 17:17:11 (9 years ago)
- Location:
- issm/trunk-jpl/examples/UncertaintyQuantification
- Files:
-
- 18 added
- 2 deleted
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/UncertaintyQuantification/runme.m
r20532 r20739 11 11 12 12 for i=1:13, 13 contour=expread(['./ MassFluxes/MassFlux' num2str(i) '.exp']);13 contour=expread(['./Exp/MassFluxes/MassFlux' num2str(i) '.exp']); 14 14 textpositions{i}=[contour.x(end) contour.y(end)]; 15 15 end 16 16 17 17 plotmodel(md,'data',md.results.StressbalanceSolution.Vel,'log',10,'expdisp',... 18 {' MassFluxes/MassFlux1.exp','MassFluxes/MassFlux2.exp',...19 ' MassFluxes/MassFlux3.exp','MassFluxes/MassFlux4.exp',...20 ' MassFluxes/MassFlux5.exp','MassFluxes/MassFlux6.exp',...21 ' MassFluxes/MassFlux7.exp','MassFluxes/MassFlux8.exp',...22 ' MassFluxes/MassFlux9.exp','MassFluxes/MassFlux10.exp',...23 ' MassFluxes/MassFlux11.exp','MassFluxes/MassFlux12.exp',...24 ' MassFluxes/MassFlux13.exp'},...18 {'Exp/MassFluxes/MassFlux1.exp','Exp/MassFluxes/MassFlux2.exp',... 19 'Exp/MassFluxes/MassFlux3.exp','Exp/MassFluxes/MassFlux4.exp',... 20 'Exp/MassFluxes/MassFlux5.exp','Exp/MassFluxes/MassFlux6.exp',... 21 'Exp/MassFluxes/MassFlux7.exp','Exp/MassFluxes/MassFlux8.exp',... 22 'Exp/MassFluxes/MassFlux9.exp','Exp/MassFluxes/MassFlux10.exp',... 23 'Exp/MassFluxes/MassFlux11.exp','Exp/MassFluxes/MassFlux12.exp',... 24 'Exp/MassFluxes/MassFlux13.exp'},... 25 25 'expstyle',{'k-','k-','k-','k-','k-','k-','k-',... 26 26 'k-','k-','k-','k-','k-','k-'},'linewidth',2,... … … 45 45 46 46 %filter out unrealistic error ranges 47 flags=ContourToNodes(md.mesh.x,md.mesh.y,'E rrorContour.exp',1);47 flags=ContourToNodes(md.mesh.x,md.mesh.y,'Exp/ErrorContour.exp',1); 48 48 pos=find(~flags); DeltaHH(pos)=0; 49 49 … … 123 123 'MassFlux12.exp',... 124 124 'MassFlux13.exp'}; 125 md.qmu.mass_flux_profile_directory='../ MassFluxes/';125 md.qmu.mass_flux_profile_directory='../Exp/MassFluxes/'; 126 126 127 127 %% sampling analysis … … 144 144 md.verbose=verbose(0); 145 145 146 % solve147 md.qmu.isdakota=1; md.inversion.iscontrol=0;146 %Here, we choose to run with 4 processors, 3 for DAKOTA 147 % while one serves as the master 148 148 md.cluster=generic('name',oshostname,'np',4); 149 149 … … 153 153 md.qmu.params.processors_per_evaluation=md.cluster.np-1; 154 154 155 %Turn dakota on 156 md.qmu.isdakota=1; md.inversion.iscontrol=0; 155 157 md=solve(md,StressbalanceSolutionEnum,'overwrite','y'); 156 158 … … 206 208 'MassFlux12.exp',... 207 209 'MassFlux13.exp'}; 208 md.qmu.mass_flux_profile_directory='../ MassFluxes/';210 md.qmu.mass_flux_profile_directory='../Exp/MassFluxes/'; 209 211 210 212 %method: local reliability … … 221 223 222 224 md.stressbalance.restol=10^-5; %tighten for qmu analyses 223 224 %solve 225 md.qmu.isdakota=1; 226 md.inversion.iscontrol=0; 225 226 %Here, we choose to run with 2 processors, 1 for DAKOTA 227 % while one serves as the master 227 228 md.cluster=generic('name',oshostname,'np',2); 228 229 229 230 %Dakota runs in parallel with a master/slave configuration. 230 231 % At least 2 cpu's are needed to run the UQ 231 232 md.qmu.params.evaluation_scheduling='master'; 232 md.qmu.params.processors_per_evaluation=md.cluster.np-1 233 233 md.qmu.params.processors_per_evaluation=md.cluster.np-1; 234 235 %Clear results and turn dakota on 236 md.results=[]; 237 md.qmu.isdakota=1; md.inversion.iscontrol=0; 234 238 md.verbose=verbose('qmu',true); 235 239 md=solve(md,StressbalanceSolutionEnum,'overwrite','y'); … … 243 247 md = loadmodel('./Models/PIG.Sampling'); 244 248 245 plotmodel(md,'data','mesh','partitionedges','on', 'meshlinewidth',1.5,...249 plotmodel(md,'data','mesh','partitionedges','on',... 246 250 'linewidth',2, 'axis#all','image','unit','km','colorbar','off',... 247 'title',' ','meshcolor','b','grid','on');251 'title','Partition Edges on ISSM mesh','grid','on'); 248 252 249 253 end … … 284 288 285 289 plotmodel(md,'data',sh,'data',sa,'data',sb,'expdisp#all',... 286 [' MassFluxes/MassFlux' num2str(index) '.exp'],...290 ['Exp/MassFluxes/MassFlux' num2str(index) '.exp'],... 287 291 'expstyle#all','b-','linewidth#all',2,... 288 292 'nlines',3,'ncols',1, 'axis#all','image',... 289 293 'colorbar#all','on','colorbarfontsize#all',10,... 290 294 'colorbartitle#1','S_{H}', 'colorbartitle#2','S_{\alpha}',... 291 'colorbartitle#3','S_{B}','unit#all','km','figure',1); 295 'colorbartitle#3','S_{B}','unit#all','km','figure',1,... 296 'title','Sensitivities: H, \alpha, B'); 292 297 293 298 %To plot importance factors … … 297 302 298 303 plotmodel(md,'data',ifh,'data',ifa,'data',ifb,'expdisp#all',... 299 [' MassFluxes/MassFlux' num2str(index) '.exp'],...304 ['Exp/MassFluxes/MassFlux' num2str(index) '.exp'],... 300 305 'expstyle#all','b-','linewidth#all',2,'log#all',10,... 301 306 'nlines',3,'ncols',1, 'axis#all','image','caxis#all',[1e-10 1],... 302 307 'colorbar#all','on','colorbarfontsize#all',10,... 303 308 'colorbartitle#1','If_{H}', 'colorbartitle#2','If_{\alpha}',... 304 'colorbartitle#3','If_{B}','unit#all','km','figure',2); 305 306 end 309 'colorbartitle#3','If_{B}','unit#all','km','figure',2,... 310 'title','Importance Factors: H, \alpha, B'); 311 312 end
Note:
See TracChangeset
for help on using the changeset viewer.