diff --git a/bfps/TEST.py b/bfps/TEST.py
index 96090076d8770bff90d382785684d5e535745363..d9b3deddedb495e0461d4ebf92c209e3df95b4e7 100644
--- a/bfps/TEST.py
+++ b/bfps/TEST.py
@@ -349,18 +349,18 @@ class TEST(_code):
                            write_to_file = False,
                            spectra_slope = 1.0,
                            amplitude = 0.05)
-                    data[:] = 0.0
-                    # ABC
-                    data[0, 0, 1, 1] = -0.5*(1j)
-                    data[0, 0, 1, 2] =  0.5*(1j)
-                    data[0, 1, 0, 0]                         = -0.5*(1j)
-                    data[0, self.parameters['nz'] - 1, 0, 0] =  0.5*(1j)
-                    data[0, 1, 0, 1]                         =  0.5
-                    data[0, self.parameters['nz'] - 1, 0, 1] =  0.5
-                    data[1, 0, 0, 0]                         =  0.5
-                    data[self.parameters['ny'] - 1, 0, 0, 0] =  0.5
-                    data[1, 0, 0, 2]                         = -0.5*(1j)
-                    data[self.parameters['ny'] - 1, 0, 0, 2] =  0.5*(1j)
+                    #data[:] = 0.0
+                    ## ABC
+                    #data[0, 0, 1, 1] = -0.5*(1j)
+                    #data[0, 0, 1, 2] =  0.5*(1j)
+                    #data[0, 1, 0, 0]                         = -0.5*(1j)
+                    #data[0, self.parameters['nz'] - 1, 0, 0] =  0.5*(1j)
+                    #data[0, 1, 0, 1]                         =  0.5
+                    #data[0, self.parameters['nz'] - 1, 0, 1] =  0.5
+                    #data[1, 0, 0, 0]                         =  0.5
+                    #data[self.parameters['ny'] - 1, 0, 0, 0] =  0.5
+                    #data[1, 0, 0, 2]                         = -0.5*(1j)
+                    #data[self.parameters['ny'] - 1, 0, 0, 2] =  0.5*(1j)
                     ofile['vorticity/complex/{0}'.format(0)] = data
                 with h5py.File(os.path.join(self.work_dir, self.simname + '_output.h5'), 'a') as ofile:
                     ofile.require_group('tracers0')
diff --git a/bfps/cpp/full_code/test_interpolation.cpp b/bfps/cpp/full_code/test_interpolation.cpp
index 4cfdca664f9d7fd03df92ec5e70afd805bfb9622..e702e7089d390168957ac440e8e723b4a1468c7a 100644
--- a/bfps/cpp/full_code/test_interpolation.cpp
+++ b/bfps/cpp/full_code/test_interpolation.cpp
@@ -108,6 +108,7 @@ int test_interpolation<rnumber>::do_work()
             this->simname + std::string("_input.h5"),
             "vorticity",
             0, true);
+    this->kk->template force_divfree<rnumber>(this->vorticity->get_cdata());
 
     // compute velocity
     invert_curl(this->kk, this->vorticity, this->velocity);