source: issm/trunk-jpl/test/NightlyRun/test476.m@ 24149

Last change on this file since 24149 was 24149, checked in by Mathieu Morlighem, 5 years ago

CHG: bed commit, reverting
File size: 3.4 KB
Line 
1%Test Name: BeckmannGoosseMeltRate_HO
2md=triangle(model(),'../Exp/Square.exp',90000.);
3md=setmask(md,'../Exp/SquareShelf.exp','');
4md=parameterize(md,'../Par/SquareSheetShelf.par');
5md.initialization.vx(:)=1.;
6md.initialization.vy(:)=1.;
7md.geometry.thickness(:)=500-md.mesh.x/10000;
8md.geometry.bed =-100-md.mesh.x/1000;
9md.geometry.base=-md.geometry.thickness*md.materials.rho_ice/md.materials.rho_water;
10md.mask.groundedice_levelset=md.geometry.thickness+md.materials.rho_water/md.materials.rho_ice*md.geometry.bed;
11pos=find(md.mask.groundedice_levelset>=0);
12md.geometry.base(pos)=md.geometry.bed(pos);
13md.geometry.surface=md.geometry.base+md.geometry.thickness;
14md=extrude(md,3,1.1);
15md=setflowequation(md,'HO','all');
16
17%Set Pico Parameters
18md.basalforcings=basalforcingsbeckmanngoosse(md.basalforcings);
19md.basalforcings.ocean_temp=-1.7*ones(md.mesh.numberofvertices,1);
20md.basalforcings.ocean_salinity=35.0*ones(md.mesh.numberofvertices,1);
21md.basalforcings.meltrate_factor=1;
22
23%Boundary conditions:
24md.mask.ice_levelset=-ones(md.mesh.numberofvertices,1);
25md.mask.ice_levelset(find(md.mesh.x==max(md.mesh.x)))=0;
26
27%Model conditions
28md.transient.isthermal=0;
29md.transient.isstressbalance=1;
30md.transient.isgroundingline=1;
31md.transient.ismasstransport=1;
32md.transient.issmb=1;
33md.transient.requested_outputs={'default','BasalforcingsFloatingiceMeltingRate'};
34md.groundingline.migration='SubelementMigration';
35md.groundingline.friction_interpolation='SubelementFriction1';
36md.groundingline.melt_interpolation='SubelementMelt1';
37md.timestepping.final_time=1.5;
38md.timestepping.time_step=0.5;
39
40md.cluster=generic('name',oshostname(),'np',3);
41md=solve(md,'Transient');
42
43field_names ={'Bed1','Surface1','Thickness1','Floatingice1','Vx1','Vy1','Pressure1','FloatingiceMeltingrate1',...
44 'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Pressure2','FloatingiceMeltingrate2',...
45 'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Pressure3','FloatingiceMeltingrate3'};
46field_tolerances={7e-09,8e-09,8e-09,7e-09,6e-08,7e-08,6e-09,8e-7,...
47 7e-09,8e-09,8e-09,7e-09,6e-08,7e-08,6e-09,8e-7,...
48 7e-09,8e-09,8e-09,7e-09,6e-08,7e-08,6e-09,8e-7};
49field_values={...
50 (md.results.TransientSolution(1).Base),...
51 (md.results.TransientSolution(1).Surface),...
52 (md.results.TransientSolution(1).Thickness),...
53 (md.results.TransientSolution(1).MaskGroundediceLevelset),...
54 (md.results.TransientSolution(1).Vx),...
55 (md.results.TransientSolution(1).Vy),...
56 (md.results.TransientSolution(1).Pressure),...
57 (md.results.TransientSolution(1).BasalforcingsFloatingiceMeltingRate),...
58 (md.results.TransientSolution(2).Base),...
59 (md.results.TransientSolution(2).Surface),...
60 (md.results.TransientSolution(2).Thickness),...
61 (md.results.TransientSolution(2).MaskGroundediceLevelset),...
62 (md.results.TransientSolution(2).Vx),...
63 (md.results.TransientSolution(2).Vy),...
64 (md.results.TransientSolution(2).Pressure),...
65 (md.results.TransientSolution(2).BasalforcingsFloatingiceMeltingRate),...
66 (md.results.TransientSolution(3).Base),...
67 (md.results.TransientSolution(3).Surface),...
68 (md.results.TransientSolution(3).Thickness),...
69 (md.results.TransientSolution(3).MaskGroundediceLevelset),...
70 (md.results.TransientSolution(3).Vx),...
71 (md.results.TransientSolution(3).Vy),...
72 (md.results.TransientSolution(3).Pressure),...
73 (md.results.TransientSolution(3).BasalforcingsFloatingiceMeltingRate),...
74 };
Note: See TracBrowser for help on using the repository browser.