Index: /issm/trunk/test/Exp/SquareShelf.exp
===================================================================
--- /issm/trunk/test/Exp/SquareShelf.exp	(revision 4996)
+++ /issm/trunk/test/Exp/SquareShelf.exp	(revision 4996)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 500000
+1000000 500000
+1000000 1000000
+0 1000000 
+0 500000
Index: /issm/trunk/test/NightlyRun/Id2Name.m
===================================================================
--- /issm/trunk/test/NightlyRun/Id2Name.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/Id2Name.m	(revision 4996)
@@ -105,4 +105,71 @@
 elseif (id==337), name='SquareSheetConstrainedCMDrag3dSerial';
 elseif (id==338), name='SquareSheetConstrainedCMDrag3dParallel';
-elseif (id==702), name='SquareShelf...';
+elseif (id==339), name='SquareSheetConstrainedSteaH3dSerial';
+elseif (id==340), name='SquareSheetConstrainedSteaH3dParallel';
+elseif (id==401), name='SquareSheetShelfDiagM2dSerial';
+elseif (id==402), name='SquareSheetShelfDiagM2dParallel';
+elseif (id==403), name='SquareSheetShelfDiagM3dSerial';
+elseif (id==404), name='SquareSheetShelfDiagM3dParallel';
+elseif (id==405), name='SquareSheetShelfDiagP3dSerial';
+elseif (id==406), name='SquareSheetShelfDiagP3dParallel';
+elseif (id==407), name='SquareSheetShelfDiagS3dSerial';
+elseif (id==408), name='SquareSheetShelfDiagS3dParallel';
+elseif (id==409), name='SquareSheetShelfDiagMP3dSerial';
+elseif (id==410), name='SquareSheetShelfDiagMP3dParallel';
+elseif (id==411), name='SquareSheetShelfTherSteaSerial';
+elseif (id==412), name='SquareSheetShelfTherSteaParallel';
+elseif (id==413), name='SquareSheetShelfTherTranSerial';
+elseif (id==414), name='SquareSheetShelfTherTranParallel';
+elseif (id==415), name='SquareSheetShelfTranM2dSerial';
+elseif (id==416), name='SquareSheetShelfTranM2dParallel';
+elseif (id==417), name='SquareSheetShelfTranMP3dSerial';
+elseif (id==418), name='SquareSheetShelfTranMP3dParallel';
+elseif (id==419), name='SquareSheetShelfCMSteaM3dSerial';
+elseif (id==420), name='SquareSheetShelfCMSteaM3dParallel';
+elseif (id==421), name='SquareSheetShelfCMSteaP3dSerial';
+elseif (id==422), name='SquareSheetShelfCMSteaP3dParallel';
+elseif (id==423), name='SquareSheetShelfSteaM3dSerial';
+elseif (id==424), name='SquareSheetShelfSteaM3dParallel';
+elseif (id==425), name='SquareSheetShelfSteaP3dSerial';
+elseif (id==426), name='SquareSheetShelfSteaP3dParallel';
+elseif (id==501), name='PigDiagM2dSerial';
+elseif (id==502), name='PigDiagM2dParallel';
+elseif (id==503), name='PigDiagP3dSerial';
+elseif (id==504), name='PigDiagP3dParallel';
+elseif (id==505), name='PigDiagS3dSerial';
+elseif (id==506), name='PigDiagS3dParallel';
+elseif (id==507), name='PigTranM2dSerial';
+elseif (id==508), name='PigTranM2dParallel';
+elseif (id==509), name='PigTranM3dSerial';
+elseif (id==510), name='PigTranM3dParallel';
+elseif (id==511), name='PigTranP3dSerial';
+elseif (id==512), name='PigTranP3dParallel';
+elseif (id==513), name='PigTranS3dSerial';
+elseif (id==514), name='PigTranS3dParallel';
+elseif (id==515), name='PigSteaM3dSerial';
+elseif (id==516), name='PigSteaM3dParallel';
+elseif (id==517), name='PigSteaP3dSerial';
+elseif (id==518), name='PigSteaP3dParallel';
+elseif (id==519), name='PigSteaS3dSerial';
+elseif (id==520), name='PigSteaS3dParallel';
+elseif (id==601), name='79NorthProg2dSerial';
+elseif (id==602), name='79NorthProg2dParallel';
+elseif (id==603), name='79NorthProg2dDGSerial';
+elseif (id==604), name='79NorthProg2dDGParallel';
+elseif (id==605), name='79NorthProg3dSerial';
+elseif (id==606), name='79NorthProg3dParallel';
+elseif (id==607), name='79NorthSurfSlop2dSerial';
+elseif (id==608), name='79NorthSurfSlop2dParallel';
+elseif (id==609), name='79NorthSurfSlop3dSerial';
+elseif (id==610), name='79NorthSurfSlop3dParallel';
+elseif (id==611), name='79NorthBedSlop2dSerial';
+elseif (id==612), name='79NorthBedSlop2dParallel';
+elseif (id==613), name='79NorthBedSlop3dSerial';
+elseif (id==614), name='79NorthBedSlop3dParallel';
+elseif (id==615), name='79NorthBalThic2dSerial';
+elseif (id==616), name='79NorthBalThic2dParallel';
+elseif (id==617), name='79NorthBalThic2dDGSerial';
+elseif (id==618), name='79NorthBalThic2dDGParallel';
+elseif (id==619), name='79NorthBalThic3dSerial';
+elseif (id==620), name='79NorthBalThic3dParallel';
 else name='N/A'; end
