diff --git a/Source/mass.f90 b/Source/mass.f90
index 0a04960569..d4e5d9fbe4 100644
--- a/Source/mass.f90
+++ b/Source/mass.f90
@@ -806,7 +806,7 @@ SUBROUTINE CHECK_MASS_DENSITY
          IF (CELL(IC)%WALL_INDEX(-3)==0) MASS_N(-3) = ABS(MIN(RHOMAX,MAX(RHOMIN,RHOP(I,J,K-1)))-RHO_CUT)*VC(-3)
          IF (CELL(IC)%WALL_INDEX( 3)==0) MASS_N( 3) = ABS(MIN(RHOMAX,MAX(RHOMIN,RHOP(I,J,K+1)))-RHO_CUT)*VC( 3)
          SUM_MASS_N = SUM(MASS_N)
-         IF (SUM_MASS_N<=TWENTY_EPSILON_EB) CYCLE
+         IF (SUM_MASS_N<=TWO_EPSILON_EB) CYCLE
          CONST = SIGN_FACTOR*MIN(1._EB,MASS_C/SUM_MASS_N)
          DELTA_RHO(I,J,K)   = DELTA_RHO(I,J,K)   + CONST*SUM_MASS_N/VC( 0)
          DELTA_RHO(I-1,J,K) = DELTA_RHO(I-1,J,K) - CONST*MASS_N(-1)/VC(-1)
@@ -882,7 +882,7 @@ SUBROUTINE CHECK_MASS_DENSITY
             IF (CELL(IC)%WALL_INDEX(-3)==0) MASS_N(-3) = ABS(MIN(RHO_ZZ_MAX,MAX(RHO_ZZ_MIN,RHO_ZZ(I,J,K-1,N)))-RHO_ZZ_CUT)*VC(-3)
             IF (CELL(IC)%WALL_INDEX( 3)==0) MASS_N( 3) = ABS(MIN(RHO_ZZ_MAX,MAX(RHO_ZZ_MIN,RHO_ZZ(I,J,K+1,N)))-RHO_ZZ_CUT)*VC( 3)
             SUM_MASS_N = SUM(MASS_N)
-            IF (SUM_MASS_N<=TWENTY_EPSILON_EB) CYCLE
+            IF (SUM_MASS_N<=TWO_EPSILON_EB) CYCLE
             CONST = SIGN_FACTOR*MIN(1._EB,MASS_C/SUM_MASS_N)
             DELTA_RHO_ZZ(I,J,K)   = DELTA_RHO_ZZ(I,J,K)   + CONST*SUM_MASS_N/VC( 0)
             DELTA_RHO_ZZ(I-1,J,K) = DELTA_RHO_ZZ(I-1,J,K) - CONST*MASS_N(-1)/VC(-1)