Index: /issm/trunk-jpl/test/NightlyRun/test201.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test201.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test201.js	(revision 20593)
@@ -0,0 +1,18 @@
+//Test Name: SquareShelfStressSSA2d
+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,StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vel','Pressure'];
+field_tolerances=[1e-13,1e-13,1e-13,1e-13];
+field_values=[
+	(md.results.StressbalanceSolution[0].Vx),
+	(md.results.StressbalanceSolution[0].Vy),
+	(md.results.StressbalanceSolution[0].Vel),
+	(md.results.StressbalanceSolution[0].Pressure),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test202.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test202.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test202.js	(revision 20593)
@@ -0,0 +1,20 @@
+//Test Name: SquareShelfStressSSA3d
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,2.);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vz','Vel','Pressure'];
+field_tolerances=[1e-13,1e-13,1e-13,1e-13,1e-13];
+field_values=[
+	(md.results.StressbalanceSolution[0].Vx),
+	(md.results.StressbalanceSolution[0].Vy),
+	(md.results.StressbalanceSolution[0].Vz),
+	(md.results.StressbalanceSolution[0].Vel),
+	(md.results.StressbalanceSolution[0].Pressure),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test203.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test203.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test203.js	(revision 20593)
@@ -0,0 +1,20 @@
+//Test Name: SquareShelfStressHO
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,2.);
+setflowequation(md,'HO','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vz','Vel','Pressure'];
+field_tolerances=[2e-09,1e-09,1e-09,1e-09,1e-09];
+field_values=[
+	(md.results.StressbalanceSolution[0].Vx),
+	(md.results.StressbalanceSolution[0].Vy),
+	(md.results.StressbalanceSolution[0].Vz),
+	(md.results.StressbalanceSolution[0].Vel),
+	(md.results.StressbalanceSolution[0].Pressure),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test204.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test204.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test204.js	(revision 20593)
@@ -0,0 +1,29 @@
+//Test Name: SquareShelfStressFS
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,2.);
+setflowequation(md,'FS','all');
+//md.cluster=generic('name',oshostname(),'np',1);
+md.stressbalance.shelf_dampening=1;
+md.timestepping.time_step=0;
+md1=solve(md,StressbalanceSolutionEnum());
+md.stressbalance.shelf_dampening=0;
+md=solve(md,StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vz','Vel','Pressure','Vx_damp','Vy_damp','Vz_damp','Vel_damp','Pressure_damp'];
+field_tolerances=[1e-08,1e-08,2e-06,1e-08,1e-08,1e-08,1e-08,2e-07,1e-08,1e-08];
+field_values=[
+	(md.results.StressbalanceSolution[0].Vx),
+	(md.results.StressbalanceSolution[0].Vy),
+	(md.results.StressbalanceSolution[0].Vz),
+	(md.results.StressbalanceSolution[0].Vel),
+	(md.results.StressbalanceSolution[0].Pressure),
+	(md1.results.StressbalanceSolution.Vx),
+	(md1.results.StressbalanceSolution.Vy),
+	(md1.results.StressbalanceSolution.Vz),
+	(md1.results.StressbalanceSolution.Vel),
+	(md1.results.StressbalanceSolution.Pressure),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test205.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test205.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test205.js	(revision 20593)
@@ -0,0 +1,20 @@
+//Test Name: SquareShelfStressMHOPenalties
+var md = new model();
+triangle(md,square[0],150000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,2.);
+setflowequation(md,'HO','../Exp/SquareHalfRight.exp','fill','SSA','coupling','penalties');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,StressbalanceSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx','Vy','Vz','Vel','Pressure'];
+field_tolerances=[2e-05,2e-05,1e-05,1e-05,1e-05];
+field_values=[
+	(md.results.StressbalanceSolution[0].Vx),
+	(md.results.StressbalanceSolution[0].Vy),
+	(md.results.StressbalanceSolution[0].Vz),
+	(md.results.StressbalanceSolution[0].Vel),
+	(md.results.StressbalanceSolution[0].Pressure),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test206.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test206.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test206.js	(revision 20593)
@@ -0,0 +1,18 @@
+//Test Name: SquareShelfTherStea
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,1.);
+setflowequation(md,'SSA','all');
+md.timestepping.time_step=0;
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,ThermalSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Temperature','BasalforcingsGroundediceMeltingRate'];
+field_tolerances=[1e-13,5e-6];
+field_values=[
+	(md.results.ThermalSolution[0].Temperature),
+	(md.results.ThermalSolution[0].BasalforcingsGroundediceMeltingRate),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test207.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test207.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test207.js	(revision 20593)
@@ -0,0 +1,26 @@
+//Test Name: SquareShelfTherTran
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,1.);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md.transient.isstressbalance=0;
+md.transient.ismasstransport=0;
+md.transient.issmb=1;
+md.transient.isthermal=1;
+md.transient.isgroundingline=0;
+md=solve(md,TransientSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Temperature1','BasalforcingsGroundediceMeltingRate1','Temperature2','BasalforcingsGroundediceMeltingRate2','Temperature3','BasalforcingsGroundediceMeltingRate3'];
+field_tolerances=[1e-13,1e-6,1e-13,1e-6,1e-13,1e-6];
+field_values=[
+	(md.results.TransientSolution[0](1).Temperature),
+	(md.results.TransientSolution[0](1).BasalforcingsGroundediceMeltingRate),
+	(md.results.TransientSolution[0](2).Temperature),
+	(md.results.TransientSolution[0](2).BasalforcingsGroundediceMeltingRate),
+	(md.results.TransientSolution[0](3).Temperature),
+	(md.results.TransientSolution[0](3).BasalforcingsGroundediceMeltingRate),
+	];
Index: /issm/trunk-jpl/test/NightlyRun/test209.js
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test209.js	(revision 20593)
+++ /issm/trunk-jpl/test/NightlyRun/test209.js	(revision 20593)
@@ -0,0 +1,50 @@
+//Test Name: SquareShelfTranSSA3d
+var md = new model();
+triangle(md,square[0],180000.);
+setmask(md,'all','');
+parameterize(md);
+md.extrude(md,3,1.);
+setflowequation(md,'SSA','all');
+//md.cluster=generic('name',oshostname(),'np',3);
+md=solve(md,TransientSolutionEnum());
+
+//Fields and tolerances to track changes
+field_names     =['Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','BasalforcingsGroundediceMeltingRate1', 
+	'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','BasalforcingsGroundediceMeltingRate2', 
+	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsGroundediceMeltingRate3'];
+field_tolerances=[
+	1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,
+	1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,2e-8,
+	5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,5e-8];
+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](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](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),
+	];
