source: issm/trunk/test/NightlyRun/test626.m@ 8295

Last change on this file since 8295 was 8295, checked in by Mathieu Morlighem, 14 years ago

no need to provide analysis_type in solve command anymore

File size: 1.5 KB
RevLine 
[6265]1md=mesh(model,'../Exp/79North.exp',10000);
2md=meshconvert(md);
3md=geography(md,'../Exp/79NorthShelf.exp','');
4md=parameterize(md,'../Par/79North.par');
5md=setelementstype(md,'macayeal','all');
6
7%Ice sheet only
8md=modelextract(md,md.elementonicesheet);
9pos=find(md.gridonboundary);
10md.spcthickness(pos,1)=1;
11md.spcthickness(:,2)=md.thickness;
12
13%control parameters
14md.thickness_obs=md.thickness;
15md.control_analysis=1;
16md.thickness_obs=md.thickness;
17md.nsteps=2;
18md.control_type=[VxEnum() VyEnum()];
19md.prognostic_DG=0;
20md.optscal=[10/md.yts*ones(md.nsteps,1) 10/md.yts*ones(md.nsteps,1)];
[8134]21md.cm_min=[-2000*ones(md.numberofgrids,1) -2000*ones(md.numberofgrids,1)];
22md.cm_max=[+2000*ones(md.numberofgrids,1) +2000*ones(md.numberofgrids,1)];
[6265]23md.cm_responses=ThicknessAbsMisfitEnum*ones(md.nsteps,1);
24md.weights=ones(md.numberofgrids,1);
25md.maxiter=4*ones(md.nsteps,1);
26md.cm_jump=0.99*ones(md.nsteps,1);
27
28md=SetParallel(md,3);
[8295]29md=solve(md,BalancethicknessSolutionEnum);
[6265]30
31%Fields and tolerances to track changes
32field_names ={'Gradient1' 'Gradient2' 'Misfits' 'Vx' 'Vy' 'Thickness'};
33field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
34field_values={...
[8291]35 PatchToVec(md.results.BalancethicknessSolution.Gradient1),...
36 PatchToVec(md.results.BalancethicknessSolution.Gradient2),...
37 md.results.BalancethicknessSolution.J,...
38 PatchToVec(md.results.BalancethicknessSolution.Vx),...
39 PatchToVec(md.results.BalancethicknessSolution.Vy),...
40 PatchToVec(md.results.BalancethicknessSolution.Thickness)
[6265]41};
Note: See TracBrowser for help on using the repository browser.