diff --git a/test/porousmediumflow/3p/implicit/infiltration3pproblem.hh b/test/porousmediumflow/3p/implicit/infiltration3pproblem.hh
index 7a73d40545438b5a058dddb3f8fcf506ef051079..4dc9f7878080c4f4d9208f7e69c811b93b8d6176 100644
--- a/test/porousmediumflow/3p/implicit/infiltration3pproblem.hh
+++ b/test/porousmediumflow/3p/implicit/infiltration3pproblem.hh
@@ -143,6 +143,8 @@ public:
         FluidSystem::init(282.15, 284.15, 3, 8e4, 3e5, 200);
 
         name_               = GET_RUNTIME_PARAM(TypeTag, std::string, Problem.Name);
+
+        this->spatialParams().plotMaterialLaw();
     }
 
     /*!
diff --git a/test/porousmediumflow/3p/implicit/infiltration3pspatialparams.hh b/test/porousmediumflow/3p/implicit/infiltration3pspatialparams.hh
index ea8b9fcc01d14286f9cc02088065c42da44d187b..c7bc3e18a863d3b6d0f2337fa4e992a66d89b201 100644
--- a/test/porousmediumflow/3p/implicit/infiltration3pspatialparams.hh
+++ b/test/porousmediumflow/3p/implicit/infiltration3pspatialparams.hh
@@ -30,6 +30,7 @@
 #include <dumux/material/fluidmatrixinteractions/3p/regularizedparkervangen3p.hh>
 #include <dumux/material/fluidmatrixinteractions/3p/regularizedparkervangen3pparams.hh>
 #include <dumux/material/fluidmatrixinteractions/3p/efftoabslaw.hh>
+#include <dumux/io/plotmateriallaw3p.hh>
 namespace Dumux
 {
 
@@ -117,11 +118,26 @@ public:
         // parameters for adsorption
         materialParams_.setKdNAPL(0.);
         materialParams_.setRhoBulk(1500.);
+
+        plotFluidMatrixInteractions_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Output,
+                                                                    PlotFluidMatrixInteractions);
     }
 
     ~InfiltrationThreePSpatialParams()
     {}
 
+     /*!
+     * \brief This is called from the problem and creates a gnuplot output
+     *        of e.g the pc-Sw curve
+     */
+    void plotMaterialLaw()
+    {
+        PlotMaterialLaw<TypeTag> plotMaterialLaw(plotFluidMatrixInteractions_);
+
+        plotMaterialLaw.plotpc(materialParams_);
+        plotMaterialLaw.plotkr(materialParams_);
+    }
+
     /*!
      * \brief Intrinsic permability
      *
@@ -188,6 +204,8 @@ private:
     Scalar porosity_;
 
     MaterialLawParams materialParams_;
+
+    bool plotFluidMatrixInteractions_;
 };
 
 }
diff --git a/test/porousmediumflow/3p/implicit/test_box3p.input b/test/porousmediumflow/3p/implicit/test_box3p.input
index edd8a86f233adbb0cbaa2229763fe06de8a7f24e..6445dff4ae034abcd0ce6f15d05dc78e0c9e4349 100644
--- a/test/porousmediumflow/3p/implicit/test_box3p.input
+++ b/test/porousmediumflow/3p/implicit/test_box3p.input
@@ -15,3 +15,6 @@ porosity = 0.40
 vanGenuchtenAlpha = 0.0005
 vanGenuchtenN = 4.0
 
+[Output]
+PlotFluidMatrixInteractions = false
+
diff --git a/test/porousmediumflow/3p/implicit/test_cc3p.input b/test/porousmediumflow/3p/implicit/test_cc3p.input
index 08652caa0f1966db45b69feded08c298bc2ec9e7..764e1a3f0c4ab0036e27c25bc2f9177e9a33673e 100644
--- a/test/porousmediumflow/3p/implicit/test_cc3p.input
+++ b/test/porousmediumflow/3p/implicit/test_cc3p.input
@@ -15,3 +15,6 @@ porosity = 0.40
 vanGenuchtenAlpha = 0.0005
 vanGenuchtenN = 4.0
 
+[Output]
+PlotFluidMatrixInteractions = false
+