Changeset 15944

Timestamp:

08/26/13 10:36:00 (12 years ago)

Author:

Mathieu Morlighem

Message:

CHG: prepare groundedice_levelset in setmask

Location:

issm/trunk-jpl/src/m/parameterization

Files:

• setmask.m (modified) (1 diff)
• setmask.py (modified) (1 diff)

issm/trunk-jpl/src/m/parameterization/setmask.m

@@ -46 +46 @@
 md.mask.elementongroundedice=elementongroundedice;
 md.mask.vertexongroundedice=vertexongroundedice;
+
+%level sets
 md.mask.ice_levelset=ones(md.mesh.numberofvertices,1);
+md.mask.groundedice_levelset=vertexongroundedice;
+md.mask.groundedice_levelset(find(vertexongroundedice==0.))=-1.;

issm/trunk-jpl/src/m/parameterization/setmask.py

@@ -46 +46 @@
         md.mask.elementongroundedice = elementongroundedice
         md.mask.vertexongroundedice = vertexongroundedice
-        md.mask.ice_levelset = numpy.ones(md.mesh.numberofvertices,bool)
+
+        #level sets
+        md.mask.ice_levelset         = numpy.ones(md.mesh.numberofvertices,bool)
+        md.mask.groundedice_levelset = -1.*numpy.ones((md.mesh.numberofvertices,1))
+        md.mask.groundedice_levelset[md.mesh.elements[numpy.nonzero(elementongroundedice),:]-1]=1.
 
         return md