Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15685)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15686)
@@ -69,4 +69,8 @@
 	case 274, name='SquareShelfDiagSSA2dDamageRift';
 	case 280, name='SquareShelfDiagSSA2dP2';
+	case 285, name='SquareShelfDiagHOP1xP2';
+	case 286, name='SquareShelfDiagHOP2xP1';
+	case 287, name='SquareShelfDiagHOP2xP2';
+	case 290, name='SquareShelfDiagFSP2P1';
 	case 301, name='SquareSheetConstrainedDiagSSA2d';
 	case 302, name='SquareSheetConstrainedDiagSIA2d';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15685)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15686)
@@ -73,4 +73,8 @@
 	274  : 'SquareShelfDiagSSA2dDamageRift',
 	280  : 'SquareShelfDiagSSA2dP2',
+	285  : 'SquareShelfDiagHOP1xP2',
+	286  : 'SquareShelfDiagHOP2xP1',
+	287  : 'SquareShelfDiagHOP2xP2',
+	290  : 'SquareShelfDiagFSP2P1',
 	301  : 'SquareSheetConstrainedDiagSSA2d',
 	302  : 'SquareSheetConstrainedDiagSIA2d',
Index: /issm/trunk-jpl/test/NightlyRun/test285.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test285.m	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test285.m	(revision 15686)
@@ -0,0 +1,19 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,2.);
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=1;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-13};
+field_values={...
+	(md.results.DiagnosticSolution.Vx),...
+	(md.results.DiagnosticSolution.Vy),...
+	(md.results.DiagnosticSolution.Vz),...
+	(md.results.DiagnosticSolution.Vel),...
+	(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test285.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test285.py	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test285.py	(revision 15686)
@@ -0,0 +1,29 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+from ContourToMesh import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.py');
+md.extrude(3,2.)
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=1;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+#Fields and tolerances to track changes
+field_names     =['Vx', 'Vy', 'Vz', 'Vel','Pressure'];
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-13];
+field_values=[\
+	md.results.DiagnosticSolution.Vx,\
+	md.results.DiagnosticSolution.Vy,\
+	md.results.DiagnosticSolution.Vz,\
+	md.results.DiagnosticSolution.Vel,\
+	md.results.DiagnosticSolution.Pressure,\
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test286.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test286.m	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test286.m	(revision 15686)
@@ -0,0 +1,19 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,2.);
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=2;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-13};
+field_values={...
+	(md.results.DiagnosticSolution.Vx),...
+	(md.results.DiagnosticSolution.Vy),...
+	(md.results.DiagnosticSolution.Vz),...
+	(md.results.DiagnosticSolution.Vel),...
+	(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test286.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test286.py	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test286.py	(revision 15686)
@@ -0,0 +1,29 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+from ContourToMesh import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.py');
+md.extrude(3,2.)
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=2;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+#Fields and tolerances to track changes
+field_names     =['Vx', 'Vy', 'Vz', 'Vel','Pressure'];
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-13];
+field_values=[\
+	md.results.DiagnosticSolution.Vx,\
+	md.results.DiagnosticSolution.Vy,\
+	md.results.DiagnosticSolution.Vz,\
+	md.results.DiagnosticSolution.Vel,\
+	md.results.DiagnosticSolution.Pressure,\
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test287.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test287.m	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test287.m	(revision 15686)
@@ -0,0 +1,19 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,2.);
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=3;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-13};
+field_values={...
+	(md.results.DiagnosticSolution.Vx),...
+	(md.results.DiagnosticSolution.Vy),...
+	(md.results.DiagnosticSolution.Vz),...
+	(md.results.DiagnosticSolution.Vel),...
+	(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test287.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test287.py	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test287.py	(revision 15686)
@@ -0,0 +1,29 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+from ContourToMesh import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.py');
+md.extrude(3,2.)
+md=setflowequation(md,'HO','all');
+md.flowequation.fe_HO=3;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+#Fields and tolerances to track changes
+field_names     =['Vx', 'Vy', 'Vz', 'Vel','Pressure'];
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-13];
+field_values=[\
+	md.results.DiagnosticSolution.Vx,\
+	md.results.DiagnosticSolution.Vy,\
+	md.results.DiagnosticSolution.Vz,\
+	md.results.DiagnosticSolution.Vel,\
+	md.results.DiagnosticSolution.Pressure,\
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test290.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test290.m	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test290.m	(revision 15686)
@@ -0,0 +1,19 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,3,2.);
+md=setflowequation(md,'FS','all');
+md.flowequation.fe_FS=4;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+%Fields and tolerances to track changes
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-12,1e-12,1e-12,1e-12,1e-13};
+field_values={...
+	(md.results.DiagnosticSolution.Vx),...
+	(md.results.DiagnosticSolution.Vy),...
+	(md.results.DiagnosticSolution.Vz),...
+	(md.results.DiagnosticSolution.Vel),...
+	(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk-jpl/test/NightlyRun/test290.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test290.py	(revision 15686)
+++ /issm/trunk-jpl/test/NightlyRun/test290.py	(revision 15686)
@@ -0,0 +1,29 @@
+from model import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+from ContourToMesh import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.py');
+md.extrude(3,2.)
+md=setflowequation(md,'FS','all');
+md.flowequation.fe_FS=4;
+md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,DiagnosticSolutionEnum());
+
+#Fields and tolerances to track changes
+field_names     =['Vx', 'Vy', 'Vz', 'Vel','Pressure'];
+field_tolerances=[1e-12,1e-12,1e-12,1e-12,1e-13];
+field_values=[\
+	md.results.DiagnosticSolution.Vx,\
+	md.results.DiagnosticSolution.Vy,\
+	md.results.DiagnosticSolution.Vz,\
+	md.results.DiagnosticSolution.Vel,\
+	md.results.DiagnosticSolution.Pressure,\
+	];