Changeset 24214 for issm/trunk-jpl/test/NightlyRun/test2101.py
- Timestamp:
- 10/11/19 00:27:00 (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2101.py
r23793 r24214 17 17 18 18 #define load 19 md.esa.deltathickness = np.zeros((md.mesh.numberofelements, ))19 md.esa.deltathickness = np.zeros((md.mesh.numberofelements, )) 20 20 pos = 449 21 md.esa.deltathickness[pos] = - 100# this is the only "icy" element21 md.esa.deltathickness[pos] = - 100 # this is the only "icy" element 22 22 23 23 #love numbers: 24 24 nlov = 10000 25 md.esa.love_h = np.array(love_numbers('h')) 25 md.esa.love_h = np.array(love_numbers('h')) #Originally had CM arg 26 26 md.esa.love_h = np.resize(md.esa.love_h, nlov + 1) 27 md.esa.love_l = np.array(love_numbers('l')) 27 md.esa.love_l = np.array(love_numbers('l')) #Originally had CM arg 28 28 md.esa.love_l = np.resize(md.esa.love_l, nlov + 1) 29 29 … … 33 33 34 34 #make sure wherever there is an ice load, that the mask is set to ice: 35 md.mask.ice_levelset = np.ones((md.mesh.numberofvertices, ))35 md.mask.ice_levelset = np.ones((md.mesh.numberofvertices, )) 36 36 pos = np.where(md.esa.deltathickness) 37 md.mask.ice_levelset[md.mesh.elements[pos, :]] = - 137 md.mask.ice_levelset[md.mesh.elements[pos, :]] = - 1 38 38 39 39 #is ice grounded? 40 md.mask.groundedice_levelset = - np.ones((md.mesh.numberofvertices,))40 md.mask.groundedice_levelset = - np.ones((md.mesh.numberofvertices, )) 41 41 pos = np.where(md.mask.ice_levelset <= 0) 42 42 md.mask.groundedice_levelset[pos] = 1 … … 52 52 #geometry: {{{ 53 53 di = md.materials.rho_ice / md.materials.rho_water 54 md.geometry.thickness = np.ones((md.mesh.numberofvertices, ))55 md.geometry.surface = (1 - di) * np.zeros((md.mesh.numberofvertices, ))54 md.geometry.thickness = np.ones((md.mesh.numberofvertices, )) 55 md.geometry.surface = (1 - di) * np.zeros((md.mesh.numberofvertices, )) 56 56 md.geometry.base = md.geometry.surface - md.geometry.thickness 57 57 md.geometry.bed = md.geometry.base 58 58 # }}} 59 59 #materials: {{{ 60 md.initialization.temperature = 273.25 * np.ones((md.mesh.numberofvertices, ))60 md.initialization.temperature = 273.25 * np.ones((md.mesh.numberofvertices, )) 61 61 md.materials.rheology_B = paterson(md.initialization.temperature) 62 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements, ))62 md.materials.rheology_n = 3 * np.ones((md.mesh.numberofelements, )) 63 63 # }}} 64 64 #Miscellaneous: {{{
Note:
See TracChangeset
for help on using the changeset viewer.