diff --git a/cpp/full_code/ornstein_uhlenbeck_process.cpp b/cpp/full_code/ornstein_uhlenbeck_process.cpp
index f3955eccd9b59067d6fdab8d08cc631313e92ff3..a73e48e465ca221855e7459f72b34c661df29ba1 100644
--- a/cpp/full_code/ornstein_uhlenbeck_process.cpp
+++ b/cpp/full_code/ornstein_uhlenbeck_process.cpp
@@ -105,9 +105,9 @@ void ornstein_uhlenbeck_process<rnumber,be>::step(
         }
     });
 
+    this->ou_field->symmetrize();
     this->calc_ou_vorticity();
 
-    // this->ou_field->symmetrize(); ???
 }
 
 // template <class rnumber,field_backend be>
diff --git a/cpp/full_code/ou_vorticity_equation.cpp b/cpp/full_code/ou_vorticity_equation.cpp
index 2b341ca7951b103c92bc974eaacc2f86b5d44f4e..f3bfe8c45db11c42eab7faf8ad6361803a1b2f39 100644
--- a/cpp/full_code/ou_vorticity_equation.cpp
+++ b/cpp/full_code/ou_vorticity_equation.cpp
@@ -38,7 +38,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[1], this->v[0]);
+    //this->impose_forcing(this->v[1], this->v[0]);
 
     this->omega_nonlin(1);
     this->kk->CLOOP_K2(
@@ -60,7 +60,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[2], this->v[0]);
+    //this->impose_forcing(this->v[2], this->v[0]);
 
     this->omega_nonlin(2);
     // store old vorticity
@@ -83,7 +83,7 @@ void ou_vorticity_equation<rnumber, be>::step(double dt)
         }
     }
     );
-    this->impose_forcing(this->v[0], this->v[1]);
+    //this->impose_forcing(this->v[0], this->v[1]);
 
     this->kk->template force_divfree<rnumber>(this->cvorticity->get_cdata());
     this->cvorticity->symmetrize();