diff --git a/dumux/porousmediumflow/richards/model.hh b/dumux/porousmediumflow/richards/model.hh
index 787d9bb9792b630385b106ca58a1f286cf4bbfa2..f08140354627b3875a6b18a311c56f245eb15dda 100644
--- a/dumux/porousmediumflow/richards/model.hh
+++ b/dumux/porousmediumflow/richards/model.hh
@@ -112,6 +112,7 @@
#include "volumevariables.hh"
#include "iofields.hh"
#include "localresidual.hh"
+#include "velocityoutput.hh"
namespace Dumux {
@@ -217,6 +218,15 @@ public:
using type = RichardsIOFields<enableWaterDiffusionInAir>;
};
+template<class TypeTag>
+struct VelocityOutput<TypeTag, TTag::Richards>
+{
+ using type = RichardsVelocityOutput<
+ GetPropType<TypeTag, Properties::GridVariables>,
+ GetPropType<TypeTag, Properties::FluxVariables>
+ >;
+};
+
//! The model traits
template<class TypeTag>
struct ModelTraits<TypeTag, TTag::Richards> { using type = RichardsModelTraits<getPropValue<TypeTag, Properties::EnableWaterDiffusionInAir>()>; };
diff --git a/dumux/porousmediumflow/richards/velocityoutput.hh b/dumux/porousmediumflow/richards/velocityoutput.hh
new file mode 100644
index 0000000000000000000000000000000000000000..fb1940837821d449787aee543b49cf543e2f6039
--- /dev/null
+++ b/dumux/porousmediumflow/richards/velocityoutput.hh
@@ -0,0 +1,51 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 3 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file
+ * \ingroup PorousmediumflowModels
+ * \brief Velocity output for the Richards model.
+ */
+
+#ifndef DUMUX_POROUSMEDIUMFLOW_RICHARDS_VELOCITYOUTPUT_HH
+#define DUMUX_POROUSMEDIUMFLOW_RICHARDS_VELOCITYOUTPUT_HH
+
+#include <dumux/porousmediumflow/velocityoutput.hh>
+
+namespace Dumux {
+
+/*!
+ * \ingroup PorousmediumflowModels
+ * \brief Velocity output policy for the Richards model.
+ */
+template<class GridVariables, class FluxVariables>
+class RichardsVelocityOutput : public PorousMediumFlowVelocityOutput<GridVariables, FluxVariables>
+{
+ using ParentType = PorousMediumFlowVelocityOutput<GridVariables, FluxVariables>;
+
+public:
+ using ParentType::ParentType;
+
+ //! Returns the number of phases.
+ int numFluidPhases() const override { return 1; }
+
+};
+
+} // end namespace Dumux
+
+#endif
diff --git a/test/porousmediumflow/richards/nonisothermal/conduction/CMakeLists.txt b/test/porousmediumflow/richards/nonisothermal/conduction/CMakeLists.txt
index 916672671fe8559b14c3cdef25e1f8eede25f106..f9b13e83e4e2fbd46cd45469736e8b7d19ce3b91 100644
--- a/test/porousmediumflow/richards/nonisothermal/conduction/CMakeLists.txt
+++ b/test/porousmediumflow/richards/nonisothermal/conduction/CMakeLists.txt
@@ -6,7 +6,7 @@ dumux_add_test(NAME test_richardsni_conduction_box
SOURCES main.cc
COMPILE_DEFINITIONS TYPETAG=RichardsNIConductionBox
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
- CMD_ARGS --script fuzzy
+ CMD_ARGS --script fuzzy --zeroThreshold {"velocity_liq \(m/s\)":1e-16}
--files ${CMAKE_SOURCE_DIR}/test/references/test_richardsni_conduction_box-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/test_richardsni_conduction_box-00024.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/test_richardsni_conduction_box params.input -Problem.Name test_richardsni_conduction_box")
@@ -16,7 +16,7 @@ dumux_add_test(NAME test_richardsni_conduction_tpfa
SOURCES main.cc
COMPILE_DEFINITIONS TYPETAG=RichardsNIConductionCC
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
- CMD_ARGS --script fuzzy
+ CMD_ARGS --script fuzzy --zeroThreshold {"velocity_liq \(m/s\)":1e-16}
--files ${CMAKE_SOURCE_DIR}/test/references/test_richardsni_conduction_tpfa-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/test_richardsni_conduction_tpfa-00024.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/test_richardsni_conduction_tpfa params.input -Problem.Name test_richardsni_conduction_tpfa")
diff --git a/test/porousmediumflow/richards/nonisothermal/conduction/params.input b/test/porousmediumflow/richards/nonisothermal/conduction/params.input
index 03be824f3d4e5aca0fac79c1d6da022cd2431b91..ca2b3f6cb66e8db2d5ac608c2120481ba392a312 100644
--- a/test/porousmediumflow/richards/nonisothermal/conduction/params.input
+++ b/test/porousmediumflow/richards/nonisothermal/conduction/params.input
@@ -25,3 +25,6 @@ Swr = 0.05
Snr = 0.0
VanGenuchtenAlpha = 0.0037
VanGenuchtenN = 4.7
+
+[Vtk]
+AddVelocity = true
diff --git a/test/references/test_richardsni_conduction_box-reference.vtu b/test/references/test_richardsni_conduction_box-reference.vtu
index 8ee8e4fcd8bcbaf476223a64faf8041f5c6be74f..a822b334c890d16e693f22fb7270a20c62879430 100644
--- a/test/references/test_richardsni_conduction_box-reference.vtu
+++ b/test/references/test_richardsni_conduction_box-reference.vtu
@@ -435,6 +435,109 @@
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="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+ -1.05591e-09 -4.25604e-19 0 -1.04981e-09 7.28418e-19 0 -1.05591e-09 -4.25604e-19 0 -1.04981e-09 7.28418e-19 0
+ -1.03178e-09 3.46407e-18 0 -1.03178e-09 3.46407e-18 0 -1.00245e-09 1.10487e-18 0 -1.00245e-09 1.10487e-18 0
+ -9.65443e-10 -2.95879e-19 0 -9.65443e-10 -2.95879e-19 0 -9.22399e-10 -1.67069e-19 0 -9.22399e-10 -1.67069e-19 0
+ -8.72297e-10 7.22805e-20 0 -8.72297e-10 7.22805e-20 0 -8.16404e-10 1.25613e-18 0 -8.16404e-10 1.25613e-18 0
+ -7.56045e-10 -6.60533e-19 0 -7.56045e-10 -6.60533e-19 0 -6.93172e-10 1.22664e-18 0 -6.93172e-10 1.22664e-18 0
+ -6.33493e-10 1.2504e-19 0 -6.33493e-10 1.2504e-19 0 -5.77581e-10 -2.02157e-18 0 -5.77581e-10 -2.02157e-18 0
+ -5.22003e-10 -4.10922e-19 0 -5.22003e-10 -4.10922e-19 0 -4.67539e-10 1.6689e-18 0 -4.67539e-10 1.6689e-18 0
+ -4.14844e-10 -7.61274e-19 0 -4.14844e-10 -7.61274e-19 0 -3.64441e-10 4.59137e-20 0 -3.64441e-10 4.59137e-20 0
+ -3.16726e-10 -3.52457e-19 0 -3.16726e-10 -3.52457e-19 0 -2.71979e-10 1.29841e-19 0 -2.71979e-10 1.29841e-19 0
+ -2.30374e-10 2.19336e-19 0 -2.30374e-10 2.19336e-19 0 -1.91993e-10 6.48811e-19 0 -1.91993e-10 6.48811e-19 0
+ -1.56842e-10 4.46935e-20 0 -1.56842e-10 4.46935e-20 0 -1.24863e-10 -1.72071e-19 0 -1.24863e-10 -1.72071e-19 0
+ -9.79422e-11 1.20041e-27 0 -9.79422e-11 1.20041e-27 0 -7.63398e-11 6.83938e-19 0 -7.63398e-11 6.83938e-19 0
+ -5.76297e-11 -5.99756e-19 0 -5.76297e-11 -5.99756e-19 0 -4.0978e-11 -1.95299e-18 0 -4.0978e-11 -1.95299e-18 0
+ -2.62252e-11 -2.20641e-18 0 -2.62252e-11 -2.20641e-18 0 -1.32097e-11 -3.38825e-18 0 -1.32097e-11 -3.38825e-18 0
+ -1.77181e-12 -2.19887e-18 0 -1.77178e-12 -2.19887e-18 0 8.24294e-12 -2.5634e-19 0 8.24297e-12 -2.5634e-19 0
+ 1.69816e-11 -8.13922e-20 0 1.69816e-11 -8.13922e-20 0 2.45825e-11 -3.37094e-19 0 2.45825e-11 -3.37094e-19 0
+ 3.11742e-11 -1.34722e-18 0 3.11742e-11 -1.34722e-18 0 3.68747e-11 5.04838e-19 0 3.68747e-11 5.04838e-19 0
+ 4.17917e-11 1.09018e-18 0 4.17917e-11 1.09018e-18 0 4.60225e-11 3.28061e-18 0 4.60225e-11 3.28061e-18 0
+ 4.96544e-11 3.02416e-18 0 4.96544e-11 3.02416e-18 0 5.27656e-11 3.35889e-18 0 5.27656e-11 3.35889e-18 0
+ 5.54254e-11 2.51835e-18 0 5.54253e-11 2.51835e-18 0 5.76948e-11 1.59157e-18 0 5.76948e-11 1.59157e-18 0
+ 5.96278e-11 1.61905e-19 0 5.96278e-11 1.61905e-19 0 6.12713e-11 1.67182e-18 0 6.12714e-11 1.67182e-18 0
+ 6.26668e-11 1.84519e-18 0 6.26668e-11 1.84519e-18 0 6.38497e-11 8.44503e-19 0 6.38497e-11 8.44503e-19 0
+ 6.48511e-11 -5.89606e-21 0 6.4851e-11 -5.89606e-21 0 6.56976e-11 -6.64845e-19 0 6.56976e-11 -6.64845e-19 0
+ 6.64124e-11 3.29447e-19 0 6.64124e-11 3.29447e-19 0 6.70153e-11 5.08873e-19 0 6.70153e-11 5.08873e-19 0
+ 6.75231e-11 -5.89423e-21 0 6.75231e-11 -5.89423e-21 0 6.79505e-11 2.34698e-18 0 6.79505e-11 2.34698e-18 0
+ 6.83098e-11 3.02332e-18 0 6.83098e-11 3.02332e-18 0 6.86116e-11 9.99881e-19 0 6.86116e-11 9.99881e-19 0
+ 6.88649e-11 6.7639e-19 0 6.88649e-11 6.7639e-19 0 6.90773e-11 1.61727e-19 0 6.90773e-11 1.61727e-19 0
+ 6.92552e-11 6.70468e-19 0 6.92552e-11 6.70468e-19 0 6.94042e-11 1.01158e-18 0 6.94042e-11 1.01158e-18 0
+ 6.95288e-11 9.99778e-19 0 6.95288e-11 9.99778e-19 0 6.96331e-11 1.17916e-18 0 6.96331e-11 1.17916e-18 0
+ 6.97201e-11 9.99757e-19 0 6.97202e-11 9.99757e-19 0 6.97929e-11 1.17325e-18 0 6.97929e-11 1.17325e-18 0
+ 6.98536e-11 -4.96925e-19 0 6.98536e-11 -4.96925e-19 0 6.99042e-11 -5.08707e-19 0 6.99042e-11 -5.08707e-19 0
+ 6.99464e-11 -9.99732e-19 0 6.99464e-11 -9.99732e-19 0 6.99816e-11 -1.34672e-18 0 6.99816e-11 -1.34672e-18 0
+ 7.00109e-11 -1.50843e-18 0 7.00109e-11 -1.50843e-18 0 7.00353e-11 -1.33493e-18 0 7.00353e-11 -1.33493e-18 0
+ 7.00556e-11 -1.17911e-18 0 7.00556e-11 -1.17911e-18 0 7.00725e-11 -8.32116e-19 0 7.00725e-11 -8.32116e-19 0
+ 7.00865e-11 1.33492e-18 0 7.00866e-11 1.33492e-18 0 7.00982e-11 1.34081e-18 0 7.00982e-11 1.34081e-18 0
+ 7.01079e-11 5.08696e-19 0 7.01079e-11 5.08696e-19 0 7.0116e-11 -3.41094e-19 0 7.0116e-11 -3.41094e-19 0
+ 7.01227e-11 -4.9691e-19 0 7.01227e-11 -4.9691e-19 0 7.01283e-11 -6.76295e-19 0 7.01283e-11 -6.76295e-19 0
+ 7.01329e-11 5.02802e-19 0 7.01329e-11 5.02802e-19 0 7.01367e-11 -8.32111e-19 0 7.01367e-11 -8.32111e-19 0
+ 7.01399e-11 -5.08694e-19 0 7.01399e-11 -5.08694e-19 0 7.01425e-11 5.89244e-21 0 7.01426e-11 5.89244e-21 0
+ 7.01447e-11 1.33491e-18 0 7.01448e-11 1.33491e-18 0 7.01466e-11 2.51401e-18 0 7.01466e-11 2.51401e-18 0
+ 7.01481e-11 2.84921e-18 0 7.01481e-11 2.84921e-18 0 7.01493e-11 3.86659e-18 0 7.01493e-11 3.86659e-18 0
+ 7.01504e-11 2.83742e-18 0 7.01504e-11 2.83742e-18 0 7.01512e-11 8.38002e-19 0 7.01512e-11 8.38002e-19 0
+ 7.01519e-11 6.70401e-19 0 7.01519e-11 6.70401e-19 0 7.01525e-11 3.35201e-19 0 7.01525e-11 3.35201e-19 0
+ 7.0153e-11 -1.49662e-18 0 7.0153e-11 -1.49662e-18 0 7.01534e-11 -1.52019e-18 0 7.01534e-11 -1.52019e-18 0
+ 7.01537e-11 5.02801e-19 0 7.01537e-11 5.02801e-19 0 7.0154e-11 -3.35201e-19 0 7.0154e-11 -3.35201e-19 0
+ 7.01542e-11 1.0056e-18 0 7.01542e-11 1.0056e-18 0 7.01544e-11 1.18499e-18 0 7.01544e-11 1.18499e-18 0
+ 7.01546e-11 4.91016e-19 0 7.01546e-11 4.91016e-19 0 7.01547e-11 5.02801e-19 0 7.01547e-11 5.02801e-19 0
+ 7.01548e-11 -6.70401e-19 0 7.01548e-11 -6.70401e-19 0 7.01549e-11 -1.55815e-19 0 7.01549e-11 -1.55815e-19 0
+ 7.0155e-11 3.23416e-19 0 7.0155e-11 3.23416e-19 0 7.0155e-11 3.35201e-19 0 7.0155e-11 3.35201e-19 0
+ 7.01551e-11 -8.38001e-19 0 7.01551e-11 -8.38001e-19 0 7.01551e-11 -1.676e-19 0 7.01551e-11 -1.676e-19 0
+ 7.01552e-11 1.79385e-19 0 7.01552e-11 1.79385e-19 0 7.01552e-11 9.93817e-19 0 7.01552e-11 9.93817e-19 0
+ 7.01552e-11 1.3408e-18 0 7.01552e-11 1.3408e-18 0 7.01552e-11 3.8944e-27 0 7.01552e-11 3.8944e-27 0
+ 7.01553e-11 3.8944e-27 0 7.01553e-11 3.8944e-27 0 7.01552e-11 1.35259e-18 0 7.01553e-11 1.35259e-18 0
+ 7.01553e-11 3.00502e-18 0 7.01553e-11 3.00502e-18 0 7.01553e-11 2.1788e-18 0 7.01553e-11 2.1788e-18 0
+ 7.01553e-11 2.6816e-18 0 7.01553e-11 2.6816e-18 0 7.01553e-11 3.51961e-18 0 7.01553e-11 3.51961e-18 0
+ 7.01553e-11 1.18499e-18 0 7.01553e-11 1.18499e-18 0 7.01553e-11 1.55816e-19 0 7.01553e-11 1.55816e-19 0
+ 7.01553e-11 -2.0112e-18 0 7.01553e-11 -2.0112e-18 0 7.01553e-11 -1.3408e-18 0 7.01553e-11 -1.3408e-18 0
+ 7.01553e-11 6.70401e-19 0 7.01553e-11 6.70401e-19 0 7.01553e-11 1.52019e-18 0 7.01553e-11 1.52019e-18 0
+ 7.01553e-11 8.26217e-19 0 7.01553e-11 8.26217e-19 0 7.01553e-11 3.8944e-27 0 7.01553e-11 3.8944e-27 0
+ 7.01553e-11 1.0056e-18 0 7.01553e-11 1.0056e-18 0 7.01553e-11 -8.38001e-19 0 7.01553e-11 -8.38001e-19 0
+ 7.01553e-11 -3.23416e-19 0 7.01553e-11 -3.23416e-19 0 7.01553e-11 -3.46985e-19 0 7.01553e-11 -3.46985e-19 0
+ 7.01553e-11 -1.1732e-18 0 7.01553e-11 -1.1732e-18 0 7.01553e-11 -1.1732e-18 0 7.01553e-11 -1.1732e-18 0
+ 7.01553e-11 1.676e-19 0 7.01553e-11 1.676e-19 0 7.01553e-11 3.46985e-19 0 7.01553e-11 3.46985e-19 0
+ 7.01553e-11 1.32902e-18 0 7.01554e-11 1.32902e-18 0 7.01553e-11 2.3464e-18 0 7.01553e-11 2.3464e-18 0
+ 7.01553e-11 1.1732e-18 0 7.01553e-11 1.1732e-18 0 7.01553e-11 2.514e-18 0 7.01553e-11 2.514e-18 0
+ 7.01553e-11 4.03419e-18 0 7.01553e-11 4.03419e-18 0 7.01553e-11 4.68102e-18 0 7.01553e-11 4.68102e-18 0
+ 7.01553e-11 2.514e-18 0 7.01553e-11 2.514e-18 0 7.01553e-11 2.1788e-18 0 7.01553e-11 2.1788e-18 0
+ 7.01553e-11 1.0056e-18 0 7.01553e-11 1.0056e-18 0 7.01553e-11 1.35259e-18 0 7.01553e-11 1.35259e-18 0
+ 7.01553e-11 -1.17846e-20 0 7.01553e-11 -1.17846e-20 0 7.01553e-11 -1.0056e-18 0 7.01553e-11 -1.0056e-18 0
+ 7.01553e-11 -1.5084e-18 0 7.01553e-11 -1.5084e-18 0 7.01553e-11 -3.1844e-18 0 7.01553e-11 -3.1844e-18 0
+ 7.01553e-11 -2.33462e-18 0 7.01553e-11 -2.33462e-18 0 7.01553e-11 -3.53139e-18 0 7.01553e-11 -3.53139e-18 0
+ 7.01553e-11 -2.6816e-18 0 7.01553e-11 -2.6816e-18 0 7.01553e-11 -3.0168e-18 0 7.01553e-11 -3.0168e-18 0
+ 7.01553e-11 -4.02241e-18 0 7.01553e-11 -4.02241e-18 0 7.01553e-11 -3.17262e-18 0 7.01553e-11 -3.17262e-18 0
+ 7.01553e-11 -2.52579e-18 0 7.01553e-11 -2.52579e-18 0 7.01553e-11 -2.514e-18 0 7.01553e-11 -2.514e-18 0
+ 7.01553e-11 -3.1844e-18 0 7.01553e-11 -3.1844e-18 0 7.01553e-11 -3.68721e-18 0 7.01553e-11 -3.68721e-18 0
+ 7.01553e-11 -3.00502e-18 0 7.01553e-11 -3.00502e-18 0 7.01553e-11 -2.02299e-18 0 7.01553e-11 -2.02299e-18 0
+ 7.01553e-11 -2.6816e-18 0 7.01553e-11 -2.6816e-18 0 7.01553e-11 -2.0112e-18 0 7.01553e-11 -2.0112e-18 0
+ 7.01553e-11 -3.0168e-18 0 7.01553e-11 -3.0168e-18 0 7.01553e-11 -1.49662e-18 0 7.01553e-11 -1.49662e-18 0
+ 7.01553e-11 -1.17846e-20 0 7.01553e-11 -1.17846e-20 0 7.01553e-11 1.676e-19 0 7.01553e-11 1.676e-19 0
+ 7.01554e-11 -2.0112e-18 0 7.01553e-11 -2.0112e-18 0 7.01553e-11 -1.676e-18 0 7.01553e-11 -1.676e-18 0
+ 7.01553e-11 -2.62268e-18 0 7.01553e-11 -2.62268e-18 0 7.01553e-11 -1.53197e-18 0 7.01553e-11 -1.53197e-18 0
+ 7.01553e-11 -1.65243e-18 0 7.01553e-11 -1.65243e-18 0 7.01553e-11 -1.1732e-18 0 7.01553e-11 -1.1732e-18 0
+ 7.01553e-11 -2.03477e-18 0 7.01553e-11 -2.03477e-18 0 7.01553e-11 -1.65243e-18 0 7.01553e-11 -1.65243e-18 0
+ 7.01553e-11 -3.54317e-18 0 7.01553e-11 -3.54317e-18 0 7.01553e-11 -5.00444e-18 0 7.01553e-11 -5.00444e-18 0
+ 7.01553e-11 -3.85481e-18 0 7.01553e-11 -3.85481e-18 0 7.01553e-11 -2.70517e-18 0 7.01553e-11 -2.70517e-18 0
+ 7.01553e-11 -2.82564e-18 0 7.01553e-11 -2.82564e-18 0 7.01553e-11 -1.02917e-18 0 7.01553e-11 -1.02917e-18 0
+ 7.01553e-11 -1.14963e-18 0 7.01553e-11 -1.14963e-18 0 7.01553e-11 -1.1732e-18 0 7.01553e-11 -1.1732e-18 0
+ 7.01553e-11 3.11632e-19 0 7.01553e-11 3.11632e-19 0 7.01553e-11 1.36437e-18 0 7.01553e-11 1.36437e-18 0
+ 7.01553e-11 -1.19677e-18 0 7.01553e-11 -1.19677e-18 0 7.01553e-11 -2.15523e-18 0 7.01553e-11 -2.15523e-18 0
+ 7.01553e-11 -2.6816e-18 0 7.01553e-11 -2.6816e-18 0 7.01553e-11 -2.87277e-18 0 7.01553e-11 -2.87277e-18 0
+ 7.01553e-11 -2.82564e-18 0 7.01553e-11 -2.82564e-18 0 7.01553e-11 -1.36437e-18 0 7.01553e-11 -1.36437e-18 0
+ 7.01553e-11 -8.14432e-19 0 7.01553e-11 -8.14432e-19 0 7.01553e-11 -3.352e-19 0 7.01553e-11 -3.352e-19 0
+ 7.01553e-11 -2.35689e-20 0 7.01553e-11 -2.35689e-20 0 7.01553e-11 -1.65243e-18 0 7.01553e-11 -1.65243e-18 0
+ 7.01553e-11 -1.91169e-19 0 7.01553e-11 -1.91169e-19 0 7.01553e-11 -3.11632e-19 0 7.01553e-11 -3.11632e-19 0
+ 7.01553e-11 1.55776e-26 0 7.01553e-11 1.55776e-26 0 7.01554e-11 -3.20797e-18 0 7.01553e-11 -3.20797e-18 0
+ 7.01553e-11 -3.32844e-18 0 7.01553e-11 -3.32844e-18 0 7.01553e-11 -2.8492e-18 0 7.01553e-11 -2.8492e-18 0
+ 7.01553e-11 -1.02917e-18 0 7.01553e-11 -1.02917e-18 0 7.01553e-11 -9.82033e-19 0 7.01553e-11 -9.82033e-19 0
+ 7.01553e-11 -1.91169e-19 0 7.01553e-11 -1.91169e-19 0 7.01553e-11 1.69957e-18 0 7.01553e-11 1.69957e-18 0
+ 7.01553e-11 1.5084e-18 0 7.01553e-11 1.5084e-18 0 7.01553e-11 4.79232e-19 0 7.01553e-11 4.79232e-19 0
+ 7.01553e-11 -1.82003e-18 0 7.01553e-11 -1.82003e-18 0 7.01553e-11 -3.58769e-19 0 7.01553e-11 -3.58769e-19 0
+ 7.01553e-11 -4.71376e-20 0 7.01553e-11 -4.71376e-20 0
+ </DataArray>
<DataArray type="Float32" Name="temperatureExact" NumberOfComponents="1" format="ascii">
300 299.373 300 299.373 298.749 298.749 298.133 298.133 297.529 297.529 296.939 296.939
296.368 296.368 295.817 295.817 295.289 295.289 294.787 294.787 294.313 294.313 293.866 293.866
diff --git a/test/references/test_richardsni_conduction_tpfa-reference.vtu b/test/references/test_richardsni_conduction_tpfa-reference.vtu
index 5187408907e2504100bd71f1f235121ecc3f33e8..2b6f63a245a0b0f50a3eb2a4db6b9c85e954b47a 100644
--- a/test/references/test_richardsni_conduction_tpfa-reference.vtu
+++ b/test/references/test_richardsni_conduction_tpfa-reference.vtu
@@ -231,6 +231,58 @@
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1
</DataArray>
+ <DataArray type="Float32" Name="velocity_liq (m/s)" NumberOfComponents="3" format="ascii">
+ -1.0497e-09 0 0 -1.0375e-09 0 0 -1.01365e-09 0 0 -9.81329e-10 0 0
+ -9.42115e-10 0 0 -8.95233e-10 0 0 -8.41905e-10 0 0 -7.83441e-10 0 0
+ -7.21184e-10 0 0 -6.61476e-10 0 0 -6.05608e-10 0 0 -5.49651e-10 0 0
+ -4.9445e-10 0 0 -4.40723e-10 0 0 -3.89059e-10 0 0 -3.39916e-10 0 0
+ -2.9363e-10 0 0 -2.50423e-10 0 0 -2.1042e-10 0 0 -1.73663e-10 0 0
+ -1.40119e-10 0 0 -1.10503e-10 0 0 -8.64195e-11 0 0 -6.66413e-11 0 0
+ -4.90002e-11 0 0 -3.33388e-11 0 0 -1.94955e-11 0 0 -7.30858e-12 0 0
+ 3.37956e-12 0 0 1.27202e-11 0 0 2.08563e-11 0 0 2.79214e-11 0 0
+ 3.4039e-11 0 0 3.93219e-11 0 0 4.38725e-11 0 0 4.7783e-11 0 0
+ 5.1136e-11 0 0 5.4005e-11 0 0 5.64552e-11 0 0 5.85437e-11 0 0
+ 6.03209e-11 0 0 6.18308e-11 0 0 6.31115e-11 0 0 6.41964e-11 0 0
+ 6.51141e-11 0 0 6.58894e-11 0 0 6.65436e-11 0 0 6.7095e-11 0 0
+ 6.75592e-11 0 0 6.79497e-11 0 0 6.82778e-11 0 0 6.85533e-11 0 0
+ 6.87844e-11 0 0 6.8978e-11 0 0 6.91402e-11 0 0 6.9276e-11 0 0
+ 6.93895e-11 0 0 6.94844e-11 0 0 6.95636e-11 0 0 6.96298e-11 0 0
+ 6.9685e-11 0 0 6.97311e-11 0 0 6.97694e-11 0 0 6.98014e-11 0 0
+ 6.98281e-11 0 0 6.98502e-11 0 0 6.98687e-11 0 0 6.9884e-11 0 0
+ 6.98967e-11 0 0 6.99074e-11 0 0 6.99162e-11 0 0 6.99235e-11 0 0
+ 6.99296e-11 0 0 6.99346e-11 0 0 6.99388e-11 0 0 6.99423e-11 0 0
+ 6.99452e-11 0 0 6.99476e-11 0 0 6.99495e-11 0 0 6.99512e-11 0 0
+ 6.99526e-11 0 0 6.99537e-11 0 0 6.99546e-11 0 0 6.99554e-11 0 0
+ 6.9956e-11 0 0 6.99566e-11 0 0 6.9957e-11 0 0 6.99574e-11 0 0
+ 6.99577e-11 0 0 6.99579e-11 0 0 6.99581e-11 0 0 6.99583e-11 0 0
+ 6.99584e-11 0 0 6.99586e-11 0 0 6.99587e-11 0 0 6.99587e-11 0 0
+ 6.99588e-11 0 0 6.99589e-11 0 0 6.99589e-11 0 0 6.99589e-11 0 0
+ 6.9959e-11 0 0 6.9959e-11 0 0 6.9959e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0 6.99591e-11 0 0
+ </DataArray>
<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