Index: /issm/trunk-jpl/test/NightlyRun/IdToName.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 18090)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.m	(revision 18091)
@@ -78,4 +78,5 @@
 	case 285, name='SquareShelfStressHOHigherOrder';
 	case 290, name='SquareShelfStressFSP2P1';
+	case 291, name='SquareShelfStressFSP4z';
 	case 301, name='SquareSheetConstrainedStressSSA2d';
 	case 302, name='SquareSheetConstrainedStressSIA2d';
Index: /issm/trunk-jpl/test/NightlyRun/IdToName.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 18090)
+++ /issm/trunk-jpl/test/NightlyRun/IdToName.py	(revision 18091)
@@ -77,4 +77,5 @@
 	285  : 'SquareShelfStressHOHigherOrder',
 	290  : 'SquareShelfStressFSP2P1',
+	291  : 'SquareShelfStressFSP4z',
 	301  : 'SquareSheetConstrainedStressSSA2d',
 	302  : 'SquareSheetConstrainedStressSIA2d',
Index: sm/trunk-jpl/test/NightlyRun/temp291.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/temp291.m	(revision 18090)
+++ 	(revision )
@@ -1,8 +1,0 @@
-md=triangle(model(),'../Exp/Square.exp',150000.);
-md=setmask(md,'all','');
-md=parameterize(md,'../Par/SquareShelf.par');
-md=extrude(md,2,1.);
-md=setflowequation(md,'FS','all');
-md.flowequation.fe_FS='OneLayerP4z';
-md.cluster=generic('name',oshostname(),'np',1);
-md=solve(md,StressbalanceSolutionEnum());
Index: sm/trunk-jpl/test/NightlyRun/temp291.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/temp291.py	(revision 18090)
+++ 	(revision )
@@ -1,18 +1,0 @@
-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(2,1.)
-md=setflowequation(md,'FS','all')
-md.flowequation.fe_FS='OneLayerP4z'
-md.cluster=generic('name',oshostname(),'np',1)
-md=solve(md,StressbalanceSolutionEnum())
Index: /issm/trunk-jpl/test/NightlyRun/test291.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test291.m	(revision 18091)
+++ /issm/trunk-jpl/test/NightlyRun/test291.m	(revision 18091)
@@ -0,0 +1,8 @@
+md=triangle(model(),'../Exp/Square.exp',150000.);
+md=setmask(md,'all','');
+md=parameterize(md,'../Par/SquareShelf.par');
+md=extrude(md,2,1.);
+md=setflowequation(md,'FS','all');
+md.flowequation.fe_FS='OneLayerP4z';
+md.cluster=generic('name',oshostname(),'np',1);
+md=solve(md,StressbalanceSolutionEnum());
Index: /issm/trunk-jpl/test/NightlyRun/test291.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test291.py	(revision 18091)
+++ /issm/trunk-jpl/test/NightlyRun/test291.py	(revision 18091)
@@ -0,0 +1,18 @@
+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(2,1.)
+md=setflowequation(md,'FS','all')
+md.flowequation.fe_FS='OneLayerP4z'
+md.cluster=generic('name',oshostname(),'np',1)
+md=solve(md,StressbalanceSolutionEnum())