source: issm/oecreview/Archive/15392-16133/ISSM-15599-15600.diff@ 16134

Index: ../trunk-jpl/src/m/classes/mask.py
===================================================================
--- ../trunk-jpl/src/m/classes/mask.py  (revision 15599)
+++ ../trunk-jpl/src/m/classes/mask.py  (revision 15600)
@@ -48,12 +48,12 @@
                md = checkfield(md,'mask.vertexongroundedice' ,'size',[md.mesh.numberofvertices],'values',[0,1])
                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
        # }}}