diff --git a/diffusivityModel/diffusivityModel/diffusivityModel.C b/diffusivityModel/diffusivityModel/diffusivityModel.C
index 94520ee..6f7ea52 100644
--- a/diffusivityModel/diffusivityModel/diffusivityModel.C
+++ b/diffusivityModel/diffusivityModel/diffusivityModel.C
@@ -173,7 +173,6 @@ void Foam::diffusivityModel::calculateMuD
  const scalar Cpi,
  const scalar pi,
  const scalar Ti,
- const scalar WbarI,
  const scalar rhoQc2i,
  const scalarField &localY,
  const scalarField &localX,
@@ -664,7 +663,7 @@ void Foam::diffusivityModel::correct()
             localX[i] = Y[i][celli] * WbarI / Wpure[i];
         }
 
-        calculateMuD ( muI, Dij, rhoi, Cpi, pi, Ti, WbarI, rhoQc2i, localY, localX, electron_.valid());
+        calculateMuD ( muI, Dij, rhoi, Cpi, pi, Ti, rhoQc2i, localY, localX, electron_.valid());
 
         forAll (Dii, i)
         {
@@ -751,7 +750,7 @@ void Foam::diffusivityModel::correct()
                 localX[i] = localY[i] * WbarI / Wpure[i];
             }
 
-            calculateMuD ( muI, Dij, rhoi, Cpi, pi, Ti, WbarI, rhoQc2i, localY, localX, electron_.valid());
+            calculateMuD ( muI, Dij, rhoi, Cpi, pi, Ti, rhoQc2i, localY, localX, electron_.valid());
 
             forAll (Dii, i)
             {
diff --git a/diffusivityModel/diffusivityModel/diffusivityModel.H b/diffusivityModel/diffusivityModel/diffusivityModel.H
index e392a8c..1a44424 100644
--- a/diffusivityModel/diffusivityModel/diffusivityModel.H
+++ b/diffusivityModel/diffusivityModel/diffusivityModel.H
@@ -131,7 +131,6 @@ class diffusivityModel
             const scalar Cpi,
             const scalar pi,
             const scalar Ti,
-            const scalar WbarI,
             const scalar rhoQc2i,
             const scalarField &localY,
             const scalarField &localX,