Changeset 15411

Timestamp:

07/03/13 13:42:34 (12 years ago)

Author:

Mathieu Morlighem

Message:

FIX: mask hopefully working in python now

File:

• issm/trunk-jpl/src/m/classes/mask.py (modified) (1 diff)

issm/trunk-jpl/src/m/classes/mask.py

@@ -56 +56 @@
                 WriteData(fid,'object',self,'fieldname','vertexongroundedice','format','DoubleMat','mattype',1)
                 WriteData(fid,'object',self,'fieldname','vertexonwater','format','DoubleMat','mattype',1)
-                icelevelset=numpy.ones(md.mesh.numberofvertices,int)
-                pos=md.diagnostic.icefront[:,:-2]
-                icelevelset[pos[:]-1]=0;
+                icelevelset=numpy.ones(md.mesh.numberofvertices,float)
+                if numpy.any(numpy.logical_not(numpy.isnan(md.diagnostic.icefront))):
+                        pos=md.diagnostic.icefront[:,:-2]
+                        icelevelset[pos[:]-1]=0;
                 WriteData(fid,'data',icelevelset,'format','DoubleMat','mattype',1,'enum',IcelevelsetEnum())
         # }}}