Index: /issm/trunk-jpl/src/m/parameterization/setmask.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setmask.py	(revision 17867)
+++ /issm/trunk-jpl/src/m/parameterization/setmask.py	(revision 17868)
@@ -55,9 +55,9 @@
 
 	#level sets
-	md.mask.groundedice_levelset = -1.*numpy.ones((md.mesh.numberofvertices,1))
+	md.mask.groundedice_levelset = -1.*numpy.ones(md.mesh.numberofvertices)
 	md.mask.groundedice_levelset[md.mesh.elements[numpy.nonzero(elementongroundedice),:]-1]=1.
 
 	if(len(args)):
-		md.mask.ice_levelset = 1.*numpy.ones((md.mesh.numberofvertices,1))
+		md.mask.ice_levelset = 1.*numpy.ones(md.mesh.numberofvertices)
 		icedomainfile = options.getfieldvalue('icedomain','none')
 		if not os.path.exists(icedomainfile):
@@ -67,5 +67,5 @@
 		md.mask.ice_levelset[numpy.nonzero(vertexinsideicedomain)[0]] = -1.
 	else:
-		md.mask.ice_levelset = -1.*numpy.ones((md.mesh.numberofvertices,1))
+		md.mask.ice_levelset = -1.*numpy.ones(md.mesh.numberofvertices)
 
 	return md