diff --git a/CMakeLists.txt b/CMakeLists.txt
index a9e435901a8e14b4c706e2005886b19c1461cdf0..7b7cd609f65e6a0afecfbb7083ad9da8edaa8135 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -214,7 +214,8 @@ set(cpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/particles/particles_inner_computer.cpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_process.cpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_test.cpp
-    ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.cpp)
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.cpp
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/NSVE_ou_forcing.cpp)
 
 set(hpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/full_code/code_base.hpp
@@ -292,6 +293,7 @@ set(hpp_for_lib
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_process.hpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ornstein_uhlenbeck_test.hpp
     ${PROJECT_SOURCE_DIR}/cpp/full_code/ou_vorticity_equation.hpp
+    ${PROJECT_SOURCE_DIR}/cpp/full_code/NSVE_ou_forcing.hpp
     )
 #file(GLOB_RECURSE hpp_for_lib ${PROJECT_SOURCE_DIR}/*.hpp)
 LIST(APPEND source_files ${hpp_for_lib} ${cpp_for_lib})
diff --git a/TurTLE/DNS.py b/TurTLE/DNS.py
index 411f2cf7fb3dc25c9436c360e74200d9841670ba..88dc3854f0985a96103ef66085952cedb310e716 100644
--- a/TurTLE/DNS.py
+++ b/TurTLE/DNS.py
@@ -655,6 +655,10 @@ class DNS(_code):
                 'NSVEp_extra_sampling',
                 help = 'plain Navier-Stokes vorticity formulation, with basic fluid tracers, that sample velocity gradient, as well as pressure and its derivatives.')
 
+        parser_NSVE_ou = subparsers.add_parser(
+                'NSVE_ou_forcing',
+                help = 'plain Navier-Stokes vorticity formulation, with ornstein-uhlenbeck forcing')
+
         for parser in ['NSVEparticles_no_output', 'NSVEp2', 'NSVEp2p', 'NSVEp_extra']:
             eval('self.simulation_parser_arguments({0})'.format('parser_' + parser))
             eval('self.job_parser_arguments({0})'.format('parser_' + parser))
@@ -1057,4 +1061,3 @@ class DNS(_code):
                 no_submit = opt.no_submit,
                 no_debug = opt.no_debug)
         return None
-
diff --git a/TurTLE/ou_quick_plot.py b/TurTLE/ou_quick_plot.py
index b42c2cd57a4312a0ad8928f752142454408cbd50..fcf5766e9d7774029caff99768d7026e3160f21a 100644
--- a/TurTLE/ou_quick_plot.py
+++ b/TurTLE/ou_quick_plot.py
@@ -56,9 +56,9 @@ class plot_field:
         
 
 
-# plot_spectrum()
-# plot_histogram()
+plot_spectrum()
+plot_histogram()
 p = plot_field()
-p.plot('ou_field.h5')
-p.plot('add_ou_test_field.h5')
+#p.plot('ou_field.h5')
+#p.plot('add_ou_test_field.h5')
 plt.show()
diff --git a/cpp/vorticity_equation.hpp b/cpp/vorticity_equation.hpp
index cbff223e2d8119b37ef30b4e6b739aa64eff7039..50db48f0cc68e7a7cea2dcfa939ef1d2bf137e3f 100644
--- a/cpp/vorticity_equation.hpp
+++ b/cpp/vorticity_equation.hpp
@@ -87,7 +87,7 @@ class vorticity_equation
         ~vorticity_equation(void);
 
         /* solver essential methods */
-        void omega_nonlin(int src);
+        virtual void omega_nonlin(int src);
         void step(double dt);
         void impose_zero_modes(void);
 
@@ -166,4 +166,3 @@ class vorticity_equation
 };
 
 #endif//VORTICITY_EQUATION
-