diff --git a/models/IonModel.cpp b/models/IonModel.cpp
index 5cf38e64..bc58a649 100644
--- a/models/IonModel.cpp
+++ b/models/IonModel.cpp
@@ -693,15 +693,9 @@ void ScaLBL_IonModel::Run(double *Velocity, double *ElectricField){
     //2. ElectricField is from Poisson model
 
     //LB-related parameter
-    vector<double> rlx(tau.begin(),tau.end());
-    for (double item : rlx){
-        item = 1.0/item; 
-    }
-    //**debug
-    if (rank==0){
-        for (unsigned int ic=0;ic<rlx.size();ic++){
-            printf("**debug: Ion %i rlx = %.5g; LB diffusivity = %.5g \n",ic+1,rlx[ic],IonDiffusivity[ic]);
-        }
+    vector<double> rlx;
+    for (unsigned int ic=0;ic<tau.size();ic++){
+        rlx.push_back(1.0/tau[ic]); 
     }
     
 	//.......create and start timer............