[test] Include plotMaterialLaw in 3p test

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();
     }
 
     /*!

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_;
 };
 
 }

test/porousmediumflow/3p/implicit/test_box3p.input

@@ -15,3 +15,6 @@ porosity = 0.40
 vanGenuchtenAlpha = 0.0005
 vanGenuchtenN = 4.0
 
+[Output]
+PlotFluidMatrixInteractions = false
+

test/porousmediumflow/3p/implicit/test_cc3p.input

@@ -15,3 +15,6 @@ porosity = 0.40
 vanGenuchtenAlpha = 0.0005
 vanGenuchtenN = 4.0
 
+[Output]
+PlotFluidMatrixInteractions = false
+