Index: /issm/trunk/test/NightlyRun/test106.m
===================================================================
--- /issm/trunk/test/NightlyRun/test106.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test106.m	(revision 4996)
@@ -3,5 +3,4 @@
 md=parameterize(md,'../Par/SquareShelfConstrained.par');
 md=extrude(md,3,2);
-md.verbose=3;
 md=setelementstype(md,'pattyn','all');
 md.cluster=oshostname;
Index: /issm/trunk/test/NightlyRun/test107.m
===================================================================
--- /issm/trunk/test/NightlyRun/test107.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test107.m	(revision 4996)
@@ -4,5 +4,4 @@
 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/test121_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test121_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test121_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test122_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test122_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test122_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test217_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test217_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test217_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test218_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test218_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test218_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test219_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test219_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test219_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test220_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test220_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test220_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test221_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test221_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test221_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test222_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test222_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test222_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test329_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test329_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test329_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test330_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test330_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test330_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test333_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test333_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test333_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test334_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test334_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test334_nightly.m	(revision 4996)
@@ -1,9 +1,10 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +16,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -24,4 +26,5 @@
 	PatchToVec(md.results.Transient3DSolution(3).Vx),...
 	PatchToVec(md.results.Transient3DSolution(3).Vy),...
+	PatchToVec(md.results.Transient3DSolution(3).Vz),...
 	PatchToVec(md.results.Transient3DSolution(3).Vel),...
 	PatchToVec(md.results.Transient3DSolution(3).Pressure),...
