![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 21408 for issm/trunk-jpl/test/NightlyRun/test329.py
- Timestamp:
- 11/22/16 02:31:19 (8 years ago)
- File:
-
- 1 edited
-
issm/trunk-jpl/test/NightlyRun/test329.py (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test329.py
r21060 r21408 1 1 #Test Name: SquareSheetConstrainedSmbGradients3d 2 import numpy3 2 import copy 4 3 from model import * 4 from socket import gethostname 5 5 from triangle import * 6 6 from setmask import * … … 8 8 from setflowequation import * 9 9 from solve import * 10 from MatlabFuncs import *11 10 12 11 md=triangle(model(),'../Exp/Square.exp',150000.) … … 15 14 md.extrude(3,1) 16 15 md=setflowequation(md,'HO','all') 17 md.smb = SMBgradients() ;16 md.smb = SMBgradients() 18 17 md.smb.b_pos=-100. + 0.00005*md.mesh.x - 0.0001*md.mesh.y 19 18 md.smb.b_neg=250. + 0.000051*md.mesh.x - 0.00011*md.mesh.y 20 md.smb.href=copy.deepcopy(md.geometry.surface) .reshape(-1)21 md.smb.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y ;19 md.smb.href=copy.deepcopy(md.geometry.surface) 20 md.smb.smbref= 1000. - 0.001*md.mesh.x - 0.005*md.mesh.y 22 21 md.transient.requested_outputs=['default','TotalSmb'] 23 md.cluster=generic('name', oshostname(),'np',3)22 md.cluster=generic('name',gethostname(),'np',3) 24 23 md=solve(md,'Transient') 25 24 26 25 #Fields and tolerances to track changes 27 field_names =['Vx1','Vy1','Vz1','Vel1','Bed1','Surface1','Thickness1','Temperature1','SMB1','TotalSmb1','Vx2','Vy2','Vz2','Vel2','Bed2','Surface2','Thickness2','Temperature2','SMB2','TotalSmb2','Vx3','Vy3','Vz3','Vel3','Bed3','Surface3','Thickness3','Temperature3','SMB3','TotalSmb3'] 28 field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ 29 1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ 30 1e-09,5e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] 31 field_values=[\ 32 md.results.TransientSolution[0].Vx,\ 33 md.results.TransientSolution[0].Vy,\ 34 md.results.TransientSolution[0].Vz,\ 35 md.results.TransientSolution[0].Vel,\ 36 md.results.TransientSolution[0].Base,\ 37 md.results.TransientSolution[0].Surface,\ 38 md.results.TransientSolution[0].Thickness,\ 39 md.results.TransientSolution[0].Temperature,\ 40 md.results.TransientSolution[0].SmbMassBalance,\ 41 md.results.TransientSolution[0].TotalSmb,\ 42 md.results.TransientSolution[1].Vx,\ 43 md.results.TransientSolution[1].Vy,\ 44 md.results.TransientSolution[1].Vz,\ 45 md.results.TransientSolution[1].Vel,\ 46 md.results.TransientSolution[1].Base,\ 47 md.results.TransientSolution[1].Surface,\ 48 md.results.TransientSolution[1].Thickness,\ 49 md.results.TransientSolution[1].Temperature,\ 50 md.results.TransientSolution[1].SmbMassBalance,\ 51 md.results.TransientSolution[1].TotalSmb,\ 52 md.results.TransientSolution[2].Vx,\ 53 md.results.TransientSolution[2].Vy,\ 54 md.results.TransientSolution[2].Vz,\ 55 md.results.TransientSolution[2].Vel,\ 56 md.results.TransientSolution[2].Base,\ 57 md.results.TransientSolution[2].Surface,\ 58 md.results.TransientSolution[2].Thickness,\ 59 md.results.TransientSolution[2].Temperature,\ 60 md.results.TransientSolution[2].SmbMassBalance,\ 61 md.results.TransientSolution[2].TotalSmb,\ 62 ] 26 field_names =['Vx1','Vy1','Vz1','Vel1', 27 'Bed1','Surface1','Thickness1', 28 'Temperature1','SMB1','TotalSmb1', 29 'Vx2','Vy2','Vz2','Vel2', 30 'Bed2','Surface2','Thickness2', 31 'Temperature2','SMB2','TotalSmb2', 32 'Vx3','Vy3','Vz3','Vel3', 33 'Bed3','Surface3','Thickness3', 34 'Temperature3','SMB3','TotalSmb3'] 35 field_tolerances=[1e-09,1e-09,1e-09,1e-09, 36 1e-10,1e-10,1e-10, 37 1e-10,1e-10,1e-10, 38 1e-09,1e-09,1e-10,1e-09, 39 1e-10,1e-10,1e-10, 40 1e-10,1e-10,1e-10, 41 1e-09,5e-09,1e-09,1e-09, 42 1e-10,1e-10,1e-10, 43 1e-10,1e-10,1e-10] 44 field_values=[md.results.TransientSolution[0].Vx, 45 md.results.TransientSolution[0].Vy, 46 md.results.TransientSolution[0].Vz, 47 md.results.TransientSolution[0].Vel, 48 md.results.TransientSolution[0].Base, 49 md.results.TransientSolution[0].Surface, 50 md.results.TransientSolution[0].Thickness, 51 md.results.TransientSolution[0].Temperature, 52 md.results.TransientSolution[0].SmbMassBalance, 53 md.results.TransientSolution[0].TotalSmb, 54 md.results.TransientSolution[1].Vx, 55 md.results.TransientSolution[1].Vy, 56 md.results.TransientSolution[1].Vz, 57 md.results.TransientSolution[1].Vel, 58 md.results.TransientSolution[1].Base, 59 md.results.TransientSolution[1].Surface, 60 md.results.TransientSolution[1].Thickness, 61 md.results.TransientSolution[1].Temperature, 62 md.results.TransientSolution[1].SmbMassBalance, 63 md.results.TransientSolution[1].TotalSmb, 64 md.results.TransientSolution[2].Vx, 65 md.results.TransientSolution[2].Vy, 66 md.results.TransientSolution[2].Vz, 67 md.results.TransientSolution[2].Vel, 68 md.results.TransientSolution[2].Base, 69 md.results.TransientSolution[2].Surface, 70 md.results.TransientSolution[2].Thickness, 71 md.results.TransientSolution[2].Temperature, 72 md.results.TransientSolution[2].SmbMassBalance, 73 md.results.TransientSolution[2].TotalSmb]
Note:
See TracChangeset
for help on using the changeset viewer.
