source: issm/trunk-jpl/test/NightlyRun/test121.py@ 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: 1.5 KB
Line 
1#Test Name: SquareShelfConstrainedEnthalpyTran
2from model import *
3from triangle import *
4from setmask import *
5from parameterize import *
6from setflowequation import *
7from EnumDefinitions import *
8from solve import *
9from MatlabFuncs import *
10
11md=triangle(model(),'../Exp/Square.exp',180000)
12md=setmask(md,'all','')
13md=parameterize(md,'../Par/SquareShelfConstrained.py')
14md.extrude(3,1.)
15md=setflowequation(md,'SSA','all')
16md.cluster=generic('name',oshostname(),'np',3);
17md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
18md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
19md.transient.isstressbalance=False
20md.transient.ismasstransport=False
21md.transient.issmb=False
22md.transient.isthermal=True
23md.transient.isgroundingline=False
24md.thermal.isenthalpy=1
25md.thermal.isdynamicbasalspc=1
26md=solve(md,TransientSolutionEnum())
27
28#Fields and tolerances to track changes
29field_names =['Enthalpy1','Waterfraction1','Temperature1',\
30 'Enthalpy2','Waterfraction2','Temperature2',\
31 'Enthalpy3','Waterfraction3','Temperature3']
32field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-9,1e-13]
33field_values=[\
34 md.results.TransientSolution[0].Enthalpy,\
35 md.results.TransientSolution[0].Waterfraction,\
36 md.results.TransientSolution[0].Temperature,\
37 md.results.TransientSolution[1].Enthalpy,\
38 md.results.TransientSolution[1].Waterfraction,\
39 md.results.TransientSolution[1].Temperature,\
40 md.results.TransientSolution[2].Enthalpy,\
41 md.results.TransientSolution[2].Waterfraction,\
42 md.results.TransientSolution[2].Temperature,\
43 ]
Note: See TracBrowser for help on using the repository browser.