Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15925)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 15926)
@@ -68,5 +68,5 @@
 	case 272, name='SquareShelfCMZSSA2dDamage';
 	case 274, name='SquareShelfStressSSA2dDamageRift';
-	case 280, name='SquareShelfStressSSA2dP2';
+	case 280, name='SquareShelfStressSSA2dHigherOrder';
 	case 285, name='SquareShelfStressHOP1xP2';
 	case 286, name='SquareShelfStressHOP2xP1';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15925)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 15926)
@@ -72,5 +72,5 @@
 	272  : 'SquareShelfCMZSSA2dDamage',
 	274  : 'SquareShelfStressSSA2dDamageRift',
-	280  : 'SquareShelfStressSSA2dP2',
+	280  : 'SquareShelfStressSSA2dHigherOrder',
 	285  : 'SquareShelfStressHOP1xP2',
 	286  : 'SquareShelfStressHOP2xP1',
Index: /issm/trunk-jpl/test/NightlyRun/test280.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test280.m	(revision 15925)
+++ /issm/trunk-jpl/test/NightlyRun/test280.m	(revision 15926)
@@ -3,15 +3,20 @@
 md=parameterize(md,'../Par/SquareShelf.par');
 md=setflowequation(md,'SSA','all');
-md.flowequation.fe_SSA='P2';
 md.cluster=generic('name',oshostname(),'np',3);
-md=solve(md,StressbalanceSolutionEnum());
 
-%Fields and tolerances to track changes
-field_names     ={'Vx','Vy','Vel','Pressure'};
-field_tolerances={1e-12,1e-13,1e-13,1e-13};
-field_values={...
-	(md.results.StressbalanceSolution.Vx),...
-	(md.results.StressbalanceSolution.Vy),...
-	(md.results.StressbalanceSolution.Vel),...
-	(md.results.StressbalanceSolution.Pressure),...
-	};
+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/test280.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test280.py	(revision 15925)
+++ /issm/trunk-jpl/test/NightlyRun/test280.py	(revision 15926)
@@ -13,15 +13,18 @@
 md=parameterize(md,'../Par/SquareShelf.py')
 md=setflowequation(md,'SSA','all')
-md.flowequation.fe_SSA='P2'
 md.cluster=generic('name',oshostname(),'np',3)
-md=solve(md,StressbalanceSolutionEnum())
 
-#Fields and tolerances to track changes
-field_names     =['Vx','Vy','Vel','Pressure']
-field_tolerances=[1e-12,1e-13,1e-13,1e-13]
-field_values=[\
-	md.results.StressbalanceSolution.Vx,\
-	md.results.StressbalanceSolution.Vy,\
-	md.results.StressbalanceSolution.Vel,\
-	md.results.StressbalanceSolution.Pressure,\
-	]
+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,\
+		]