1 | //Test Name: SquareShelfTranFS
|
---|
2 | var md = new model();
|
---|
3 | triangle(md,square[0],200000.);
|
---|
4 | setmask(md,'all','');
|
---|
5 | parameterize(md);
|
---|
6 | md.extrude(md,3,1.);
|
---|
7 | setflowequation(md,'FS','all');
|
---|
8 | md.stressbalance.reltol=null;
|
---|
9 | //md.cluster=generic('name',oshostname(),'np',3);
|
---|
10 | md=solve(md,'Transient');
|
---|
11 |
|
---|
12 | //Fields and tolerances to track changes
|
---|
13 | field_names=[
|
---|
14 | 'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1',
|
---|
15 | 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2',
|
---|
16 | 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3'];
|
---|
17 | field_tolerances=[
|
---|
18 | 2e-08,2e-08,5e-05,2e-08,1e-08,1e-08,1e-08,1e-08,1e-08,1e-08,
|
---|
19 | 5e-06,5e-06,8e-05,5e-06,5e-07,5e-07,5e-07,5e-07,3e-06,5e-05,
|
---|
20 | 8e-06,8e-06,8e-05,8e-06,8e-07,8e-07,8e-07,8e-07,5e-06,8e-05];
|
---|
21 | field_values=[
|
---|
22 | (md.results.TransientSolution[0](1).Vx),
|
---|
23 | (md.results.TransientSolution[0](1).Vy),
|
---|
24 | (md.results.TransientSolution[0](1).Vz),
|
---|
25 | (md.results.TransientSolution[0](1).Vel),
|
---|
26 | (md.results.TransientSolution[0](1).Pressure),
|
---|
27 | (md.results.TransientSolution[0](1).Base),
|
---|
28 | (md.results.TransientSolution[0](1).Surface),
|
---|
29 | (md.results.TransientSolution[0](1).Thickness),
|
---|
30 | (md.results.TransientSolution[0](1).Temperature),
|
---|
31 | (md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate),
|
---|
32 | (md.results.TransientSolution[0](2).Vx),
|
---|
33 | (md.results.TransientSolution[0](2).Vy),
|
---|
34 | (md.results.TransientSolution[0](2).Vz),
|
---|
35 | (md.results.TransientSolution[0](2).Vel),
|
---|
36 | (md.results.TransientSolution[0](2).Pressure),
|
---|
37 | (md.results.TransientSolution[0](2).Base),
|
---|
38 | (md.results.TransientSolution[0](2).Surface),
|
---|
39 | (md.results.TransientSolution[0](2).Thickness),
|
---|
40 | (md.results.TransientSolution[0](2).Temperature),
|
---|
41 | (md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate),
|
---|
42 | (md.results.TransientSolution[0](3).Vx),
|
---|
43 | (md.results.TransientSolution[0](3).Vy),
|
---|
44 | (md.results.TransientSolution[0](3).Vz),
|
---|
45 | (md.results.TransientSolution[0](3).Vel),
|
---|
46 | (md.results.TransientSolution[0](3).Pressure),
|
---|
47 | (md.results.TransientSolution[0](3).Base),
|
---|
48 | (md.results.TransientSolution[0](3).Surface),
|
---|
49 | (md.results.TransientSolution[0](3).Thickness),
|
---|
50 | (md.results.TransientSolution[0](3).Temperature),
|
---|
51 | (md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate),
|
---|
52 | ];
|
---|