diff --git a/dumux/porousmediumflow/2pdfm/implicit/localresidual.hh b/dumux/porousmediumflow/2pdfm/implicit/localresidual.hh
index 37ba13300baf2c219c94ef325b0c0b62de2f5d23..80d018919e40d915fd614428dcc6b1e1d0aea653 100644
--- a/dumux/porousmediumflow/2pdfm/implicit/localresidual.hh
+++ b/dumux/porousmediumflow/2pdfm/implicit/localresidual.hh
@@ -168,12 +168,6 @@ public:
         storageMatrix[nPhaseIdx]    = 0.0;
         //        const GlobalPosition &globalPos = geometry.corner(scvIdx);
 
-        Scalar dsm_dsf = volVars.dsm_dsf();
-        if (!this->problem_().useInterfaceCondition())
-        {
-            dsm_dsf = 1.0;
-        }
-
         if (isFracture)
         {
             for (int phaseIdx = 0; phaseIdx<2; phaseIdx++)
@@ -185,7 +179,6 @@ public:
                 storageMatrix[phaseIdx] = volVars.density(phaseIdx)
                                         * volVars.porosity()
                                         * wm
-                                        * dsm_dsf
                                         * volVars.saturationMatrix(phaseIdx);
             }
         }
diff --git a/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh b/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
index d12c7a7a6c4c213f0331917935c7fa941b4a83af..4cc216aea564cbef11ee05e4036727725997302a 100644
--- a/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
+++ b/dumux/porousmediumflow/2pdfm/implicit/volumevariables.hh
@@ -189,13 +189,6 @@ public:
             mobilityFracture_[nPhaseIdx] =
                     MaterialLaw::krn(materialParamsFracture, fluidStateFracture_.saturation(wPhaseIdx))
                         / fluidStateFracture_.viscosity(nPhaseIdx);
-
-            // derivative resulted from BrooksCorey pc_Sw formulation
-            dsm_dsf_ = (1 - problem.spatialParams().swrm_) / (1 - problem.spatialParams().swrf_)
-                    * pow((problem.spatialParams().pdm_/ problem.spatialParams().pdf_),problem.spatialParams().lambdaM_)
-                    * (problem.spatialParams().lambdaM_ / problem.spatialParams().lambdaF_)
-                    * pow((satWFracture_ - problem.spatialParams().swrf_ ) / (1 - problem.spatialParams().swrf_),
-                            (problem.spatialParams().lambdaM_ / problem.spatialParams().lambdaF_) - 1);
         }// end if (node)
         ///////////////////////////////////////////////////////////////////////////////
         else
@@ -343,12 +336,6 @@ public:
     Scalar permeabilityFracture() const
     { return permeabilityFracture_; }
 
-    /*!
-     * \brief Returns the derivative dsm/dsf
-     */
-    Scalar dsm_dsf() const
-    { return dsm_dsf_;}
-
 protected:
     FluidState fluidState_;
     FluidState fluidStateFracture_;
@@ -370,7 +357,6 @@ protected:
     Scalar pcFracture_;
     Scalar pcMatrix_;
     Scalar pEntryMatrix_;
-    Scalar dsm_dsf_;
 
     bool isNodeOnFracture_;