Index: /issm/trunk/src/m/classes/public/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 5618)
+++ /issm/trunk/src/m/classes/public/ismodelselfconsistent.m	(revision 5619)
@@ -320,4 +320,9 @@
 			error(['model not consistent: model ' md.name ' should have all the nodes on boundary constrained in field spcthickness']);			 
 		end 
+	end
+
+	%Triangle with zero velocity
+	if any(sum(abs(md.vx(md.elements)),2)==0 & sum(abs(md.vy(md.elements)),2)==0)
+		error('model not consistent: at least one triangle has all its vertices with a zero velocity');
 	end
 end