diff --git a/examples/2pinfiltration/.doc_config b/examples/2pinfiltration/.doc_config
index 4151fcbee0b494affb2ee3939d9f42481c16eb03..a8f7fb29c68864f042611612790f533869144d5e 100644
--- a/examples/2pinfiltration/.doc_config
+++ b/examples/2pinfiltration/.doc_config
@@ -4,7 +4,6 @@
         "spatialparams.hh",
         "problem.hh",
         "properties.hh",
-        "main.cc",
-        "doc/results.md"
+        "main.cc"
     ]
 }
diff --git a/examples/2pinfiltration/README.md b/examples/2pinfiltration/README.md
index 724c139acfd42c8d65fa286c235bcbbc40aa18cf..a4d4e86babc2205a337de424cf4534634bf21a27 100644
--- a/examples/2pinfiltration/README.md
+++ b/examples/2pinfiltration/README.md
@@ -1,13 +1,38 @@
 <!-- Important: This file has been automatically generated by generate_example_docs.py. Do not edit this file directly! -->
 
-This tutorial is similar to tests/porousmediumflow/2p/adaptive and restricted to the cell-centered finite volume TPFA discretization scheme.
-You need [ALUGrid][0] in order to compile and run it.
+# Two-phase flow infiltration problem with adaptive grid
 
-# Two-phase flow with infiltration and adaptive grid
+In this example, you will learn how to
+* solve a two-phase flow in porous media problem with two immiscible phases
+* set boundary conditions and a simple injection well
+* specify a lens with different porous material parameters
+* use adaptive grid refinement around the saturation front
+* specify a point source
+* read the initial solution from a text file
 
