diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C b/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C
index 0dca7e3b..98b3f05e 100644
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/MeiLift/MeiLift.C
@@ -190,12 +190,8 @@ void MeiLift::setForce() const
                     else
                     {
                         J_star = 0.3
-                                *(     1.0
-                                      +tanh(  2.5 * log10(epsilon+0.191)  )
-                                 )
-                                *(    0.667
-                                     +tanh(  6.0 * (epsilon-0.32)  )
-                                  );
+                                 * (1.0   + tanh(2.5 * (log10(epsilon) + 0.191)))
+                                 * (0.667 + tanh(6.0 * (      epsilon  - 0.32 )));
                     }
                     Cl = J_star * 4.11 * epsilon; //multiply McLaughlin's correction to the basic Saffman model