diff --git a/src/testutil/compare_hdf5.py b/src/testutil/compare_hdf5.py
index e382bfe..006a151 100644
--- a/src/testutil/compare_hdf5.py
+++ b/src/testutil/compare_hdf5.py
@@ -26,9 +26,9 @@
 if ds1.shape!=ds2.shape:
     print "\tDatasets have different size!"
 else:
-    ds=ds1-ds2
+    ds=(ds1-ds2)*(ds1-ds2)
     rmsd=math.sqrt(ds.mean())
-    print "\tRMSD=%.3e"%(rmsd,)
+    print "\tRMSD=%.3e Mean1=%.3e RMSD/Mean1=%.3e"%(rmsd,ds1.mean(),rmsd/ds1.mean())
 
 #compare Mash
 print "Comparing Mesh/Connections"