Changeset 17283

Timestamp:

02/14/14 14:15:24 (11 years ago)

Author:

jbondzio

Message:

BUG: in lsm, an element has ice only if at least one of its vertices has a negative levelset function value

File:

• issm/trunk-jpl/src/m/classes/mask.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/mask.py

@@ -46 +46 @@
 
                 # get mask of vertices of elements with ice
-                isice=numpy.array(md.mask.ice_levelset<=0.,int)
+                isice=numpy.array(md.mask.ice_levelset<0.,int)
                 vlist = numpy.zeros((md.mesh.numberofvertices,1), dtype=int)
                 pos=numpy.nonzero(numpy.sum(isice[md.mesh.elements-1],axis=1))[0]