source: issm/trunk-jpl/test/NightlyRun/test122.m@ 17555

Last change on this file since 17555 was 17555, checked in by Mathieu Morlighem, 11 years ago

NEW: renamed Bed -> Base in src/c (matlab to be done)

File size: 2.6 KB
Line 
1md=triangle(model(),'../Exp/Square.exp',200000.);
2md=setmask(md,'all','');
3md=parameterize(md,'../Par/SquareShelfConstrained.par');
4md=extrude(md,3,1.);
5md=setflowequation(md,'HO','all');
6md.initialization.waterfraction=zeros(md.mesh.numberofvertices,1);
7md.initialization.watercolumn=zeros(md.mesh.numberofvertices,1);
8md.thermal.isenthalpy=1;
9md.thermal.isdynamicbasalspc=1;
10md.thermal.stabilization=2;
11md.cluster=generic('name',oshostname(),'np',3);
12md=solve(md,TransientSolutionEnum());
13
14%Fields and tolerances to track changes
15field_names ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','Enthalpy1','Waterfraction1', ...
16 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','Enthalpy2','Waterfraction2', ...
17 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','Enthalpy3','Waterfraction3'};
18field_tolerances={...
19 1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
20 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,...
21 1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09};
22field_values={...
23 (md.results.TransientSolution(1).Vx),...
24 (md.results.TransientSolution(1).Vy),...
25 (md.results.TransientSolution(1).Vz),...
26 (md.results.TransientSolution(1).Vel),...
27 (md.results.TransientSolution(1).Pressure),...
28 (md.results.TransientSolution(1).Base),...
29 (md.results.TransientSolution(1).Surface),...
30 (md.results.TransientSolution(1).Thickness),...
31 (md.results.TransientSolution(1).Temperature),...
32 (md.results.TransientSolution(1).Enthalpy),...
33 (md.results.TransientSolution(1).Waterfraction),...
34 (md.results.TransientSolution(2).Vx),...
35 (md.results.TransientSolution(2).Vy),...
36 (md.results.TransientSolution(2).Vz),...
37 (md.results.TransientSolution(2).Vel),...
38 (md.results.TransientSolution(2).Pressure),...
39 (md.results.TransientSolution(2).Base),...
40 (md.results.TransientSolution(2).Surface),...
41 (md.results.TransientSolution(2).Thickness),...
42 (md.results.TransientSolution(2).Temperature),...
43 (md.results.TransientSolution(2).Enthalpy),...
44 (md.results.TransientSolution(2).Waterfraction),...
45 (md.results.TransientSolution(3).Vx),...
46 (md.results.TransientSolution(3).Vy),...
47 (md.results.TransientSolution(3).Vz),...
48 (md.results.TransientSolution(3).Vel),...
49 (md.results.TransientSolution(3).Pressure),...
50 (md.results.TransientSolution(3).Base),...
51 (md.results.TransientSolution(3).Surface),...
52 (md.results.TransientSolution(3).Thickness),...
53 (md.results.TransientSolution(3).Temperature),...
54 (md.results.TransientSolution(3).Enthalpy),...
55 (md.results.TransientSolution(3).Waterfraction),...
56 };
Note: See TracBrowser for help on using the repository browser.