- Timestamp:
- 08/26/13 10:41:10 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/src/m/boundaryconditions/SetIceShelfBC.py
r15816 r15946 40 40 #Icefront position 41 41 pos=numpy.nonzero(nodeonicefront)[0] 42 md.mask.ice levelset[pos]=042 md.mask.ice_levelset[pos]=0 43 43 44 44 #First find segments that are not completely on the front … … 47 47 else: 48 48 numbernodesfront=4 49 if any(md.mask.ice levelset<=0):50 values=md.mask.ice levelset[md.mesh.segments[:,0:-1]-1]49 if any(md.mask.ice_levelset<=0): 50 values=md.mask.ice_levelset[md.mesh.segments[:,0:-1]-1] 51 51 segmentsfront=1-values 52 52 numpy.sum(segmentsfront,axis=1)!=numbernodesfront
Note:
See TracChangeset
for help on using the changeset viewer.