Index: /issm/trunk-jpl/src/m/parameterization/setmask.m
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setmask.m	(revision 15943)
+++ /issm/trunk-jpl/src/m/parameterization/setmask.m	(revision 15944)
@@ -46,3 +46,7 @@
 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.;
Index: /issm/trunk-jpl/src/m/parameterization/setmask.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setmask.py	(revision 15943)
+++ /issm/trunk-jpl/src/m/parameterization/setmask.py	(revision 15944)
@@ -46,5 +46,9 @@
 	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