diff --git a/test/porousmediumflow/2pncmin/implicit/nonisothermal/CMakeLists.txt b/test/porousmediumflow/2pncmin/implicit/nonisothermal/CMakeLists.txt
index 082a0bbdb4c4b6ef36d12ccebb05deb233a15301..94ebabe05e17bdef8b6784f263eb96d04b4ca9fa 100644
--- a/test/porousmediumflow/2pncmin/implicit/nonisothermal/CMakeLists.txt
+++ b/test/porousmediumflow/2pncmin/implicit/nonisothermal/CMakeLists.txt
@@ -1,12 +1,22 @@
 dune_symlink_to_source_files(FILES params.input)
 
 # isothermal tests
-dune_add_test(NAME test_2pncminni_salinization
+dune_add_test(NAME test_2pncminni_salinization_box
               LABELS porousmediumflow 2pncmin
               SOURCES main.cc
               COMPILE_DEFINITIONS TYPETAG=SalinizationBox
               COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
               CMD_ARGS      --script fuzzy
-                            --files ${CMAKE_SOURCE_DIR}/test/references/test_2pncminni_salinization-reference.vtu
-                               ${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization-00034.vtu
-                            --command "${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization -ParameterFile params.input -Problem.Name test_2pncminni_salinization")
+                            --files ${CMAKE_SOURCE_DIR}/test/references/test_2pncminni_salinization_box-reference.vtu
+                               ${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization_box-00035.vtu
+                               --command "${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization_box -ParameterFile params.input -Problem.Name test_2pncminni_salinization_box")
+
+dune_add_test(NAME test_2pncminni_salinization_tpfa
+              LABELS porousmediumflow 2pncmin
+              SOURCES main.cc
+              COMPILE_DEFINITIONS TYPETAG=SalinizationCCTpfa
+              COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+              CMD_ARGS      --script fuzzy
+                            --files ${CMAKE_SOURCE_DIR}/test/references/test_2pncminni_salinization_tpfa-reference.vtu
+                               ${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization_tpfa-00035.vtu
+                            --command "${CMAKE_CURRENT_BINARY_DIR}/test_2pncminni_salinization_tpfa -ParameterFile params.input -Problem.Name test_2pncminni_salinization_tpfa")
diff --git a/test/porousmediumflow/2pncmin/implicit/nonisothermal/problem.hh b/test/porousmediumflow/2pncmin/implicit/nonisothermal/problem.hh
index 8756092803ca45855f217ec4f27838ddbb018349..51f6b08f29901c3bb4a60cc0b59ac1034778eea1 100644
--- a/test/porousmediumflow/2pncmin/implicit/nonisothermal/problem.hh
+++ b/test/porousmediumflow/2pncmin/implicit/nonisothermal/problem.hh
@@ -43,7 +43,7 @@
 namespace Dumux {
 /*!
  * \ingroup TwoPNCMinTests
- * \brief Problem where brine is evaporating at the top boundary. The system is closed the remaining boundaries.
+ * \brief Problem where brine is evaporating at the top boundary. The system is closed at the remaining boundaries.
  */
 template <class TypeTag>
 class SalinizationProblem;
diff --git a/test/references/test_2pncminni_salinization-reference.vtu b/test/references/test_2pncminni_salinization_box-reference.vtu
similarity index 55%
rename from test/references/test_2pncminni_salinization-reference.vtu
rename to test/references/test_2pncminni_salinization_box-reference.vtu
index 1d14a543b1d604fcc81529dec7c6120bcfbecb57..9a18a0de3113b90cd90ece975d7b37f8b5a6502e 100644
--- a/test/references/test_2pncminni_salinization-reference.vtu
+++ b/test/references/test_2pncminni_salinization_box-reference.vtu
@@ -4,76 +4,76 @@
     <Piece NumberOfCells="10" NumberOfPoints="18">
       <PointData Scalars="S_liq" Vectors="velocity_liq (m/s)">
         <DataArray type="Float32" Name="S_liq" NumberOfComponents="1" format="ascii">
-          0.371972 0.371972 0.354829 0.354829 0.371972 0.354829 0.328148 0.328148 0.328148 0.315167 0.315167 0.315167
-          0.300425 0.300425 0.300425 0.283783 0.283783 0.283783
+          0.371984 0.371984 0.354842 0.354842 0.371984 0.354842 0.32817 0.32817 0.32817 0.315198 0.315198 0.315198
+          0.300464 0.300464 0.300464 0.283892 0.283892 0.283892
         </DataArray>
         <DataArray type="Float32" Name="p_liq" NumberOfComponents="1" format="ascii">
-          92441.8 92441.8 92339.1 92339.1 92441.8 92339.1 92142.7 92142.7 92142.7 92022.1 92022.1 92022.1
-          91851.9 91851.9 91851.9 91583.7 91583.7 91583.7
+          92441.6 92441.6 92339 92339 92441.6 92339 92142.7 92142.7 92142.7 92022.2 92022.2 92022.2
+          91852.1 91852.1 91852.1 91585.9 91585.9 91585.9
         </DataArray>
         <DataArray type="Float32" Name="rho_liq" NumberOfComponents="1" format="ascii">
-          1012.93 1012.93 1012.93 1012.93 1012.93 1012.93 1012.94 1012.94 1012.94 1013.29 1013.29 1013.29
-          1021.51 1021.51 1021.51 1200.29 1200.29 1200.29
+          1012.93 1012.93 1012.93 1012.93 1012.93 1012.93 1012.94 1012.94 1012.94 1013.27 1013.27 1013.27
+          1021.11 1021.11 1021.11 1200.29 1200.29 1200.29
         </DataArray>
         <DataArray type="Float32" Name="mob_liq" NumberOfComponents="1" format="ascii">
-          6.3638 6.3638 4.25016 4.25016 6.3638 4.25016 1.94427 1.94427 1.94427 1.19814 1.19814 1.19814
-          0.594915 0.594915 0.594915 0.111125 0.111125 0.111125
+          6.36616 6.36616 4.25205 4.25205 6.36616 4.25205 1.94592 1.94592 1.94592 1.1998 1.1998 1.1998
+          0.596671 0.596671 0.596671 0.112089 0.112089 0.112089
         </DataArray>
         <DataArray type="Float32" Name="S_gas" NumberOfComponents="1" format="ascii">
-          0.628028 0.628028 0.645171 0.645171 0.628028 0.645171 0.671852 0.671852 0.671852 0.684833 0.684833 0.684833
-          0.699575 0.699575 0.699575 0.716217 0.716217 0.716217
+          0.628016 0.628016 0.645158 0.645158 0.628016 0.645158 0.67183 0.67183 0.67183 0.684802 0.684802 0.684802
+          0.699536 0.699536 0.699536 0.716108 0.716108 0.716108
         </DataArray>
         <DataArray type="Float32" Name="p_gas" NumberOfComponents="1" format="ascii">
           100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002 100002
           100002 100002 100002 100001 100001 100001
         </DataArray>
         <DataArray type="Float32" Name="rho_gas" NumberOfComponents="1" format="ascii">
-          1.21014 1.21014 1.21013 1.21013 1.21014 1.21013 1.21012 1.21012 1.21012 1.21012 1.21012 1.21012
-          1.21014 1.21014 1.21014 1.21073 1.21073 1.21073
+          1.21011 1.21011 1.21011 1.21011 1.21011 1.21011 1.2101 1.2101 1.2101 1.21009 1.21009 1.21009
+          1.21011 1.21011 1.21011 1.21071 1.21071 1.21071
         </DataArray>
         <DataArray type="Float32" Name="mob_gas" NumberOfComponents="1" format="ascii">
-          39281.7 39281.7 41576.5 41576.5 39281.7 41576.5 45236.5 45236.5 45236.5 47150.2 47150.2 47150.2
-          49764.6 49764.6 49764.6 52753.5 52753.5 52753.5
+          39279.6 39279.6 41574.1 41574.1 39279.6 41574.1 45232.9 45232.9 45232.9 47144.6 47144.6 47144.6
+          49756.9 49756.9 49756.9 52734.6 52734.6 52734.6
         </DataArray>
         <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
-          7560.19 7560.19 7662.75 7662.75 7560.19 7662.75 7859.09 7859.09 7859.09 7979.69 7979.69 7979.69
-          8149.9 8149.9 8149.9 8417.42 8417.42 8417.42
+          7560.12 7560.12 7662.67 7662.67 7560.12 7662.67 7858.9 7858.9 7858.9 7979.38 7979.38 7979.38
+          8149.4 8149.4 8149.4 8415.26 8415.26 8415.26
         </DataArray>
         <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
           0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
-          0.4 0.4 0.4 0.399824 0.399824 0.399824
+          0.4 0.4 0.4 0.399753 0.399753 0.399753
         </DataArray>
         <DataArray type="Float32" Name="x^H2O_liq" NumberOfComponents="1" format="ascii">
-          0.993731 0.993731 0.993731 0.993731 0.993731 0.993731 0.993726 0.993726 0.993726 0.993565 0.993565 0.993565
-          0.989756 0.989756 0.989756 0.902258 0.902258 0.902258
+          0.993731 0.993731 0.993731 0.993731 0.993731 0.993731 0.993726 0.993726 0.993726 0.993572 0.993572 0.993572
+          0.98994 0.98994 0.98994 0.902257 0.902257 0.902257
         </DataArray>
         <DataArray type="Float32" Name="x^Air_liq" NumberOfComponents="1" format="ascii">
-          1.69669e-05 1.69669e-05 1.69666e-05 1.69666e-05 1.69669e-05 1.69666e-05 1.69657e-05 1.69657e-05 1.69657e-05 1.69654e-05 1.69654e-05 1.69654e-05
-          1.69661e-05 1.69661e-05 1.69661e-05 1.69884e-05 1.69884e-05 1.69884e-05
+          1.6965e-05 1.6965e-05 1.69647e-05 1.69647e-05 1.6965e-05 1.69647e-05 1.69639e-05 1.69639e-05 1.69639e-05 1.69636e-05 1.69636e-05 1.69636e-05
+          1.69642e-05 1.69642e-05 1.69642e-05 1.69866e-05 1.69866e-05 1.69866e-05
         </DataArray>
         <DataArray type="Float32" Name="x^NaCl_liq" NumberOfComponents="1" format="ascii">
-          0.00625221 0.00625221 0.00625235 0.00625235 0.00625221 0.00625235 0.00625709 0.00625709 0.00625709 0.0064178 0.0064178 0.0064178
-          0.0102271 0.0102271 0.0102271 0.0977251 0.0977251 0.0977251
+          0.00625221 0.00625221 0.00625235 0.00625235 0.00625221 0.00625235 0.00625695 0.00625695 0.00625695 0.0064111 0.0064111 0.0064111
+          0.0100428 0.0100428 0.0100428 0.097726 0.097726 0.097726
         </DataArray>
         <DataArray type="Float32" Name="rhoMolar_liq" NumberOfComponents="1" format="ascii">
-          55447.9 55447.9 55447.9 55447.9 55447.9 55447.9 55447.8 55447.8 55447.8 55447.1 55447.1 55447.1
-          55430 55430 55430 54642.7 54642.7 54642.7
+          55447.9 55447.9 55447.8 55447.8 55447.9 55447.8 55447.8 55447.8 55447.8 55447.1 55447.1 55447.1
+          55430.9 55430.9 55430.9 54642.6 54642.6 54642.6
         </DataArray>
         <DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
-          0.0149431 0.0149431 0.0149438 0.0149438 0.0149431 0.0149438 0.014946 0.014946 0.014946 0.0149445 0.0149445 0.0149445
-          0.014888 0.014888 0.014888 0.0135726 0.0135726 0.0135726
+          0.0149479 0.0149479 0.0149487 0.0149487 0.0149479 0.0149487 0.0149508 0.0149508 0.0149508 0.0149493 0.0149493 0.0149493
+          0.0148954 0.0148954 0.0148954 0.0135767 0.0135767 0.0135767
         </DataArray>
         <DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
-          0.985057 0.985057 0.985056 0.985056 0.985057 0.985056 0.985054 0.985054 0.985054 0.985056 0.985056 0.985056
-          0.985112 0.985112 0.985112 0.986427 0.986427 0.986427
+          0.985052 0.985052 0.985051 0.985051 0.985052 0.985051 0.985049 0.985049 0.985049 0.985051 0.985051 0.985051
+          0.985105 0.985105 0.985105 0.986423 0.986423 0.986423
         </DataArray>
         <DataArray type="Float32" Name="x^NaCl_gas" NumberOfComponents="1" format="ascii">
           -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0 -0
           -0 -0 -0 -0 -0 -0
         </DataArray>
         <DataArray type="Float32" Name="rhoMolar_gas" NumberOfComponents="1" format="ascii">
-          42.0241 42.0241 42.024 42.024 42.0241 42.024 42.0236 42.0236 42.0236 42.0234 42.0234 42.0234
-          42.0233 42.0233 42.0233 42.0228 42.0228 42.0228
+          42.0233 42.0233 42.0231 42.0231 42.0233 42.0231 42.0228 42.0228 42.0228 42.0227 42.0227 42.0227
+          42.0225 42.0225 42.0225 42.0222 42.0222 42.0222
         </DataArray>
         <DataArray type="Float32" Name="phase presence" NumberOfComponents="1" format="ascii">
           3 3 3 3 3 3 3 3 3 3 3 3
@@ -81,29 +81,29 @@
         </DataArray>
         <DataArray type="Float32" Name="precipitateVolumeFraction^NaCl" NumberOfComponents="1" format="ascii">
           0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0.00017601 0.00017601 0.00017601
+          0 0 0 0.000246746 0.000246746 0.000246746
         </DataArray>
         <DataArray type="Float32" Name="T" NumberOfComponents="1" format="ascii">
-          286.208 286.208 286.209 286.209 286.208 286.209 286.211 286.211 286.211 286.212 286.212 286.212
-          286.213 286.213 286.213 286.214 286.214 286.214
+          286.213 286.213 286.214 286.214 286.213 286.214 286.216 286.216 286.216 286.217 286.217 286.217
+          286.218 286.218 286.218 286.218 286.218 286.218
         </DataArray>
         <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
-          -6.7602e-22 2.38612e-09 0 -3.79683e-22 2.38612e-09 0 -6.3085e-22 4.64661e-09 0 -3.23156e-22 4.64661e-09 0
-          -8.3345e-23 2.38612e-09 0 -1.5462e-23 4.64661e-09 0 -3.42343e-22 8.41592e-09 0 -2.7727e-22 8.41592e-09 0
-          -2.12196e-22 8.41592e-09 0 -4.95383e-22 1.05518e-08 0 -3.18852e-22 1.05518e-08 0 -1.42321e-22 1.05518e-08 0
-          -1.73772e-22 1.1066e-08 0 -5.82074e-23 1.1066e-08 0 5.73568e-23 1.1066e-08 0 1.94852e-22 1.0953e-08 0
-          1.33006e-22 1.0953e-08 0 7.11605e-23 1.0953e-08 0
+          -2.96447e-22 2.38548e-09 0 -5.92894e-22 2.38548e-09 0 -4.95004e-22 4.64498e-09 0 -7.42505e-22 4.64498e-09 0
+          -8.89342e-22 2.38548e-09 0 -9.90007e-22 4.64498e-09 0 -2.26535e-22 8.4126e-09 0 -6.56951e-22 8.4126e-09 0
+          -1.08737e-21 8.4126e-09 0 -2.23481e-22 1.05503e-08 0 -6.56475e-22 1.05503e-08 0 -1.08947e-21 1.05503e-08 0
+          -1.66671e-22 1.104e-08 0 -2.01397e-22 1.104e-08 0 -2.36123e-22 1.104e-08 0 5.20911e-23 1.09002e-08 0
+          2.60456e-23 1.09002e-08 0 0 1.09002e-08 0
         </DataArray>
         <DataArray type="Float32" Name="velocity_gas (m/s)" NumberOfComponents="3" format="ascii">
-          -9.77477e-18 -2.47197e-09 0 -5.91631e-18 -2.47197e-09 0 -6.62494e-18 -4.83592e-09 0 -6.20143e-18 -4.83592e-09 0
-          -2.05785e-18 -2.47197e-09 0 -5.77792e-18 -4.83592e-09 0 -5.75993e-18 -7.66282e-09 0 -7.17523e-18 -7.66282e-09 0
-          -8.59053e-18 -7.66282e-09 0 -6.91496e-18 2.53119e-08 0 -1.30386e-17 2.53119e-08 0 -1.91623e-17 2.53119e-08 0
-          5.96753e-18 1.0712e-06 0 -8.72983e-18 1.0712e-06 0 -2.34272e-17 1.0712e-06 0 -3.41432e-19 2.08366e-06 0
-          -2.79445e-18 2.08366e-06 0 -5.24746e-18 2.08366e-06 0
+          0 -2.41187e-09 0 -4.57275e-18 -2.41187e-09 0 -4.83987e-18 -4.70772e-09 0 -4.83987e-18 -4.70772e-09 0
+          -9.14549e-18 -2.41187e-09 0 -4.83987e-18 -4.70772e-09 0 -6.31896e-18 -7.86737e-09 0 -1.21113e-17 -7.86737e-09 0
+          -1.79037e-17 -7.86737e-09 0 -1.20744e-17 1.09526e-08 0 -3.12836e-17 1.09526e-08 0 -5.04928e-17 1.09526e-08 0
+          -2.08467e-17 6.65984e-07 0 -2.08475e-17 6.65984e-07 0 -2.08482e-17 6.65984e-07 0 -9.80292e-18 1.30133e-06 0
+          -7.35219e-18 1.30133e-06 0 -4.90146e-18 1.30133e-06 0
         </DataArray>
         <DataArray type="Float32" Name="Permeability" NumberOfComponents="1" format="ascii">
           1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13
-          9.99523e-14 9.99523e-14 9.99523e-14 9.98571e-14 9.98571e-14 9.98571e-14
+          9.99332e-14 9.99332e-14 9.99332e-14 9.97997e-14 9.97997e-14 9.97997e-14
         </DataArray>
       </PointData>
       <CellData Scalars="process rank">
diff --git a/test/references/test_2pncminni_salinization_tpfa-reference.vtu b/test/references/test_2pncminni_salinization_tpfa-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..400add18b90324c20e55ee425019fd88147db3a3
--- /dev/null
+++ b/test/references/test_2pncminni_salinization_tpfa-reference.vtu
@@ -0,0 +1,111 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfCells="10" NumberOfPoints="18">
+      <CellData Scalars="S_liq" Vectors="velocity_liq (m/s)">
+        <DataArray type="Float32" Name="S_liq" NumberOfComponents="1" format="ascii">
+          0.362084 0.362084 0.340382 0.340382 0.316467 0.316467 0.301085 0.301085 0.282783 0.282783
+        </DataArray>
+        <DataArray type="Float32" Name="p_liq" NumberOfComponents="1" format="ascii">
+          92385.2 92385.2 92240.5 92240.5 92036 92036 91861.3 91861.3 91564.8 91564.8
+        </DataArray>
+        <DataArray type="Float32" Name="rho_liq" NumberOfComponents="1" format="ascii">
+          1012.97 1012.97 1012.98 1012.98 1013.13 1013.13 1016.92 1016.92 1116.91 1116.91
+        </DataArray>
+        <DataArray type="Float32" Name="mob_liq" NumberOfComponents="1" format="ascii">
+          5.03768 5.03768 2.83968 2.83968 1.25285 1.25285 0.61688 0.61688 0.148255 0.148255
+        </DataArray>
+        <DataArray type="Float32" Name="S_gas" NumberOfComponents="1" format="ascii">
+          0.637916 0.637916 0.659618 0.659618 0.683533 0.683533 0.698915 0.698915 0.717217 0.717217
+        </DataArray>
+        <DataArray type="Float32" Name="p_gas" NumberOfComponents="1" format="ascii">
+          100003 100003 100003 100003 100003 100003 100003 100003 100002 100002
+        </DataArray>
+        <DataArray type="Float32" Name="rho_gas" NumberOfComponents="1" format="ascii">
+          1.21187 1.21187 1.21186 1.21186 1.21186 1.21186 1.21187 1.21187 1.21218 1.21218
+        </DataArray>
+        <DataArray type="Float32" Name="mob_gas" NumberOfComponents="1" format="ascii">
+          40641 40641 43590.7 43590.7 46987.1 46987.1 49693.9 49693.9 52974.1 52974.1
+        </DataArray>
+        <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+          7617.59 7617.59 7762.18 7762.18 7966.6 7966.6 8141.26 8141.26 8437.63 8437.63
+        </DataArray>
+        <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+          0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+        </DataArray>
+        <DataArray type="Float32" Name="x^H2O_liq" NumberOfComponents="1" format="ascii">
+          0.993731 0.993731 0.99373 0.99373 0.993658 0.993658 0.991904 0.991904 0.944436 0.944436
+        </DataArray>
+        <DataArray type="Float32" Name="x^Air_liq" NumberOfComponents="1" format="ascii">
+          1.71059e-05 1.71059e-05 1.71056e-05 1.71056e-05 1.71052e-05 1.71052e-05 1.71055e-05 1.71055e-05 1.71174e-05 1.71174e-05
+        </DataArray>
+        <DataArray type="Float32" Name="x^NaCl_liq" NumberOfComponents="1" format="ascii">
+          0.00625227 0.00625227 0.00625305 0.00625305 0.00632463 0.00632463 0.00807839 0.00807839 0.0555472 0.0555472
+        </DataArray>
+        <DataArray type="Float32" Name="rhoMolar_liq" NumberOfComponents="1" format="ascii">
+          55450.4 55450.4 55450.4 55450.4 55450 55450 55442.3 55442.3 55126.2 55126.2
+        </DataArray>
+        <DataArray type="Float32" Name="x^H2O_gas" NumberOfComponents="1" format="ascii">
+          0.0145878 0.0145878 0.0145886 0.0145886 0.0145886 0.0145886 0.0145633 0.0145633 0.0138667 0.0138667
+        </DataArray>
+        <DataArray type="Float32" Name="x^Air_gas" NumberOfComponents="1" format="ascii">
+          0.985412 0.985412 0.985411 0.985411 0.985411 0.985411 0.985437 0.985437 0.986133 0.986133
+        </DataArray>
+        <DataArray type="Float32" Name="x^NaCl_gas" NumberOfComponents="1" format="ascii">
+          -0 -0 -0 -0 -0 -0 -0 -0 -0 -0
+        </DataArray>
+        <DataArray type="Float32" Name="rhoMolar_gas" NumberOfComponents="1" format="ascii">
+          42.0784 42.0784 42.0782 42.0782 42.0781 42.0781 42.078 42.078 42.0778 42.0778
+        </DataArray>
+        <DataArray type="Float32" Name="phase presence" NumberOfComponents="1" format="ascii">
+          3 3 3 3 3 3 3 3 3 3
+        </DataArray>
+        <DataArray type="Float32" Name="precipitateVolumeFraction^NaCl" NumberOfComponents="1" format="ascii">
+          0 0 0 0 0 0 0 0 0 0
+        </DataArray>
+        <DataArray type="Float32" Name="T" NumberOfComponents="1" format="ascii">
+          285.841 285.841 285.842 285.842 285.843 285.843 285.843 285.843 285.844 285.844
+        </DataArray>
+        <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+          0 2.35559e-09 0 0 2.35559e-09 0 0 7.07053e-09 0 0 7.07053e-09 0
+          9.11568e-23 1.01265e-08 0 9.11568e-23 1.01265e-08 0 0 1.11427e-08 0 0 1.11427e-08 0
+          2.15739e-23 1.14623e-08 0 2.15739e-23 1.14623e-08 0
+        </DataArray>
+        <DataArray type="Float32" Name="velocity_gas (m/s)" NumberOfComponents="3" format="ascii">
+          0 -2.37122e-09 0 0 -2.37122e-09 0 0 -7.00433e-09 0 0 -7.00433e-09 0
+          3.41876e-18 -1.19557e-09 0 3.41876e-18 -1.19557e-09 0 0 3.1111e-07 0 0 3.1111e-07 0
+          3.85438e-18 6.15346e-07 0 3.85438e-18 6.15346e-07 0
+        </DataArray>
+        <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+          0 0 0 0 0 0 0 0 0 0
+        </DataArray>
+        <DataArray type="Float32" Name="Permeability" NumberOfComponents="1" format="ascii">
+          1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13 1e-13
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+          0 0 0 0.01 0 0 0 0.0075 0 0.01 0.0075 0
+          0.02 0 0 0.02 0.0075 0 0 0.015 0 0.01 0.015 0
+          0.02 0.015 0 0 0.016978 0 0.01 0.016978 0 0.02 0.016978 0
+          0 0.0186264 0 0.01 0.0186264 0 0.02 0.0186264 0 0 0.02 0
+          0.01 0.02 0 0.02 0.02 0
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+          0 1 3 2 1 4 5 3 2 3 7 6
+          3 5 8 7 6 7 10 9 7 8 11 10
+          9 10 13 12 10 11 14 13 12 13 16 15
+          13 14 17 16
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+          4 8 12 16 20 24 28 32 36 40
+        </DataArray>
+        <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+          9 9 9 9 9 9 9 9 9 9
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>