Index: /issm/trunk-jpl/src/m/classes/mask.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/mask.py	(revision 15599)
+++ /issm/trunk-jpl/src/m/classes/mask.py	(revision 15600)
@@ -49,10 +49,10 @@
 		md = checkfield(md,'mask.vertexonwater'       ,'size',[md.mesh.numberofvertices],'values',[0,1])
 		md = checkfield(md,'mask.icelevelset'         ,'size',[md.mesh.numberofvertices])
-			isice=numpy.zeros_like(md.mesh.numberofvertices,int)
-			pos=npy.nonzero(md.mask.icelevelset>0)
-			isice[pos]=1
-			totalyicefree=(numpy.sum(isice[md.mesh.elements-1],axis=1)==0).astype(int)
-			if any(totalyicefree)
-				raise TypeError("elements with no ice not implemented yet, each element should have at least one vertex with md.mask.icelevelset > 0")
+		isice=numpy.zeros(md.mesh.numberofvertices,int)
+		pos=npy.nonzero(md.mask.icelevelset>0)
+		isice[pos]=1
+		totallyicefree=(numpy.sum(isice[md.mesh.elements-1],axis=1)==0).astype(int)
+		if any(totallyicefree):
+			raise TypeError("elements with no ice not implemented yet, each element should have at least one vertex with md.mask.icelevelset > 0")
 
 		return md