Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test460.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test460.py
r24025 r24214 15 15 md = md.extrude(3, 1.) 16 16 md.materials = matestar() 17 md.materials.rheology_B = 3.15e8 * np.ones((md.mesh.numberofvertices, ))18 md.materials.rheology_Ec = np.ones((md.mesh.numberofvertices, ))19 md.materials.rheology_Es = 3 * np.ones((md.mesh.numberofvertices, ))17 md.materials.rheology_B = 3.15e8 * np.ones((md.mesh.numberofvertices, )) 18 md.materials.rheology_Ec = np.ones((md.mesh.numberofvertices, )) 19 md.materials.rheology_Es = 3 * np.ones((md.mesh.numberofvertices, )) 20 20 md.cluster = generic('name', gethostname(), 'np', 3) 21 21 … … 24 24 field_tolerances = [] 25 25 field_values = [] 26 #md.initialization.pressure = md.constants.g *md.materials.rho_ice*(md.geometry.surface-md.mesh.y);26 #md.initialization.pressure = md.constants.g * md.materials.rho_ice * (md.geometry.surface - md.mesh.y) 27 27 for i in ['SSA', 'HO', 'FS']: 28 28 md = setflowequation(md, i, 'all') 29 29 md = solve(md, 'Stressbalance') 30 30 field_names = field_names + ['Vx' + i, 'Vy' + i, 'Vz' + i, 'Vel' + i, 'Pressure' + i] 31 field_tolerances = field_tolerances + [7e-06, 2e-05,2e-06,5e-06,8e-07]31 field_tolerances = field_tolerances + [7e-06, 2e-05, 2e-06, 5e-06, 8e-07] 32 32 field_values = field_values + [md.results.StressbalanceSolution.Vx, 33 33 md.results.StressbalanceSolution.Vy,
Note:
See TracChangeset
for help on using the changeset viewer.