![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test231.py
- Timestamp:
- 03/13/19 03:17:46 (6 years ago)
- File:
-
- 1 edited
-
issm/trunk-jpl/test/NightlyRun/test231.py (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test231.py
r21408 r23793 10 10 11 11 12 md =triangle(model(),'../Exp/Square.exp',350000.)13 md =setmask(md,'all','')14 md =parameterize(md,'../Par/SquareShelf.py')15 md.extrude(3, 1.)16 md =setflowequation(md,'SSA','all')17 md.cluster =generic('name',gethostname(),'np',3)12 md = triangle(model(), '../Exp/Square.exp', 350000.) 13 md = setmask(md, 'all', '') 14 md = parameterize(md, '../Par/SquareShelf.py') 15 md.extrude(3, 1.) 16 md = setflowequation(md, 'SSA', 'all') 17 md.cluster = generic('name', gethostname(), 'np', 3) 18 18 19 md.timestepping.time_step =1.20 md.settings.output_frequency =121 md.timestepping.final_time =4.19 md.timestepping.time_step = 1. 20 md.settings.output_frequency = 1 21 md.timestepping.final_time = 4. 22 22 23 23 #Set up transient 24 smb =np.ones((md.mesh.numberofvertices))*3.625 smb =np.vstack((smb,smb*2.)).T24 smb = np.ones((md.mesh.numberofvertices)) * 3.6 25 smb = np.vstack((smb, smb * 2.)).T 26 26 27 md.smb.mass_balance =np.vstack((smb,[1.5,3.]))28 md.transient.isthermal =False27 md.smb.mass_balance = np.vstack((smb, [1.5, 3.])) 28 md.transient.isthermal = False 29 29 30 md =solve(md,'Transient')30 md = solve(md, 'Transient') 31 31 32 32 #Fields and tolerances to track changes 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 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 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.
