Changeset 20739


Ignore:
Timestamp:
06/15/16 17:17:11 (9 years ago)
Author:
schlegel
Message:

CHG: move from workshop directory

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  
    1111
    1212        for i=1:13,
    13                 contour=expread(['./MassFluxes/MassFlux' num2str(i) '.exp']);
     13                contour=expread(['./Exp/MassFluxes/MassFlux' num2str(i) '.exp']);
    1414                textpositions{i}=[contour.x(end) contour.y(end)];
    1515        end
    1616
    1717        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'},...
    2525                'expstyle',{'k-','k-','k-','k-','k-','k-','k-',...
    2626                'k-','k-','k-','k-','k-','k-'},'linewidth',2,...
     
    4545
    4646        %filter out unrealistic error ranges
    47         flags=ContourToNodes(md.mesh.x,md.mesh.y,'ErrorContour.exp',1);
     47        flags=ContourToNodes(md.mesh.x,md.mesh.y,'Exp/ErrorContour.exp',1);
    4848        pos=find(~flags); DeltaHH(pos)=0;
    4949
     
    123123                                   'MassFlux12.exp',...
    124124                                   'MassFlux13.exp'};
    125         md.qmu.mass_flux_profile_directory='../MassFluxes/';
     125        md.qmu.mass_flux_profile_directory='../Exp/MassFluxes/';
    126126
    127127        %%  sampling analysis
     
    144144        md.verbose=verbose(0);
    145145       
    146         %solve
    147         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
    148148        md.cluster=generic('name',oshostname,'np',4);
    149149
     
    153153        md.qmu.params.processors_per_evaluation=md.cluster.np-1;
    154154
     155        %Turn dakota on
     156        md.qmu.isdakota=1; md.inversion.iscontrol=0;
    155157        md=solve(md,StressbalanceSolutionEnum,'overwrite','y');
    156158
     
    206208                                   'MassFlux12.exp',...
    207209                                   'MassFlux13.exp'};
    208         md.qmu.mass_flux_profile_directory='../MassFluxes/';
     210        md.qmu.mass_flux_profile_directory='../Exp/MassFluxes/';
    209211
    210212        %method: local reliability
     
    221223
    222224        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
    227228        md.cluster=generic('name',oshostname,'np',2);
    228 
     229       
    229230        %Dakota runs in parallel with a master/slave configuration.
    230231        % At least 2 cpu's are needed to run the UQ
    231232        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;
    234238        md.verbose=verbose('qmu',true);
    235239        md=solve(md,StressbalanceSolutionEnum,'overwrite','y');
     
    243247        md = loadmodel('./Models/PIG.Sampling');
    244248       
    245         plotmodel(md,'data','mesh','partitionedges','on','meshlinewidth',1.5,...
     249        plotmodel(md,'data','mesh','partitionedges','on',...
    246250        'linewidth',2, 'axis#all','image','unit','km','colorbar','off',...
    247         'title','','meshcolor','b','grid','on');
     251        'title','Partition Edges on ISSM mesh','grid','on');
    248252
    249253end
     
    284288
    285289        plotmodel(md,'data',sh,'data',sa,'data',sb,'expdisp#all',...
    286                 ['MassFluxes/MassFlux' num2str(index) '.exp'],...
     290                ['Exp/MassFluxes/MassFlux' num2str(index) '.exp'],...
    287291                'expstyle#all','b-','linewidth#all',2,...
    288292                'nlines',3,'ncols',1, 'axis#all','image',...
    289293                'colorbar#all','on','colorbarfontsize#all',10,...
    290294                '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');
    292297
    293298        %To plot importance factors
     
    297302
    298303        plotmodel(md,'data',ifh,'data',ifa,'data',ifb,'expdisp#all',...
    299                 ['MassFluxes/MassFlux' num2str(index) '.exp'],...
     304                ['Exp/MassFluxes/MassFlux' num2str(index) '.exp'],...
    300305                'expstyle#all','b-','linewidth#all',2,'log#all',10,...
    301306                'nlines',3,'ncols',1, 'axis#all','image','caxis#all',[1e-10 1],...
    302307                'colorbar#all','on','colorbarfontsize#all',10,...
    303308                '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
     312end
Note: See TracChangeset for help on using the changeset viewer.