diff --git a/CMakeLists.txt b/CMakeLists.txt
index 19ac5373fad8b3ea1cbd24e02211c86b2f78c45d..a781214e4373602352394f97e3e7d70b9d412dfe 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -392,7 +392,11 @@ file(COPY test/boxmodels/mpnc/forchheimer1p-reference.vtp DESTINATION references
 file(COPY test/boxmodels/mpnc/test_forchheimer1p.input DESTINATION test/boxmodels/mpnc)
 file(COPY test/boxmodels/mpnc/forchheimer2p-reference.vtu DESTINATION references)
 file(COPY test/boxmodels/mpnc/test_forchheimer2p.input DESTINATION test/boxmodels/mpnc)
-file(COPY test/boxmodels/richards/richardslens-reference.vtu DESTINATION references)
+if(MPI_FOUND)
+  file(COPY test/boxmodels/richards/richardslens-reference.vtu DESTINATION references)
+else(MPI_FOUND)
+  file(COPY test/boxmodels/richards/richardslens-reference-parallel.vtu DESTINATION references)
+endif(MPI_FOUND)
 file(COPY test/boxmodels/richards/test_richards.input DESTINATION test/boxmodels/richards)
 file(COPY test/decoupled/1p/diffusion-reference.vtu DESTINATION references)
 file(COPY test/decoupled/1p/test_1p-reference.vtu DESTINATION references)
@@ -449,7 +453,11 @@ add_test(test_3p3cni   bin/runTest.sh references/kuevette3p3cni-reference.vtu ku
 add_test(test_mpnc     bin/runTest.sh references/obstacle-reference.vtu       obstacle-00010.vtu       test/boxmodels/mpnc/test_mpnc         -Grid.File test/boxmodels/mpnc/grids/obstacle_24x16.dgf         -TimeManager.TEnd 1e4     -TimeManager.DtInitial 250)
 add_test(test_forchheimer1p    bin/runTest.sh references/forchheimer1p-reference.vtp      test_forchheimer1p-00002.vtp       test/boxmodels/mpnc/test_forchheimer1p  -Grid.File test/boxmodels/mpnc/grids/forchheimer1d.dgf)
 add_test(test_forchheimer2p    bin/runTest.sh references/forchheimer2p-reference.vtu      test_forchheimer2p-00009.vtu       test/boxmodels/mpnc/test_forchheimer2p  -Grid.File test/boxmodels/mpnc/grids/obstacle_24x16.dgf)
-add_test(test_richards bin/runTest.sh references/richardslens-reference.vtu   richardslens-00008.vtu   test/boxmodels/richards/test_richards -Grid.File test/boxmodels/richards/grids/richardslens-24x16.dgf -TimeManager.TEnd 3000    -TimeManager.DtInitial 100)
+if(MPI_FOUND)
+  add_test(test_richards bin/runTest.sh references/richardslens-reference.vtu   richardslens-00008.vtu   test/boxmodels/richards/test_richards -Grid.File test/boxmodels/richards/grids/richardslens-24x16.dgf -TimeManager.TEnd 3000    -TimeManager.DtInitial 100)
+else(MPI_FOUND)
+  add_test(test_richards bin/runTest.sh references/richardslens-reference-parallel.vtu s0002-p0000-richardslens-00008.vtu   'mpirun -np 2 test/boxmodels/richards/test_richards' -Grid.File test/boxmodels/richards/grids/richardslens-24x16.dgf -TimeManager.TEnd 3000    -TimeManager.DtInitial 100)
+endif(MPI_FOUND)
 
 add_test(test_diffusion        bin/runTest.sh references/diffusion-reference.vtu             mimeticdiffusion-00001.vtu      test/decoupled/1p/test_diffusion 3)
 add_test(test_dec1p            bin/runTest.sh references/test_1p-reference.vtu               test_1p-00001.vtu               test/decoupled/1p/test_dec1p -ParameterFile test/decoupled/1p/test_1p.input)
diff --git a/test/boxmodels/richards/richardslens-reference-parallel.vtu b/test/boxmodels/richards/richardslens-reference-parallel.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ca49f203be3f60ad9279937754c94c2f91a1588b
--- /dev/null
+++ b/test/boxmodels/richards/richardslens-reference-parallel.vtu
@@ -0,0 +1,424 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+  <UnstructuredGrid>
+    <Piece NumberOfCells="192" NumberOfPoints="221">
+      <PointData Scalars="Sn">
+        <DataArray type="Float32" Name="Sn" NumberOfComponents="1" format="ascii">
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 0.999976 0.999964
+          0.999976 1 1 1 1 1 1 1 1 1 0.997534 0.909917
+          0.905361 0.909917 0.997534 1 1 1 1 1 1 1 0.999 0.771847
+          0.384776 0.36421 0.384776 0.771847 1 1 1 1 1 1 0.999996 0.918527
+          0.286966 0.0671314 0.0634913 0.0671314 0.286966
+        </DataArray>
+        <DataArray type="Float32" Name="Sw" NumberOfComponents="1" format="ascii">
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 9.02056e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 9.02056e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 9.02056e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17
+          9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17 -5.55112e-17
+          -5.55112e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 2.42731e-05 3.56558e-05
+          2.42731e-05 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 0.0024663 0.0900827
+          0.0946392 0.0900827 0.0024663 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 0.00099992 0.228153
+          0.615224 0.63579 0.615224 0.228153 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 9.02056e-17 3.93768e-06 0.0814728
+          0.713034 0.932869 0.936509 0.932869 0.713034
+        </DataArray>
+        <DataArray type="Float32" Name="pn" NumberOfComponents="1" format="ascii">
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+          100000 100000 100000 100000 100000
+        </DataArray>
+        <DataArray type="Float32" Name="pw" NumberOfComponents="1" format="ascii">
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4
+          78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 97470.4 97470.4 97470.4
+          97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 97470.4 97470.4
+          97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 97470.4
+          97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+          97470.4 97470.4 97470.4 97470.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4 78501.4
+          78501.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97471 97471.3
+          97471 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97536.4 99366.6
+          99385.1 99366.6 97536.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97497.1 99586.5
+          99733.6 99739 99733.6 99586.5 97470.4 97470.4 97470.4 97470.4 97470.4 97470.4 97470.5 99323.1
+          99759.4 99835.3 99837.4 99835.3 99759.4
+        </DataArray>
+        <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62
+          21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 2529.62 2529.62 2529.62
+          2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 2529.62 2529.62
+          2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 2529.62
+          2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+          2529.62 2529.62 2529.62 2529.62 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6 21498.6
+          21498.6 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2528.97 2528.67
+          2528.97 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2463.64 633.409
+          614.893 633.409 2463.64 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2502.87 413.55
+          266.37 260.962 266.37 413.55 2529.62 2529.62 2529.62 2529.62 2529.62 2529.62 2529.52 676.884
+          240.611 164.722 162.618 164.722 240.611
+        </DataArray>
+        <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+          1000 1000 1000 1000 1000
+        </DataArray>
+        <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+          1e-10 1e-10 1e-10 1e-10 1e-10
+        </DataArray>
+        <DataArray type="Float32" Name="mobW" NumberOfComponents="1" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0.0411001
+          0.0570506 0.0411001 0 0 0 0 0 0 0 0 0 3.92159
+          146.749 164.984 146.749 3.92159 0 0 0 0 0 0 0 0.0196823
+          249.082 733.751 753.199 733.751 249.082
+        </DataArray>
+        <DataArray type="Float32" Name="mobN" NumberOfComponents="1" format="ascii">
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+          1e+10 1e+10 1e+10 1e+10 1e+10
+        </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.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.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.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.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.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.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.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.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.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.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.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.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.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.4 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="temperature" NumberOfComponents="1" format="ascii">
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+          283.15 283.15 283.15 283.15 283.15
+        </DataArray>
+      </PointData>
+      <CellData Scalars="process rank">
+        <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 0 0 0 0 0 0
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+          0 0 0 0.25 0 0 0 0.25 0 0.25 0.25 0
+          0.5 0 0 0.5 0.25 0 0.75 0 0 0.75 0.25 0
+          1 0 0 1 0.25 0 1.25 0 0 1.25 0.25 0
+          1.5 0 0 1.5 0.25 0 1.75 0 0 1.75 0.25 0
+          2 0 0 2 0.25 0 2.25 0 0 2.25 0.25 0
+          2.5 0 0 2.5 0.25 0 2.75 0 0 2.75 0.25 0
+          3 0 0 3 0.25 0 0 0.5 0 0.25 0.5 0
+          0.5 0.5 0 0.75 0.5 0 1 0.5 0 1.25 0.5 0
+          1.5 0.5 0 1.75 0.5 0 2 0.5 0 2.25 0.5 0
+          2.5 0.5 0 2.75 0.5 0 3 0.5 0 0 0.75 0
+          0.25 0.75 0 0.5 0.75 0 0.75 0.75 0 1 0.75 0
+          1.25 0.75 0 1.5 0.75 0 1.75 0.75 0 2 0.75 0
+          2.25 0.75 0 2.5 0.75 0 2.75 0.75 0 3 0.75 0
+          0 1 0 0.25 1 0 0.5 1 0 0.75 1 0
+          1 1 0 1.25 1 0 1.5 1 0 1.75 1 0
+          2 1 0 2.25 1 0 2.5 1 0 2.75 1 0
+          3 1 0 0 1.25 0 0.25 1.25 0 0.5 1.25 0
+          0.75 1.25 0 1 1.25 0 1.25 1.25 0 1.5 1.25 0
+          1.75 1.25 0 2 1.25 0 2.25 1.25 0 2.5 1.25 0
+          2.75 1.25 0 3 1.25 0 0 1.5 0 0.25 1.5 0
+          0.5 1.5 0 0.75 1.5 0 1 1.5 0 1.25 1.5 0
+          1.5 1.5 0 1.75 1.5 0 2 1.5 0 2.25 1.5 0
+          2.5 1.5 0 2.75 1.5 0 3 1.5 0 0 1.75 0
+          0.25 1.75 0 0.5 1.75 0 0.75 1.75 0 1 1.75 0
+          1.25 1.75 0 1.5 1.75 0 1.75 1.75 0 2 1.75 0
+          2.25 1.75 0 2.5 1.75 0 2.75 1.75 0 3 1.75 0
+          0 2 0 0.25 2 0 0.5 2 0 0.75 2 0
+          1 2 0 1.25 2 0 1.5 2 0 1.75 2 0
+          2 2 0 2.25 2 0 2.5 2 0 2.75 2 0
+          3 2 0 0 2.25 0 0.25 2.25 0 0.5 2.25 0
+          0.75 2.25 0 1 2.25 0 1.25 2.25 0 1.5 2.25 0
+          1.75 2.25 0 2 2.25 0 2.25 2.25 0 2.5 2.25 0
+          2.75 2.25 0 3 2.25 0 0 2.5 0 0.25 2.5 0
+          0.5 2.5 0 0.75 2.5 0 1 2.5 0 1.25 2.5 0
+          1.5 2.5 0 1.75 2.5 0 2 2.5 0 2.25 2.5 0
+          2.5 2.5 0 2.75 2.5 0 3 2.5 0 0 2.75 0
+          0.25 2.75 0 0.5 2.75 0 0.75 2.75 0 1 2.75 0
+          1.25 2.75 0 1.5 2.75 0 1.75 2.75 0 2 2.75 0
+          2.25 2.75 0 2.5 2.75 0 2.75 2.75 0 3 2.75 0
+          0 3 0 0.25 3 0 0.5 3 0 0.75 3 0
+          1 3 0 1.25 3 0 1.5 3 0 1.75 3 0
+          2 3 0 2.25 3 0 2.5 3 0 2.75 3 0
+          3 3 0 0 3.25 0 0.25 3.25 0 0.5 3.25 0
+          0.75 3.25 0 1 3.25 0 1.25 3.25 0 1.5 3.25 0
+          1.75 3.25 0 2 3.25 0 2.25 3.25 0 2.5 3.25 0
+          2.75 3.25 0 3 3.25 0 0 3.5 0 0.25 3.5 0
+          0.5 3.5 0 0.75 3.5 0 1 3.5 0 1.25 3.5 0
+          1.5 3.5 0 1.75 3.5 0 2 3.5 0 2.25 3.5 0
+          2.5 3.5 0 2.75 3.5 0 3 3.5 0 0 3.75 0
+          0.25 3.75 0 0.5 3.75 0 0.75 3.75 0 1 3.75 0
+          1.25 3.75 0 1.5 3.75 0 1.75 3.75 0 2 3.75 0
+          2.25 3.75 0 2.5 3.75 0 2.75 3.75 0 3 3.75 0
+          0 4 0 0.25 4 0 0.5 4 0 0.75 4 0
+          1 4 0 1.25 4 0 1.5 4 0 1.75 4 0
+          2 4 0 2.25 4 0 2.5 4 0 2.75 4 0
+          3 4 0
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+          0 1 3 2 1 4 5 3 4 6 7 5
+          6 8 9 7 8 10 11 9 10 12 13 11
+          12 14 15 13 14 16 17 15 16 18 19 17
+          18 20 21 19 20 22 23 21 22 24 25 23
+          2 3 27 26 3 5 28 27 5 7 29 28
+          7 9 30 29 9 11 31 30 11 13 32 31
+          13 15 33 32 15 17 34 33 17 19 35 34
+          19 21 36 35 21 23 37 36 23 25 38 37
+          26 27 40 39 27 28 41 40 28 29 42 41
+          29 30 43 42 30 31 44 43 31 32 45 44
+          32 33 46 45 33 34 47 46 34 35 48 47
+          35 36 49 48 36 37 50 49 37 38 51 50
+          39 40 53 52 40 41 54 53 41 42 55 54
+          42 43 56 55 43 44 57 56 44 45 58 57
+          45 46 59 58 46 47 60 59 47 48 61 60
+          48 49 62 61 49 50 63 62 50 51 64 63
+          52 53 66 65 53 54 67 66 54 55 68 67
+          55 56 69 68 56 57 70 69 57 58 71 70
+          58 59 72 71 59 60 73 72 60 61 74 73
+          61 62 75 74 62 63 76 75 63 64 77 76
+          65 66 79 78 66 67 80 79 67 68 81 80
+          68 69 82 81 69 70 83 82 70 71 84 83
+          71 72 85 84 72 73 86 85 73 74 87 86
+          74 75 88 87 75 76 89 88 76 77 90 89
+          78 79 92 91 79 80 93 92 80 81 94 93
+          81 82 95 94 82 83 96 95 83 84 97 96
+          84 85 98 97 85 86 99 98 86 87 100 99
+          87 88 101 100 88 89 102 101 89 90 103 102
+          91 92 105 104 92 93 106 105 93 94 107 106
+          94 95 108 107 95 96 109 108 96 97 110 109
+          97 98 111 110 98 99 112 111 99 100 113 112
+          100 101 114 113 101 102 115 114 102 103 116 115
+          104 105 118 117 105 106 119 118 106 107 120 119
+          107 108 121 120 108 109 122 121 109 110 123 122
+          110 111 124 123 111 112 125 124 112 113 126 125
+          113 114 127 126 114 115 128 127 115 116 129 128
+          117 118 131 130 118 119 132 131 119 120 133 132
+          120 121 134 133 121 122 135 134 122 123 136 135
+          123 124 137 136 124 125 138 137 125 126 139 138
+          126 127 140 139 127 128 141 140 128 129 142 141
+          130 131 144 143 131 132 145 144 132 133 146 145
+          133 134 147 146 134 135 148 147 135 136 149 148
+          136 137 150 149 137 138 151 150 138 139 152 151
+          139 140 153 152 140 141 154 153 141 142 155 154
+          143 144 157 156 144 145 158 157 145 146 159 158
+          146 147 160 159 147 148 161 160 148 149 162 161
+          149 150 163 162 150 151 164 163 151 152 165 164
+          152 153 166 165 153 154 167 166 154 155 168 167
+          156 157 170 169 157 158 171 170 158 159 172 171
+          159 160 173 172 160 161 174 173 161 162 175 174
+          162 163 176 175 163 164 177 176 164 165 178 177
+          165 166 179 178 166 167 180 179 167 168 181 180
+          169 170 183 182 170 171 184 183 171 172 185 184
+          172 173 186 185 173 174 187 186 174 175 188 187
+          175 176 189 188 176 177 190 189 177 178 191 190
+          178 179 192 191 179 180 193 192 180 181 194 193
+          182 183 196 195 183 184 197 196 184 185 198 197
+          185 186 199 198 186 187 200 199 187 188 201 200
+          188 189 202 201 189 190 203 202 190 191 204 203
+          191 192 205 204 192 193 206 205 193 194 207 206
+          195 196 209 208 196 197 210 209 197 198 211 210
+          198 199 212 211 199 200 213 212 200 201 214 213
+          201 202 215 214 202 203 216 215 203 204 217 216
+          204 205 218 217 205 206 219 218 206 207 220 219
+        </DataArray>
+        <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+          4 8 12 16 20 24 28 32 36 40 44 48
+          52 56 60 64 68 72 76 80 84 88 92 96
+          100 104 108 112 116 120 124 128 132 136 140 144
+          148 152 156 160 164 168 172 176 180 184 188 192
+          196 200 204 208 212 216 220 224 228 232 236 240
+          244 248 252 256 260 264 268 272 276 280 284 288
+          292 296 300 304 308 312 316 320 324 328 332 336
+          340 344 348 352 356 360 364 368 372 376 380 384
+          388 392 396 400 404 408 412 416 420 424 428 432
+          436 440 444 448 452 456 460 464 468 472 476 480
+          484 488 492 496 500 504 508 512 516 520 524 528
+          532 536 540 544 548 552 556 560 564 568 572 576
+          580 584 588 592 596 600 604 608 612 616 620 624
+          628 632 636 640 644 648 652 656 660 664 668 672
+          676 680 684 688 692 696 700 704 708 712 716 720
+          724 728 732 736 740 744 748 752 756 760 764 768
+        </DataArray>
+        <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+          9 9 9 9 9 9 9 9 9 9 9 9
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>