Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test344.py
- Timestamp:
- 03/13/19 03:17:46 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test344.py
r22267 r23793 10 10 from SMBgradientsela import * 11 11 12 md = triangle(model(), '../Exp/Square.exp',150000.)13 md = setmask(md, '','')14 md = parameterize(md, '../Par/SquareSheetConstrained.py')12 md = triangle(model(), '../Exp/Square.exp', 150000.) 13 md = setmask(md, '', '') 14 md = parameterize(md, '../Par/SquareSheetConstrained.py') 15 15 16 16 #Change geometry … … 18 18 md.geometry.surface = md.geometry.base + md.geometry.thickness 19 19 20 md = md.extrude(3, 1.)21 md = setflowequation(md, 'HO','all')20 md = md.extrude(3, 1.) 21 md = setflowequation(md, 'HO', 'all') 22 22 md.smb = SMBgradientsela() 23 md.smb.ela = 1500. * np.ones((md.mesh.numberofvertices +1,))24 md.smb.b_pos = 0.002 * np.ones((md.mesh.numberofvertices +1,))25 md.smb.b_neg = 0.005 * np.ones((md.mesh.numberofvertices +1,))26 md.smb.b_max = 4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices +1,))27 md.smb.b_min = -4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices +1,))23 md.smb.ela = 1500. * np.ones((md.mesh.numberofvertices + 1,)) 24 md.smb.b_pos = 0.002 * np.ones((md.mesh.numberofvertices + 1,)) 25 md.smb.b_neg = 0.005 * np.ones((md.mesh.numberofvertices + 1,)) 26 md.smb.b_max = 4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices + 1,)) 27 md.smb.b_min = -4. * (md.materials.rho_freshwater / md.materials.rho_ice) * np.ones((md.mesh.numberofvertices + 1,)) 28 28 29 29 30 30 #Transient options 31 md.transient.requested_outputs = ['default', 'TotalSmb']32 md.cluster = generic('name', gethostname(),'np',3)33 md = solve(md, 'Transient')31 md.transient.requested_outputs = ['default', 'TotalSmb'] 32 md.cluster = generic('name', gethostname(), 'np', 3) 33 md = solve(md, 'Transient') 34 34 35 35 #Fields and tolerances to track changes 36 field_names = ['Vx1','Vy1','Vz1','Vel1','Bed1','Surface1','Thickness1','Temperature1','SMB1','TotalSmb1', 37 'Vx2','Vy2','Vz2','Vel2','Bed2','Surface2','Thickness2','Temperature2','SMB2','TotalSmb2', 38 'Vx3','Vy3','Vz3','Vel3','Bed3','Surface3','Thickness3','Temperature3','SMB3','TotalSmb3'] 39 field_tolerances = [1e-09,1e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10, 40 1e-09,1e-09,1e-10,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10, 41 1e-09,5e-09,1e-09,1e-09,1e-10,1e-10,1e-10,1e-10,1e-10,1e-10] 42 field_values = [ 43 md.results.TransientSolution[0].Vx, 44 md.results.TransientSolution[0].Vy, 45 md.results.TransientSolution[0].Vz, 46 md.results.TransientSolution[0].Vel, 47 md.results.TransientSolution[0].Base, 48 md.results.TransientSolution[0].Surface, 49 md.results.TransientSolution[0].Thickness, 50 md.results.TransientSolution[0].Temperature, 51 md.results.TransientSolution[0].SmbMassBalance, 52 md.results.TransientSolution[0].TotalSmb, 53 md.results.TransientSolution[1].Vx, 54 md.results.TransientSolution[1].Vy, 55 md.results.TransientSolution[1].Vz, 56 md.results.TransientSolution[1].Vel, 57 md.results.TransientSolution[1].Base, 58 md.results.TransientSolution[1].Surface, 59 md.results.TransientSolution[1].Thickness, 60 md.results.TransientSolution[1].Temperature, 61 md.results.TransientSolution[1].SmbMassBalance, 62 md.results.TransientSolution[1].TotalSmb, 63 md.results.TransientSolution[2].Vx, 64 md.results.TransientSolution[2].Vy, 65 md.results.TransientSolution[2].Vz, 66 md.results.TransientSolution[2].Vel, 67 md.results.TransientSolution[2].Base, 68 md.results.TransientSolution[2].Surface, 69 md.results.TransientSolution[2].Thickness, 70 md.results.TransientSolution[2].Temperature, 71 md.results.TransientSolution[2].SmbMassBalance, 72 md.results.TransientSolution[2].TotalSmb, 73 ] 36 field_names = ['Vx1', 'Vy1', 'Vz1', 'Vel1', 'Bed1', 'Surface1', 'Thickness1', 'Temperature1', 'SMB1', 'TotalSmb1', 37 'Vx2', 'Vy2', 'Vz2', 'Vel2', 'Bed2', 'Surface2', 'Thickness2', 'Temperature2', 'SMB2', 'TotalSmb2', 38 'Vx3', 'Vy3', 'Vz3', 'Vel3', 'Bed3', 'Surface3', 'Thickness3', 'Temperature3', 'SMB3', 'TotalSmb3'] 39 field_tolerances = [1e-09, 1e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 40 1e-09, 1e-09, 1e-10, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 41 1e-09, 5e-09, 1e-09, 1e-09, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10, 1e-10] 42 field_values = [md.results.TransientSolution[0].Vx, 43 md.results.TransientSolution[0].Vy, 44 md.results.TransientSolution[0].Vz, 45 md.results.TransientSolution[0].Vel, 46 md.results.TransientSolution[0].Base, 47 md.results.TransientSolution[0].Surface, 48 md.results.TransientSolution[0].Thickness, 49 md.results.TransientSolution[0].Temperature, 50 md.results.TransientSolution[0].SmbMassBalance, 51 md.results.TransientSolution[0].TotalSmb, 52 md.results.TransientSolution[1].Vx, 53 md.results.TransientSolution[1].Vy, 54 md.results.TransientSolution[1].Vz, 55 md.results.TransientSolution[1].Vel, 56 md.results.TransientSolution[1].Base, 57 md.results.TransientSolution[1].Surface, 58 md.results.TransientSolution[1].Thickness, 59 md.results.TransientSolution[1].Temperature, 60 md.results.TransientSolution[1].SmbMassBalance, 61 md.results.TransientSolution[1].TotalSmb, 62 md.results.TransientSolution[2].Vx, 63 md.results.TransientSolution[2].Vy, 64 md.results.TransientSolution[2].Vz, 65 md.results.TransientSolution[2].Vel, 66 md.results.TransientSolution[2].Base, 67 md.results.TransientSolution[2].Surface, 68 md.results.TransientSolution[2].Thickness, 69 md.results.TransientSolution[2].Temperature, 70 md.results.TransientSolution[2].SmbMassBalance, 71 md.results.TransientSolution[2].TotalSmb]
Note:
See TracChangeset
for help on using the changeset viewer.