source: issm/trunk-jpl/test/NightlyRun/test246.m@ 10976

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

removed all PatchToVec not needed anymore as default is md.settings.results_as_patches = 0

File size: 1.7 KB
Line 
1md=setmesh(model,'../Exp/Square.exp',150000);
2md=setmask(md,'all','');
3md=parameterize(md,'../Par/SquareShelf.par');
4md=setflowequation(md,'macayeal','all');
5md.geometry.bed=md.geometry.bed+50; md.geometry.surface=md.geometry.surface+50;
6md.cluster=generic('name',oshostname(),'np',3);
7md.prognostic.hydrostatic_adjustment='Incremental';
8md=solve(md,TransientSolutionEnum);
9
10%Fields and tolerances to track changes
11field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
12field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
13field_values={...
14 (md.results.TransientSolution(1).Vx),...
15 (md.results.TransientSolution(1).Vy),...
16 (md.results.TransientSolution(1).Vel),...
17 (md.results.TransientSolution(1).Pressure),...
18 (md.results.TransientSolution(1).Bed),...
19 (md.results.TransientSolution(1).Surface),...
20 (md.results.TransientSolution(1).Thickness),...
21 (md.results.TransientSolution(2).Vx),...
22 (md.results.TransientSolution(2).Vy),...
23 (md.results.TransientSolution(2).Vel),...
24 (md.results.TransientSolution(2).Pressure),...
25 (md.results.TransientSolution(2).Bed),...
26 (md.results.TransientSolution(2).Surface),...
27 (md.results.TransientSolution(2).Thickness),...
28 (md.results.TransientSolution(3).Vx),...
29 (md.results.TransientSolution(3).Vy),...
30 (md.results.TransientSolution(3).Vel),...
31 (md.results.TransientSolution(3).Pressure),...
32 (md.results.TransientSolution(3).Bed),...
33 (md.results.TransientSolution(3).Surface),...
34 (md.results.TransientSolution(3).Thickness),...
35 };
Note: See TracBrowser for help on using the repository browser.