Index: ../trunk-jpl/test/NightlyRun/test410.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test410.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test410.m	(revision 17747)
@@ -9,7 +9,7 @@
 
 %Fields and tolerances to track changes
 field_names     ={'Vx','Vy','Vz','Vel','Pressure','Temperature','BasalforcingsMeltingRate'};
-field_tolerances={1e-09,1e-09,1e-08,3e-10,1e-13,2e-10,1e-06};
+field_tolerances={1e-09,1e-09,1e-08,3e-10,1e-13,3e-10,1e-06};
 field_values={...
 	(md.results.SteadystateSolution.Vx),...
 	(md.results.SteadystateSolution.Vy),...
Index: ../trunk-jpl/test/NightlyRun/test435.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test435.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test435.m	(revision 17747)
@@ -47,7 +47,7 @@
 	'Bed2','Surface2','Thickness2','Floatingice2','Vx2','Vy2','Vz2','Pressure2',...
 	'Bed3','Surface3','Thickness3','Floatingice3','Vx3','Vy3','Vz3','Pressure3'};
 field_tolerances={2e-11,3e-12,2e-11,1e-11,4e-10,1e-08,6e-10,1e-13,...
-	2e-11,3e-11,9e-10,7e-11,1e-09,1e-08,1e-09,1e-10,...
+	3e-11,3e-11,9e-10,7e-11,1e-09,1e-08,1e-09,1e-10,...
 	1e-10,3e-11,1e-10,7e-11,1e-09,1e-08,1e-09,1e-10};
 field_values={...
 	(md.results.TransientSolution(1).Base),...
Index: ../trunk-jpl/test/NightlyRun/test206.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test206.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test206.m	(revision 17747)
@@ -9,7 +9,7 @@
 
 %Fields and tolerances to track changes
 field_names     ={'Temperature','BasalforcingsMeltingRate'};
-field_tolerances={1e-13,2.6e-6};
+field_tolerances={1e-13,5e-6};
 field_values={...
 	(md.results.ThermalSolution.Temperature),...
 	(md.results.ThermalSolution.BasalforcingsMeltingRate),...
Index: ../trunk-jpl/test/NightlyRun/test322.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test322.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test322.m	(revision 17747)
@@ -22,7 +22,7 @@
 
 %Fields and tolerances to track changes
 field_names     ={'Gradient','Misfits','FrictionCoefficient','Pressure','Vel','Vx','Vy'};
-field_tolerances={1e-06,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,2e-09,2e-09,2e-09};
+field_tolerances={2e-06,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,2e-09,2e-09,2e-09};
 field_values={...
 	(md.results.StressbalanceSolution.Gradient1),...
 	(md.results.StressbalanceSolution.J),...
Index: ../trunk-jpl/test/NightlyRun/test221.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test221.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test221.py	(revision 17747)
@@ -18,7 +18,7 @@
 
 #Fields and tolerances to track changes
 field_names     =['Vx','Vy','Vz','Vel','Pressure']
-field_tolerances=[1e-09,1e-09,1e-06,1e-09,1e-09]
+field_tolerances=[1e-09,1e-09,5e-06,1e-09,1e-09]
 field_values=[\
 	md.results.StressbalanceSolution.Vx,\
 	md.results.StressbalanceSolution.Vy,\
Index: ../trunk-jpl/test/NightlyRun/test209.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test209.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test209.py	(revision 17747)
@@ -25,8 +25,8 @@
 	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3']
 field_tolerances=[\
 	1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,\
-	1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,1.9e-8,\
-	5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,1.85e-8]
+	1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,2e-8,\
+	5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,2e-8]
 
 field_values=[\
 	md.results.TransientSolution[0].Vx,\
Index: ../trunk-jpl/test/NightlyRun/test233.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test233.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test233.m	(revision 17747)
@@ -15,8 +15,8 @@
 	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3', ...
 	'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsMeltingRate4'};
 field_tolerances={1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, ...
-	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-07, ...
 	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, ...
+	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, ...
 	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06};
 field_values={...
 	(md.results.TransientSolution(1).Vx),...
Index: ../trunk-jpl/test/NightlyRun/test455.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test455.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test455.py	(revision 17747)
@@ -22,7 +22,7 @@
 	md.flowequation.fe_HO=i
 	md=solve(md,StressbalanceSolutionEnum())
 	field_names     =field_names+['Vx'+i,'Vy'+i,'Vz'+i,'Vel'+i,'Pressure'+i]
-	field_tolerances=field_tolerances+[4e-08,4e-08,1e-08,4e-08,3e-13]
+	field_tolerances=field_tolerances+[5e-08,4e-08,1e-08,4e-08,3e-13]
 	field_values=field_values+[\
 			md.results.StressbalanceSolution.Vx,\
 			md.results.StressbalanceSolution.Vy,\
Index: ../trunk-jpl/test/NightlyRun/test455.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test455.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test455.m	(revision 17747)
@@ -12,7 +12,7 @@
 	md.flowequation.fe_HO=i{1};
 	md=solve(md,StressbalanceSolutionEnum());
 	field_names     ={field_names{:},['Vx' i{1}],['Vy' i{1}],['Vz' i{1}],['Vel' i{1}],['Pressure' i{1}]};
-	field_tolerances={field_tolerances{:},4e-08,4e-08,1e-08,4e-08,3e-13};
+	field_tolerances={field_tolerances{:},5e-08,4e-08,1e-08,4e-08,3e-13};
 	field_values={field_values{:},...
 	(md.results.StressbalanceSolution.Vx),...
 	(md.results.StressbalanceSolution.Vy),...
Index: ../trunk-jpl/test/NightlyRun/test410.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test410.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test410.py	(revision 17747)
@@ -19,7 +19,7 @@
 
 #Fields and tolerances to track changes
 field_names     =['Vx','Vy','Vz','Vel','Pressure','Temperature','BasalforcingsMeltingRate']
-field_tolerances=[1e-09,1e-09,1e-08,1e-09,1e-13,2e-10,1e-06]
+field_tolerances=[1e-09,1e-09,1e-08,1e-09,1e-13,3e-10,1e-06]
 field_values=[\
 	md.results.SteadystateSolution.Vx,\
 	md.results.SteadystateSolution.Vy,\
Index: ../trunk-jpl/test/NightlyRun/test209.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test209.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test209.m	(revision 17747)
@@ -12,8 +12,8 @@
 	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3'};
 field_tolerances={...
 	1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-8,...
-	1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,1.9e-8,...
-	5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,1.85e-8};
+	1e-13,1e-13,3e-9,1e-13,3.8e-11,3.85e-11,3.85e-11,3.85e-11,1e-13,2e-8,...
+	5e-11,6e-12,1e-08,6e-12,3.9e-11,1e-10,1e-10,1e-10,9e-12,2-8};
 field_values={...
 	(md.results.TransientSolution(1).Vx),...
 	(md.results.TransientSolution(1).Vy),...
Index: ../trunk-jpl/test/NightlyRun/test285.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test285.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test285.m	(revision 17747)
@@ -12,7 +12,7 @@
 	md.flowequation.fe_HO=i{1};
 	md=solve(md,StressbalanceSolutionEnum());
 	field_names     ={field_names{:},['Vx' i{1}],['Vy' i{1}],['Vz' i{1}],['Vel' i{1}],['Pressure' i{1}]};
-	field_tolerances={field_tolerances{:},6e-08,2e-08,2e-08,2e-08,1e-13};
+	field_tolerances={field_tolerances{:},6e-08,3e-08,2e-08,3e-08,1e-13};
 	field_values={field_values{:},...
 	(md.results.StressbalanceSolution.Vx),...
 	(md.results.StressbalanceSolution.Vy),...
Index: ../trunk-jpl/test/NightlyRun/test285.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test285.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test285.py	(revision 17747)
@@ -22,7 +22,7 @@
 	md.flowequation.fe_HO=i
 	md=solve(md,StressbalanceSolutionEnum())
 	field_names     =field_names+['Vx'+i,'Vy'+i,'Vz'+i,'Vel'+i,'Pressure'+i]
-	field_tolerances=field_tolerances+[6e-08,2e-08,2e-08,2e-08,1e-13]
+	field_tolerances=field_tolerances+[6e-08,3e-08,2e-08,3e-08,1e-13]
 	field_values=field_values+[\
 		md.results.StressbalanceSolution.Vx,\
 		md.results.StressbalanceSolution.Vy,\
Index: ../trunk-jpl/test/NightlyRun/test206.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test206.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test206.py	(revision 17747)
@@ -22,7 +22,7 @@
 # Fields and tolerances to track changes
 
 field_names     =['Temperature','BasalforcingsMeltingRate']
-field_tolerances=[1e-13,1.35e-6]
+field_tolerances=[1e-13,5e-6]
 field_values=[\
 	md.results.ThermalSolution.Temperature,\
 	md.results.ThermalSolution.BasalforcingsMeltingRate,\
Index: ../trunk-jpl/test/NightlyRun/test221.m
===================================================================
--- ../trunk-jpl/test/NightlyRun/test221.m	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test221.m	(revision 17747)
@@ -9,7 +9,7 @@
 
 %Fields and tolerances to track changes
 field_names     ={'Vx','Vy','Vz','Vel','Pressure'};
-field_tolerances={1e-09,1e-09,1e-06,1e-09,1e-09};
+field_tolerances={1e-09,1e-09,5e-06,1e-09,1e-09};
 field_values={...
 	(md.results.StressbalanceSolution.Vx),...
 	(md.results.StressbalanceSolution.Vy),...
Index: ../trunk-jpl/test/NightlyRun/test233.py
===================================================================
--- ../trunk-jpl/test/NightlyRun/test233.py	(revision 17746)
+++ ../trunk-jpl/test/NightlyRun/test233.py	(revision 17747)
@@ -25,8 +25,8 @@
 	'Vx3','Vy3','Vz3','Vel3','Pressure3','Bed3','Surface3','Thickness3','Temperature3','BasalforcingsMeltingRate3', \
 	'Vx4','Vy4','Vz4','Vel4','Pressure4','Bed4','Surface4','Thickness4','Temperature4','BasalforcingsMeltingRate4']
 field_tolerances=[1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09, \
-	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-07, \
 	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, \
+	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06, \
 	1e-09,1e-09,1e-08,1e-09,1e-09,1e-09,1e-09,1e-09,1e-09,1e-06]
 field_values=[\
 	md.results.TransientSolution[0].Vx,\