![(please configure the [header_logo] section in trac.ini)](/trac/issm/chrome/common/trac_banner.png)
Changeset 24313 for issm/trunk/test/NightlyRun/test241.py
- Timestamp:
- 11/01/19 12:01:57 (5 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/test241.py
r21729 r24313 10 10 11 11 12 md =triangle(model(),'../Exp/Square.exp',150000.)13 md =setmask(md,'all','')14 md =parameterize(md,'../Par/SquareShelf.py')15 md =setflowequation(md,'SSA','all')16 md.cluster =generic('name',gethostname(),'np',3)12 md = triangle(model(), '../Exp/Square.exp', 150000.) 13 md = setmask(md, 'all', '') 14 md = parameterize(md, '../Par/SquareShelf.py') 15 md = setflowequation(md, 'SSA', 'all') 16 md.cluster = generic('name', gethostname(), 'np', 3) 17 17 18 md.timestepping.time_step =1.19 md.settings.output_frequency =120 md.timestepping.final_time =4.21 md.timestepping.interp_forcings =False18 md.timestepping.time_step = 1. 19 md.settings.output_frequency = 1 20 md.timestepping.final_time = 4. 21 md.timestepping.interp_forcings = False 22 22 23 23 #Set up transient 24 smb =np.ones((md.mesh.numberofvertices))*3.625 smb =np.vstack((smb,smb*-1.)).T24 smb = np.ones((md.mesh.numberofvertices)) * 3.6 25 smb = np.vstack((smb, smb * - 1.)).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','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 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 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.