Index: /issm/trunk/test/NightlyRun/test339.m
===================================================================
--- /issm/trunk/test/NightlyRun/test339.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test339.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'hutter','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test339_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test339_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test339_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test340.m
===================================================================
--- /issm/trunk/test/NightlyRun/test340.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test340.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'','');
+md=parameterize(md,'../Par/SquareSheetConstrained.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'hutter','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test340_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test340_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test340_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test401.m
===================================================================
--- /issm/trunk/test/NightlyRun/test401.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test401.m	(revision 4996)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test401_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test401_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test401_nightly.m	(revision 4996)
@@ -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/test402.m
===================================================================
--- /issm/trunk/test/NightlyRun/test402.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test402.m	(revision 4996)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test402_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test402_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test402_nightly.m	(revision 4996)
@@ -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/test403.m
===================================================================
--- /issm/trunk/test/NightlyRun/test403.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test403.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test403_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test403_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test403_nightly.m	(revision 4996)
@@ -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/test404.m
===================================================================
--- /issm/trunk/test/NightlyRun/test404.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test404.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test404_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test404_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test404_nightly.m	(revision 4996)
@@ -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/test405.m
===================================================================
--- /issm/trunk/test/NightlyRun/test405.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test405.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test405_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test405_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test405_nightly.m	(revision 4996)
@@ -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/test406.m
===================================================================
--- /issm/trunk/test/NightlyRun/test406.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test406.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test406_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test406_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test406_nightly.m	(revision 4996)
@@ -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/test407.m
===================================================================
--- /issm/trunk/test/NightlyRun/test407.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test407.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.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/test407_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test407_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test407_nightly.m	(revision 4996)
@@ -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/test408.m
===================================================================
--- /issm/trunk/test/NightlyRun/test408.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test408.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.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/test408_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test408_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test408_nightly.m	(revision 4996)
@@ -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/test409.m
===================================================================
--- /issm/trunk/test/NightlyRun/test409.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test409.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
+md.cluster='none';
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test409_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test409_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test409_nightly.m	(revision 4996)
@@ -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/test410.m
===================================================================
--- /issm/trunk/test/NightlyRun/test410.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test410.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,5,1);
+md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',DiagnosticSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test410_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test410_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test410_nightly.m	(revision 4996)
@@ -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/test411.m
===================================================================
--- /issm/trunk/test/NightlyRun/test411.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test411.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test411_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test411_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test411_nightly.m	(revision 4996)
@@ -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/test412.m
===================================================================
--- /issm/trunk/test/NightlyRun/test412.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test412.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test412_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test412_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test412_nightly.m	(revision 4996)
@@ -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/test413.m
===================================================================
--- /issm/trunk/test/NightlyRun/test413.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test413.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test413_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test413_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test413_nightly.m	(revision 4996)
@@ -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/test414.m
===================================================================
--- /issm/trunk/test/NightlyRun/test414.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test414.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,4,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',ThermalSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test414_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test414_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test414_nightly.m	(revision 4996)
@@ -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/test415.m
===================================================================
--- /issm/trunk/test/NightlyRun/test415.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test415.m	(revision 4996)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test415_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test415_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test415_nightly.m	(revision 4996)
@@ -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/test416.m
===================================================================
--- /issm/trunk/test/NightlyRun/test416.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test416.m	(revision 4996)
@@ -0,0 +1,6 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient2DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test416_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test416_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test416_nightly.m	(revision 4996)
@@ -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/test417.m
===================================================================
--- /issm/trunk/test/NightlyRun/test417.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test417.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
+md.cluster='none';
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test417_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test417_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test417_nightly.m	(revision 4996)
@@ -0,0 +1,36 @@
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
+	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).Vz),...
+	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).Vz),...
+	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/test418.m
===================================================================
--- /issm/trunk/test/NightlyRun/test418.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test418.m	(revision 4996)
@@ -0,0 +1,7 @@
+md=mesh(model,'../Exp/Square.exp',180000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'macayeal','../Exp/SquareHalfRight.exp','fill','pattyn');
+md.cluster=oshostname();
+md=solve(md,'analysis_type',Transient3DSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test418_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test418_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test418_nightly.m	(revision 4996)
@@ -0,0 +1,36 @@
+field_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2', ...
+					   'Vx3','Vy3','Vz3','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,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).Vz),...
+	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).Vz),...
+	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).Vz),...
+	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/test423.m
===================================================================
--- /issm/trunk/test/NightlyRun/test423.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test423.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test423_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test423_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test423_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test424.m
===================================================================
--- /issm/trunk/test/NightlyRun/test424.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test424.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test424_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test424_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test424_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test425.m
===================================================================
--- /issm/trunk/test/NightlyRun/test425.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test425.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test425_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test425_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test425_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test426.m
===================================================================
--- /issm/trunk/test/NightlyRun/test426.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test426.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Square.exp',150000);
+md=geography(md,'../Exp/SquareShelf.exp','');
+md=parameterize(md,'../Par/SquareSheetShelf.par');
+md=extrude(md,3,2);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test426_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test426_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test426_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test507_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test507_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test507_nightly.m	(revision 4996)
@@ -1,4 +1,4 @@
-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_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient2DSolution(1).Vx),...
@@ -16,10 +16,3 @@
 	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/test508_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test508_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test508_nightly.m	(revision 4996)
