source: issm/trunk-jpl/test/NightlyRun/test121.py@ 15767

Last change on this file since 15767 was 15767, checked in by Mathieu Morlighem, 12 years ago

CHG: Prognostic is now Masstransport

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