Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp	(revision 15643)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateConstraintsDiagnosticHoriz.cpp	(revision 15644)
@@ -86,7 +86,7 @@
 			switch(temp){
 				case 0 : finiteelement = P1Enum;    break;
-				case 1 : finiteelement = P2Enum;    break;
+				case 1 : finiteelement = P1xP2Enum; break;
 				case 2 : finiteelement = P2xP1Enum; break;
-				case 3 : finiteelement = P1xP2Enum; break;
+				case 3 : finiteelement = P2Enum;    break;
 				default: _error_("finite element "<<temp<<" not supported");
 			}
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateNodesDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateNodesDiagnosticHoriz.cpp	(revision 15643)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/CreateNodesDiagnosticHoriz.cpp	(revision 15644)
@@ -52,7 +52,7 @@
 			switch(temp){
 				case 0 : finiteelement = P1Enum;    break;
-				case 1 : finiteelement = P2Enum;    break;
+				case 1 : finiteelement = P1xP2Enum; break;
 				case 2 : finiteelement = P2xP1Enum; break;
-				case 3 : finiteelement = P1xP2Enum; break;
+				case 3 : finiteelement = P2Enum;    break;
 				default: _error_("finite element "<<temp<<" not supported");
 			}
Index: /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/UpdateElementsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/UpdateElementsDiagnosticHoriz.cpp	(revision 15643)
+++ /issm/trunk-jpl/src/c/modules/ModelProcessorx/DiagnosticHoriz/UpdateElementsDiagnosticHoriz.cpp	(revision 15644)
@@ -57,7 +57,7 @@
 			switch(temp){
 				case 0 : finiteelement = P1Enum;    break;
-				case 1 : finiteelement = P2Enum;    break;
+				case 1 : finiteelement = P1xP2Enum; break;
 				case 2 : finiteelement = P2xP1Enum; break;
-				case 3 : finiteelement = P1xP2Enum; break;
+				case 3 : finiteelement = P2Enum;    break;
 				default: _error_("finite element "<<temp<<" not supported");
 			}
Index: /issm/trunk-jpl/src/m/classes/flowequation.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/flowequation.m	(revision 15643)
+++ /issm/trunk-jpl/src/m/classes/flowequation.m	(revision 15644)
@@ -113,5 +113,6 @@
 			fielddisplay(obj,'isHO','is the Higher-Order (HO) approximation used ?');
 			fielddisplay(obj,'isFS','are the Full-FS (FS) equations used ?');
-			fielddisplay(obj,'fe_SSA','Finite Element for SSA   0: Lagrange P1 (linear), 1: Lagrange P2 (quadratic)');
+			fielddisplay(obj,'fe_SSA','Finite Element for SSA  0: Lagrange P1 (linear), 1: Lagrange P2 (quadratic)');
+			fielddisplay(obj,'fe_HO', 'Finite Element for HO   0: P1xP1, 1: P1xP2, 2: P2xP1, 3: P2xP2');
 			fielddisplay(obj,'vertex_equation','flow equation for each vertex');
 			fielddisplay(obj,'element_equation','flow equation for each element');
Index: /issm/trunk-jpl/src/m/classes/flowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15643)
+++ /issm/trunk-jpl/src/m/classes/flowequation.py	(revision 15644)
@@ -43,4 +43,5 @@
 		string="%s\n%s"%(string,fielddisplay(self,'isFS',"are the Full-FS (FS) equations used ?"))
 		string="%s\n%s"%(string,fielddisplay(self,'fe_SSA',"Finite Element for SSA   0: Lagrange P1 (linear), 1: Lagrange P2 (quadratic)"))
+		string="%s\n%s"%(string,fielddisplay(self,'fe_HO' ,"Finite Element for HO   0: P1xP1, 1: P1xP2, 2: P2xP1, 3: P2xP2"))
 		string="%s\n%s"%(string,fielddisplay(self,'vertex_equation',"flow equation for each vertex"))
 		string="%s\n%s"%(string,fielddisplay(self,'element_equation',"flow equation for each element"))
