source: issm/oecreview/Archive/15392-16133/ISSM-15526-15527.diff

../trunk-jpl/src/m/classes/flowequation.py

@@ -66 +66 @@
                                 md = checkfield(md,'flowequation.vertex_equation','size',[md.mesh.numberofvertices],'values',[1,2])
                                 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")
 
                 if DiagnosticHutterAnalysisEnum() in analyses:

../trunk-jpl/src/m/classes/flowequation.m

@@ -49 +49 @@
                                         md = checkfield(md,'flowequation.vertex_equation','size',[md.mesh.numberofvertices 1],'values',[0:8]);
                                         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
                         end