Changeset 24758

Timestamp:

04/28/20 16:40:09 (5 years ago)

Author:

jdquinn

Message:

CHG: Arbitrary geometry

File:

• issm/trunk-jpl/test/NightlyRun/test2002.py (modified) (2 diffs)

issm/trunk-jpl/test/NightlyRun/test2002.py

@@ -9 +9 @@
 from paterson import *
 from love_numbers import *
+from materials import *
 
 #mesh earth:
@@ -66 +67 @@
 md.slr.ocean_area_scaling = 0
 
-#geometry
-di = md.materials.rho_ice / md.materials.rho_water
+#geometry for the bed, arbitrary
 md.geometry.thickness = np.ones((md.mesh.numberofvertices))
-md.geometry.surface = (1 - di) * np.zeros((md.mesh.numberofvertices))
+md.geometry.surface = np.zeros((md.mesh.numberofvertices))
 md.geometry.base = md.geometry.surface - md.geometry.thickness
 md.geometry.bed = md.geometry.base
 
 #materials
-md.initialization.temperature = 273.25 * np.ones((md.mesh.numberofvertices))
-md.materials.rheology_B = paterson(md.initialization.temperature)
-md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements))
+md.materials=materials('hydro')
 
 #New stuff