Changeset 23793 for issm/trunk-jpl/test/NightlyRun/test442.py
- Timestamp:
- 03/13/19 03:17:46 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test442.py
r23004 r23793 9 9 from solve import * 10 10 11 md = triangle(model(), '../Exp/Square.exp',100000.)12 md = setmask(md, '../Exp/SquareShelf.exp','')13 md = parameterize(md, '../Par/SquareSheetShelf.py')11 md = triangle(model(), '../Exp/Square.exp', 100000.) 12 md = setmask(md, '../Exp/SquareShelf.exp', '') 13 md = parameterize(md, '../Par/SquareSheetShelf.py') 14 14 md.initialization.vx[:] = 1. 15 15 md.initialization.vy[:] = 1. … … 21 21 md.geometry.base[pos] = md.geometry.bed[pos] 22 22 md.geometry.surface = md.geometry.base + md.geometry.thickness 23 md = md.extrude(4, 1.)24 md = setflowequation(md, 'HO','all')23 md = md.extrude(4, 1.) 24 md = setflowequation(md, 'HO', 'all') 25 25 26 26 #Boundary conditions: … … 30 30 md.stressbalance.spcvy[:] = float('Nan') 31 31 md.stressbalance.spcvz[:] = float('Nan') 32 posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)), np.array(np.where(md.mesh.y > 999999.9)))33 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > -0.1)))34 pos = np.unique(np.concatenate((posA, posB)))32 posA = np.intersect1d(np.array(np.where(md.mesh.y < 1000000.1)), np.array(np.where(md.mesh.y > 999999.9))) 33 posB = np.intersect1d(np.array(np.where(md.mesh.y < 0.1)), np.array(np.where(md.mesh.y > -0.1))) 34 pos = np.unique(np.concatenate((posA, posB))) 35 35 md.stressbalance.spcvy[pos] = 0. 36 36 pos2 = np.intersect1d(np.array(np.where(md.mesh.x < 0.1)), np.array(np.where(md.mesh.x > -0.1))) … … 38 38 md.stressbalance.spcvy[pos2] = 0. 39 39 40 md.materials.rheology_B = 1. / ((10**-25)**(1. /3.)) * np.ones((md.mesh.numberofvertices,))40 md.materials.rheology_B = 1. / ((10**-25)**(1. / 3.)) * np.ones((md.mesh.numberofvertices,)) 41 41 md.materials.rheology_law = 'None' 42 42 md.friction.coefficient[:] = np.sqrt(1e7) * np.ones((md.mesh.numberofvertices,)) … … 50 50 md.transient.ismasstransport = 1 51 51 md.transient.issmb = 1 52 md.transient.requested_outputs = ['default', 'BasalforcingsFloatingiceMeltingRate']52 md.transient.requested_outputs = ['default', 'BasalforcingsFloatingiceMeltingRate'] 53 53 md.groundingline.migration = 'SubelementMigration' 54 md.groundingline.friction_interpolation = 'SubelementFriction2'55 md.groundingline.melt_interpolation = 'SubelementMelt2'54 md.groundingline.friction_interpolation = 'SubelementFriction2' 55 md.groundingline.melt_interpolation = 'SubelementMelt2' 56 56 md.timestepping.final_time = 30 57 57 md.timestepping.time_step = 10 58 58 59 md.cluster = generic('name', gethostname(),'np',3)60 md = solve(md, 'Transient')59 md.cluster = generic('name', gethostname(), 'np', 3) 60 md = solve(md, 'Transient') 61 61 62 62 #Fields and tolerances to track changes 63 field_names = [ 64 'Bed1','Surface1','Thickness1','Floatingice1','Vx1','Vy1','Vz1','Pressure1','FloatingiceMeltingrate1', 65 'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Vz2','Pressure2','FloatingiceMeltingrate2', 66 'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Vz3','Pressure3','FloatingiceMeltingrate3'] 67 field_tolerances = [ 68 2e-11,5e-12,2e-11,1e-11,5e-10,3e-08,6e-10,1e-13,1e-13, 69 3e-11,3e-11,9e-10,7e-11,7e-09,1e-07,1e-09,1e-10,1e-13, 70 1e-8,2e-08,7e-09,2e-7 ,1e-03,8e-04,2e-09,1e-10,1e-13] 71 field_values = [ 72 md.results.TransientSolution[0].Base, 73 md.results.TransientSolution[0].Surface, 74 md.results.TransientSolution[0].Thickness, 75 md.results.TransientSolution[0].MaskGroundediceLevelset, 76 md.results.TransientSolution[0].Vx, 77 md.results.TransientSolution[0].Vy, 78 md.results.TransientSolution[0].Vz, 79 md.results.TransientSolution[0].Pressure, 80 md.results.TransientSolution[0].BasalforcingsFloatingiceMeltingRate, 81 md.results.TransientSolution[1].Base, 82 md.results.TransientSolution[1].Surface, 83 md.results.TransientSolution[1].Thickness, 84 md.results.TransientSolution[1].MaskGroundediceLevelset, 85 md.results.TransientSolution[1].Vx, 86 md.results.TransientSolution[1].Vy, 87 md.results.TransientSolution[1].Vz, 88 md.results.TransientSolution[1].Pressure, 89 md.results.TransientSolution[1].BasalforcingsFloatingiceMeltingRate, 90 md.results.TransientSolution[2].Base, 91 md.results.TransientSolution[2].Surface, 92 md.results.TransientSolution[2].Thickness, 93 md.results.TransientSolution[2].MaskGroundediceLevelset, 94 md.results.TransientSolution[2].Vx, 95 md.results.TransientSolution[2].Vy, 96 md.results.TransientSolution[2].Vz, 97 md.results.TransientSolution[2].Pressure, 98 md.results.TransientSolution[2].BasalforcingsFloatingiceMeltingRate, 99 ] 63 field_names = ['Bed1', 'Surface1', 'Thickness1', 'Floatingice1', 'Vx1', 'Vy1', 'Vz1', 'Pressure1', 'FloatingiceMeltingrate1', 64 'Bed2', 'Surface2', 'Thickness2', 'Floatingice2', 'Vx2', 'Vy2', 'Vz2', 'Pressure2', 'FloatingiceMeltingrate2', 65 'Bed3', 'Surface3', 'Thickness3', 'Floatingice3', 'Vx3', 'Vy3', 'Vz3', 'Pressure3', 'FloatingiceMeltingrate3'] 66 field_tolerances = [2e-11, 5e-12, 2e-11, 1e-11, 5e-10, 3e-08, 6e-10, 1e-13, 1e-13, 67 3e-11, 3e-11, 9e-10, 7e-11, 7e-09, 1e-07, 1e-09, 1e-10, 1e-13, 68 1e-8, 2e-08, 7e-09, 2e-7, 1e-03, 8e-04, 2e-09, 1e-10, 1e-13] 69 field_values = [md.results.TransientSolution[0].Base, 70 md.results.TransientSolution[0].Surface, 71 md.results.TransientSolution[0].Thickness, 72 md.results.TransientSolution[0].MaskGroundediceLevelset, 73 md.results.TransientSolution[0].Vx, 74 md.results.TransientSolution[0].Vy, 75 md.results.TransientSolution[0].Vz, 76 md.results.TransientSolution[0].Pressure, 77 md.results.TransientSolution[0].BasalforcingsFloatingiceMeltingRate, 78 md.results.TransientSolution[1].Base, 79 md.results.TransientSolution[1].Surface, 80 md.results.TransientSolution[1].Thickness, 81 md.results.TransientSolution[1].MaskGroundediceLevelset, 82 md.results.TransientSolution[1].Vx, 83 md.results.TransientSolution[1].Vy, 84 md.results.TransientSolution[1].Vz, 85 md.results.TransientSolution[1].Pressure, 86 md.results.TransientSolution[1].BasalforcingsFloatingiceMeltingRate, 87 md.results.TransientSolution[2].Base, 88 md.results.TransientSolution[2].Surface, 89 md.results.TransientSolution[2].Thickness, 90 md.results.TransientSolution[2].MaskGroundediceLevelset, 91 md.results.TransientSolution[2].Vx, 92 md.results.TransientSolution[2].Vy, 93 md.results.TransientSolution[2].Vz, 94 md.results.TransientSolution[2].Pressure, 95 md.results.TransientSolution[2].BasalforcingsFloatingiceMeltingRate]
Note:
See TracChangeset
for help on using the changeset viewer.