diff --git a/dumux/discretization/cellcentered/mpfa/fickslaw.hh b/dumux/discretization/cellcentered/mpfa/fickslaw.hh
index dc79901c82530a324bd50298f72a411ca16419fa..be41e0885a68bfb0b2aed9d3e498df2afe76ea52 100644
--- a/dumux/discretization/cellcentered/mpfa/fickslaw.hh
+++ b/dumux/discretization/cellcentered/mpfa/fickslaw.hh
@@ -219,7 +219,7 @@ public:
             return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
 
         // Handle interior boundaries
-        flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, fluxVarsCache, getX, phaseIdx, compIdx);
+        flux += Implementation::computeInteriorBoundaryContribution(fvGeometry, elemVolVars, fluxVarsCache, getX, phaseIdx, compIdx);
 
         // return overall resulting flux
         return useTpfaBoundary ? flux*rho*effFactor : flux*rho*effFactor + fluxVarsCache.componentNeumannFlux(compIdx);
@@ -314,6 +314,7 @@ public:
 
     template<typename GetXFunction>
     static Scalar computeInteriorBoundaryContribution(const FVElementGeometry& fvGeometry,
+                                                      const ElementVolumeVariables& elemVolVars,
                                                       const FluxVariablesCache& fluxVarsCache,
                                                       const GetXFunction& getX,
                                                       unsigned int phaseIdx, unsigned int compIdx)