source: issm/trunk/test/NightlyRun/test622.m@ 9681

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

Added inversion

File size: 1.3 KB
RevLine 
[5668]1md=mesh(model,'../Exp/79North.exp',10000);
[9641]2md=setmask(md,'../Exp/79NorthShelf.exp','');
[5668]3md=parameterize(md,'../Par/79North.par');
[9664]4md=setflowequation(md,'macayeal','all');
[5668]5
6%control parameters
[9681]7md.inversion.nsteps=2;
[9642]8md.prognostic.stabilization=1;
[9681]9md.inversion.iscontrol=1;
10md.inversion.control_parameters={'BalancethicknessThickeningRate'};
11md.inversion.thickness_obs=md.thickness;
12md.inversion.min_parameters=-50*ones(md.numberofnodes,1);
13md.inversion.max_parameters=50*ones(md.numberofnodes,1);
14md.inversion.cost_functions=201*ones(md.inversion.nsteps,1);
15md.inversion.cost_functions_coefficients=ones(md.numberofnodes,1);
16md.inversion.gradient_scaling=10/md.constants.yts*ones(md.inversion.nsteps,1);
17md.inversion.maxiter_per_step=4*ones(md.inversion.nsteps,1);
18md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
[5668]19
[8589]20md.cluster=generic('name',oshostname(),'np',3);
[8295]21md=solve(md,BalancethicknessSolutionEnum);
[5668]22
23%Fields and tolerances to track changes
[9660]24field_names ={'Gradient' 'Misfits' 'BalancethicknessThickeningRate' 'Thickness'};
[5668]25field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-12,1e-12};
26field_values={...
[8291]27 PatchToVec(md.results.BalancethicknessSolution.Gradient1),...
28 md.results.BalancethicknessSolution.J,...
[9660]29 PatchToVec(md.results.BalancethicknessSolution.BalancethicknessThickeningRate),...
[8291]30 PatchToVec(md.results.BalancethicknessSolution.Thickness)
[5668]31};
Note: See TracBrowser for help on using the repository browser.