Changeset 23039 for issm/trunk-jpl/test/NightlyRun/test2002.py
- Timestamp:
- 08/01/18 14:40:07 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/test/NightlyRun/test2002.py
r22981 r23039 55 55 md.mask.groundedice_levelset=-np.ones((md.mesh.numberofvertices)) 56 56 57 #make sure that the elements that have loads are fully grounded: 58 pos=np.nonzero(md.slr.deltathickness)[0] 59 md.mask.groundedice_levelset[md.mesh.elements[pos,:]-1]=1; 60 57 61 #make sure wherever there is an ice load, that the mask is set to ice: 58 pos=np.nonzero(md.slr.deltathickness)[0]59 62 icemask[md.mesh.elements[pos,:]-1]=-1 63 md.mask.ice_levelset=icemask 60 64 65 md.slr.ocean_area_scaling=0 61 66 62 67 #geometry … … 71 76 md.materials.rheology_B=paterson(md.initialization.temperature) 72 77 md.materials.rheology_n=3*np.ones((md.mesh.numberofelements)) 78 79 #New stuff 80 md.slr.spcthickness = np.nan((md.mesh.numberofvertices,)); 81 md.slr.Ngia = np.zeros((md.mesh.numberofvertices,)) 82 md.slr.Ugia = np.zeros((md.mesh.numberofvertices,)) 73 83 74 84 #Miscellaneous
Note:
See TracChangeset
for help on using the changeset viewer.