source: issm/trunk-jpl/test/NightlyRun/test445.m@ 24870

Last change on this file since 24870 was 24870, checked in by jdquinn, 5 years ago

CHG: Translated MATLAB changes to Python that allow for distributing variable over partition.

File size: 3.5 KB
RevLine 
[23438]1%Test Name: SquareSheetShelfSteaEnthalpyHO3dDakotaSampNeff
2md=triangle(model(),'../Exp/Square.exp',150000.);
3md=setmask(md,'../Exp/SquareShelf.exp','');
4md=parameterize(md,'../Par/SquareSheetShelf.par');
5md=extrude(md,3,2.);
6md=setflowequation(md,'HO','all');
7md.cluster=generic('name',oshostname(),'np',3);
8md.timestepping.time_step=0.;
9md.thermal.isenthalpy=1;
10md.thermal.isdynamicbasalspc=1;
11md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
12md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
13
14md.friction.coupling=3;
15md.friction.effective_pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness+md.materials.rho_water*md.constants.g*md.geometry.base;
16
17%dakota version
18version=IssmConfig('_DAKOTA_VERSION_'); version=version(1:3); version=str2num(version);
19
[24836]20%partitioning
21md.qmu.numberofpartitions=10;
22md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
23md.qmu.vpartition=md.qmu.vpartition-1;
24md.qmu.isdakota=1;
25
[23438]26%variables
[24836]27md.qmu.variables.neff=normal_uncertain('descriptor','scaled_FrictionEffectivePressure',...
28 'mean',ones(md.qmu.numberofpartitions,1),...
29 'stddev',.05*ones(md.qmu.numberofpartitions,1),...
30 'partition',md.qmu.vpartition);
[24870]31md.qmu.variables.geoflux=normal_uncertain('descriptor','scaled_BasalforcingsGeothermalflux',...
[24836]32 'mean',ones(md.qmu.numberofpartitions,1),...
33 'stddev',.05*ones(md.qmu.numberofpartitions,1),...
34 'partition',md.qmu.vpartition);
[23438]35
36%responses
[24839]37md.qmu.responses.MaxVel=response_function('descriptor','MaxVel');
38md.qmu.responses.MassFlux1=response_function('descriptor','indexed_MassFlux_1');
39md.qmu.responses.MassFlux2=response_function('descriptor','indexed_MassFlux_2');
40md.qmu.responses.MassFlux3=response_function('descriptor','indexed_MassFlux_3');
41md.qmu.responses.MassFlux4=response_function('descriptor','indexed_MassFlux_4');
42md.qmu.responses.MassFlux5=response_function('descriptor','indexed_MassFlux_5');
43md.qmu.responses.MassFlux6=response_function('descriptor','indexed_MassFlux_6');
44md.qmu.responses.MassFlux7=response_function('descriptor','indexed_MassFlux_7');
[23438]45
46%mass flux profiles
47md.qmu.mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp','../Exp/Square.exp'};
48md.qmu.mass_flux_profile_directory=pwd;
49
50md.qmu.method =dakota_method('nond_samp');
51md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
52 'seed',1234,...
53 'samples',20,...
54 'sample_type','random');
55
56%% a variety of parameters
57md.qmu.params.direct=true;
58md.qmu.params.analysis_components='';
59md.qmu.params.tabular_graphics_data=true;
60
61if version>=6,
62 md.qmu.params.analysis_driver='matlab';
63 md.qmu.params.evaluation_scheduling='master';
64 md.qmu.params.processors_per_evaluation=2;
65else
[23865]66 md.qmu.params.analysis_driver='stressbalance';
[23438]67 md.qmu.params.evaluation_concurrency=1;
68end
69
70
71md.stressbalance.reltol=10^-5; %tighten for qmu analyses
72
73md=solve(md,'Steadystate','overwrite','y');
74
75%Fields and tolerances to track changes
76md.qmu.results=md.results.dakota;
77
78%we put all the mean and stdev data in the montecarlo field, which we will use to test for success.
79md.results.dakota.montecarlo=[];
80for i=1:8,
81 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).mean];
82end
83for i=1:8,
84 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).stddev];
85end
86field_names ={'montecarlo'};
[23445]87field_tolerances={2e-10};
[23438]88field_values={...
89 md.results.dakota.montecarlo,...
90 };
91
Note: See TracBrowser for help on using the repository browser.