diff --git a/bfps/cpp/full_code/NSVE.cpp b/bfps/cpp/full_code/NSVE.cpp
index ecec7db31235bc827b17b55f0c733f305e488761..da26ad22b0afba8d60ed2a3b45568e3eec723ced 100644
--- a/bfps/cpp/full_code/NSVE.cpp
+++ b/bfps/cpp/full_code/NSVE.cpp
@@ -156,6 +156,7 @@ int NSVE<rnumber>::read_parameters(void)
     this->fmode = hdf5_tools::read_value<int>(parameter_file, "parameters/fmode");
     this->famplitude = hdf5_tools::read_value<double>(parameter_file, "parameters/famplitude");
     this->friction_coefficient = hdf5_tools::read_value<double>(parameter_file, "parameters/friction_coefficient");
+    this->injection_rate = hdf5_tools::read_value<double>(parameter_file, "parameters/injection_rate");
     this->fk0 = hdf5_tools::read_value<double>(parameter_file, "parameters/fk0");
     this->fk1 = hdf5_tools::read_value<double>(parameter_file, "parameters/fk1");
     this->energy = hdf5_tools::read_value<double>(parameter_file, "parameters/energy");
diff --git a/bfps/cpp/full_code/NSVEcomplex_particles.cpp b/bfps/cpp/full_code/NSVEcomplex_particles.cpp
index 02a199317ae6bdf294adae1b8805b89df7f276b9..487fe4bb436268d9a2a305dab0c077a572b84cce 100644
--- a/bfps/cpp/full_code/NSVEcomplex_particles.cpp
+++ b/bfps/cpp/full_code/NSVEcomplex_particles.cpp
@@ -24,7 +24,7 @@
 
 
 
-#define NDEBUG
+//#define NDEBUG
 
 #include <string>
 #include <cmath>
diff --git a/bfps/cpp/full_code/code_base.cpp b/bfps/cpp/full_code/code_base.cpp
index a6487c726de44b018392128f955ccebf7e7100a1..210f2c963089751920788f8f2004c0331a565f6c 100644
--- a/bfps/cpp/full_code/code_base.cpp
+++ b/bfps/cpp/full_code/code_base.cpp
@@ -24,7 +24,7 @@
 
 
 
-#define NDEBUG
+//#define NDEBUG
 
 #include "code_base.hpp"
 #include "scope_timer.hpp"
diff --git a/bfps/cpp/particles/particles_distr_mpi.hpp b/bfps/cpp/particles/particles_distr_mpi.hpp
index 251119be950ad8219ee611d858ca8ac6c9986dbe..dbf15b6ba3ccb5f22ca9d297e549ee2135336253 100644
--- a/bfps/cpp/particles/particles_distr_mpi.hpp
+++ b/bfps/cpp/particles/particles_distr_mpi.hpp
@@ -141,7 +141,7 @@ public:
                        const real_number particles_positions[],
                        real_number particles_current_rhs[],
                        const int interpolation_size){
-        TIMEZONE("compute_distr");
+        TIMEZONE("particle_distr_mpi::compute_distr");
 
         // Some processes might not be involved
         if(nb_processes_involved <= my_rank){
diff --git a/bfps/cpp/vorticity_equation.cpp b/bfps/cpp/vorticity_equation.cpp
index ead9345af5a2f0555e7fa6e2b6ee45cecd9f3624..320c51d80d9bde3f45259f62880dbaa1c4ebeaee 100644
--- a/bfps/cpp/vorticity_equation.cpp
+++ b/bfps/cpp/vorticity_equation.cpp
@@ -24,7 +24,7 @@
 
 
 
-#define NDEBUG
+//#define NDEBUG
 
 #include <limits>
 #include <cassert>