Changeset 18649

Timestamp:

10/15/14 14:40:07 (11 years ago)

Author:

seroussi

Message:

NEW: added check consistency of front (not all vertices on the front)

Location:

issm/trunk-jpl/src/m/classes

Files:

• mask.m (modified) (1 diff)
• mask.py (modified) (1 diff)

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

@@ -53 +53 @@
                                 warning('no ice front provided');
                         end
+                        icefront=sum(md.mask.ice_levelset(md.mesh.elements)==0,2);
+                        if (max(icefront)==3 & strcmp(elementtype(md.mesh),'Tria')) | (max(icefront==6) & strcmp(elementtype(md.mesh),'Penta')),
+                                error('At least one element has all nodes on ice front, change md.mask.ice_levelset to fix it')
+                        end
                 end % }}}
                 function disp(obj) % {{{

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

@@ -38 +38 @@
                         raise TypeError("no ice present in the domain")
 
+                icefront=numpy.sum(md.mask.ice_levelset[md.mesh.elements]==0,2)
+                if (max(icefront)==3 and m.strcmp(md.mesh.elementtype(),'Tria')) or (max(icefront==6) and m.strcmp(md.mesh.elementtype(),'Penta')):
+                        raise TypeError("At least one element has all nodes on ice front, change md.mask.ice_levelset to fix it")
+
                 return md
         # }}}