diff --git a/appl/coupled/stokes_darcy/beaversjoseph/params.input b/appl/coupled/stokes_darcy/beaversjoseph/params.input
index 0219f08f968d3816c2a1845d72ff97f6fb7c1957..1f0bfcae9dcc6d0e3c78a90e59a1707282ee32ef 100644
--- a/appl/coupled/stokes_darcy/beaversjoseph/params.input
+++ b/appl/coupled/stokes_darcy/beaversjoseph/params.input
@@ -29,6 +29,7 @@ NumericDifference.BaseEpsilon = 1e-8
 [Darcy]
 Problem.Name = darcy
 SpatialParams.Permeability = 1e-4
+SpatialParams.Porosity = 1.0
 SpatialParams.AlphaBeaversJoseph = 1.0
 
 [Vtk]
diff --git a/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh b/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
index 18113c4df250ea5b77dbf0da44e3c24f0cff4060..18bd803740fe490ab4317751d92c23cf7f52cac8 100644
--- a/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
+++ b/appl/coupled/stokes_darcy/beaversjoseph/spatialparams.hh
@@ -56,6 +56,7 @@ public:
     {
         permeability_ = getParam<Scalar>("Darcy.SpatialParams.Permeability");
         alphaBJ_ = getParam<Scalar>("Darcy.SpatialParams.AlphaBeaversJoseph");
+        porosity_ = getParam<Scalar>("Darcy.SpatialParams.Porosity");
     }
 
     /*!
@@ -72,7 +73,7 @@ public:
      * \param globalPos The global position
      */
     Scalar porosityAtPos(const GlobalPosition& globalPos) const
-    { return 0.4; }
+    { return porosity_; }
 
     /*! \brief Defines the Beavers-Joseph coefficient in [-].
      *