Index: /issm/trunk-jpl/src/m/classes/flowequation.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/flowequation.m	(revision 15526)
+++ /issm/trunk-jpl/src/m/classes/flowequation.m	(revision 15527)
@@ -50,5 +50,5 @@
 					md = checkfield(md,'flowequation.element_equation','size',[md.mesh.numberofelements 1],'values',[0:8]);
 				end
-				if ~(obj.ismacayeal || obj.ispattyn || obj.ishutter || obj.isstokes || obj.isl1l2),
+				if ~(obj.ishutter || obj.ismacayeal || obj.isl1l2 || obj.ispattyn || obj.isstokes),
 					md = checkmessage(md,['no element types set for this model']);
 				end
Index: /issm/trunk-jpl/src/m/classes/flowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15526)
+++ /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15527)
@@ -67,7 +67,7 @@
 				md = checkfield(md,'flowequation.element_equation','size',[md.mesh.numberofelements],'values',[1,2])
 			else:
-				md = checkfield(md,'flowequation.vertex_equation','size',[md.mesh.numberofvertices],'values',numpy.arange(0,7+1))
-				md = checkfield(md,'flowequation.element_equation','size',[md.mesh.numberofelements],'values',numpy.arange(0,7+1))
-			if not (self.ismacayeal or self.ispattyn or self.ishutter or self.isstokes or self.isl1l2):
+				md = checkfield(md,'flowequation.vertex_equation','size',[md.mesh.numberofvertices],'values',numpy.arange(0,8))
+				md = checkfield(md,'flowequation.element_equation','size',[md.mesh.numberofelements],'values',numpy.arange(0,8))
+			if not (self.ishutter or self.ismacayeal or self.isl1l2 or self.ispattyn or self.isstokes):
 				md.checkmessage("no element types set for this model")