Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15927)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15928)
@@ -136,4 +136,6 @@
 	case 433, name='RoundSheetShelfGLMigrationSSA3d';
 	case 434, name='SquareSheetShelfL1L2';
+	case 450, name='SquareSheetShelfStressSSAHigherOrder';
+	case 455, name='SquareSheetShelfStressHOHigherOrder';
 	case 501, name='PigStressSSA2d';
 	case 502, name='PigStressHO';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15927)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15928)
@@ -140,4 +140,6 @@
 	433  : 'RoundSheetShelfGLMigrationSSA3d',
 	434  : 'SquareSheetShelfL1L2',
+	450  : 'SquareSheetShelfStressSSAHigherOrder',
+	455  : 'SquareSheetShelfStressHOHigherOrder',
 	501  : 'PigStressSSA2d',
 	502  : 'PigStressHO',
Index: /issm/trunk-jpl/test/NightlyRun/test280.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test280.m	(revision 15927)
+++ /issm/trunk-jpl/test/NightlyRun/test280.m	(revision 15928)
@@ -20,3 +20,2 @@
 		};
 end
-
Index: /issm/trunk-jpl/test/NightlyRun/test285.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test285.m	(revision 15927)
+++ /issm/trunk-jpl/test/NightlyRun/test285.m	(revision 15928)
@@ -4,5 +4,4 @@
 md=extrude(md,3,2.);
 md=setflowequation(md,'HO','all');
-md.flowequation.fe_HO='P1xP2';
 md.cluster=generic('name',oshostname(),'np',3);
 
Index: /issm/trunk-jpl/test/NightlyRun/test450.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test450.m	(revision 15928)
+++ /issm/trunk-jpl/test/NightlyRun/test450.m	(revision 15928)
@@ -0,0 +1,21 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setflowequation(md,'SSA','all');
+md.cluster=generic('name',oshostname(),'np',3);
+
+field_names={};
+field_tolerances={};
+field_values={};
+for i={'P1bubble','P1bubblecondensed','P2'}
+	md.flowequation.fe_SSA=i{1};
+	md=solve(md,StressbalanceSolutionEnum());
+	field_names     ={field_names{:},['Vx' i{1}],['Vy' i{1}],['Vel' i{1}],['Pressure' i{1}]};
+	field_tolerances={field_tolerances{:},1e-12,1e-13,1e-13,1e-13};
+	field_values={field_values{:},...
+		(md.results.StressbalanceSolution.Vx),...
+		(md.results.StressbalanceSolution.Vy),...
+		(md.results.StressbalanceSolution.Vel),...
+		(md.results.StressbalanceSolution.Pressure),...
+		};
+end
Index: /issm/trunk-jpl/test/NightlyRun/test450.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test450.py	(revision 15928)
+++ /issm/trunk-jpl/test/NightlyRun/test450.py	(revision 15928)
@@ -0,0 +1,30 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/Square.exp',150000.)
+md=setmask(md,'../Exp/SquareShelf.exp','')
+md=parameterize(md,'../Par/SquareSheetShelf.py')
+md=setflowequation(md,'SSA','all')
+md.cluster=generic('name',oshostname(),'np',3)
+
+field_names=[]
+field_tolerances=[]
+field_values=[]
+for i in ['P1bubble','P1bubblecondensed','P2']:
+	md.flowequation.fe_SSA=i
+	md=solve(md,StressbalanceSolutionEnum())
+	field_names     =field_names+['Vx'+i,'Vy'+i,'Vel'+i,'Pressure'+i]
+	field_tolerances=field_tolerances+[1e-12,1e-13,1e-13,1e-13]
+	field_values=field_values+[\
+		md.results.StressbalanceSolution.Vx,\
+		md.results.StressbalanceSolution.Vy,\
+		md.results.StressbalanceSolution.Vel,\
+		md.results.StressbalanceSolution.Pressure,\
+		]
Index: /issm/trunk-jpl/test/NightlyRun/test455.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test455.m	(revision 15928)
+++ /issm/trunk-jpl/test/NightlyRun/test455.m	(revision 15928)
@@ -0,0 +1,23 @@
+md=triangle(model(),'../Exp/Square.exp',180000.);
+md=setmask(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,5,1.);
+md=setflowequation(md,'HO','all');
+md.cluster=generic('name',oshostname(),'np',3);
+
+field_names={};
+field_tolerances={};
+field_values={};
+for i={'P1bubble','P1bubblecondensed','P1xP2','P2xP1','P2'}
+	md.flowequation.fe_HO=i{1};
+	md=solve(md,StressbalanceSolutionEnum());
+	field_names     ={field_names{:},['Vx' i{1}],['Vy' i{1}],['Vz' i{1}],['Vel' i{1}],['Pressure' i{1}]};
+	field_tolerances={field_tolerances{:},1e-08,1e-08,1e-08,1e-08,1e-13};
+	field_values={field_values{:},...
+	(md.results.StressbalanceSolution.Vx),...
+	(md.results.StressbalanceSolution.Vy),...
+	(md.results.StressbalanceSolution.Vz),...
+	(md.results.StressbalanceSolution.Vel),...
+	(md.results.StressbalanceSolution.Pressure),...
+	};
+end
Index: /issm/trunk-jpl/test/NightlyRun/test455.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test455.py	(revision 15928)
+++ /issm/trunk-jpl/test/NightlyRun/test455.py	(revision 15928)
@@ -0,0 +1,32 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/Square.exp',180000.)
+md=setmask(md,'../Exp/SquareShelf.exp','')
+md=parameterize(md,'../Par/SquareSheetShelf.py')
+md.extrude(5,1.)
+md=setflowequation(md,'HO','all')
+md.cluster=generic('name',oshostname(),'np',3)
+
+field_names=[]
+field_tolerances=[]
+field_values=[]
+for i in ['P1bubble','P1bubblecondensed','P1xP2','P2xP1','P2']:
+	md.flowequation.fe_HO=i
+	md=solve(md,StressbalanceSolutionEnum())
+	field_names     =field_names+['Vx'+i,'Vy'+i,'Vz'+i,'Vel'+i,'Pressure'+i]
+	field_tolerances=field_tolerances+[1e-08,1e-08,1e-08,1e-08,1e-13]
+	field_values=field_values+[\
+			md.results.StressbalanceSolution.Vx,\
+			md.results.StressbalanceSolution.Vy,\
+			md.results.StressbalanceSolution.Vz,\
+			md.results.StressbalanceSolution.Vel,\
+			md.results.StressbalanceSolution.Pressure,\
+			]