Index: /issm/trunk-jpl/test/NightlyRun/test111.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test111.js	(revision 20594)
+++ /issm/trunk-jpl/test/NightlyRun/test111.js	(revision 20594)
@@ -0,0 +1,55 @@
+//Test Name: SquareShelfConstrainedTranHO
+var md = new model();
+triangle(md,square[0],200000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,1.);
+setflowequation(md,'HO','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md.transient.requested_outputs=['IceVolume'];
+md=solve(md,TransientSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names=[
+	'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface2','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1','Volume1', 
+	'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2','Volume2', 
+	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3','Volume3'];
+field_tolerances=[
+	1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-13,
+	1e-09,1e-08,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,2e-11,
+	1e-09,1e-08,1e-08,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-05,1e-11];
+field_values=[
+	(md.results.TransientSolution[0](1).Vx),
+	(md.results.TransientSolution[0](1).Vy),
+	(md.results.TransientSolution[0](1).Vz),
+	(md.results.TransientSolution[0](1).Vel),
+	(md.results.TransientSolution[0](1).Pressure),
+	(md.results.TransientSolution[0](1).Base),
+	(md.results.TransientSolution[0](1).Surface),
+	(md.results.TransientSolution[0](1).Thickness),
+	(md.results.TransientSolution[0](1).Temperature),
+	(md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate),
+	(md.results.TransientSolution[0](1).IceVolume),
+	(md.results.TransientSolution[0](2).Vx),
+	(md.results.TransientSolution[0](2).Vy),
+	(md.results.TransientSolution[0](2).Vz),
+	(md.results.TransientSolution[0](2).Vel),
+	(md.results.TransientSolution[0](2).Pressure),
+	(md.results.TransientSolution[0](2).Base),
+	(md.results.TransientSolution[0](2).Surface),
+	(md.results.TransientSolution[0](2).Thickness),
+	(md.results.TransientSolution[0](2).Temperature),
+	(md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate),
+	(md.results.TransientSolution[0](2).IceVolume),
+	(md.results.TransientSolution[0](3).Vx),
+	(md.results.TransientSolution[0](3).Vy),
+	(md.results.TransientSolution[0](3).Vz),
+	(md.results.TransientSolution[0](3).Vel),
+	(md.results.TransientSolution[0](3).Pressure),
+	(md.results.TransientSolution[0](3).Base),
+	(md.results.TransientSolution[0](3).Surface),
+	(md.results.TransientSolution[0](3).Thickness),
+	(md.results.TransientSolution[0](3).Temperature),
+	(md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate),
+	(md.results.TransientSolution[0](3).IceVolume),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test112.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test112.js	(revision 20594)
+++ /issm/trunk-jpl/test/NightlyRun/test112.js	(revision 20594)
@@ -0,0 +1,16 @@
+//Test Name: SquareShelfConstrainedSurfSlop2d
+var md = new model();
+triangle(md,square[0],150000.);
+setmask(md,'all','');
+parameterize(md);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,SurfaceSlopeSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['SurfaceSlopeX','SurfaceSlopeY'];
+field_tolerances=[1e-13,1e-13];
+field_values=[
+	(md.results.SurfaceSlopeSolution[0].SurfaceSlopeX),
+	(md.results.SurfaceSlopeSolution[0].SurfaceSlopeY),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test113.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test113.js	(revision 20594)
+++ /issm/trunk-jpl/test/NightlyRun/test113.js	(revision 20594)
@@ -0,0 +1,17 @@
+//Test Name: SquareShelfConstrainedSurfSlope3d
+var md = new model();
+triangle(md,square[0],150000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,5,1.);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,SurfaceSlopeSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['SurfaceSlopeX','SurfaceSlopeY'];
+field_tolerances=[1e-13,1e-13];
+field_values=[
+	(md.results.SurfaceSlopeSolution[0].SurfaceSlopeX),
+	(md.results.SurfaceSlopeSolution[0].SurfaceSlopeY),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test114.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test114.js	(revision 20594)
+++ /issm/trunk-jpl/test/NightlyRun/test114.js	(revision 20594)
@@ -0,0 +1,16 @@
+//Test Name: SquareShelfConstrainedBedSlop2d
+var md = new model();
+triangle(md,square[0],150000.);
+setmask(md,'all','');
+parameterize(md);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,BedSlopeSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['BedSlopeX','BedSlopeY'];
+field_tolerances=[1e-13,1e-13];
+field_values=[
+	(md.results.BedSlopeSolution[0].BedSlopeX),
+	(md.results.BedSlopeSolution[0].BedSlopeY),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test115.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test115.js	(revision 20594)
+++ /issm/trunk-jpl/test/NightlyRun/test115.js	(revision 20594)
@@ -0,0 +1,17 @@
+//Test Name: SquareShelfConstrainedBedSlop3d
+var md = new model();
+triangle(md,square[0],150000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,5,1.);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,BedSlopeSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['BedSlopeX','BedSlopeY'];
+field_tolerances=[1e-13,1e-13];
+field_values=[
+	(md.results.BedSlopeSolution[0].BedSlopeX),
+	(md.results.BedSlopeSolution[0].BedSlopeY),
+	];