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

Last change on this file since 19049 was 19049, checked in by Mathieu Morlighem, 10 years ago

NEW: test name is now included in the file itself so that people don't forget to include it

File size: 1.5 KB
RevLine 
[19049]1#Test Name: SquareShelfConstrainedEnthalpyTran
[13335]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')
[13672]14md.extrude(3,1.)
[15565]15md=setflowequation(md,'SSA','all')
[13335]16md.cluster=generic('name',oshostname(),'np',3);
[13463]17md.initialization.waterfraction=numpy.zeros((md.mesh.numberofvertices,1))
[16041]18md.initialization.watercolumn=numpy.zeros((md.mesh.numberofvertices,1))
[15771]19md.transient.isstressbalance=False
[15767]20md.transient.ismasstransport=False
[14021]21md.transient.isthermal=True
22md.transient.isgroundingline=False
[13335]23md.thermal.isenthalpy=1
[16322]24md.thermal.isdynamicbasalspc=1
[13335]25md=solve(md,TransientSolutionEnum())
26
27#Fields and tolerances to track changes
28field_names =['Enthalpy1','Waterfraction1','Temperature1',\
29 'Enthalpy2','Waterfraction2','Temperature2',\
30 'Enthalpy3','Waterfraction3','Temperature3']
[15460]31field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-9,1e-13]
[13335]32field_values=[\
[13967]33 md.results.TransientSolution[0].Enthalpy,\
34 md.results.TransientSolution[0].Waterfraction,\
35 md.results.TransientSolution[0].Temperature,\
36 md.results.TransientSolution[1].Enthalpy,\
37 md.results.TransientSolution[1].Waterfraction,\
38 md.results.TransientSolution[1].Temperature,\
39 md.results.TransientSolution[2].Enthalpy,\
40 md.results.TransientSolution[2].Waterfraction,\
41 md.results.TransientSolution[2].Temperature,\
[13335]42 ]
Note: See TracBrowser for help on using the repository browser.