Index: /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/Square.par
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/Square.par	(revision 1159)
+++ /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/Square.par	(revision 1160)
@@ -6,5 +6,6 @@
 md.ndt=md.dt*3; 
 md.artificial_diffusivity=1;
-md.eps_abs=5000;
+md.eps_abs=NaN;
+md.eps_rel=0.1;
 
 hmin=300;
@@ -26,4 +27,5 @@
 
 md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+md.min_thermal_constraints=3;
 
 md.B=paterson(md.observed_temperature);
Index: /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m	(revision 1159)
+++ /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m	(revision 1160)
@@ -66,6 +66,6 @@
 		md=extrude(md,8,4);
 		md=setelementstype(md,'macayeal','all','stokes','all');
-		if md.numberofgrids==832
-			load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+		if md.numberofgrids==728
+			load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature; md.pressure=pressure;
 		end
 
Index: /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/updatearchive.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/updatearchive.m	(revision 1159)
+++ /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/updatearchive.m	(revision 1160)
@@ -60,5 +60,5 @@
 		md=extrude(md,8,4);
 		md=setelementstype(md,'macayeal','all','stokes','all');
-		if md.numberofgrids==832
+		if md.numberofgrids==728
 			load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
 		end