Index: /issm/trunk/test/NightlyRun/test301.m
===================================================================
--- /issm/trunk/test/NightlyRun/test301.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test301.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test301_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test301_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test301_nightly.m	(revision 4975)
@@ -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/test302.m
===================================================================
--- /issm/trunk/test/NightlyRun/test302.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test302.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test302_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test302_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test302_nightly.m	(revision 4975)
@@ -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/test303.m
===================================================================
--- /issm/trunk/test/NightlyRun/test303.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test303.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test303_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test303_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test303_nightly.m	(revision 4975)
@@ -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/test304.m
===================================================================
--- /issm/trunk/test/NightlyRun/test304.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test304.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test304_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test304_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test304_nightly.m	(revision 4975)
@@ -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/test305.m
===================================================================
--- /issm/trunk/test/NightlyRun/test305.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test305.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test305_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test305_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test305_nightly.m	(revision 4975)
@@ -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/test306.m
===================================================================
--- /issm/trunk/test/NightlyRun/test306.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test306.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test306_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test306_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test306_nightly.m	(revision 4975)
@@ -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/test307.m
===================================================================
--- /issm/trunk/test/NightlyRun/test307.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test307.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test307_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test307_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test307_nightly.m	(revision 4975)
@@ -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/test308.m
===================================================================
--- /issm/trunk/test/NightlyRun/test308.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test308.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test308_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test308_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test308_nightly.m	(revision 4975)
@@ -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/test309.m
===================================================================
--- /issm/trunk/test/NightlyRun/test309.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test309.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,5,2);
+md=setelementstype(md,'hutter','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test309_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test309_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test309_nightly.m	(revision 4975)
@@ -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/test310.m
===================================================================
--- /issm/trunk/test/NightlyRun/test310.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test310.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,5,2);
+md=setelementstype(md,'hutter','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test310_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test310_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test310_nightly.m	(revision 4975)
@@ -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/test311.m
===================================================================
--- /issm/trunk/test/NightlyRun/test311.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test311.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,4,2);
+md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test311_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test311_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test311_nightly.m	(revision 4975)
@@ -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/test312.m
===================================================================
--- /issm/trunk/test/NightlyRun/test312.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test312.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,4,2);
+md=setelementstype(md,'hutter','../Exp/SquareHalfRight.exp','fill','macayeal');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test312_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test312_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test312_nightly.m	(revision 4975)
@@ -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/test313.m
===================================================================
--- /issm/trunk/test/NightlyRun/test313.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test313.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test313_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test313_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test313_nightly.m	(revision 4975)
@@ -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/test314.m
===================================================================
--- /issm/trunk/test/NightlyRun/test314.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test314.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test314_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test314_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test314_nightly.m	(revision 4975)
@@ -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/test315.m
===================================================================
--- /issm/trunk/test/NightlyRun/test315.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test315.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all','stokes','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test315_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test315_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test315_nightly.m	(revision 4975)
@@ -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/test316.m
===================================================================
--- /issm/trunk/test/NightlyRun/test316.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test316.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all','stokes','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test316_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test316_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test316_nightly.m	(revision 4975)
@@ -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/test317.m
===================================================================
--- /issm/trunk/test/NightlyRun/test317.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test317.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test317_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test317_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test317_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test318.m
===================================================================
--- /issm/trunk/test/NightlyRun/test318.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test318.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test318_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test318_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test318_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test319.m
===================================================================
--- /issm/trunk/test/NightlyRun/test319.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test319.m	(revision 4975)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.prognostic_DG=1;
+md.cluster='none';
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test319_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test319_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test319_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test320.m
===================================================================
--- /issm/trunk/test/NightlyRun/test320.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test320.m	(revision 4975)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=meshconvert(md);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.prognostic_DG=1;
+md.cluster=oshostname();
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test320_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test320_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test320_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test321.m
===================================================================
--- /issm/trunk/test/NightlyRun/test321.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test321.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md=extrude(md,5,0.5);
+md.cluster='none';
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test321_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test321_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test321_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test322.m
===================================================================
--- /issm/trunk/test/NightlyRun/test322.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test322.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md=extrude(md,5,0.5);
+md.cluster=oshostname();
+md=solve(md,'analysis_type',PrognosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test322_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test322_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test322_nightly.m	(revision 4975)
@@ -0,0 +1,5 @@
+field_names     ={'Thickness'};
+field_tolerances={1e-13};
+field_values={...
+	PatchToVec(md.results.PrognosticSolution.Thickness),...
+	};
Index: /issm/trunk/test/NightlyRun/test323.m
===================================================================
--- /issm/trunk/test/NightlyRun/test323.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test323.m	(revision 4975)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test323_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test323_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test323_nightly.m	(revision 4975)
@@ -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/test324.m
===================================================================
--- /issm/trunk/test/NightlyRun/test324.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test324.m	(revision 4975)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test324_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test324_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test324_nightly.m	(revision 4975)
@@ -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/test325.m
===================================================================
--- /issm/trunk/test/NightlyRun/test325.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test325.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test325_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test325_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test325_nightly.m	(revision 4975)
@@ -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/test326.m
===================================================================
--- /issm/trunk/test/NightlyRun/test326.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test326.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test326_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test326_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test326_nightly.m	(revision 4975)
@@ -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/test327.m
===================================================================
--- /issm/trunk/test/NightlyRun/test327.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test327.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test327_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test327_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test327_nightly.m	(revision 4975)
@@ -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/test328.m
===================================================================
--- /issm/trunk/test/NightlyRun/test328.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test328.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'hutter','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test328_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test328_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test328_nightly.m	(revision 4975)
@@ -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/test329.m
===================================================================
--- /issm/trunk/test/NightlyRun/test329.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test329.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,5,1.2);
+md=setelementstype(md,'hutter','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test329_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test329_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test329_nightly.m	(revision 4975)
@@ -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/test330.m
===================================================================
--- /issm/trunk/test/NightlyRun/test330.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test330.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,5,1.2);
+md=setelementstype(md,'hutter','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test330_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test330_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test330_nightly.m	(revision 4975)
@@ -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/test331.m
===================================================================
--- /issm/trunk/test/NightlyRun/test331.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test331.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test331_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test331_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test331_nightly.m	(revision 4975)
@@ -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/test332.m
===================================================================
--- /issm/trunk/test/NightlyRun/test332.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test332.m	(revision 4975)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test332_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test332_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test332_nightly.m	(revision 4975)
@@ -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/test333.m
===================================================================
--- /issm/trunk/test/NightlyRun/test333.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test333.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test333_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test333_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test333_nightly.m	(revision 4975)
@@ -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/test334.m
===================================================================
--- /issm/trunk/test/NightlyRun/test334.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test334.m	(revision 4975)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',200000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.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/test334_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test334_nightly.m	(revision 4975)
+++ /issm/trunk/test/NightlyRun/test334_nightly.m	(revision 4975)
@@ -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/Par/SquareSheetConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 4974)
+++ /issm/trunk/test/Par/SquareSheetConstrained.par	(revision 4975)
@@ -33,5 +33,5 @@
 
 %Numerical parameters
-md.viscosity_overshoot=0.3;
+md.viscosity_overshoot=0.0;
 md.artificial_diffusivity=1;
 md.verbose=0;
@@ -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;
Index: /issm/trunk/test/Par/SquareShelfConstrained.par
===================================================================
--- /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4974)
+++ /issm/trunk/test/Par/SquareShelfConstrained.par	(revision 4975)
@@ -24,4 +24,8 @@
 md.temperature=md.observed_temperature;
 
+%Accumulation and melting
+md.accumulation_rate=10*ones(md.numberofgrids,1);
+md.melting_rate=5*ones(md.numberofgrids,1);
+
 %Friction
 pos=find(md.elementoniceshelf);
@@ -33,5 +37,5 @@
 
 %Numerical parameters
-md.viscosity_overshoot=0.3;
+md.viscosity_overshoot=0.0;
 md.artificial_diffusivity=1;
 md.verbose=0;
