source: issm/trunk/test/NightlyRun/test408.m@ 22758

Last change on this file since 22758 was 22758, checked in by Mathieu Morlighem, 7 years ago

merged trunk-jpl and trunk for revision 22757

File size: 3.3 KB
Line 
1%Test Name: SquareSheetShelfTranSSA3d
2md=triangle(model(),'../Exp/Square.exp',150000.);
3md=setmask(md,'../Exp/SquareShelf.exp','');
4md=parameterize(md,'../Par/SquareSheetShelf.par');
5md.geometry.bed=md.geometry.base;
6pos=find(md.mask.groundedice_levelset<0);
7md.geometry.bed(pos)=md.geometry.base(pos)-10;
8md.friction=frictioncoulomb();
9md.friction.coefficient=20.*ones(md.mesh.numberofvertices,1);
10md.friction.p=ones(md.mesh.numberofelements,1);
11md.friction.q=ones(md.mesh.numberofelements,1);
12md.friction.coefficientcoulomb=0.02*ones(md.mesh.numberofvertices,1);
13md.transient.isthermal=0;
14md.transient.isgroundingline=1;
15md=setflowequation(md,'SSA','all');
16md.cluster=generic('name',oshostname(),'np',3);
17md.transient.requested_outputs={'default','GroundedArea','FloatingArea','TotalFloatingBmb','TotalGroundedBmb','TotalSmb'};
18md=extrude(md,3,1);
19md=solve(md,'Transient');
20
21%Fields and tolerances to track changes
22field_names ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','GroundedArea1','FloatingArea1','TotalFloatingBmb1','TotalGroundedBmb1','TotalSmb1',...
23 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','GroundedArea2','FloatingArea2','TotalFloatingBmb2','TotalGroundedBmb2','TotalSmb2',...
24 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','GroundedArea3','FloatingArea3','TotalFloatingBmb3','TotalGroundedBmb3','TotalSmb3'};
25field_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,...
26 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,...
27 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};
28field_values={...
29 (md.results.TransientSolution(1).Vx),...
30 (md.results.TransientSolution(1).Vy),...
31 (md.results.TransientSolution(1).Vel),...
32 (md.results.TransientSolution(1).Pressure),...
33 (md.results.TransientSolution(1).Base),...
34 (md.results.TransientSolution(1).Surface),...
35 (md.results.TransientSolution(1).Thickness),...
36 (md.results.TransientSolution(1).GroundedArea),...
37 (md.results.TransientSolution(1).FloatingArea),...
38 (md.results.TransientSolution(1).TotalFloatingBmb),...
39 (md.results.TransientSolution(1).TotalGroundedBmb),...
40 (md.results.TransientSolution(1).TotalSmb),...
41 (md.results.TransientSolution(2).Vx),...
42 (md.results.TransientSolution(2).Vy),...
43 (md.results.TransientSolution(2).Vel),...
44 (md.results.TransientSolution(2).Pressure),...
45 (md.results.TransientSolution(2).Base),...
46 (md.results.TransientSolution(2).Surface),...
47 (md.results.TransientSolution(2).Thickness),...
48 (md.results.TransientSolution(2).GroundedArea),...
49 (md.results.TransientSolution(2).FloatingArea),...
50 (md.results.TransientSolution(2).TotalFloatingBmb),...
51 (md.results.TransientSolution(2).TotalGroundedBmb),...
52 (md.results.TransientSolution(2).TotalSmb),...
53 (md.results.TransientSolution(3).Vx),...
54 (md.results.TransientSolution(3).Vy),...
55 (md.results.TransientSolution(3).Vel),...
56 (md.results.TransientSolution(3).Pressure),...
57 (md.results.TransientSolution(3).Base),...
58 (md.results.TransientSolution(3).Surface),...
59 (md.results.TransientSolution(3).Thickness),...
60 (md.results.TransientSolution(3).GroundedArea),...
61 (md.results.TransientSolution(3).FloatingArea),...
62 (md.results.TransientSolution(3).TotalFloatingBmb),...
63 (md.results.TransientSolution(3).TotalGroundedBmb),...
64 (md.results.TransientSolution(3).TotalSmb),...
65 };
Note: See TracBrowser for help on using the repository browser.