source: issm/trunk-jpl/test/NightlyRun/test803.py@ 21049

Last change on this file since 21049 was 21049, checked in by agscott1, 9 years ago

CHG: Replaced Enums with Strings in matlab and python. Updated corresponding cpp code.

File size: 2.6 KB
RevLine 
[19049]1#Test Name: ValleyGlacierLevelsetEnthalpyHO3d
[17452]2from model import *
3from triangle import *
4from squaremesh import *
5from setmask import *
6from parameterize import *
7from setflowequation import *
8from EnumDefinitions import *
9from solve import *
10from MatlabFuncs import *
11
12md=triangle(model(),'../Exp/Square.exp',50000)
13md=setmask(md,'','')
14md=parameterize(md,'../Par/ValleyGlacierShelf.py')
15md.extrude(3,2.)
16md=setflowequation(md,'HO','all')
17md.cluster=generic('name',oshostname(),'np',3)
18
19#Thermal model
20pos_surf=numpy.nonzero(md.mesh.vertexonsurface)[0]
21md.thermal.spctemperature[pos_surf]=md.initialization.temperature[pos_surf]
22md.thermal.isenthalpy=True
[17470]23md.thermal.isdynamicbasalspc=True
[17452]24
25#Transient
26md.transient.isstressbalance=True
[20452]27md.transient.ismovingfront=True
[17452]28md.transient.ismasstransport=True
[19527]29md.transient.issmb=True
[17452]30md.transient.isthermal=True
31md.transient.isgroundingline=True
32md.transient.isgia=False
33
[21049]34md=solve(md,'TransientSolution')
[17452]35
36#Fields and tolerances to track changes
37field_names =['Vx1','Vy1','Vel1','Pressure1','Thickness1','Surface1','MaskIceLevelset1','Enthalpy1','Watercolumn1',\
38 'Vx2','Vy2','Vel2','Pressure2','Thickness2','Surface2','MaskIceLevelset2','Enthalpy2','Watercolumn2',\
39 'Vx3','Vy3','Vel3','Pressure3','Thickness3','Surface3','MaskIceLevelset3','Enthalpy3','Watercolumn3']
[17478]40field_tolerances=[1e-10,1e-10,1e-10,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\
[20119]41 1e-9,1e-9,1e-10,1e-10,1e-10,1e-10,1e-10,1e-9,1e-10,\
[17452]42 1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10]
43field_values=[\
44 md.results.TransientSolution[0].Vx,\
45 md.results.TransientSolution[0].Vy,\
46 md.results.TransientSolution[0].Vel,\
47 md.results.TransientSolution[0].Pressure,\
48 md.results.TransientSolution[0].Thickness,\
49 md.results.TransientSolution[0].Surface,\
50 md.results.TransientSolution[0].MaskIceLevelset,\
51 md.results.TransientSolution[0].Enthalpy,\
52 md.results.TransientSolution[0].Watercolumn,\
53 md.results.TransientSolution[1].Vx,\
54 md.results.TransientSolution[1].Vy,\
55 md.results.TransientSolution[1].Vel,\
56 md.results.TransientSolution[1].Pressure,\
57 md.results.TransientSolution[1].Thickness,\
58 md.results.TransientSolution[1].Surface,\
59 md.results.TransientSolution[1].MaskIceLevelset,\
60 md.results.TransientSolution[1].Enthalpy,\
61 md.results.TransientSolution[1].Watercolumn,\
62 md.results.TransientSolution[2].Vx,\
63 md.results.TransientSolution[2].Vy,\
64 md.results.TransientSolution[2].Vel,\
65 md.results.TransientSolution[2].Pressure,\
66 md.results.TransientSolution[2].Thickness,\
67 md.results.TransientSolution[2].Surface,\
68 md.results.TransientSolution[2].MaskIceLevelset,\
69 md.results.TransientSolution[2].Enthalpy,\
70 md.results.TransientSolution[2].Watercolumn,\
71 ]
Note: See TracBrowser for help on using the repository browser.