diff --git a/diffusivityModel/diffusivityModel/diffusivityModel.C b/diffusivityModel/diffusivityModel/diffusivityModel.C
index 83135e6..9105991 100644
--- a/diffusivityModel/diffusivityModel/diffusivityModel.C
+++ b/diffusivityModel/diffusivityModel/diffusivityModel.C
@@ -74,9 +74,9 @@ Foam::diffusivityModel::mixAvgDi(UList<scalar>& Di, const scalarSymmetricSquareM
     // Pure Condition Test ( Xi = 1 )
     Switch pure = false;
     label pureSpecieI = -1;
-    forAll (X, i)
+    forAll (Y, i)
     {
-        if (1. - X[i] < 1e-12 || std::isinf(X[i]/(1.-Y[i])))
+        if (1. - Y[i] < 1e-12)
             
         {
             pureSpecieI = i;