diff --git a/geiger/results/mfd/mfd_geiger_blocking.mat b/geiger/results/mfd/mfd_geiger_blocking.mat
index 427dfae7c0612ba83b615a8a208a43c7b2fe30d2..d2c1cbb851f98105f77542b3cd6f2a1d78d8bddd 100644
Binary files a/geiger/results/mfd/mfd_geiger_blocking.mat and b/geiger/results/mfd/mfd_geiger_blocking.mat differ
diff --git a/geiger/results/mfd/mfd_geiger_conductive.mat b/geiger/results/mfd/mfd_geiger_conductive.mat
index 02ffa519d10f80235ddc3cd76279e8120ae3dc91..2fdef9cfeddcf423fd697cfe9fc2642626b41b6e 100644
Binary files a/geiger/results/mfd/mfd_geiger_conductive.mat and b/geiger/results/mfd/mfd_geiger_conductive.mat differ
diff --git a/scripts/README.md b/scripts/README.md
index 709e971df35191875e430f3a09398fc47ec0a14f..a5c2073c11aac140e1928f74cc2107c78706158e 100644
--- a/scripts/README.md
+++ b/scripts/README.md
@@ -4,13 +4,13 @@ This folder contains several Matlab scripts which can be used to evaluate a solu
 
 #### 1.1 Preparation of the solution data
 
-The solution data has to be provided in form of a Matlab `mat`-file which has to contain the following variables:
+The matrix solution has to be provided in form of a Matlab `mat`-file which has to contain the following variables:
 
-* `Points`: an `np x 2` matrix that contains the x and y coordinates of the `np` vertices of the employed computational grid in columns `1` and `2`, respectively. 
-* `T`: an `nt x nvelem` connectivity matrix which contains in each row the `nvelem` vertex indices with respect to the numbering induced by `Points`.
-* `P`: an `nt`-size vector containing the solution value for each cell with respect to the numbering induced by `T`.
+* `Points`: an `np x 2` matrix that contains the x and y coordinates of the `np` vertices of the computational grid in columns `1` and `2`, respectively. 
+* `T`: an `nt x nvcell` cell connectivity matrix which contains in each row a cell's `nvcell` vertex indices with respect to the numbering induced by `Points`.
+* `P`: an `nt`-size vector containing the cell-wise solution values with respect to the numbering induced by `T`.
 
-The original solution data is available as a VTK file? A helper function `hdf5_to_mat.m` is provided which extracts the required `mat`-file from an HDF5-file, which in turn can be exported from Paraview: select `File` -> `Save Data...` -> `Files of type: Xdmf Data File`, use the produced `h5`-file.
+Is your original solution data available as a VTK file? A helper function `hdf5_to_mat.m` is provided which extracts the required `mat`-file from an HDF5-file, which in turn can be exported from Paraview: select `File` -> `Save Data...` -> `Files of type: Xdmf Data File`, use the produced `h5`-file.
 
 #### 1.2 Performing the comparison
 
@@ -19,3 +19,10 @@ The function [compute_2d_error](compute_2d_error.m) can be used to calculate a r
 function [relativeMatrixError] = compute_2d_error(coarseFile, referenceFile)
 ```
 `coarseFile` is the `mat`-file explained above, `referenceFile` is one of the `mat`-files found in the `results/mfd`-subfolders, for example `../hydrocoin/results/mfd/mfd_hydrocoin.mat`.
+
+### 2 Comparison of the solution in the fractures
+
+#### 2.1 Preparation of the solution data
+
+#### 2.2 Performing the comparison
+