![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 21729 for issm/trunk/test/NightlyRun/test231.py
- Timestamp:
- 05/19/17 14:48:02 (8 years ago)
- Location:
- issm/trunk
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk
-
issm/trunk/test
- Property svn:mergeinfo changed
-
issm/trunk/test/NightlyRun/test231.py
r21341 r21729 1 1 #Test Name: SquareShelfTranForcePos3d 2 import numpy 2 import numpy as np 3 3 from model import * 4 from socket import gethostname 4 5 from triangle import * 5 6 from setmask import * … … 7 8 from setflowequation import * 8 9 from solve import * 9 from MatlabFuncs import * 10 10 11 11 12 md=triangle(model(),'../Exp/Square.exp',350000.) … … 14 15 md.extrude(3,1.) 15 16 md=setflowequation(md,'SSA','all') 16 md.cluster=generic('name', oshostname(),'np',3)17 md.cluster=generic('name',gethostname(),'np',3) 17 18 18 19 md.timestepping.time_step=1. … … 21 22 22 23 #Set up transient 23 smb=n umpy.ones((md.mesh.numberofvertices,1))*3.624 smb=n umpy.hstack((smb,smb*2.))24 smb=np.ones((md.mesh.numberofvertices))*3.6 25 smb=np.vstack((smb,smb*2.)).T 25 26 26 md.smb.mass_balance=n umpy.vstack((smb,[1.5,3.]))27 md.smb.mass_balance=np.vstack((smb,[1.5,3.])) 27 28 md.transient.isthermal=False 28 29 … … 30 31 31 32 #Fields and tolerances to track changes 32 field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', \ 33 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', \ 34 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', \ 35 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassbalance4'] 36 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 37 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 38 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,\ 39 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13] 40 field_values=[\ 41 md.results.TransientSolution[0].Vx,\ 42 md.results.TransientSolution[0].Vy,\ 43 md.results.TransientSolution[0].Vz,\ 44 md.results.TransientSolution[0].Vel,\ 45 md.results.TransientSolution[0].Pressure,\ 46 md.results.TransientSolution[0].Base,\ 47 md.results.TransientSolution[0].Surface,\ 48 md.results.TransientSolution[0].Thickness,\ 49 md.results.TransientSolution[0].SmbMassBalance,\ 50 md.results.TransientSolution[1].Vx,\ 51 md.results.TransientSolution[1].Vy,\ 52 md.results.TransientSolution[1].Vz,\ 53 md.results.TransientSolution[1].Vel,\ 54 md.results.TransientSolution[1].Pressure,\ 55 md.results.TransientSolution[1].Base,\ 56 md.results.TransientSolution[1].Surface,\ 57 md.results.TransientSolution[1].Thickness,\ 58 md.results.TransientSolution[1].SmbMassBalance,\ 59 md.results.TransientSolution[2].Vx,\ 60 md.results.TransientSolution[2].Vy,\ 61 md.results.TransientSolution[2].Vz,\ 62 md.results.TransientSolution[2].Vel,\ 63 md.results.TransientSolution[2].Pressure,\ 64 md.results.TransientSolution[2].Base,\ 65 md.results.TransientSolution[2].Surface,\ 66 md.results.TransientSolution[2].Thickness,\ 67 md.results.TransientSolution[2].SmbMassBalance,\ 68 md.results.TransientSolution[3].Vx,\ 69 md.results.TransientSolution[3].Vy,\ 70 md.results.TransientSolution[3].Vz,\ 71 md.results.TransientSolution[3].Vel,\ 72 md.results.TransientSolution[3].Pressure,\ 73 md.results.TransientSolution[3].Base,\ 74 md.results.TransientSolution[3].Surface,\ 75 md.results.TransientSolution[3].Thickness,\ 76 md.results.TransientSolution[3].SmbMassBalance,\ 77 ] 33 field_names =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', 34 'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', 35 'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', 36 'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassbalance4'] 37 field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13, 38 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13, 39 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13, 40 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13] 41 field_values=[md.results.TransientSolution[0].Vx, 42 md.results.TransientSolution[0].Vy, 43 md.results.TransientSolution[0].Vz, 44 md.results.TransientSolution[0].Vel, 45 md.results.TransientSolution[0].Pressure, 46 md.results.TransientSolution[0].Base, 47 md.results.TransientSolution[0].Surface, 48 md.results.TransientSolution[0].Thickness, 49 md.results.TransientSolution[0].SmbMassBalance, 50 md.results.TransientSolution[1].Vx, 51 md.results.TransientSolution[1].Vy, 52 md.results.TransientSolution[1].Vz, 53 md.results.TransientSolution[1].Vel, 54 md.results.TransientSolution[1].Pressure, 55 md.results.TransientSolution[1].Base, 56 md.results.TransientSolution[1].Surface, 57 md.results.TransientSolution[1].Thickness, 58 md.results.TransientSolution[1].SmbMassBalance, 59 md.results.TransientSolution[2].Vx, 60 md.results.TransientSolution[2].Vy, 61 md.results.TransientSolution[2].Vz, 62 md.results.TransientSolution[2].Vel, 63 md.results.TransientSolution[2].Pressure, 64 md.results.TransientSolution[2].Base, 65 md.results.TransientSolution[2].Surface, 66 md.results.TransientSolution[2].Thickness, 67 md.results.TransientSolution[2].SmbMassBalance, 68 md.results.TransientSolution[3].Vx, 69 md.results.TransientSolution[3].Vy, 70 md.results.TransientSolution[3].Vz, 71 md.results.TransientSolution[3].Vel, 72 md.results.TransientSolution[3].Pressure, 73 md.results.TransientSolution[3].Base, 74 md.results.TransientSolution[3].Surface, 75 md.results.TransientSolution[3].Thickness, 76 md.results.TransientSolution[3].SmbMassBalance]
Note:
See TracChangeset
for help on using the changeset viewer.
