diff --git a/TurTLE/test/test_turtle_NSVE_Stokes_particles.py b/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
index 825ad1c8890aa699611c72638bfdafb11e549d42..9428b3f8dbfe74647a1b7399ec70e6e5ff94bf41 100644
--- a/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
+++ b/TurTLE/test/test_turtle_NSVE_Stokes_particles.py
@@ -63,7 +63,8 @@ def main():
              sys.argv[1:])
     f = h5py.File('quiescent_nsve_stokes_particles_particles.h5', 'r')
     print(f['tracers0/momentum/0'][:])
-    print(f['tracers0/momentum/100'][:])
+    print(f['tracers0/momentum/0'][:]-f['tracers0/momentum/100'][:])
+    print(f['tracers0/position/0'][:]-f['tracers0/position/100'][:])
     return None
 
 if __name__ == '__main__':
diff --git a/cpp/full_code/NSVE_Stokes_particles.cpp b/cpp/full_code/NSVE_Stokes_particles.cpp
index eb9b4084eb04583a6b6a28bb50432c93f0e213bc..4c90dbc226e2bbc7a92ff43de75d6eea1f5ca36b 100644
--- a/cpp/full_code/NSVE_Stokes_particles.cpp
+++ b/cpp/full_code/NSVE_Stokes_particles.cpp
@@ -88,7 +88,7 @@ int NSVE_Stokes_particles<rnumber>::step(void)
     TIMEZONE("NSVE_Stokes_particles::step");
     this->fs->compute_velocity(this->fs->cvorticity);
     this->fs->cvelocity->ift();
-    this->ps->completeLoop(this->dt);
+    this->ps->complete2ndOrderLoop(this->dt);
     this->NSVE<rnumber>::step();
     return EXIT_SUCCESS;
 }