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

Last change on this file since 25005 was 25005, checked in by Eric.Larour, 5 years ago

CHG: updating tests to not rely on md.qmu.vpartition and md.qmu.epartition.

File size: 3.3 KB
Line 
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
20%partitioning
21npart=10;
22partition=partitioner(md,'package','chaco','npart',npart,'weighting','on')-1;
23md.qmu.isdakota=1;
24
25%variables
26md.qmu.variables.neff=normal_uncertain('descriptor','scaled_FrictionEffectivePressure',...
27 'mean',ones(npart,1),...
28 'stddev',.05*ones(npart,1),...
29 'partition',partition);
30md.qmu.variables.geoflux=normal_uncertain('descriptor','scaled_BasalforcingsGeothermalflux',...
31 'mean',ones(npart,1),...
32 'stddev',.05*ones(npart,1),...
33 'partition',partition);
34
35%responses
36md.qmu.responses.MaxVel=response_function('descriptor','MaxVel');
37md.qmu.responses.MassFlux1=response_function('descriptor','indexed_MassFlux_1');
38md.qmu.responses.MassFlux2=response_function('descriptor','indexed_MassFlux_2');
39md.qmu.responses.MassFlux3=response_function('descriptor','indexed_MassFlux_3');
40md.qmu.responses.MassFlux4=response_function('descriptor','indexed_MassFlux_4');
41md.qmu.responses.MassFlux5=response_function('descriptor','indexed_MassFlux_5');
42md.qmu.responses.MassFlux6=response_function('descriptor','indexed_MassFlux_6');
43md.qmu.responses.MassFlux7=response_function('descriptor','indexed_MassFlux_7');
44
45%mass flux profiles
46md.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'};
47md.qmu.mass_flux_profile_directory=pwd;
48
49md.qmu.method =dakota_method('nond_samp');
50md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
51 'seed',1234,...
52 'samples',20,...
53 'sample_type','random');
54
55%% a variety of parameters
56md.qmu.params.direct=true;
57md.qmu.params.analysis_components='';
58md.qmu.params.tabular_graphics_data=true;
59
60if version>=6,
61 md.qmu.params.analysis_driver='matlab';
62 md.qmu.params.evaluation_scheduling='master';
63 md.qmu.params.processors_per_evaluation=2;
64else
65 md.qmu.params.analysis_driver='stressbalance';
66 md.qmu.params.evaluation_concurrency=1;
67end
68
69
70md.stressbalance.reltol=10^-5; %tighten for qmu analyses
71
72md=solve(md,'Steadystate','overwrite','y');
73
74%Fields and tolerances to track changes
75md.qmu.results=md.results.dakota;
76
77%we put all the mean and stdev data in the montecarlo field, which we will use to test for success.
78md.results.dakota.montecarlo=[];
79for i=1:8,
80 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).mean];
81end
82for i=1:8,
83 md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).stddev];
84end
85field_names ={'montecarlo'};
86field_tolerances={2e-10};
87field_values={...
88 md.results.dakota.montecarlo,...
89 };
90
Note: See TracBrowser for help on using the repository browser.