Index: /issm/trunk/test/Verification/IceSheetIceFrontM2d_7/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontM2d_7/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetIceFrontM2d_7/runme.m	(revision 707)
@@ -63,5 +63,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetIceFrontM3d_8/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontM3d_8/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetIceFrontM3d_8/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetIceFrontMP3d_9/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontMP3d_9/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetIceFrontMP3d_9/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetIceFrontP3d_10/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontP3d_10/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetIceFrontP3d_10/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetIceFrontS3d_11/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontH2d_12/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontH2d_12/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontH2d_12/runme.m	(revision 707)
@@ -62,5 +62,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontH3d_13/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontH3d_13/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontH3d_13/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontHM2d_14/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontHM2d_14/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontHM2d_14/runme.m	(revision 707)
@@ -62,5 +62,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontHM3d_15/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontHM3d_15/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontHM3d_15/runme.m	(revision 707)
@@ -63,5 +63,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontM2d_16/runme.m	(revision 707)
@@ -63,5 +63,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontM3d_17/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontM3d_17/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontM3d_17/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontP3d_18/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontP3d_18/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontP3d_18/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceSheetNoIceFrontS3d_19/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceSheetNoIceFrontS3d_19/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceSheetNoIceFrontS3d_19/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontM2d_1/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontM2d_1/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontM2d_1/runme.m	(revision 707)
@@ -63,5 +63,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontM3d_2/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontM3d_2/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontM3d_2/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontMP3d_3/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontMP3d_3/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontMP3d_3/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontMPS3d_4/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontMPS3d_4/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontMPS3d_4/runme.m	(revision 707)
@@ -63,5 +63,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontP3d_5/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontP3d_5/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontP3d_5/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/IceShelfIceFrontS3d_6/runme.m
===================================================================
--- /issm/trunk/test/Verification/IceShelfIceFrontS3d_6/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/IceShelfIceFrontS3d_6/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package solution '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/PigControlMethodBM2d_21/runme.m	(revision 707)
@@ -67,5 +67,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package misfit '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result
Index: /issm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m
===================================================================
--- /issm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m	(revision 706)
+++ /issm/trunk/test/Verification/PigControlMethodDragM2d_20/runme.m	(revision 707)
@@ -65,5 +65,5 @@
 			%compare to archive
 			eval(['Archive=Archive' package misfit '_field' num2str(k) ';']);
-			eval(['error_diff=abs(norm((Archive(find(Archive>tolerance))-md.results.' field '(find(Archive>tolerance)))./Archive(find(Archive>tolerance)),2));']);
+			eval(['error_diff=max(abs(Archive-md.results.' field '))/max(abs(Archive));']);
 
 			%disp test result