diff --git a/dependencies/reactor_class.js b/dependencies/reactor_class.js
index 76a6340..c1e8989 100644
--- a/dependencies/reactor_class.js
+++ b/dependencies/reactor_class.js
@@ -111,10 +111,10 @@ class Reactor_PFR {
         let index_in = input.payload.inlet;
         this.Fs[index_in] = input.payload.F;
         this.Cs_in[index_in] = input.payload.C;
-        console.log("Pe total " + this.length*math.sum(this.Fs)/(this.D*this.A));
-        console.log("Pe local " + this.d_x*math.sum(this.Fs)/(this.D*this.A));
-        console.log("Co ad " + math.sum(this.Fs)*this.timeStep/(this.A*this.d_x));
-        console.log("Co D " + this.D*this.timeStep/(this.d_x*this.d_x));
+        // console.log("Pe total " + this.length*math.sum(this.Fs)/(this.D*this.A));
+        // console.log("Pe local " + this.d_x*math.sum(this.Fs)/(this.D*this.A));
+        // console.log("Co ad " + math.sum(this.Fs)*this.timeStep/(this.A*this.d_x));
+        // console.log("Co D " + this.D*this.timeStep/(this.d_x*this.d_x));
     }
 
     set setOTR(input) { // setter for OTR (WIP) [g O2 d-1]
@@ -185,7 +185,6 @@ class Reactor_PFR {
             BC_gradient[1] = 1;
             let Pe = this.length*math.sum(this.Fs)/(this.D*this.A)
             const BC_dispersion = math.multiply((1-(1+4*this.volume/math.sum(this.Fs)/Pe)^0.5)/Pe, [BC_gradient], new_state)[0];
-            console.log(math.add(BC_C_in, BC_dispersion));
             new_state[0] = math.add(BC_C_in, BC_dispersion).map(val => val < 0 ? 0 : val);
             
         } else { // Neumann BC (no flux)