source: issm/trunk/test/NightlyRun/test339.m@ 5279

Last change on this file since 5279 was 5279, checked in by Mathieu Morlighem, 15 years ago

moved md.fit to md.cm_responses and 0-4 to enums

File size: 1.2 KB
Line 
1md=mesh(model,'../Exp/Square.exp',200000);
2md=geography(md,'','');
3md=parameterize(md,'../Par/SquareSheetConstrained.par');
4md=extrude(md,3,1);
5md=setelementstype(md,'macayeal','all');
6
7%control parameters
8md.control_analysis=1;
9md.control_type='drag_coefficient';
10md.cm_min=1;
11md.cm_max=200;
12md.nsteps=2;
13md.cm_responses=SurfaceLogVelMisfitEnum()*ones(md.nsteps,1);
14md.weights=ones(md.numberofgrids,1);
15md.optscal=3*ones(md.nsteps,1);
16md.maxiter=2*ones(md.nsteps,1);
17md.cm_jump=0.3*ones(md.nsteps,1);
18md.cm_noisedmp=2*10^-7;
19md.vx_obs=md.vx; md.vy_obs=md.vy;
20
21md.cluster='none';
22md=solve(md,'analysis_type',DiagnosticSolutionEnum);
23
24%Fields and tolerances to track changes
25field_names ={'Gradient' 'Misfits' 'DragCoefficient' 'Pressure' 'Vel' 'Vx' 'Vy'};
26field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
27field_values={...
28 PatchToVec(md.results.DiagnosticSolution.Gradient),...
29 md.results.DiagnosticSolution.J,...
30 PatchToVec(md.results.DiagnosticSolution.DragCoefficient),...
31 PatchToVec(md.results.DiagnosticSolution.Pressure),...
32 PatchToVec(md.results.DiagnosticSolution.Vel),...
33 PatchToVec(md.results.DiagnosticSolution.Vx),...
34 PatchToVec(md.results.DiagnosticSolution.Vy)
35};
Note: See TracBrowser for help on using the repository browser.