[6265] | 1 | md=mesh(model,'../Exp/79North.exp',10000);
|
---|
| 2 | md=meshconvert(md);
|
---|
[9641] | 3 | md=setmask(md,'../Exp/79NorthShelf.exp','');
|
---|
[6265] | 4 | md=parameterize(md,'../Par/79North.par');
|
---|
[9664] | 5 | md=setflowequation(md,'macayeal','all');
|
---|
[6265] | 6 |
|
---|
| 7 | %Ice sheet only
|
---|
[9641] | 8 | md=modelextract(md,md.mask.elementongroundedice);
|
---|
[8304] | 9 | pos=find(md.nodeonboundary);
|
---|
[9645] | 10 | md.balancethickness.spcthickness(pos)=md.thickness(pos);
|
---|
[6265] | 11 |
|
---|
| 12 | %control parameters
|
---|
| 13 | md.thickness_obs=md.thickness;
|
---|
| 14 | md.control_analysis=1;
|
---|
| 15 | md.thickness_obs=md.thickness;
|
---|
| 16 | md.nsteps=2;
|
---|
[9535] | 17 | md.control_type={'Vx' 'Vy'};
|
---|
[9642] | 18 | md.prognostic.stabilization=1;
|
---|
[9597] | 19 | md.optscal=[10/md.constants.yts*ones(md.nsteps,1) 10/md.constants.yts*ones(md.nsteps,1)];
|
---|
[8304] | 20 | md.cm_min=[-2000*ones(md.numberofnodes,1) -2000*ones(md.numberofnodes,1)];
|
---|
| 21 | md.cm_max=[+2000*ones(md.numberofnodes,1) +2000*ones(md.numberofnodes,1)];
|
---|
[9541] | 22 | md.cm_responses=201*ones(md.nsteps,1);
|
---|
[8304] | 23 | md.weights=ones(md.numberofnodes,1);
|
---|
[6265] | 24 | md.maxiter=4*ones(md.nsteps,1);
|
---|
| 25 | md.cm_jump=0.99*ones(md.nsteps,1);
|
---|
| 26 |
|
---|
[8589] | 27 | md.cluster=generic('name',oshostname(),'np',3);
|
---|
[8295] | 28 | md=solve(md,BalancethicknessSolutionEnum);
|
---|
[6265] | 29 |
|
---|
| 30 | %Fields and tolerances to track changes
|
---|
| 31 | field_names ={'Gradient1' 'Gradient2' 'Misfits' 'Vx' 'Vy' 'Thickness'};
|
---|
| 32 | field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
|
---|
| 33 | field_values={...
|
---|
[8291] | 34 | PatchToVec(md.results.BalancethicknessSolution.Gradient1),...
|
---|
| 35 | PatchToVec(md.results.BalancethicknessSolution.Gradient2),...
|
---|
| 36 | md.results.BalancethicknessSolution.J,...
|
---|
| 37 | PatchToVec(md.results.BalancethicknessSolution.Vx),...
|
---|
| 38 | PatchToVec(md.results.BalancethicknessSolution.Vy),...
|
---|
| 39 | PatchToVec(md.results.BalancethicknessSolution.Thickness)
|
---|
[6265] | 40 | };
|
---|