diff --git a/bfps/NavierStokes.py b/bfps/NavierStokes.py
index 229f1480afd05068ff293d364c553cd14e8f6175..ffdbd43751633106b4daa533499ca7ab9a80c4a7 100644
--- a/bfps/NavierStokes.py
+++ b/bfps/NavierStokes.py
@@ -156,7 +156,7 @@ class NavierStokes(bfps.fluid_base.fluid_particle_base):
                         hist_gradu,
                         hist_trS2_Q_R,
                         hist_QR2D,
-                        max_estimates,
+                        trS2QR_max_estimates,
                         gradu_max_estimates,
                         histogram_bins,
                         QR2D_histogram_bins);
diff --git a/tests/test_plain.py b/tests/test_plain.py
index 9e25cf029f294f804656f18afe88e25594c6667b..7437f1f61c9694c0e8828b4946c54e840a22461e 100755
--- a/tests/test_plain.py
+++ b/tests/test_plain.py
@@ -42,7 +42,7 @@ def plain(opt):
     print ('Rlambda = {0:.0f}'.format(c0.statistics['Rlambda']))
     print ('Lint = {0:.4e}, etaK = {1:.4e}'.format(c0.statistics['Lint'], c0.statistics['etaK']))
     print ('Tint = {0:.4e}, tauK = {1:.4e}'.format(c0.statistics['Tint'], c0.statistics['tauK']))
-    print ('kMetaK = {0:.4e}'.format(c0.statistics['kM']*c0.statistics['etaK']))
+    print ('kMetaK = {0:.4e}'.format(c0.statistics['kMeta']))
     for s in range(c0.particle_species):
         acceleration_test(c0, species = s, m = 1)
     if not opt.multiplejob: