source: issm/trunk-jpl/test/NightlyRun/test227.m@ 13670

Last change on this file since 13670 was 13670, checked in by Mathieu Morlighem, 12 years ago

CHG: added point at the end of triangle command, resolution is a float

File size: 2.3 KB
Line 
1md=triangle(model(),'../Exp/Square.exp',200000.);
2md=setmask(md,'all','');
3md=parameterize(md,'../Par/SquareShelf.par');
4md=extrude(md,3,1);
5md=setflowequation(md,'pattyn','all');
6md.cluster=generic('name',oshostname(),'np',3);
7md.timestepping.time_adapt=1.;
8md.timestepping.final_time=10.;
9md=solve(md,TransientSolutionEnum());
10
11%Fields and tolerances to track changes
12field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsMeltingRate1', ...
13 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsMeltingRate2', ...
14 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3'};
15field_tolerances={...
16 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
17 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06,...
18 1e-09,1e-09,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05};
19field_values={...
20 (md.results.TransientSolution(1).Vx),...
21 (md.results.TransientSolution(1).Vy),...
22 (md.results.TransientSolution(1).Vz),...
23 (md.results.TransientSolution(1).Vel),...
24 (md.results.TransientSolution(1).Pressure),...
25 (md.results.TransientSolution(1).Bed),...
26 (md.results.TransientSolution(1).Surface),...
27 (md.results.TransientSolution(1).Thickness),...
28 (md.results.TransientSolution(1).Temperature),...
29 (md.results.TransientSolution(1).BasalforcingsMeltingRate),...
30 (md.results.TransientSolution(2).Vx),...
31 (md.results.TransientSolution(2).Vy),...
32 (md.results.TransientSolution(2).Vz),...
33 (md.results.TransientSolution(2).Vel),...
34 (md.results.TransientSolution(2).Pressure),...
35 (md.results.TransientSolution(2).Bed),...
36 (md.results.TransientSolution(2).Surface),...
37 (md.results.TransientSolution(2).Thickness),...
38 (md.results.TransientSolution(2).Temperature),...
39 (md.results.TransientSolution(2).BasalforcingsMeltingRate),...
40 (md.results.TransientSolution(3).Vx),...
41 (md.results.TransientSolution(3).Vy),...
42 (md.results.TransientSolution(3).Vz),...
43 (md.results.TransientSolution(3).Vel),...
44 (md.results.TransientSolution(3).Pressure),...
45 (md.results.TransientSolution(3).Bed),...
46 (md.results.TransientSolution(3).Surface),...
47 (md.results.TransientSolution(3).Thickness),...
48 (md.results.TransientSolution(3).Temperature),...
49 (md.results.TransientSolution(3).BasalforcingsMeltingRate),...
50 };
Note: See TracBrowser for help on using the repository browser.