diff --git a/doc/handbook/tutorial-decoupled.tex b/doc/handbook/tutorial-decoupled.tex
index d0ca25af9..21dd7fe07 100644
--- a/doc/handbook/tutorial-decoupled.tex
+++ b/doc/handbook/tutorial-decoupled.tex
@@ -166,7 +166,7 @@ Hence, they all feature a common argument list:
 \end{itemize}
 In the following, there are the methods for general parameters, source- or
 sinkterms, boundary conditions (lines \ref{tutorial-decoupled:bctypePress} to
-\ref{tutorial-decoupled:neumannSat}) and initial values for the transported
+\ref{tutorial-decoupled:neumann}) and initial values for the transported
 quantity in line \label{tutorial-decoupled:initSat}. For more information
 on the functions, it is referred to the documentation in the code.
 
diff --git a/examples/tutorialproblem_decoupled.hh b/examples/tutorialproblem_decoupled.hh
index afd29ec3f..3c77df509 100644
--- a/examples/tutorialproblem_decoupled.hh
+++ b/examples/tutorialproblem_decoupled.hh
@@ -253,11 +253,11 @@ public:
     {
         return 1;
     }
-    //! Value for pressure neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
+    //! Value for neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
     /*! In case of a neumann boundary condition, the flux of matter
      *  is returned as a vector.
      */
-    std::vector<Scalar> neumannPress(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:neumannPress}@*/
+    std::vector<Scalar> neumann(const GlobalPosition& globalPos, const Intersection& intersection) const /*@\label{tutorial-decoupled:neumann}@*/
     {
         std::vector<Scalar> neumannFlux(2,0.0);
         if (globalPos[0] > this->bboxMax()[0] - eps_)
@@ -266,14 +266,6 @@ public:
         }
         return neumannFlux;
     }
-    //! Value for transport neumann boundary condition \f$ [\frac{kg}{m^3 \cdot s}] \f$.
-    /*! In case of a neumann boundary condition for the transport equation
-     *  the flux of matter for the primary variable is returned as a scalar.
-     */
-    Scalar neumannSat(const GlobalPosition& globalPos, const Intersection& intersection, Scalar factor) const /*@\label{tutorial-decoupled:neumannSat}@*/
-    {
-        return 0;
-    }
     //! Saturation initial condition (dimensionless)
     /*! The problem is initialized with the following saturation.
      */