-## Problem set-up
-In this example we model a soil contamination problem where DNAPL infiltrates a porous medium. The initial distribution of DNAPL is known and we can read it from a txt-file.
+__Prerequisites:__ You need [dune-alugrid](https://gitlab.dune-project.org/extensions/dune-alugrid) in order to compile and run this example.
 
+__Result:__ The resulting saturation distribution in this example will look like this:
+
+![](./img/test_2p_pointsource_adaptive.png)
+
+
+## Folder layout and files
+
+```
+└── 2pinfiltration/
+    ├── CMakeLists.txt          -> build system file
+    ├── main.cc                 -> main program flow
+    ├── params.input            -> runtime parameters
+    ├── properties.hh           -> compile time configuration
+    ├── problem.hh              -> boundary & initial conditions
+    ├── spatialparams.hh        -> spatial parameter fields
+    └── initialsolutioncc.txt   -> text file with initial solution
+```
+
+## Scenario and mathematical model
+
+We model a soil contamination problem where DNAPL infiltrates a porous medium. The initial distribution of DNAPL is known and we can read it from a txt-file.
 To describe that problem we use a two phase model of two immiscible fluids with the multiphase Darcy's law as the description of momentum, i.e.:
 
 ```math
@@ -33,19 +58,15 @@ The lens and the initial saturation can be seen in Figures 1 and 2.
 
 ![](./img/test_2p_pointsource_lens.png)
 
-![](./img/test_2p_pointsource_initial.png)
-
 At the left and the right boundary of the domain we set a linear pressure gradient as a Dirichlet boundary condition. On the upper and lower boundary we set Neumann boundary conditions.
 DNAPL enters the model domain at the upper boundary between 1.75m ≤ x ≤ 2m with a rate of 0.04 kg/ms. That means that we set the value for the Neumann boundary flux to that rate in between 1.75m and 2m. On the rest of the Neumann boundary we set the flux to zero, which means we define it as a no-flow boundary.
-In addition, the DNAPL is injected at a point source at x = 0.502 and y = 3.02 with a rate of 0.1 kg/s.
+In addition, the DNAPL is injected at a point source at x = 0.502m and y = 3.02m with a rate of 0.1 kg/s.
 
-## Discretization
-We discretize the equations with a cell-centered finite volume TPFA scheme in space and an implicit Euler scheme in time. We use Newton's method to solve the system of nonlinear equations. For more information about the discretization please have a look at the [handbook](https://dumux.org/handbook).
-
-## Adaptive grid
+We discretize the equations with a cell-centered finite volume TPFA scheme in space and an implicit Euler scheme in time. We use Newton's method to solve the system of nonlinear equations.
 The grid is adapitvely refined around the injection. The adaptive behaviour can be changed with input parameters in the `params.input` file.
+For more information about the discretization please have a look at the [handbook](https://dumux.org/handbook).
 
-[0]: https://gitlab.dune-project.org/extensions/dune-alugrid
+# Implementation
 
 
 ## The file `spatialparams.hh`
@@ -1060,7 +1081,3 @@ catch (...)
 ```
 
 
-## Results
-The 2p-infiltration model computes the water saturation distribution after 500 seconds according tthe following figure:
-
-![](./img/test_2p_pointsource_adaptive.png)
diff --git a/examples/2pinfiltration/doc/intro.md b/examples/2pinfiltration/doc/intro.md
index 6829a1eebbb23bf92fd9b3b3e1ccd7fc20c7d36f..cd6657b35a1a8248fce55b1ad696e5e9616fa73a 100644
--- a/examples/2pinfiltration/doc/intro.md
+++ b/examples/2pinfiltration/doc/intro.md
@@ -1,11 +1,36 @@
-This tutorial is similar to tests/porousmediumflow/2p/adaptive and restricted to the cell-centered finite volume TPFA discretization scheme.
-You need [ALUGrid][0] in order to compile and run it.
+# Two-phase flow infiltration problem with adaptive grid
 
-# Two-phase flow with infiltration and adaptive grid
+In this example, you will learn how to
+* solve a two-phase flow in porous media problem with two immiscible phases
+* set boundary conditions and a simple injection well
+* specify a lens with different porous material parameters
+* use adaptive grid refinement around the saturation front
+* specify a point source
+* read the initial solution from a text file
 
-## Problem set-up
-In this example we model a soil contamination problem where DNAPL infiltrates a porous medium. The initial distribution of DNAPL is known and we can read it from a txt-file.
+__Prerequisites:__ You need [dune-alugrid](https://gitlab.dune-project.org/extensions/dune-alugrid) in order to compile and run this example.
 
+__Result:__ The resulting saturation distribution in this example will look like this:
+
+![](./img/test_2p_pointsource_adaptive.png)
+
+
+## Folder layout and files
+
+```
+└── 2pinfiltration/
+    ├── CMakeLists.txt          -> build system file
+    ├── main.cc                 -> main program flow
+    ├── params.input            -> runtime parameters
+    ├── properties.hh           -> compile time configuration
+    ├── problem.hh              -> boundary & initial conditions
+    ├── spatialparams.hh        -> spatial parameter fields
+    └── initialsolutioncc.txt   -> text file with initial solution
+```
+
+## Scenario and mathematical model
+
+We model a soil contamination problem where DNAPL infiltrates a porous medium. The initial distribution of DNAPL is known and we can read it from a txt-file.
 To describe that problem we use a two phase model of two immiscible fluids with the multiphase Darcy's law as the description of momentum, i.e.:
 
 ```math
@@ -31,16 +56,12 @@ The lens and the initial saturation can be seen in Figures 1 and 2.
 
 ![](./img/test_2p_pointsource_lens.png)
 
-![](./img/test_2p_pointsource_initial.png)
-
 At the left and the right boundary of the domain we set a linear pressure gradient as a Dirichlet boundary condition. On the upper and lower boundary we set Neumann boundary conditions.
 DNAPL enters the model domain at the upper boundary between 1.75m ≤ x ≤ 2m with a rate of 0.04 kg/ms. That means that we set the value for the Neumann boundary flux to that rate in between 1.75m and 2m. On the rest of the Neumann boundary we set the flux to zero, which means we define it as a no-flow boundary.
-In addition, the DNAPL is injected at a point source at x = 0.502 and y = 3.02 with a rate of 0.1 kg/s.
-
-## Discretization
-We discretize the equations with a cell-centered finite volume TPFA scheme in space and an implicit Euler scheme in time. We use Newton's method to solve the system of nonlinear equations. For more information about the discretization please have a look at the [handbook](https://dumux.org/handbook).
+In addition, the DNAPL is injected at a point source at x = 0.502m and y = 3.02m with a rate of 0.1 kg/s.
 
-## Adaptive grid
+We discretize the equations with a cell-centered finite volume TPFA scheme in space and an implicit Euler scheme in time. We use Newton's method to solve the system of nonlinear equations.
 The grid is adapitvely refined around the injection. The adaptive behaviour can be changed with input parameters in the `params.input` file.
+For more information about the discretization please have a look at the [handbook](https://dumux.org/handbook).
 
-[0]: https://gitlab.dune-project.org/extensions/dune-alugrid
+# Implementation
diff --git a/examples/2pinfiltration/doc/results.md b/examples/2pinfiltration/doc/results.md
deleted file mode 100644
index 0429a0d52ea3909cbbc279138164142ea1313df7..0000000000000000000000000000000000000000
--- a/examples/2pinfiltration/doc/results.md
+++ /dev/null
@@ -1,4 +0,0 @@
-## Results
-The 2p-infiltration model computes the water saturation distribution after 500 seconds according tthe following figure:
-
-![](./img/test_2p_pointsource_adaptive.png)
diff --git a/examples/2pinfiltration/img/test_2p_pointsource_initial.png b/examples/2pinfiltration/img/test_2p_pointsource_initial.png
deleted file mode 100644
index e9bb688c32dc26ef851a5887aa390c50c971971c..0000000000000000000000000000000000000000
Binary files a/examples/2pinfiltration/img/test_2p_pointsource_initial.png and /dev/null differ