diff --git a/src/fluid_solver.cpp b/src/fluid_solver.cpp
index ab91edafa27b45b788f649d4b1b15c94054d24e8..c6fde9c58f89bf077921845f45e1342907279e05 100644
--- a/src/fluid_solver.cpp
+++ b/src/fluid_solver.cpp
@@ -171,7 +171,7 @@ fluid_solver<R>::fluid_solver( \
  \
     /* ``physical'' parameters etc */ \
  \
-    this->nu = 0.01; \
+    this->nu = 0.1; \
     this->dkx = DKX; \
     this->dky = DKY; \
     this->dkz = DKZ; \
diff --git a/test.py b/test.py
index 0f19c7276c501ddfb180ba13ddee4614a5e2068d..b2bdd40f29b451ffb6b72a8dbb24e3e6f91c8056 100755
--- a/test.py
+++ b/test.py
@@ -61,10 +61,11 @@ Kdata0[..., 1] = Kdata01
 Kdata0[..., 2] = Kdata02
 Kdata0.tofile("Kdata0")
 run_test('test_step')
-#Kdata1 = np.fromfile('Kdata1', dtype = np.complex64).reshape(Kdata0.shape)
-#
+Kdata1 = np.fromfile('Kdata1', dtype = np.complex64).reshape(Kdata0.shape)
+
 #print np.max(np.abs(Kdata0 - Kdata1))
-#print np.max(np.abs(Kdata0))
+print np.max(np.abs(Kdata0))
+print np.max(np.abs(Kdata1))
 #
 #fig = plt.figure(figsize=(12, 6))
 #a = fig.add_subplot(121)