![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
- Timestamp:
- 06/07/17 10:50:54 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/branches/trunk-larour-NatGeoScience2016/test/NightlyRun/test241.py
r21060 r21759 1 1 #Test Name: SquareShelfTranForceNoInterp2d 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',150000.) … … 13 14 md=parameterize(md,'../Par/SquareShelf.py') 14 15 md=setflowequation(md,'SSA','all') 15 md.cluster=generic('name', oshostname(),'np',3)16 md.cluster=generic('name',gethostname(),'np',3) 16 17 17 18 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*-1.))24 smb=np.ones((md.mesh.numberofvertices))*3.6 25 smb=np.vstack((smb,smb*-1.)).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','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', \ 33 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', \ 34 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', \ 35 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassBalance4'] 36 field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ 37 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ 38 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10,\ 39 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] 40 field_values=[\ 41 md.results.TransientSolution[0].Vx,\ 42 md.results.TransientSolution[0].Vy,\ 43 md.results.TransientSolution[0].Vel,\ 44 md.results.TransientSolution[0].Pressure,\ 45 md.results.TransientSolution[0].Base,\ 46 md.results.TransientSolution[0].Surface,\ 47 md.results.TransientSolution[0].Thickness,\ 48 md.results.TransientSolution[0].SmbMassBalance,\ 49 md.results.TransientSolution[1].Vx,\ 50 md.results.TransientSolution[1].Vy,\ 51 md.results.TransientSolution[1].Vel,\ 52 md.results.TransientSolution[1].Pressure,\ 53 md.results.TransientSolution[1].Base,\ 54 md.results.TransientSolution[1].Surface,\ 55 md.results.TransientSolution[1].Thickness,\ 56 md.results.TransientSolution[1].SmbMassBalance,\ 57 md.results.TransientSolution[2].Vx,\ 58 md.results.TransientSolution[2].Vy,\ 59 md.results.TransientSolution[2].Vel,\ 60 md.results.TransientSolution[2].Pressure,\ 61 md.results.TransientSolution[2].Base,\ 62 md.results.TransientSolution[2].Surface,\ 63 md.results.TransientSolution[2].Thickness,\ 64 md.results.TransientSolution[2].SmbMassBalance,\ 65 md.results.TransientSolution[3].Vx,\ 66 md.results.TransientSolution[3].Vy,\ 67 md.results.TransientSolution[3].Vel,\ 68 md.results.TransientSolution[3].Pressure,\ 69 md.results.TransientSolution[3].Base,\ 70 md.results.TransientSolution[3].Surface,\ 71 md.results.TransientSolution[3].Thickness,\ 72 md.results.TransientSolution[3].SmbMassBalance,\ 73 ] 33 field_names=['Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','SmbMassBalance1', 34 'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','SmbMassBalance2', 35 'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','SmbMassBalance3', 36 'Vx4','Vy4','Vel4','Pressure4','Bed4','Surface4','Thickness4','SmbMassBalance4'] 37 field_tolerances=[1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10, 38 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10, 39 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10, 40 1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] 41 field_values=[md.results.TransientSolution[0].Vx, 42 md.results.TransientSolution[0].Vy, 43 md.results.TransientSolution[0].Vel, 44 md.results.TransientSolution[0].Pressure, 45 md.results.TransientSolution[0].Base, 46 md.results.TransientSolution[0].Surface, 47 md.results.TransientSolution[0].Thickness, 48 md.results.TransientSolution[0].SmbMassBalance, 49 md.results.TransientSolution[1].Vx, 50 md.results.TransientSolution[1].Vy, 51 md.results.TransientSolution[1].Vel, 52 md.results.TransientSolution[1].Pressure, 53 md.results.TransientSolution[1].Base, 54 md.results.TransientSolution[1].Surface, 55 md.results.TransientSolution[1].Thickness, 56 md.results.TransientSolution[1].SmbMassBalance, 57 md.results.TransientSolution[2].Vx, 58 md.results.TransientSolution[2].Vy, 59 md.results.TransientSolution[2].Vel, 60 md.results.TransientSolution[2].Pressure, 61 md.results.TransientSolution[2].Base, 62 md.results.TransientSolution[2].Surface, 63 md.results.TransientSolution[2].Thickness, 64 md.results.TransientSolution[2].SmbMassBalance, 65 md.results.TransientSolution[3].Vx, 66 md.results.TransientSolution[3].Vy, 67 md.results.TransientSolution[3].Vel, 68 md.results.TransientSolution[3].Pressure, 69 md.results.TransientSolution[3].Base, 70 md.results.TransientSolution[3].Surface, 71 md.results.TransientSolution[3].Thickness, 72 md.results.TransientSolution[3].SmbMassBalance]
Note:
See TracChangeset
for help on using the changeset viewer.
