force divfree on test vorticity

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')

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);