source: issm/oecreview/Archive/24684-25833/ISSM-24757-24758.diff@ 25834

Last change on this file since 25834 was 25834, checked in by Mathieu Morlighem, 4 years ago

CHG: added 24684-25833
File size: 1.1 KB
RevLine 
[25834]1Index: ../trunk-jpl/test/NightlyRun/test2002.py
2===================================================================
3--- ../trunk-jpl/test/NightlyRun/test2002.py (revision 24757)
4+++ ../trunk-jpl/test/NightlyRun/test2002.py (revision 24758)
5@@ -8,6 +8,7 @@
6 from gmtmask import *
7 from paterson import *
8 from love_numbers import *
9+from materials import *
10
11 #mesh earth:
12 md = model()
13@@ -65,17 +66,14 @@
14
15 md.slr.ocean_area_scaling = 0
16
17-#geometry
18-di = md.materials.rho_ice / md.materials.rho_water
19+#geometry for the bed, arbitrary
20 md.geometry.thickness = np.ones((md.mesh.numberofvertices))
21-md.geometry.surface = (1 - di) * np.zeros((md.mesh.numberofvertices))
22+md.geometry.surface = np.zeros((md.mesh.numberofvertices))
23 md.geometry.base = md.geometry.surface - md.geometry.thickness
24 md.geometry.bed = md.geometry.base
25
26 #materials
27-md.initialization.temperature = 273.25 * np.ones((md.mesh.numberofvertices))
28-md.materials.rheology_B = paterson(md.initialization.temperature)
29-md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements))
30+md.materials=materials('hydro')
31
32 #New stuff
33 md.slr.spcthickness = np.nan * np.ones((md.mesh.numberofvertices, ))
Note: See TracBrowser for help on using the repository browser.