Ignore:
Timestamp:
12/30/20 13:58:25 (4 years ago)
Author:
jdquinn
Message:

CHG: Committing changes to UQ example (still need to modify plotting of partitions); minor adjustment to output assignment

File:
1 edited

Legend:

Unmodified
Added
Removed

  • TabularUnified issm/trunk-jpl/examples/UncertaintyQuantification/runme.m

    r25915 r25923  
    6969        %partition the mesh
    7070        npart=50;
    71         partition=partitioner(md,'package','chaco','npart',npart,'weighting','on')-1;
     71        [partition,md]=partitioner(md,'package','chaco','npart',npart,'weighting','on');
     72        partition=partition-1;
    7273
    7374        %make DeltaHH into our 3 sigma deviation
    … …  
    8283
    8384        %responses
    84         md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],...
    85                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    86         md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],...
    87                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]); %grounding line
    88         md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],...
    89                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    90         md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],...
    91                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    92         md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],...
    93                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    94         md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],...
    95                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    96         md.qmu.responses.MassFlux7=response_function('indexed_MassFlux_7',[],...
    97                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    98         md.qmu.responses.MassFlux8=response_function('indexed_MassFlux_8',[],...
    99                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    100         md.qmu.responses.MassFlux9=response_function('indexed_MassFlux_9',[],...
    101                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    102         md.qmu.responses.MassFlux10=response_function('indexed_MassFlux_10',[],...
    103                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    104         md.qmu.responses.MassFlux11=response_function('indexed_MassFlux_11',[],...
    105                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    106         md.qmu.responses.MassFlux12=response_function('indexed_MassFlux_12',[],...
    107                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    108         md.qmu.responses.MassFlux13=response_function('indexed_MassFlux_13',[],...
    109                 [0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     85        md.qmu.responses.MassFlux1=response_function('descriptor','indexed_MassFlux_1');
     86        md.qmu.responses.MassFlux2=response_function('descriptor','indexed_MassFlux_2'); %grounding line
     87        md.qmu.responses.MassFlux3=response_function('descriptor','indexed_MassFlux_3');
     88        md.qmu.responses.MassFlux4=response_function('descriptor','indexed_MassFlux_4');
     89        md.qmu.responses.MassFlux5=response_function('descriptor','indexed_MassFlux_5');
     90        md.qmu.responses.MassFlux6=response_function('descriptor','indexed_MassFlux_6');
     91        md.qmu.responses.MassFlux7=response_function('descriptor','indexed_MassFlux_7');
     92        md.qmu.responses.MassFlux8=response_function('descriptor','indexed_MassFlux_8');
     93        md.qmu.responses.MassFlux9=response_function('descriptor','indexed_MassFlux_9');
     94        md.qmu.responses.MassFlux10=response_function('descriptor','indexed_MassFlux_10');
     95        md.qmu.responses.MassFlux11=response_function('descriptor','indexed_MassFlux_11');
     96        md.qmu.responses.MassFlux12=response_function('descriptor','indexed_MassFlux_12');
     97        md.qmu.responses.MassFlux13=response_function('descriptor','indexed_MassFlux_13');
    11098
    11199        %mass flux profiles
    … …  
    125113                'MassFlux13.exp'...
    126114        };
    127         md.qmu.mass_flux_profile_directory='../MassFluxes/';
     115        md.qmu.mass_flux_profile_directory='./MassFluxes/';
    128116
    129117        %%  sampling analysis
    … …  
    171159
    172160        %partition the mesh
    173         md.qmu.numberofpartitions=10;
    174         md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
    175         md.qmu.partition=md.qmu.partition-1; %switch partition to c-indexing
     161        npart=10;
     162        [partition,md]=partitioner(md,'package','chaco','weighting','on');
     163        partition=partition-1; %switch partition to c-indexing
    176164
    177165        %all types of variables and responses: scaled_Thickness, indexed_MassFlux_i,MaxVel,nodal_DragCoefficient_i. scaled variables are expanded.
    178166
    179167        %variables
    180         md.qmu.variables.DragCoefficient=normal_uncertain(...
    181         'scaled_FrictionCoefficient',1,0.05);
    182         md.qmu.variables.rheology_B=normal_uncertain(...
    183         'scaled_MaterialsRheologyB',1,0.05);
    184         md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);
     168        md.qmu.variables.drag_coefficient=normal_uncertain('descriptor','scaled_FrictionCoefficient',...
     169                'mean',ones(npart,1),...
     170                'stddev',0.05*ones(npart,1),...
     171                'partition',partition);
     172        md.qmu.variables.rheology_B=normal_uncertain('descriptor','scaled_MaterialsRheologyB',...
     173                'mean',ones(npart,1),...
     174                'stddev',0.05*ones(npart,1),...
     175                'partition',partition);
     176        md.qmu.variables.thickness=normal_uncertain('descriptor','scaled_Thickness',...
     177                'mean',ones(npart,1),...
     178                'stddev',0.05*ones(npart,1),...
     179                'partition',partition);
    185180
    186181        %responses
    187         md.qmu.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    188         md.qmu.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    189         md.qmu.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    190         md.qmu.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    191         md.qmu.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    192         md.qmu.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    193         md.qmu.responses.MassFlux7=response_function('indexed_MassFlux_7',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    194         md.qmu.responses.MassFlux8=response_function('indexed_MassFlux_8',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    195         md.qmu.responses.MassFlux9=response_function('indexed_MassFlux_9',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    196         md.qmu.responses.MassFlux10=response_function('indexed_MassFlux_10',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    197         md.qmu.responses.MassFlux11=response_function('indexed_MassFlux_11',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    198         md.qmu.responses.MassFlux12=response_function('indexed_MassFlux_12',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
    199         md.qmu.responses.MassFlux13=response_function('indexed_MassFlux_13',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
     182        md.qmu.responses.MassFlux1=response_function('descriptor','indexed_MassFlux_1');
     183        md.qmu.responses.MassFlux2=response_function('descriptor','indexed_MassFlux_2');
     184        md.qmu.responses.MassFlux3=response_function('descriptor','indexed_MassFlux_3');
     185        md.qmu.responses.MassFlux4=response_function('descriptor','indexed_MassFlux_4');
     186        md.qmu.responses.MassFlux5=response_function('descriptor','indexed_MassFlux_5');
     187        md.qmu.responses.MassFlux6=response_function('descriptor','indexed_MassFlux_6');
     188        md.qmu.responses.MassFlux7=response_function('descriptor','indexed_MassFlux_7');
     189        md.qmu.responses.MassFlux8=response_function('descriptor','indexed_MassFlux_8');
     190        md.qmu.responses.MassFlux9=response_function('descriptor','indexed_MassFlux_9');
     191        md.qmu.responses.MassFlux10=response_function('descriptor','indexed_MassFlux_10');
     192        md.qmu.responses.MassFlux11=response_function('descriptor','indexed_MassFlux_11');
     193        md.qmu.responses.MassFlux12=response_function('descriptor','indexed_MassFlux_12');
     194        md.qmu.responses.MassFlux13=response_function('descriptor','indexed_MassFlux_13');
    200195
    201196        %mass flux profiles
    … …  
    215210                'MassFlux13.exp'...
    216211        };
    217         md.qmu.mass_flux_profile_directory='../MassFluxes/';
     212        md.qmu.mass_flux_profile_directory='./MassFluxes/';
    218213
    219214        %method: local reliability
    … …  
    267262
    268263        %which profile are we looking at?
    269         index=1;       
     264        index=1;
    270265
    271266        %retrieve results for the specific profile, mass flux in m^3 water equiv/s
    272         result=md.results.dakota.dresp_dat(md.qmu.numberofpartitions+index);
     267        result=md.results.dakota.dresp_dat(npart+index);
    273268        result.sample=result.sample/1e12*60*60*24*365;
    274                
     269
    275270        %plot histogram
    276271        plot_hist_norm(result,'cdfleg','off','cdfplt','off','nrmplt','off',...
    … …  
    292287
    293288        %To plot sensitivities
    294         sa=md.results.dakota.dresp_out(index).sens(1:10); sa=sa(md.qmu.partition+1)/1e12*60*60*24*365;
    295         sb=md.results.dakota.dresp_out(index).sens(11:20); sb=sb(md.qmu.partition+1)/1e12*60*60*24*365;
    296         sh=md.results.dakota.dresp_out(index).sens(21:30); sh=sh(md.qmu.partition+1)/1e12*60*60*24*365;
     289        sa=md.results.dakota.dresp_out(index).sens(1:10); sa=sa(partition+1)/1e12*60*60*24*365;
     290        sb=md.results.dakota.dresp_out(index).sens(11:20); sb=sb(partition+1)/1e12*60*60*24*365;
     291        sh=md.results.dakota.dresp_out(index).sens(21:30); sh=sh(partition+1)/1e12*60*60*24*365;
    297292
    298293        plotmodel(md,'data',sh,'data',sa,'data',sb,'expdisp#all',...
Note: See TracChangeset for help on using the changeset viewer.