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 # }}}