@@ -1,4 +1,4 @@
-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_names     ={'Vx1','Vy1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Vx2','Vy2','Vel2','Pressure2','Bed2','Surface2','Thickness2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient2DSolution(1).Vx),...
@@ -16,10 +16,3 @@
 	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/test509_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test509_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test509_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test510_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test510_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test510_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test511_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test511_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test511_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test512_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test512_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test512_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test513_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test513_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test513_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test514_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test514_nightly.m	(revision 4995)
+++ /issm/trunk/test/NightlyRun/test514_nightly.m	(revision 4996)
@@ -1,9 +1,9 @@
-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_names     ={'Vx1','Vy1','Vz1','Vel1','Pressure1','Bed1','Surface1','Thickness1','Temperature1','MeltingRate1', ...
+				      'Vx2','Vy2','Vz2','Vel2','Pressure2','Bed2','Surface2','Thickness2','Temperature2','MeltingRate2'};
+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};
 field_values={...
 	PatchToVec(md.results.Transient3DSolution(1).Vx),...
 	PatchToVec(md.results.Transient3DSolution(1).Vy),...
+	PatchToVec(md.results.Transient3DSolution(1).Vz),...
 	PatchToVec(md.results.Transient3DSolution(1).Vel),...
 	PatchToVec(md.results.Transient3DSolution(1).Pressure),...
@@ -15,4 +15,5 @@
 	PatchToVec(md.results.Transient3DSolution(2).Vx),...
 	PatchToVec(md.results.Transient3DSolution(2).Vy),...
+	PatchToVec(md.results.Transient3DSolution(2).Vz),...
 	PatchToVec(md.results.Transient3DSolution(2).Vel),...
 	PatchToVec(md.results.Transient3DSolution(2).Pressure),...
@@ -22,12 +23,3 @@
 	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/test515.m
===================================================================
--- /issm/trunk/test/NightlyRun/test515.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test515.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,4,1.1);
+md=setelementstype(md,'macayeal','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test515_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test515_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test515_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test516.m
===================================================================
--- /issm/trunk/test/NightlyRun/test516.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test516.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,4,1.1);
+md=setelementstype(md,'macayeal','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test516_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test516_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test516_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test517.m
===================================================================
--- /issm/trunk/test/NightlyRun/test517.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test517.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test517_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test517_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test517_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test518.m
===================================================================
--- /issm/trunk/test/NightlyRun/test518.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test518.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,3,1);
+md=setelementstype(md,'pattyn','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test518_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test518_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test518_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test519.m
===================================================================
--- /issm/trunk/test/NightlyRun/test519.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test519.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,2,1);
+md=setelementstype(md,'macayeal','all','stokes','all');
+md.cluster='none';
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test519_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test519_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test519_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/NightlyRun/test520.m
===================================================================
--- /issm/trunk/test/NightlyRun/test520.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test520.m	(revision 4996)
@@ -0,0 +1,8 @@
+md=mesh(model,'../Exp/Pig.exp',80000);
+md=geography(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp');
+md=parameterize(md,'../Par/Pig.par');
+md=extrude(md,2,1);
+md=setelementstype(md,'macayeal','all','stokes','all');
+md.cluster=oshostname();
+md.dt=0;
+md=solve(md,'analysis_type',SteadystateSolutionEnum);
Index: /issm/trunk/test/NightlyRun/test520_nightly.m
===================================================================
--- /issm/trunk/test/NightlyRun/test520_nightly.m	(revision 4996)
+++ /issm/trunk/test/NightlyRun/test520_nightly.m	(revision 4996)
@@ -0,0 +1,12 @@
+field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','MeltingRate'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13
+};
+field_values={...
+	PatchToVec(md.results.SteadystateSolution.Vx),...
+	PatchToVec(md.results.SteadystateSolution.Vy),...
+	PatchToVec(md.results.SteadystateSolution.Vz),...
+	PatchToVec(md.results.SteadystateSolution.Vel),...
+	PatchToVec(md.results.SteadystateSolution.Pressure),...
+	PatchToVec(md.results.SteadystateSolution.Temperature),...
+	PatchToVec(md.results.SteadystateSolution.MeltingRate),...
+	};
Index: /issm/trunk/test/Par/Pig.par
===================================================================
--- /issm/trunk/test/Par/Pig.par	(revision 4995)
+++ /issm/trunk/test/Par/Pig.par	(revision 4996)
@@ -37,5 +37,5 @@
 md.ndt=2;
 md.eps_res=0.05;
-md.eps_rel=0.05;
+md.eps_rel=1;
 md.eps_abs=NaN;
 
Index: /issm/trunk/test/Par/SquareSheetShelf.par
===================================================================
--- /issm/trunk/test/Par/SquareSheetShelf.par	(revision 4996)
+++ /issm/trunk/test/Par/SquareSheetShelf.par	(revision 4996)
@@ -0,0 +1,54 @@
+%Start defining model parameters here
+
+%Geometry
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+bed_sheet=-md.rho_ice/md.rho_water*(hmax+(hmin-hmax)*(ymax/2-ymin)/(ymax-ymin));
+pos=find(md.y<=ymax/2);
+md.bed(pos)=bed_sheet;
+md.surface=md.bed+md.thickness;
+
+%Initial velocity 
+load('./../Data/SquareSheetShelf.data','-mat');
+md.vx=InterpFromMeshToMesh2d(index,x,y,vx,md.x,md.y,0);
+md.vy=InterpFromMeshToMesh2d(index,x,y,vy,md.x,md.y,0);
+clear vx vy x y index;
+md.vz=zeros(md.numberofgrids,1);
+md.pressure=zeros(md.numberofgrids,1);
+
+%Materials
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+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);
+md.drag_type=2;
+md.drag_coefficient=20*ones(md.numberofgrids,1);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+
+%Numerical parameters
+md.viscosity_overshoot=0.0;
+md.artificial_diffusivity=1;
+md.verbose=0;
+md.waitonlock=30;
+md.np=3;
+md.eps_res=0.05;
+md.eps_rel=0.05;
+md.eps_abs=NaN;
+md.dt=1;
+md.ndt=3;
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'./../Exp/SquareFront.exp');
