source: issm/trunk-jpl/test/NightlyRun/test235.m@ 24840

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

CHG: fixing bad typo.

File size: 2.9 KB
Line 
1%Test Name: SquareShelfTranForceNeg2dDakotaLocal
2md=triangle(model(),'../Exp/Square.exp',180000.);
3md=setmask(md,'all','');
4md=parameterize(md,'../Par/SquareShelf.par');
5md=setflowequation(md,'SSA','all');
6md.cluster=generic('name',oshostname(),'np',3);
7
8md.timestepping.time_step=1;
9md.settings.output_frequency=1;
10md.timestepping.final_time=4;
11
12smb = ones(md.mesh.numberofvertices,1)*3.6;
13smb=[ smb smb*-1 ];
14
15md.smb.mass_balance= smb;
16md.smb.mass_balance(end+1,:)=[1.5 3];
17md.transient.isthermal=0;
18%Dakota options
19
20%dakota version
21version=IssmConfig('_DAKOTA_VERSION_'); version=version(1:3); version=str2num(version);
22
23%partitioning
24md.qmu.numberofpartitions=20;
25md=partitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
26md.qmu.vpartition=md.qmu.vpartition-1;
27
28%variables
29md.qmu.variables.surface_mass_balance=normal_uncertain('descriptor','scaled_SmbMassBalance',...
30 'mean',ones(md.qmu.numberofpartitions,1),...
31 'stddev',100*ones(md.qmu.numberofpartitions,1),...
32 'partition',md.qmu.vpartition);
33
34%responses
35md.qmu.responses.MaxVel=response_function('descriptor','MaxVel');
36md.qmu.responses.IceVolume=response_function('descriptor','IceVolume');
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');
43
44%mass flux profiles
45md.qmu.mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
46md.qmu.mass_flux_profile_directory=pwd;
47
48%method
49md.qmu.method =dakota_method('nond_l');
50
51%parameters
52md.qmu.params.direct=true;
53md.qmu.params.analysis_components='';
54md.qmu.params.interval_type='forward';
55md.qmu.params.fd_gradient_step_size='0.1';
56md.qmu.isdakota=1;
57
58if version>=6,
59 md.qmu.params.analysis_driver='matlab';
60 md.qmu.params.evaluation_scheduling='master';
61 md.qmu.params.processors_per_evaluation=2;
62else
63 md.qmu.params.analysis_driver='stressbalance';
64 md.qmu.params.evaluation_concurrency=1;
65end
66
67
68md.stressbalance.reltol=10^-5; %tighten for qmu analyses
69md.transient.requested_outputs={'IceVolume'};
70
71%solve
72md=solve(md,'Transient','overwrite','y');
73md.qmu.results=md.results.dakota;
74
75%Fields and tolerances to track changes
76md.results.dakota.moments=[];
77for i=1:8,
78 md.results.dakota.moments=[md.results.dakota.moments md.results.dakota.dresp_out(i).mean];
79end
80for i=1:8,
81 md.results.dakota.moments=[md.results.dakota.moments md.results.dakota.dresp_out(i).stddev];
82end
83field_names ={'moments'};
84field_tolerances={1e-11};
85field_values={...
86 md.results.dakota.moments,...
87 };
Note: See TracBrowser for help on using the repository browser.