source: issm/trunk-jpl/test/NightlyRun/test338.m@ 19527

Last change on this file since 19527 was 19527, checked in by Eric.Larour, 10 years ago

CHG: moved md.surfaceforcings to md.smb.
By doing so, had to rename the SMB class to SMBforcing class (it's just that, a mass_balance forcing inside
a SMB class, hence the name).
We also now have an smb_core solution, taken out of the mass transport core. Makes more sense long term.
Synced all enums according to the new changes, and operated the adjustments in all the test decks.

In addition, progressing in terms of GEMB integration into ISSM, specifically at the SMBgemb level (which
is spurring all the changes described above). Brought the class up to the level of the GEMB.m call in Alex's
code. Starting the C integration now.

File size: 2.7 KB
Line 
1%Test Name: SquareSheetConstrainedSmbMeltComponents2d
2md=triangle(model(),'../Exp/Square.exp',150000.);
3md=setmask(md,'all','');
4md=parameterize(md,'../Par/SquareShelf.par');
5md=setflowequation(md,'SSA','all');
6md.cluster=generic('name',oshostname(),'np',3);
7
8md.timestepping.time_step=1.;
9md.settings.output_frequency=1;
10md.timestepping.final_time=4.;
11
12%Set up transient
13smb=ones(md.mesh.numberofvertices,1)*3.6;
14smb=[ smb smb*-1. ];
15
16md.smb=SMBmeltcomponents();
17md.smb.accumulation=[smb; [1.5 3.]];
18md.smb.melt=[smb/2; [1.5 3.]];
19md.smb.refreeze=[smb; [1.5 3.]];
20md.smb.evaporation=[smb/2; [1.5 3.]];
21md.transient.isthermal=0;
22
23md=solve(md,TransientSolutionEnum());
24
25%Fields and tolerances to track changes
26field_names={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', ...
27 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', ...
28 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', ...
29 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassBalance4'};
30field_tolerances={1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
31 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
32 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,...
33 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10};
34field_values={...
35 (md.results.TransientSolution(1).Vx),...
36 (md.results.TransientSolution(1).Vy),...
37 (md.results.TransientSolution(1).Vel),...
38 (md.results.TransientSolution(1).Pressure),...
39 (md.results.TransientSolution(1).Base),...
40 (md.results.TransientSolution(1).Surface),...
41 (md.results.TransientSolution(1).Thickness),...
42 (md.results.TransientSolution(1).SmbMassBalance),...
43 (md.results.TransientSolution(2).Vx),...
44 (md.results.TransientSolution(2).Vy),...
45 (md.results.TransientSolution(2).Vel),...
46 (md.results.TransientSolution(2).Pressure),...
47 (md.results.TransientSolution(2).Base),...
48 (md.results.TransientSolution(2).Surface),...
49 (md.results.TransientSolution(2).Thickness),...
50 (md.results.TransientSolution(2).SmbMassBalance),...
51 (md.results.TransientSolution(3).Vx),...
52 (md.results.TransientSolution(3).Vy),...
53 (md.results.TransientSolution(3).Vel),...
54 (md.results.TransientSolution(3).Pressure),...
55 (md.results.TransientSolution(3).Base),...
56 (md.results.TransientSolution(3).Surface),...
57 (md.results.TransientSolution(3).Thickness),...
58 (md.results.TransientSolution(3).SmbMassBalance),...
59 (md.results.TransientSolution(4).Vx),...
60 (md.results.TransientSolution(4).Vy),...
61 (md.results.TransientSolution(4).Vel),...
62 (md.results.TransientSolution(4).Pressure),...
63 (md.results.TransientSolution(4).Base),...
64 (md.results.TransientSolution(4).Surface),...
65 (md.results.TransientSolution(4).Thickness),...
66 (md.results.TransientSolution(4).SmbMassBalance),...
67 };
Note: See TracBrowser for help on using the repository browser.