bugfix: NSVE_Stokes_particles uses correct Loop in step()

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__':

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