source: issm/trunk/test/NightlyRun/test812.m@ 27035

Last change on this file since 27035 was 27035, checked in by Mathieu Morlighem, 3 years ago

merged trunk-jpl and trunk for revision 27033
File size: 3.3 KB
Line 
1%Test Name: SquareShelfLevelsetCalvingMOLHO2dLevermann
2md=triangle(model(),'../Exp/Square.exp',50000.);
3md=setmask(md,'all','');
4md=parameterize(md,'../Par/SquareShelf.par');
5md=setflowequation(md,'MOLHO','all');
6md.cluster=generic('name',oshostname(),'np',3);
7
8%Do not kill ice bergs as all is floating
9md.levelset.kill_icebergs=0;
10
11x = md.mesh.x;
12xmin = min(x);
13xmax = max(x);
14Lx = (xmax-xmin);
15alpha = 2./3.;
16md.mask.ice_levelset = ((x - alpha*Lx)>0) - ((x - alpha*Lx)<0);
17
18md.timestepping.time_step=10;
19md.timestepping.final_time=30;
20
21%Transient
22md.transient.isstressbalance=1;
23md.transient.ismasstransport=1;
24md.transient.issmb=1;
25md.transient.isthermal=0;
26md.transient.isgroundingline=0;
27md.transient.ismovingfront=1;
28
29md.calving=calvinglevermann();
30md.calving.coeff=4.89e13*ones(md.mesh.numberofvertices,1);
31md.frontalforcings.meltingrate=zeros(md.mesh.numberofvertices,1);
32md.levelset.spclevelset=NaN(md.mesh.numberofvertices,1);
33md.levelset.migration_max = 1e8;
34
35md.transient.requested_outputs={'default','StrainRateparallel','StrainRateperpendicular','Calvingratex','Calvingratey','CalvingCalvingrate'};
36
37md=SetMOLHOBC(md);
38md=solve(md,'Transient');
39
40%Fields and tolerances to track changes
41field_names ={'Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1','StrainRateparallel1','StrainRateperpendicular1','CalvingCalvingrate1'...
42 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2','StrainRateparallel2','StrainRateperpendicular2','CalvingCalvingrate2'...
43 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','StrainRateparallel3','StrainRateperpendicular3','CalvingCalvingrate3'};
44field_tolerances={1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,...
45 2e-11,2e-11,2e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,1e-11,...
46 2e-11,2e-11,2e-11,1e-11,1e-11,1e-11,1e-11,5e-11,5e-11,5e-11};
47field_values={...
48 md.results.TransientSolution(1).Vx,...
49 md.results.TransientSolution(1).Vy,...
50 md.results.TransientSolution(1).Vel,...
51 md.results.TransientSolution(1).Pressure,...
52 md.results.TransientSolution(1).Thickness,...
53 md.results.TransientSolution(1).Surface,...
54 md.results.TransientSolution(1).MaskIceLevelset,...
55 md.results.TransientSolution(1).StrainRateparallel,...
56 md.results.TransientSolution(1).StrainRateperpendicular,...
57 md.results.TransientSolution(1).CalvingCalvingrate,...
58 md.results.TransientSolution(2).Vx,...
59 md.results.TransientSolution(2).Vy,...
60 md.results.TransientSolution(2).Vel,...
61 md.results.TransientSolution(2).Pressure,...
62 md.results.TransientSolution(2).Thickness,...
63 md.results.TransientSolution(2).Surface,...
64 md.results.TransientSolution(2).MaskIceLevelset,...
65 md.results.TransientSolution(2).StrainRateparallel,...
66 md.results.TransientSolution(2).StrainRateperpendicular,...
67 md.results.TransientSolution(2).CalvingCalvingrate,...
68 md.results.TransientSolution(3).Vx,...
69 md.results.TransientSolution(3).Vy,...
70 md.results.TransientSolution(3).Vel,...
71 md.results.TransientSolution(3).Pressure,...
72 md.results.TransientSolution(3).Thickness,...
73 md.results.TransientSolution(3).Surface,...
74 md.results.TransientSolution(3).MaskIceLevelset,...
75 md.results.TransientSolution(3).StrainRateparallel,...
76 md.results.TransientSolution(3).StrainRateperpendicular,...
77 md.results.TransientSolution(3).CalvingCalvingrate,...
78 };
Note: See TracBrowser for help on using the repository browser.