Changeset 19607

Timestamp:

09/30/15 13:00:11 (10 years ago)

Author:

seroussi

Message:

NEW: improved error message to figure out how many elements have all their nodes on ice front

File:

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

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

@@ -56 +56 @@
                         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')
+                                if strcmp(elementtype(md.mesh),'Tria'),
+                                        numberemptyelements=length(find(icefront==3));
+                                elseif strcmp(elementtype(md.mesh),'Penta')
+                                        numberemptyelements=length(find(icefront==6));
+                                end
+                                error([' ' num2str(numberemptyelements) ' have all nodes on ice front, change md.mask.ice_levelset to fix it'])
                         end
                 end % }}}