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

Last change on this file since 16041 was 16041, checked in by seroussi, 12 years ago

BUG: missing watercolumn initialization

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.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
18md.transient.isstressbalance=False
19md.transient.ismasstransport=False
20md.transient.isthermal=True
21md.transient.isgroundingline=False
22md.thermal.isenthalpy=1
23md=solve(md,TransientSolutionEnum())
24
25#Fields and tolerances to track changes
26field_names =['Enthalpy1','Waterfraction1','Temperature1',\
27 'Enthalpy2','Waterfraction2','Temperature2',\
28 'Enthalpy3','Waterfraction3','Temperature3']
29field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-9,1e-13]
30field_values=[\
31 md.results.TransientSolution[0].Enthalpy,\
32 md.results.TransientSolution[0].Waterfraction,\
33 md.results.TransientSolution[0].Temperature,\
34 md.results.TransientSolution[1].Enthalpy,\
35 md.results.TransientSolution[1].Waterfraction,\
36 md.results.TransientSolution[1].Temperature,\
37 md.results.TransientSolution[2].Enthalpy,\
38 md.results.TransientSolution[2].Waterfraction,\
39 md.results.TransientSolution[2].Temperature,\
40 ]
Note: See TracBrowser for help on using the repository browser.