diff --git a/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/model.c b/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/model.c
index 5463a3a4..28e2b1fe 100644
--- a/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/model.c
+++ b/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/model.c
@@ -169,7 +169,7 @@ namespace Spatialpy{
             //fflush(stdout);
         }
         //printf("pairwiseForce(id=%i) num_chem_rxns=%i\n",me->id,system->num_chem_rxns);
-        printf("pairwiseForce(id=%i) F=[%e %e %e] Fbp=[%e %e %e]\n",me->id,me->F[0],me->F[1],me->F[2],me->Fbp[0],me->Fbp[1],me->Fbp[2]);
+        //printf("pairwiseForce(id=%i) F=[%e %e %e] Fbp=[%e %e %e]\n",me->id,me->F[0],me->F[1],me->F[2],me->Fbp[0],me->Fbp[1],me->Fbp[2]);
         //fflush(stdout);
 
         // after processing all neighbors
diff --git a/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/simulate.cpp b/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/simulate.cpp
index 353daef2..4362982a 100644
--- a/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/simulate.cpp
+++ b/spatialpy/solvers/c_base/ssa_sdpd-c-simulation-engine/src/simulate.cpp
@@ -156,7 +156,7 @@ namespace Spatialpy{
         // Step 3.3: Apply BVF
         if (me->solidTag == 0 && system->static_domain == 0) {
             computeBoundaryVolumeFraction(me,system);
-            //applyBoundaryVolumeFraction(me, system);
+            applyBoundaryVolumeFraction(me, system);
         }
 
         //  Solve deterministic/stochastic reaction-diffusion system