source: issm/trunk-jpl/test/NightlyRun/test525.m@ 11014

Last change on this file since 11014 was 11014, checked in by Mathieu Morlighem, 13 years ago

moved setmesh to triangle in all nightly runs

File size: 1.7 KB
RevLine 
[11014]1md=triangle(model,'../Exp/Pig.exp',30000);
[9641]2md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
[5052]3md=parameterize(md,'../Par/Pig.par');
4md=extrude(md,3,1);
[9664]5md=setflowequation(md,'macayeal','all');
[5052]6
7%control parameters
[9681]8md.inversion.iscontrol=1;
9md.inversion.control_parameters={'FrictionCoefficient'};
[9725]10md.inversion.min_parameters=1*ones(md.mesh.numberofvertices,1);
11md.inversion.max_parameters=200*ones(md.mesh.numberofvertices,1);
[9681]12md.inversion.nsteps=2;
13md.inversion.cost_functions=[103*ones(md.inversion.nsteps,1) 501*ones(md.inversion.nsteps,1)];
[9725]14md.inversion.cost_functions_coefficients=ones(md.mesh.numberofvertices,2); md.inversion.cost_functions_coefficients(:,2)=2*10^-7;
[9681]15md.inversion.gradient_scaling=3*ones(md.inversion.nsteps,1);
16md.inversion.maxiter_per_step=2*ones(md.inversion.nsteps,1);
17md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
[9684]18md.inversion.vx_obs=md.initialization.vx; md.inversion.vy_obs=md.initialization.vy;
[9628]19md.timestepping.time_step=0;
[5052]20
[5955]21md.cluster=none;
[8295]22md=solve(md,SteadystateSolutionEnum);
[5098]23
24%Fields and tolerances to track changes
[9611]25field_names ={'Gradient' 'Misfits' 'FrictionCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy' 'Vz' 'Temperature' 'BasalforcingsMeltingRate'};
[5098]26field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
27field_values={...
[10976]28 (md.results.SteadystateSolution.Gradient1),...
[5098]29 md.results.SteadystateSolution.J,...
[10976]30 (md.results.SteadystateSolution.FrictionCoefficient),...
31 (md.results.SteadystateSolution.Pressure),...
32 (md.results.SteadystateSolution.Vel),...
33 (md.results.SteadystateSolution.Vx),...
34 (md.results.SteadystateSolution.Vy),...
35 (md.results.SteadystateSolution.Vz),...
36 (md.results.SteadystateSolution.Temperature),...
37 (md.results.SteadystateSolution.BasalforcingsMeltingRate)
[5098]38};
Note: See TracBrowser for help on using the repository browser.