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

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

CHG: transferred Smb code from the Masstransport side of things into a full fledged SmbAnalysis
object.

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=True
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.