Index: /issm/trunk/test/NightlyRun/test101.m
===================================================================
--- /issm/trunk/test/NightlyRun/test101.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test101.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test101_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test101_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test101_nightly.m	(revision 4963)
@@ -0,0 +1,8 @@
+field_names     ={'Vx','Vy','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test102.m
===================================================================
--- /issm/trunk/test/NightlyRun/test102.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test102.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test102_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test102_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test102_nightly.m	(revision 4963)
@@ -0,0 +1,8 @@
+field_names     ={'Vx','Vy','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test103.m
===================================================================
--- /issm/trunk/test/NightlyRun/test103.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test103.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test103_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test103_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test103_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test104.m
===================================================================
--- /issm/trunk/test/NightlyRun/test104.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test104.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test104_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test104_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test104_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test105.m
===================================================================
--- /issm/trunk/test/NightlyRun/test105.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test105.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test105_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test105_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test105_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test106.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test106.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md.verbose=3;
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test106_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test106_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test106_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test107.m
===================================================================
--- /issm/trunk/test/NightlyRun/test107.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test107.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all','stokes','all');
+md.verbose=3;
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test107_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test107_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test107_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test108.m
===================================================================
--- /issm/trunk/test/NightlyRun/test108.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test108.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all','stokes','all');
+md.cluster=oshostname
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test108_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test108_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test108_nightly.m	(revision 4963)
@@ -0,0 +1,9 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.DiagnosticSolution.Vx),...
+	PatchToVec(md.results.DiagnosticSolution.Vy),...
+	PatchToVec(md.results.DiagnosticSolution.Vz),...
+	PatchToVec(md.results.DiagnosticSolution.Vel),...
+	PatchToVec(md.results.DiagnosticSolution.Pressure),...
+	};
Index: /issm/trunk/test/NightlyRun/test109.m
===================================================================
--- /issm/trunk/test/NightlyRun/test109.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test109.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test109_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test109_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test109_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test110.m
===================================================================
--- /issm/trunk/test/NightlyRun/test110.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test110.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test110_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test110_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test110_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test111.m
===================================================================
--- /issm/trunk/test/NightlyRun/test111.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test111.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md.prognostic_DG=1
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test111_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test111_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test111_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test112.m
===================================================================
--- /issm/trunk/test/NightlyRun/test112.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test112.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md.prognostic_DG=1
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test112_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test112_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test112_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test113.m
===================================================================
--- /issm/trunk/test/NightlyRun/test113.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test113.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md=extrude(md,5,3);
+md.cluster='none';
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test113_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test113_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test113_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test114.m
===================================================================
--- /issm/trunk/test/NightlyRun/test114.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test114.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md=extrude(md,5,3);
+md.cluster=oshostname;
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test114_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test114_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test114_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test115.m
===================================================================
--- /issm/trunk/test/NightlyRun/test115.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test115.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+md.dt=0;
+md.cluster='none';
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test115_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test115_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test115_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test116.m
===================================================================
--- /issm/trunk/test/NightlyRun/test116.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test116.m	(revision 4963)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+md.dt=0;
+md.cluster=oshostname();
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test116_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test116_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test116_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test117.m
===================================================================
--- /issm/trunk/test/NightlyRun/test117.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test117.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test117_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test117_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test117_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test118.m
===================================================================
--- /issm/trunk/test/NightlyRun/test118.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test118.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test118_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test118_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test118_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test119.m
===================================================================
--- /issm/trunk/test/NightlyRun/test119.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test119.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test119_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test119_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test119_nightly.m	(revision 4963)
@@ -0,0 +1,25 @@
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test120.m
===================================================================
--- /issm/trunk/test/NightlyRun/test120.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test120.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test120_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test120_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test120_nightly.m	(revision 4963)
@@ -0,0 +1,25 @@
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient2DSolution(1).Vx),...
+	PatchToVec(md.results.Transient2DSolution(1).Vy),...
+	PatchToVec(md.results.Transient2DSolution(1).Vel),...
+	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(1).Bed),...
+	PatchToVec(md.results.Transient2DSolution(1).Surface),...
+	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(2).Vx),...
+	PatchToVec(md.results.Transient2DSolution(2).Vy),...
+	PatchToVec(md.results.Transient2DSolution(2).Vel),...
+	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(2).Bed),...
+	PatchToVec(md.results.Transient2DSolution(2).Surface),...
+	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient2DSolution(3).Vx),...
+	PatchToVec(md.results.Transient2DSolution(3).Vy),...
+	PatchToVec(md.results.Transient2DSolution(3).Vel),...
+	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient2DSolution(3).Bed),...
+	PatchToVec(md.results.Transient2DSolution(3).Surface),...
+	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test121.m
===================================================================
--- /issm/trunk/test/NightlyRun/test121.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test121.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test121_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test121_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test121_nightly.m	(revision 4963)
@@ -0,0 +1,33 @@
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test122.m
===================================================================
--- /issm/trunk/test/NightlyRun/test122.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test122.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test122_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test122_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test122_nightly.m	(revision 4963)
@@ -0,0 +1,33 @@
+field_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','MeltingRate3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test123.m
===================================================================
--- /issm/trunk/test/NightlyRun/test123.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test123.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test123_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test123_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test123_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'SurfaceSlopeX','SurfaceSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeX),...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test124.m
===================================================================
--- /issm/trunk/test/NightlyRun/test124.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test124.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test124_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test124_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test124_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'SurfaceSlopeX','SurfaceSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeX),...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test125.m
===================================================================
--- /issm/trunk/test/NightlyRun/test125.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test125.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test125_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test125_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test125_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'SurfaceSlopeX','SurfaceSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeX),...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test126.m
===================================================================
--- /issm/trunk/test/NightlyRun/test126.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test126.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',SurfaceSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test126_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test126_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test126_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'SurfaceSlopeX','SurfaceSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeX),...
+	PatchToVec(md.results.SurfaceSlopeSolution.SurfaceSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test127.m
===================================================================
--- /issm/trunk/test/NightlyRun/test127.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test127.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test127_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test127_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test127_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'BedSlopeX','BedSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeX),...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test128.m
===================================================================
--- /issm/trunk/test/NightlyRun/test128.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test128.m	(revision 4963)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test128_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test128_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test128_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'BedSlopeX','BedSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeX),...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test129.m
===================================================================
--- /issm/trunk/test/NightlyRun/test129.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test129.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test129_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test129_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test129_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'BedSlopeX','BedSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeX),...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test130.m
===================================================================
--- /issm/trunk/test/NightlyRun/test130.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test130.m	(revision 4963)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',BedSlopeSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test130_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test130_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test130_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'BedSlopeX','BedSlopeY'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeX),...
+	PatchToVec(md.results.BedSlopeSolution.BedSlopeY),...
+	};
Index: /issm/trunk/test/NightlyRun/test131.m
===================================================================
--- /issm/trunk/test/NightlyRun/test131.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test131.m	(revision 4963)
@@ -0,0 +1,10 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test131_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test131_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test131_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test132.m
===================================================================
--- /issm/trunk/test/NightlyRun/test132.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test132.m	(revision 4963)
@@ -0,0 +1,10 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test132_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test132_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test132_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test133.m
===================================================================
--- /issm/trunk/test/NightlyRun/test133.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test133.m	(revision 4963)
@@ -0,0 +1,11 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test133_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test133_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test133_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test134.m
===================================================================
--- /issm/trunk/test/NightlyRun/test134.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test134.m	(revision 4963)
@@ -0,0 +1,11 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+md=extrude(md,5,1);
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test134_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test134_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test134_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test135.m
===================================================================
--- /issm/trunk/test/NightlyRun/test135.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test135.m	(revision 4963)
@@ -0,0 +1,12 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md.prognostic_DG=1;
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test135_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test135_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test135_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test136.m
===================================================================
--- /issm/trunk/test/NightlyRun/test136.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test136.m	(revision 4963)
@@ -0,0 +1,12 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'all','');
+md=parameterize(md,'../Par/SquareShelfConstrained.par');
+%Add boundary conditions on thickness on the border
+pos=find(md.gridonboundary);
+md.spcthickness(pos,1)=1;
+md.spcthickness(pos,2)=md.thickness(pos);
+md.prognostic_DG=1;
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname;
+md=solve(md,'analysis_type',BalancedthicknessSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test136_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test136_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test136_nightly.m	(revision 4963)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.BalancedthicknessSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test205.m
===================================================================
--- /issm/trunk/test/NightlyRun/test205.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test205.m	(revision 4963)
@@ -4,5 +4,4 @@
 md=extrude(md,3,2);
 md=setelementstype(md,'pattyn','all');
-md.verbose=2;
 md.cluster='none';
 md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: sm/trunk/test/NightlyRun/test211_nighlty.m
===================================================================
--- /issm/trunk/test/NightlyRun/test211_nighlty.m	(revision 4962)
+++ 	(revision )
@@ -1,6 +1,0 @@
-field_names     ={'Temperature','MeltingRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-	PatchToVec(md.results.ThermalSolution.Temperature),...
-	PatchToVec(md.results.ThermalSolution.MeltingRate),...
-	};
Index: /issm/trunk/test/NightlyRun/test211_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test211_nightly.m	(revision 4963)
+++ /issm/trunk/test/NightlyRun/test211_nightly.m	(revision 4963)
@@ -0,0 +1,6 @@
+field_names     ={'Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13};
+field_values={...
+	PatchToVec(md.results.ThermalSolution.Temperature),...
+	PatchToVec(md.results.ThermalSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test217_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test217_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test217_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/NightlyRun/test218_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test218_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test218_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/NightlyRun/test219_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test219_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test219_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/NightlyRun/test220_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test220_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test220_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/NightlyRun/test221_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test221_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test221_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/NightlyRun/test222_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test222_nightly.m	(revision 4962)
+++ /issm/trunk/test/NightlyRun/test222_nightly.m	(revision 4963)
@@ -4,30 +4,30 @@
 field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-	PatchToVec(md.results.Transient2DSolution(1).Vx),...
-	PatchToVec(md.results.Transient2DSolution(1).Vy),...
-	PatchToVec(md.results.Transient2DSolution(1).Vel),...
-	PatchToVec(md.results.Transient2DSolution(1).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(1).Bed),...
-	PatchToVec(md.results.Transient2DSolution(1).Surface),...
-	PatchToVec(md.results.Transient2DSolution(1).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(1).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(1).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(2).Vx),...
-	PatchToVec(md.results.Transient2DSolution(2).Vy),...
-	PatchToVec(md.results.Transient2DSolution(2).Vel),...
-	PatchToVec(md.results.Transient2DSolution(2).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(2).Bed),...
-	PatchToVec(md.results.Transient2DSolution(2).Surface),...
-	PatchToVec(md.results.Transient2DSolution(2).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(2).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(2).MeltingRate),...
-	PatchToVec(md.results.Transient2DSolution(3).Vx),...
-	PatchToVec(md.results.Transient2DSolution(3).Vy),...
-	PatchToVec(md.results.Transient2DSolution(3).Vel),...
-	PatchToVec(md.results.Transient2DSolution(3).Pressure),...
-	PatchToVec(md.results.Transient2DSolution(3).Bed),...
-	PatchToVec(md.results.Transient2DSolution(3).Surface),...
-	PatchToVec(md.results.Transient2DSolution(3).Thickness),...
-	PatchToVec(md.results.Transient2DSolution(3).Temperature),...
-	PatchToVec(md.results.Transient2DSolution(3).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(1).Vx),...
+	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vel),...
+	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(1).Bed),...
+	PatchToVec(md.results.Transient3DSolution(1).Surface),...
+	PatchToVec(md.results.Transient3DSolution(1).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(1).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(1).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(2).Vx),...
+	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vel),...
+	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(2).Bed),...
+	PatchToVec(md.results.Transient3DSolution(2).Surface),...
+	PatchToVec(md.results.Transient3DSolution(2).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(2).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(2).MeltingRate),...
+	PatchToVec(md.results.Transient3DSolution(3).Vx),...
+	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vel),...
+	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
+	PatchToVec(md.results.Transient3DSolution(3).Bed),...
+	PatchToVec(md.results.Transient3DSolution(3).Surface),...
+	PatchToVec(md.results.Transient3DSolution(3).Thickness),...
+	PatchToVec(md.results.Transient3DSolution(3).Temperature),...
+	PatchToVec(md.results.Transient3DSolution(3).MeltingRate),...
 	};
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4962)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4963)
@@ -39,5 +39,5 @@
 md.np=3;
 md.eps_res=0.05;
-md.eps_rel=0.005;
+md.eps_rel=0.05;
 md.eps_abs=NaN;
 md.dt=1;
