diff --git a/NumLib/DOF/MeshComponentMap.cpp b/NumLib/DOF/MeshComponentMap.cpp
index d3a08afbfb3aad535393608f569c2dc953c4d616..89730484aeea87048d0b6db74414d233147d4fcf 100644
--- a/NumLib/DOF/MeshComponentMap.cpp
+++ b/NumLib/DOF/MeshComponentMap.cpp
@@ -24,9 +24,6 @@ namespace NumLib
{
using namespace detail;
-GlobalIndexType const MeshComponentMap::nop =
- std::numeric_limits<GlobalIndexType>::max();
-
#ifdef USE_PETSC
MeshComponentMap::MeshComponentMap(
std::vector<MeshLib::MeshSubset> const& components, ComponentOrder order)
diff --git a/NumLib/DOF/MeshComponentMap.h b/NumLib/DOF/MeshComponentMap.h
index c7c2cb619e29f1e59976e454cd6acd4a5dbd7151..02d647f50ffcebdf092bd00f3b1479a6a8b9e766 100644
--- a/NumLib/DOF/MeshComponentMap.h
+++ b/NumLib/DOF/MeshComponentMap.h
@@ -137,7 +137,8 @@ public:
/// A value returned if no global index was found for the requested
/// location/component. The value is implementation dependent.
- static NUMLIB_EXPORT GlobalIndexType const nop;
+ static constexpr NUMLIB_EXPORT GlobalIndexType const nop =
+ std::numeric_limits<GlobalIndexType>::max();
#ifndef NDEBUG
const detail::ComponentGlobalIndexDict& getDictionary() const
diff --git a/ProcessLib/Output/AddProcessDataToMesh.cpp b/ProcessLib/Output/AddProcessDataToMesh.cpp
index 13a57345f440c4227a50cfd9d77b52b313f97b1f..d9e23519af667c2ce3c66bd0265ddcf590fdec7b 100644
--- a/ProcessLib/Output/AddProcessDataToMesh.cpp
+++ b/ProcessLib/Output/AddProcessDataToMesh.cpp
@@ -326,6 +326,14 @@ static std::set<std::string> addPrimaryVariablesToMesh(
// per node ordering of components
auto const out_index = node_id * n_components + component_id;
+ // request for index of linear quantities at higher order nodes
+ // results in returning NumLib::MeshComponentMap::nop
+ if (in_index == NumLib::MeshComponentMap::nop)
+ {
+ output_data[out_index] = 0;
+ continue;
+ }
+
output_data[out_index] = x_copy[in_index];
}
}
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3cfb02b2615e4a5dc2f9bc4be24cdee948f15f6a
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="152" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="212" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="272" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="336" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="516" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="576" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="640" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="704" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="0" offset="796" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="856" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="916" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="1008" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="1068" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="1160" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1220" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15" offset="1284" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15" offset="1372" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1440" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1504" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="1568" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="1660" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="1720" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1900" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2024" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2116" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAHwAAAAAAAAA=eF5LSwOCRUUOaVCaAQ2kockTUkeufnT16HwA5KYkhw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAEgAAAAAAAAA=eF5LSwOCRUUOaYOcBgAG90W7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..75daf8f00382c4fe70d7e1577768ab78273fd306
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5.000000031e-05" RangeMax="5.000000031e-05" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.03056133107" RangeMax="0.028248973988" offset="212" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="352" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190561" RangeMax="0.050476190468" offset="468" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="840" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-5.9428570479e-10" RangeMax="5.942856941e-10" offset="960" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-5.9428570479e-10" RangeMax="5.942856941e-10" offset="1132" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142719e-12" offset="1380" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="2.7527042325e-26" RangeMax="3.7142857143e-11" offset="1604" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2328" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2420" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2588" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2844" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2936" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="3112" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3204" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="5.4782669632e-16" RangeMax="0.58466805513" offset="3264" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15583048" offset="3956" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15583048" offset="4116" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.00099999999703" RangeMax="0.0010000000045" offset="4352" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4724" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="5096" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="5188" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="5248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5428" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5552" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5644" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5rt7/69pEZlz0DGPDtYcAKLuxuhKjbD+H/sMau7sVu7OIPdqPKv8GhDg4aULkXdgMAKW8UaA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAASAAAAAAAAAA=eF5LNt3Q6/N2jj0DEDR80rBlwAa0mG27Xt0zrfCevx+sTlYLu7odHHuxil/Ug6g/828PmHYWw64fBj5csEHhh3DbAgA88RREAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANgAAAAAAAAA=eF4zdpZ5/ciMy54BDPj2MGAFF3brQtTth/B/WGNX92I3dvEHu1Hl3+BQBwcNqNwbuwH3kRJaAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9gAAAAAAAAA=eF7jEQSBxv1Gay9Krb3YZc8ABk17lha98eR4PQ3Cb6jdo59Tupbn32KovMUe82pD111ci6D8CpuEp4GVySyT7fkg5tlLXnGUveKYC5WP2PPo04yHW79PgZpXtedLduRLLY85UPlmm1ubDonGf1wKU2/zh2P+/b72JijfwKbnp3y5zZQJUH6cTZLqpGfMN1ZC+UV7pMH2ToTyA2xY4rM3GOe1Qflle9rmBm2KujcDan+zzeNlH+6kiy+Aykft4TRh9XhYBXWfQ4ONfId83cYny6HyDTaZcUKNN0Jg4ROw59C9fRJtiTPg7o2rXBBYa91nDwDqkWKHAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYICAl4oVBSnZLPsY0EDtZJuKlOwve2H8Jo/tj7/fZcRQx/Jn2nKxhf9s0cXPuWQdqW/N2jNjuW8W0Hy7ucZ7VgPV7WXAAQKbkm4BzbeD8bWLOOOA9tsCAFsPJC4=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAlwAAAAAAAAA=eF5jYICAl4oVBSnZLPscJl2VXyX8Z298cH1rSvaTvbWTbSpSsr/shSpjaPLY/vj7XcZ9adbpnEB1tmZbxbXWqTDvY/kzbbnYwn+2EPqd7VzjPauB9N5zLllH6luz9sxY7psFNN8OIv5vr+emYwFA823PTuK5DzTP7o6sNtjcwKakWyD+F5EIJaC5dtpFnHFA+20BIP9Kmw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhQAAAAAAAAA=eF5jYEAFO9K3f+JXnrUXxldnyA17O2URnG98oitUznwJnL/r9v2v3rMWw/noYFOfYTmT+kK4vBqaeVOLIzaILZwH52fzuX4sTO+F8939pn1YxFkP58ei8WNy8kpTshHu2Y3m/nw089DdrydYs/XD5SlwviYafyOa+yeiuXcnmv8BPjVObw==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA/gEAAAAAAAA=eF7bvuNDaUr2EVux33W5cq927mKAAusE+U0FLces/AxkP99r2GE7adcXO8avK+Hyf3MXV3203WyVznT+072GFbYL/1z3MCmbaAWTL5yrdpDz1jKrBh/l7y9DJ9j27GM3k3jVAJfPi9m/Su920a7dYPkFtgkH93f5W2yGmy/m8HjaW6Nlu65A3Lf35pynOg/zY+HyXdf+zDhzp2TXy2hboP4Ntmuensh7t/Q03PwJP50vX7ONsvKHuG+vw9yCjJ7W2XD9nWD7m3bJuOkCzd9i27REQj94+XS4/k9vMqxag7utjGcuBZp/YK+Pc03rsi0FcHn21EfPTtZN28WgedpsecKhvfl8fHeyi6Ph8juvH2p9kZNpxQSRt20P19V/V70Rbn8jT8r+3NedYH7D9ElWf2L09h5ZOQuu/4m6pvaO0pZdupDw31sRorfyu8EEuH6mrNC8b7ELrTQg7rON2skgof9uAVz/rE96/llT91nthYT/3jO7xepOHkqG678RPE3dNWy9laehSWHktD17rzcmngzkr4TrT9mRXeXAkLjLBiJvayrwIcLbtwcu7/e3I+ZsQeUulhhw+O/dt3ragWeKbnDzuU/YsT09NHGXGiR89249dOpTiEUdXP7oHQ8JW7O5uxh8wfG/9238mlTH6/lwedGLJwt6xCZbAQDvKgHBAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXQAAAAAAAAA=eF5jYACBD/bsb9xnzpz50p4BDWTFcBkbGz+Gi7e/k0tLS3uGoY5ldfPsmTN3YogzgskH9vf9H8yaOfOk/YKu/jNnzrzBUAcDH5UOAe27DJfXdn8NdNdNewC2pCLTAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAngAAAAAAAAA=eF5jYACBD/bsb9xnzpz50r5i6ZQzZ848sTfTawXyH9pnxXAZGxs/tmeAgvZ3cmlpac/s+bZeA9LX7N1zN/////+5Pcvq5tkzZ+6E0ivtA466A+Xf2TOCdT2wv+//YNbMmSftF3T1A81/Yz9t8Vcgfcd+g1CgibHxYfuQ3a+NjI1f239UOgS077J9gepMoL4L9trur4HuuGkPAPHLS0A=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYQAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL6Dqd0LplVD5d/aMYPkH9jPA4ieh5r+xPwem79gbgc0/DLXnNZR/2f7ffxA4DzXvpj0AbEpLzA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA5gEAAAAAAAA=eF47+R8E7tt3bcu9vS33tP1iKC1aG84tpHXFRuZGv5Gx8WV7f6aX/uylB+yjobR7he4RQ+uDNp+mBgHlD9tPbtgCFN9g3wWlxbJl7HaWbbHhehw0Y+bMlfa/Q26KxH+cav8KSi9P5v104FHDHm+w/E777WDxpfZ/oPSz2G31wg837/kDcd/+mxB3wWn/SsY3m9Y17VHeVQ3Uf9L+S9gfoL6t9hNg9GL+4uv16TaPIO7b3w1x1/46KL0HbP+EPSpbz/77//+8/Rzd9FvbcnfbT4TSPsHn6hsrZth4icrPnDnz5v6k6A7R+I9H4fRVjYv5TKpL9szr6j9z5syd/dKc1T+61p/YLwWljSWtQ/cdKbOZBZG3Z4KI2wtC6QiXa9+4zKbvYQCChkPbbUD0AfPdYJqhbJ1NvctWlpSgSXs4IOG/PwQS7vsDoPTH1ZO27ZVab2MFcZ99AsRd9pFQuuDfvwNOyadt0iDhv98qFBzu+02g9N4lV8u/8e+xWVZ2PTUt7dr+G+t2gMMJRnvxHZv2vChvz1yIvP01iLj9JSj9IcmY0fpE2x41SPjvz4KEO5zm6+Y+3HR+/h59SPjuh4bv/ulQ+pJz5r/jGqv2nIfE//7JYHct3d8PpcPXP5r2SnyhDQBdKF+tAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVgAAAAAAAAA=eF5LSwOCRUUOh7Zz9lUAaQY0cHytpEcdUDwNqk716be9NVjUcU9qLcnBIj5T0XdmLVDcWf6aRhGQvnlmk0shFnUw8Nz5u0Y5kvz1G/skqoF8AEpEKWo=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAkAAAAAAAAAA=eF5LSwOCRUUOh7Zz9lUA6TJPjuZaIL0hffMHEH18raRHHZBOg6pTffptbw2Qjue7rVAJpEv8GJZVAWnuSa0lOUDa4k9pbiaQXjCz52o2kJ6p6DsTZI6z/DWNIiB988wml0KQur+skSBzlt9Y65cPpP2c4l6WAOnnzt81yoH01I3PVoD412/sk6gG0gCeiU+CAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9gAAAAAAAAA=eF475Kd9xTcpYL+Xd+8n0wcKe3jPcl8G8Q9e5Jrych3jntnXQy6B+JP/5lissrKwKZmoAOa3PX1++k+n2p7OLVZg/teyq76Fz79a74GaxwAFu2Q7wfKLjnoa3E9R2zMDah7v3W+bbUtZ9lxJ3gvmz3v05nD6DdU91YcOge1/bF74v8X2p/UuNPP2Q/mrIeptqqHuYWC5d8evr82aG+r+/H1eXl9KXu5uhpqnfuXd94k/rG0aoOp/rfwIpA1sLMyiwPIRLx5dnmj9fbcTlL8Jav5uqPsLTObxLXz9xfo21L3rwPJfd/dB/R/3cEXax3QBGwCDyKo7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9AAAAAAAAAA=eF5jXvXxkm9SwP7vk+0+mz5Q2PN3JYTvxZD8fVcg054fUL6/9o2lnn/Nbb5A+YZyj40uzVHb8xPKdznWkXxB8os1TD8Dw4F1Vo6Cu2Dyq5zuxa5VUIebt7nsnjyzI/MemHmlrO/ZVnqp7fkN5e/TjK8zXvfT+g+U/3jp7CMLtHusGaDuXbKi/nbGZlWbXzD73C5fe7pmJVz9r8qWG65Zr3azQNVP+ZAgOPmytc1/qPzSkDmC9asMbGDu0Us886vP4ftuJqj6fUUZG/nmq9rA3JOxaCnfurCP1jD+1lWrVpk++L4bZt+erZe2ZzQK2QAAOJO7aA==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..05e064742b13e942fdcbd99bc7ea356d9be45141
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.029894096593" RangeMax="0.029723375654" offset="208" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="348" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190476" RangeMax="0.050476190476" offset="464" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="824" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-3.3646614802e-18" RangeMax="3.9588653941e-18" offset="944" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-3.3646614802e-18" RangeMax="3.9588653941e-18" offset="1116" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142857e-12" offset="1376" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="7.505834107e-26" RangeMax="3.7142857143e-11" offset="1600" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2308" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2400" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2532" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2704" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2796" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2968" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3060" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.0356994338e-15" RangeMax="0.58466805513" offset="3120" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15053296" offset="3776" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15053296" offset="3928" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4152" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4492" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="4860" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="4952" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="5012" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5192" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5316" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5408" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANAAAAAAAAAA=eF6TcJZ5/ciMy54BAvYwYAUfrPUg6vZDBRqwq7thjV38wW4oowGNxgHg6qHghzUAw9gPAA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAARgAAAAAAAAA=eF5bHFP882z+PHsGEGCUt2XABuRlbfnkb97pmzNvP4jboKe4F6u6w7rYxU+JQcQTeCC0FxN2dTAQ8sEGhX9Oci8AxoET8A==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANAAAAAAAAAA=eF6TcJZ5/ciMy54BDP7tZsAKvljrQdTthwo0YFd3wxq7+AOYuQ1oNA7wAM0d36wBXh4P/w==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA7QAAAAAAAAA=eF7jEASBxv16ay9Krb3YZc8ABgk2rdZyXhyvp0H5VTbuy6av5fm3GC5f8NbWdRfXIig/Zc+5k1mVySyT7fkg5tmLXXGUveKYC5V32LOq6cjDrd+nQPgNtXv2fRd5qeUxBypfvEc3/qNo/MelUD5Dw/LokPt97U1w/sNFv8pspkyA8mNstl5Y8oz5xkqY/TYSYHsnQvkZNo5n2TcY57XB+WpfyjdF3ZsB5SfvWW7y5E66+AK4+TUi69wfVkHdx5Cz57G4d93GJ8uh/AabB3a7G26EwMIndY9EWyIQwcxrtKnZkBdYa91nDwAz8GGoAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYAAAAAAAAAA=eF5jYICAAxJtj69cd7VhQAMOF24sPesYBBfnUJnNLWFgiqFu2k5xy1X7vDDEpV7OFZzK4b8nRPfd9dbdZjbnJCW43q10wlAHAz/FSwMXywTC5bUmepv/53DcAwCnGx6BAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAogAAAAAAAAA=eF5jYICAAxJtj69cd7Vxmrl/xeVed5sZd/1eSLwRtXG4cGPpWccgG6gyBg6V2dwSBqY2J+bNy3Y3MLRJurqy5IO5g820neKWq/Z5QWkrm3OSElzvVhrZSL2cKziVw39PiO676627zaDiTjZKe0tUTnF47FHIaJ9kxuds4y55vuyhg4vNT/HSwMUygTZ+iwuOt5x2s9Ga6G3+n8NxDwDO00FrAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhwAAAAAAAAA=eF5jYEAFr5bIf+ZXnrUXxneuvxH2dsoiOH9CJVeYnPkSOH/z7vBv3rMWw/nowO5ITzmT+kK4vAOaeR6iLzeILZwH57fObP5YmN4L51+ZvPXDIs56OP8tGv/Wjg+lKdkI97xBc38/mnn9aO7PNt+79cPlKXB+BhrfBs39bmju3YjmfwAOJVQ2AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA8QEAAAAAAAA=eF7bvuNDaUr2EduoT2stV83YtIsBCqbn1VfMcTthtWHV28/3GnbYnlY7Pill2SK4vMnhrg97WbdbdYLlV9gq/CiudYpcaAWTlyp9f1eKfYKV9PHP31+GTrD9V+zU9doPIf9HdGf/zjnTdhmB5RfYimWo/bGxWAY333KK8oPCRUutrkPct3eRgPWvWReT4PLNpksYm/+k7OIC699gu8ypdopH6Cm4+TXOPwMEJRutlkDct5fpfnOQfkw/XH+07RFVN+cpu+6Czd9ia6J4X/dO3AS4/ic6S5j/qU6zkgSbf2Dvis0riw58jYHLT1lzckGL1NRdn+w/mS1POLT3Rnonn5dhOFy+5xKfgcnOBquvEHlbN6Fk7mTmNXD7D7QLNjZxLAfzDxxabPVt6+P10gyL4foZvfx+qx7cZNUHCf+9DesPr97fjLD/r2ejtif/UisRiPtsszQf3OPvWASX/yfjpqhiv8uKDRL+e2cZ2B6P/lsMlz8WmBPw4PlCq+g5rwojp+3Zuy71uxa7bS9c/qN1fdq8BWlWgRB5203MIs5MMp1weS2WkHUn73RZKUDCf+9lhfuylxfUw+VzRTznFF9bb/UYEr57zyq8eMzzuwGh32d+x7O9LVaskPjfe6Dm6TJ7Vxe4/JHtphZOmvOtAOXxBgk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAQQAAAAAAAAA=eF5jYACBD/YzweClPQMaMAaDx3Dx9DQQeIahbgZY/04McUYw+QBq/kn7s2dA4A2GOhiA2HcZLg/Rd9MeAPoQIv0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXgAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL+xngNXvhNIrofLv7BnB8g+g5p+Emv/G/hyYvmNvBDb/MNSe11D6MtTcC1B9N+0BGpVF0Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXgAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL+xngNXvhNIrofLv7BnB8g+g5p+Emv/G/hyYvmNvBDb/MNSe11D6MtTcC1B9N+0BGpVF0Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAAyQEAAAAAAAA=eF678B8E7tvP2pZ7e1vuafs1UHpxsqei8d6rNs7GIHDZPoi9NIC99IB9FJRmW+vsdWfuERs9sPxh+2Sw+Ab7QCid880zY92WOTaHZoLASvsJ8R9F4z9OtQ+B0k2/vJwv2i/ZcxAsv9P+XxxIfKm9ajyEtue48WJD/GYbVgYQeLD/AcRdcNrhruIvidNFezaA9Z+03wTWt9U+FEobLFzXGzm/z8YW4r79KRB37Y+B0pzFrNXmqov2cILNv2D/BWzubvs3UNonecHP2quLbZaBzb+5XwFs7lE43VeY/vj700V7Tp8BgTv7f3StB6ITcHrzbk/XFZ97bc5B5O2fQ8Ttv0LpAx07zoct2LIHbP2lozYgKmEKhGb4vtNmUjbPloLTB2y0IeG/3xMS7vs9oHSceN+LF7c32SyCuM9eE+IueyUorVC17Q1j43GbVkj47/8JDt+p+xkg4b+//9jhmPscG2xS00Dg2n5oPO43gNLTXJzf+zqX2MRB5O31ofGsBaVj9yib6fHNsIGG/341SLjv14DSZ7oZn+5v2mPDDQnf/a8g4br/LZQWCVlRwD55kk0HJP73s0HifT8HlL6+aNuy5c9X2wAAhqxO9g==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAUQAAAAAAAAA=eF5LSwOCRUUOc3+vOpoOpBnQQKnV/tYMoHgaVN2L2Rn2GVjU/Xf5zIlN/21P+6kg9e/f75oPkpfeYKeETR0M1NQfkUI2f8rzhfkgPgD1QypeAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAhQAAAAAAAAA=eF5LSwOCRUUOc3+vOpoOpPVeT07KANIaTvW9ILrUan8riE6DqnsxO8MexO8wNnYF0YdOdzKB6P8unzlB+jdd1dkBUufQHXQERN/2tJ8Kkn//ftd8kLz0BjslEL1jRkETSLxc6noIiH/n6pdiEF1Tf0QKJL5tu/1dEH/K84X5ID4AFNpQaQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA3QAAAAAAAAA=eF5jXPXxkm9SwH5+RyuGyZ/09nxfCeHP4Lf7mpyjZvMXyg+RfrGl0c/a5guUv/HRm8PpN0T3/IHyM66vqyg8ZmTzG8pnmLkhMGrRlt0/oPzL9q0v7uY423yF8g95fG91/8e75z2UvxJing3MvG/3lQzOn3y4+xeUDxH/sPsflM8KdreBDUw+SV9vj/JdNbj9IT8mFLLG/LBmgPpvGfeKpZoFDntg+g+B7WOFu//T2xuPd897sRvGh/hPdc9PKF+kXyLqxWZOGxifAQw+7oCpz+HidpmwWc0GAOpfv6c=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF5jXPXxkm9SwH4BRyuGyZ/09nxfCeFbMvQIiexXs/kL5UtnBCmuvmdl8wXK72d0VC34KLrnD5RvZzXNpPOzoc1vKJ9j5obAqEVbdv+A8o1u/y6rTnO2+Qrl+wmbGCpU8+x5D+WfcBCV7hMUs4GZd3L/PsNctee7f0H5HxdGBJfu/Lj7H5Rv4qJw6Fi5gQ1M/oYX9+PXmRpw+8+8mabNa/3HmgHqv7LCym0bohz2wPQ7qV1INtzBAHd/6ssFWsx5z3fD+EzbXlyTbVHd8xPKv/Du5r1XazhsYPwlNt0WZwTWWsPUm/va+D7epm4DAAtitE0=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..72639127d7e75eeb0bc99bdc02296fec3bca6739
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="144" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="204" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="264" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="328" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="0" offset="448" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="508" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="568" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="632" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="696" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="0" offset="792" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="852" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="912" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="1008" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="1068" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="1160" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1220" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15" offset="1284" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15" offset="1372" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1440" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1504" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="1568" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="1660" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="1720" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1900" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2024" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2116" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PszQ3SU1JTAUALuIFKA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAHwAAAAAAAAA=eF5LSwOCRUUOaVCaAQ2kockTUkeufnT16HwA5KYkhw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAEgAAAAAAAAA=eF5LSwOCRUUOaYOcBgAG90W7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7c69115d040786574795e48bc71558d48138e572
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5.000000031e-05" RangeMax="5.000000031e-05" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.03056133107" RangeMax="0.028248973988" offset="204" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="344" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190561" RangeMax="0.050476190468" offset="464" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="836" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-5.9428564254e-10" RangeMax="5.9428574849e-10" offset="956" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-5.9428564254e-10" RangeMax="5.9428574849e-10" offset="1128" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142719e-12" offset="1384" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="1.0423184196e-25" RangeMax="3.7142857143e-11" offset="1608" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2328" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2424" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2592" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2848" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2944" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="3116" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3208" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.1047392924e-15" RangeMax="0.58466805513" offset="3268" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15583048" offset="3968" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15583048" offset="4128" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.00099999999703" RangeMax="0.0010000000045" offset="4368" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4736" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="5104" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="5196" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="5256" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5436" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5560" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5652" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PszQ3SU1JTAUALuIFKA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOQAAAAAAAAA=eF5rt7/69pEZlz0DGLDtYcAKPlg3QtTth/B/WGNXd2E3dvEPUPELUH1PcOiHgQNo8hd2AwC86xNTAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAARwAAAAAAAAA=eF4LD93Q6/N2jj0DGNy2YcAGjOX3Ttt1z7TCe/5+MJ+H3Rarui2f9mAVVxYFq3f4a74XzFf5gd0eGIhUQDX/gNFeAN8AE70=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOAAAAAAAAAA=eF4zcZZ5/ciMy54BDFj3MGAFH6y1Ier2Q/i/rLGru7Abu/gHmHgDhHqBQz8MHEAz58FuAHSjEUY=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9QAAAAAAAAA=eF7jEASBxv16ay9Krb3YZc8ABil7Fhe98eR4PQ3Cb6jdo5dTupbn32KofMEes2pD111ci6D8GpuEp4GVySyT7QUg5tnLXHGUveKYC5XPsHn0acbDrd+nQM3r3/M5O/KllsccqLyBzY1Nh0TjPy6F8hP2/OWYf7+vvQnKD7Dp+SlfbjNlApQfsydeddIz5hsrofw6GwmwvRNh9u1hi8/eYJzXBtO/p2Fu0KaoezOg/OI9D5Z9uJMuvgDKz7VhN2H1eFgFdR9Dso10h3zdxifLofxmm5w4ocYbIbDwybc5cG+fRFsizLwWm+TKBYG11n32ANpRYfk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYAAAAAAAAAA=eF5jYIAAmetTfFKyWfYxoIFZS3hZUrK/7IXxL14o2PL9LiOGOustO8+JLfxniy7eKCNy62PJbJsXS1d0As23u3ujvw2obi8DDpD+7hDHj7uMdjD+Hdv0M0D7bQF9zyZFAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAnQAAAAAAAAA=eF5jYIAAmetTfFKyWfbJLS7ftlL4z16ju5bTkrOf7J21hJclJfvLXqgyhosXCrZ8v8u4L7T3W94q4T+23y43MaxTYd5nvWXnObGF/2wh9Dvbuzf624D03kYZkVsfS2bbvFi6ohNovh1E/N9eF1Zrq9TsJ7bMB5p+AM2zO5r47jjI3PR3hzh+APkap0w9geba3bFNPwO03xYAf8hPgw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhgAAAAAAAAA=eF5jYEAFO9O3f+JXnrUXxtdkyA17O2URnG92oitUznwJnL/n9v2v3rMWw/noYEufYTmT+kK4vDqaedOKIzaILZwH5+fyuX4sTO+F8338pn1YxFkP5yeh8eNy8kpTshHu2Yvm/hI080zQ3K8vWLP1w+UpcL42Gn8TmvunoLl3F5r/AVkVTpk=AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA+QEAAAAAAAA=eF47uONDaUr2Edtd7hdvJF4+sIsBClaU+25sPHncKsZA9vO9hh22004HagZlLofLrxI5FcHTu9cqi+n8p3sNK2x/JbrLp5ROsILJ/8u5f6bsR4NVvY/y95ehE2y/bc7wdfs+ES4PARXmm8DyC2yZrhjUK+1fDzff4deLzISflVa3Ie7byym1fcMmzmlw/c0vmjsKLPN2ccbYAvVvsF2+fakxZ8dpuPwDN88pYqlJVm4Q9+3tqs1pfnOqH27+nV//7tZqx+1SctMFmr/FtmjH+0YO9ilw/U+953c829tiZThzKdD8A3vLLdYJuIq1wuWLrWKtFhbO3sWsedpsecKhvds3u4bdOd4Pl59XEVeiLtFtxQ6Rt80Ns1K8sHcb3P6Vm5beihddBOY7fNm563hGcejMrQvh+me52qsc/rPBygUS/nu9b/gxC1wrhsvP9J3h7/dymZUxxH22fm1skSLls+Dy7NE2ng08u60eQcJ/b/mTQqM9J7IR4S+2s3/nnGVWHoYmhZHT9uz97zntyWmXqXB5zw+fdlyZUWZlB5G3fbVEceXLAz1w+X7BV4GHtvRafY8Gh//eeSe951i/6IXLL/97PNQ/cNIuWUj47t06O3dJdQIi/CYzCN+W2Vayawsk/vfO41JREM8xgMuLvfp4pqarYRcAKDL5yw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXQAAAAAAAAA=eF5jYACBD/bsb9xnzpz50p4BDWTGcBkbGz+Gi7e/k0tLS3uGoY55dfPsmTN3YogzgskH9g/8H8yaOfOk/fyu/jNnzrzBUAcDH5UOAe27DJfXdn8NdNdNewC1xCLRAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAngAAAAAAAAA=eF5jYACBD/bsb9xnzpz50r5i6ZQzZ848sTfVawXyH9pnxnAZGxs/tmeAgvZ3cmlpac/s+bZeA9LX7N1zN/////+5PfPq5tkzZ+6E0ivtA466A+Xf2TOCdT2wf+D/YNbMmSft53f1A81/Yz9t8Vcgfcd+o1CgibHxYfvg3a+NjI1f239UOgS077J9gepMoL4L9trur4HuuGkPAPALSzw=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYACBD/YzZoLAS/szYPAEyn9ob2gMAo/tGaAgDQye2aeA6Wv2/8Hguf00sPqdUHqlPUTdO6i+B/bTweInoea/gdJ3oOYftjcC06+h/Mv2f8Dmnoe646Y9AMFVUag=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA6gEAAAAAAAA=eF67+R8E7tv3bMu9vS33tP1qKL1poqaQtu9VG/Ub/UbGxpfto5he+rOXHrBPhdJh+hN+tJSdsvkyNQgof9h+UsMWoPgG+24o/ebZk2dfUvpshB8HzZg5c6W9ROhNkfiPU+2/h0BoBjBosXQCy++0PwAWX2r/E0p7vNNo7HBqt2EBq3uw/w7EXfsfQmmrDYwlTwWq9wTtqgbqP2kvE/4HqG+r/e4wCJ0hbcfJUFlgcxPivv1VEHftL4bSsWzrX6azZO/R2nr23///5+3n66bf2pa7234alL4evKKAffIkG1dR+ZkzZ97cHxPdIRr/8SicdnhuzVExecWexV39Z86cubNfnLP6R9f6E/ulofRT9p+3wp1n2CyEyNszQsTtRaH0oYuJh+cabtgD8p1CykUwzVA9G0K/3LZnaVTCvOU5+2zEIOG/Hxr++6Oh9JfVIbGpLVtsPCHus8+BuMs+DUqbyq/JfmFxwmYeJPz3R0HCf388lO745eV80X6LzbKy66lpadf2X1i3AxxOV6D0ug1B00NsWmzmQeTtr0LE7S9D6W95LAtEZs2y4YeE/34fSLjvD4bS+/bnP/ZMX7BHDhK++5sg4bq/DUo7HdguwdfTtCcEEv/71cDuWrpfBUrPcf/y6c/83j0ABH5TbQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVgAAAAAAAAA=eF5LSwOCRUUOh7dz9lUAaQY0cHStpEcdUDwNqk716be9NVjUcU9qLcnBIj5D0XdmLVDcWf6aRhGQvnlmk0shFnUw8Nz5u0Y5kvy1G/skqoF8AEmEKWc=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAkQAAAAAAAAA=eF5LSwOCRUUOh7dz9lUA6VJPjuZaIL0uffMHEH10raRHHZBOg6pTffptbw2QjuO7rVAJpEv8GJZVAWnuSa0lOUDa4k9pbiaQXjCz52o2kJ6h6DsTZI6z/DWNIiB988wml0Igbf6XNRJkzvIba/3ygbSfU9zLEiD93Pm7RjmQnrrx2QoQ/9qNfRLVQBoAm3FPeg==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF4746d9xTcpYP8EV7t7Lwy89wic5b4M4nMJ/Jmi2yGxZ+H1kEsg/oZHbw6n32DdUzRRAcxngIK+LVZgfpmA+Os1ri579kDNg8nvl+0Ey3s/UveW3aNvMwtq3ofIhP8LjI5ZX0/eC+avAZsvalN36BDYfruOe2/cHzDZ7EIz7xiUX3SdYYmwkPCeOqh7nBM6ZoivtrHhg7r/zdHrD6UTFG3aoOZ5HHQN6lSzsmmAqveSfrGl0U/bxsosCix/7/7H3ZfUWPfYQfkbwe5RtdkHdT+D7sSoLoWzu+9A3QuR/7p7EtT/+6TDnOdOUNwDABS+mjw=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8wAAAAAAAAA=eF7jW/Xxkm9SwP7tmfn3Xhh472GC8qfsv/vyHYvknv8rIfzGWw7pdwQZ93yA8s+2LuHMnfHQ+g+Uv5zhxJS8N857fkH5DAwMf8z0VHcxQM1juvftvEmsgc13qLz8j/XsVlM+W/+A8ntUZ+QHzxCz+Qfls3z/657JyGoDk3/2e+VHLs6D1hxQ82Iur7o6467gHpj9Wspn9ioK2drA3FtV5FU6O0bFhhOq/oVQ6e5jB6zg8l4fEirUNuvawNx7JbDrTKvVn90w964xv3SjJFgN7h6HN6dXf5Xbtxum/ponD9NGrkdw9TEspgIXtKT3AABXfLS2AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..21d6d0e2dd58b981d5c0fefb921760ffe69357f0
--- /dev/null
+++ b/Tests/Data/TH2M/HM/empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.029894096594" RangeMax="0.029723375655" offset="204" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="352" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190476" RangeMax="0.050476190476" offset="468" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="836" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="2.6164164567e-17" offset="956" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="2.6164164567e-17" offset="1128" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142857e-12" offset="1388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="5.1677569392e-26" RangeMax="3.7142857143e-11" offset="1612" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2308" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2404" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2532" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2692" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2788" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2948" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3040" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="3.9998130556e-16" RangeMax="0.58466805513" offset="3100" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15053296" offset="3740" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15053296" offset="3892" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4116" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4436" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="4792" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="4884" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="4944" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5124" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5248" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5340" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PszQ3SU1JTAUALuIFKA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOQAAAAAAAAA=eF4zcJZ5/ciMy54BDDj2MGAFH3ZrQdTth/Af7MaujqEBu/APawj9AqrvgzV2dTBwAF2+AQCdchA3AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAATQAAAAAAAAA=eF5zry77eTZ/nj0DEDQYhO1lwAZWBtim9t260zdn3n4Qd0GAJ1Z1B74722ITL3iQAhafYZ4Kpi+0xGK3Bwoc9piiyL+4m2ELABr8Gsk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANQAAAAAAAAA=eF7Td5Z5/ciMy54BDNj2MGAFH3ZrQ9Tth/Gxq7tgjV38D1T8BVTfJxzqYADdnAvWANhAEXE=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF7jEASBxv1Gay9Krb3YZc8ABiU2bdZyXhyvp0H5rXvclk1fy/NvMZSfZJP/1tZ1F9ciKL/O5tzJrMpklsn2bBDz7AWuOMpeccyFyjfYrGg68nDr9ylQfrnNvu8iL7U85sDs26MV/1E0/uNSKD9hz6rokPt97U1Qfr7Ng0W/ymymTIDy0222XVjyjPnGSii/eY8E2N6JUH7GHoez7BuM89pg5tlofinfFHVvBoSvMMlmhcmTO+niC+D+qRFZ5/6wCua+mj2PxL3rNj5ZDjfvpt3uhhshsPDJt5FpS5RoS4Sax9C5p25DXmCtdZ89AItoYvI=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYICANQdKQ78dybZhQANKzHd8WsJL4OIb4h7fqHb0wFBn94ZNcl1OJYb4mQeyF24x1NpobZH+82NCrQ3T14dbOrgyMdTBgI5SU6KPVh1cfsox/k6rA1U2AGjeILk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAogAAAAAAAAA=eF5jYICANQdKQ78dybb5kffj+Ob4dJvvH2v69NdU2Cgx3/FpCS+xgSpj2BD3+Ea1o4dNwhePUufSGhsdn+tLV71KsLF7wya5LqcSSmfaMH19uKWDK9LmzAPZC7cYam20tkj/+TGhFiqeabPh8rQ1XQ+qbUwCJnZdrqu2uX5Lu3V/epaNjlJToo9WnY3263PrFt2ttZlyjL/T6kCVDQAMNEn4AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhwAAAAAAAAA=eF5jYEAFb5bIf+ZXnrUXxneuvxH2dsoiOH9iJVeYnPkSOH/L7vBv3rMWw/nowP5ITzmT+kK4vCOaed6iLzeILZwH53fNbP5YmN4L59+evPXDIs56OP8TGv/Wjg+lKdkI97xFc/9ENPMmoLk/z3zv1g+Xp8D5mWh8GzT3u6K5dyOa/wEllVRXAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA6QEAAAAAAAA=eF47vONDaUr2Eduqi0ZtMh8O7GKAAiufBmlb2eNWa1e9/XyvYYctS4yi0221tXB5/hStGUKle63mgeVX2B5ZYb8inWmBFUze1md+x7O9U6xeHfv8/WXoBNtLtRzfViQg5J3B8i1WPMdB8gtsHaZEHVyVsRpuPovify7JJRus7kDctzf98i67c5oT4PrjQsr/LJYqsfoLNn+D7a6tyv5CP07C5XWqb4WGPi60mgFx3163vmmS2++Wws3XkTu5IHb11F03wOZvsT3c861r9ukmuH4waJi3XQLsvgN7l4hEFah/74fLn5qwPanhV+uur/afzJYnHNqbqBIXFNHdDZd/nXP/TNmPgF0/IPK2HpVGPgtDd8Htj1AROCEYuADMd/i8aFfaw6BPyhxL4PpfrFugdMd0gdUaSPjvVREPt2vRb4Hrn3HO8NzX5MVWQhD32eoFvH0mapcCl68I27ZwZf9+Kw2w/IS9n1SMxYP6euHmx/z6d7dWu84qbc6rwshpe/au2Flh2vMWIQ8BBeYhEHlbt/YiTV3ZVrj5fJZfn/tFV1tBwmfD3vnyrzmi+ifC9T/5bTTt9P8lVvch4bt31uUQw8W+9XD5E4E5AQ+eL7Rih8T/XvvHHH83/euAy98r9H+9tr9hFwDCfglUAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAPQAAAAAAAAA=eF5jYACBD/YzweClPQMaMAaDx3DxNDB4hqFuFlj/TgxxCHhgD5E/aX8GDN7gUAez7zJcHuKum/YA+hgi/Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVwAAAAAAAAA=eF5jYACBD/YzweCl/RkweGI/C8x/aG8MBo/tGaAgDQye2UPoa1DxF1D1O6H0Sqj8O6j8A6j4Saj5b6D0HXsTsPmH7Y3A9GuofZeh+i5A3XXTHgAahUXRAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVwAAAAAAAAA=eF5jYACBD/YzweCl/RkweGI/C8x/aG8MBo/tGaAgDQye2UPoa1DxF/YzwOp3QumVUPl3UPkHUPNPQs1/A6Xv2BuBzT8MpV9D7bsM1XcBqu+mPQAYhUXKAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAAvwEAAAAAAAA=eF67/x8E7ttP2pZ7e1vuafvNUFrS0OXxs01XbOyMQeCyvRt7aQB76QH7YCgd+dp4Jk/iKRtbsPxh+2yw+Ab7OCitFbKigH3yIpvSmSCw0t4z/qNo/Mep9nxQWgcsP8mmCSy/035XHEh8qf0bKH1yl47Z7bN7bVgYQODB/kcQd+1/AqXfTtkfeIWhyWYeWP9J+/lgc7faO0Hp4oLyl2dZGmxMIO7bHwJx1/5AKB3R9Jd5XfriPYxg8y/YPwWbu9v+AZSGWLtm53Kw+Tf3y4DNPbpfDkrXzag60b1u8p4LZ0Dgzv5vXet/dK0/sf8HlD737MmzLylxe85B5O1vQcTtP0BpXq0H79Qs1+4B2/N9PYQW2mwDpitybbYtkcnN5lhr4wIJ//1JkHDfHw+l9+/bImYbu9FmFsR99mIQd9mLQ2nvn/Lzop3P2GyHhP9+Tki471eG0o9Z179MZ+m2yUoDgWv7DaHhZAylGSAOsUyAyNurQ+NZG0over7velV9h81iSPjvF4KE+35pKK2+g2UH57lNNhyQ8N3/GBKu+59B6UnHDsfc59gAi//9XyDxvv8vlD7LY5rhubJ3DwBNXUclAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAUQAAAAAAAAA=eF5LSwOCRUUOc36vOpoOpBnQQJHV/tYMoHgaVN2z2Rn2GVjU/XX5zIlN/01P+6kg9e/e75oPkpfeYKeETR0MVNcfkUI2f9LzhfkgPgDwYypQAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAhQAAAAAAAAA=eF5LSwOCRUUOc36vOpoOpHVeT07KANJqTvW9ILrIan8riE6Dqns2O8MexG8zNnYF0QdPdzKB6L8unzlB+jdd1dkBUufQHXQERN/0tJ8Kkn/3ftd8kLz0BjslEL1tRkETSLxM6noIiH/n6pdiEF1df0QKJL5lu/1dEH/S84X5ID4AC7pQTw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAzQAAAAAAAAA=eF7jW/Xxkm9SwH4V6565H/fy2fxdCeFHZM37rrpBcs9PKP/b2xuPd8/jsIHxGcH6FPb8g/I9IyUKk9+LwPkMUPAHylfJyEqVXv5t93co/4TH91b3f7x7vkD5yx69OZx+w9SGAeqeX8vPbWYP/Gz9B808jlXI/B87YPaZrYiO8BI22fMfym9yW9qaxypiwwZVv6vg97qrz533/ILKrwLbx7oHZj4rWB0DnL8FLC8K50v1S0S92MwJ9z/E/g0WMPvXvuk/vmO+sg0AmeGyAg==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA6gAAAAAAAAA=eF7jW/Xxkm9SwH5V6565H/fy2fxdCeFzLJoYW10svec3lP8mcdqlDyoCNj+h/CzmdTkaXkp7GKD6p1lK93JkSuz5A5VnAIMGCxh/RcEHfoE3X62/Q/kur9+7HKrn2PMFyi8Nulk0jcXc5j+Uf0Fm1RW7Z8+tf8HUpz1TXVh+djcH1D4tL2Yvs+/ndv+Dynd4iq/eLWO6B6af55bp0lfFAjZsUPV+5f3SR1Rd9sDMK1VOjmYuYoW7917N/qPepxH8azPMNM9sFIfzJQ6alc6fwQ33/6XZTNOdfl6B2/+ez2fh2aeqNgDdrLMOAQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj b/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
index b5dbc9a9266ec53f93d1fae7e342281348d841d7..c5d87805a7fb721aaa148ab45ed599352386eab2 100644
--- a/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
@@ -228,39 +228,95 @@
</time_stepping>
</process>
</processes>
- <output>
- <type>VTK</type>
- <prefix>flow_fully_saturated_gas</prefix>
- <timesteps>
- <pair>
- <repeat>10000</repeat>
- <each_steps>1</each_steps>
- </pair>
- </timesteps>
- <variables>
- <variable>gas_pressure</variable>
- <variable>gas_pressure_interpolated</variable>
- <variable>capillary_pressure</variable>
- <variable>capillary_pressure_interpolated</variable>
- <variable>temperature</variable>
- <variable>temperature_interpolated</variable>
- <variable>displacement</variable>
- <variable>sigma</variable>
- <variable>epsilon</variable>
- <variable>velocity_gas</variable>
- <variable>velocity_liquid</variable>
- <variable>liquid_density</variable>
- <variable>gas_density</variable>
- <variable>porosity</variable>
- <variable>saturation</variable>
- <variable>GasMassFlowRate</variable>
- <variable>HeatFlowRate</variable>
- <variable>LiquidMassFlowRate</variable>
- <variable>NodalForces</variable>
- <variable>liquid_pressure_interpolated</variable>
- </variables>
- <suffix>_ts_{:timestep}_t_{:time}</suffix>
- </output>
+ <outputs>
+ <output>
+ <type>VTK</type>
+ <prefix>flow_fully_saturated_gas</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <timesteps>
+ <pair>
+ <repeat>10000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>gas_pressure</variable>
+ <variable>gas_pressure_interpolated</variable>
+ <variable>capillary_pressure</variable>
+ <variable>capillary_pressure_interpolated</variable>
+ <variable>temperature</variable>
+ <variable>temperature_interpolated</variable>
+ <variable>displacement</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity_gas</variable>
+ <variable>velocity_liquid</variable>
+ <variable>liquid_density</variable>
+ <variable>gas_density</variable>
+ <variable>porosity</variable>
+ <variable>saturation</variable>
+ <variable>GasMassFlowRate</variable>
+ <variable>HeatFlowRate</variable>
+ <variable>LiquidMassFlowRate</variable>
+ <variable>NodalForces</variable>
+ <variable>liquid_pressure_interpolated</variable>
+ </variables>
+ </output>
+ <!-- test output configuration with sub-domain mesh and output of linear primary variables -->
+ <output>
+ <type>VTK</type>
+ <prefix>flow_fully_saturated_gas_{:meshname}</prefix>
+ <timesteps>
+ <pair>
+ <repeat>10000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <meshes>
+ <mesh>square_1x1_quad8_1e2_bottom</mesh>
+ </meshes>
+ <variables>
+ <variable>gas_pressure</variable>
+ <variable>gas_pressure_interpolated</variable>
+ <variable>capillary_pressure</variable>
+ <variable>capillary_pressure_interpolated</variable>
+ <variable>temperature</variable>
+ <variable>temperature_interpolated</variable>
+ <variable>displacement</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity_gas</variable>
+ <variable>velocity_liquid</variable>
+ <variable>liquid_density</variable>
+ <variable>gas_density</variable>
+ <variable>porosity</variable>
+ <variable>saturation</variable>
+ <variable>GasMassFlowRate</variable>
+ <variable>HeatFlowRate</variable>
+ <variable>LiquidMassFlowRate</variable>
+ <variable>NodalForces</variable>
+ <variable>liquid_pressure_interpolated</variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ <!-- test output configuration with sub-domain mesh and empty variables section -->
+ <output>
+ <type>VTK</type>
+ <prefix>empty_variables_flow_fully_saturated_gas_{:meshname}</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <timesteps>
+ <pair>
+ <repeat>10000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <meshes>
+ <mesh>square_1x1_quad8_1e2_bottom</mesh>
+ </meshes>
+ <variables>
+ </variables>
+ </output>
+ </outputs>
</time_loop>
<parameters>
<!-- Mechanics -->
@@ -502,5 +558,29 @@
<absolute_tolerance>1e-15</absolute_tolerance>
<relative_tolerance>0</relative_tolerance>
</vtkdiff>
+ <vtkdiff>
+ <regex>flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_.*.vtu</regex>
+ <field>gas_pressure</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>2e-15</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_.*.vtu</regex>
+ <field>capillary_pressure</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>1e-15</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_.*.vtu</regex>
+ <field>gas_pressure</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>2e-15</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_.*.vtu</regex>
+ <field>capillary_pressure</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>1e-15</relative_tolerance>
+ </vtkdiff>
</test_definition>
</OpenGeoSysProject>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_newton.xml b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_newton.xml
index 756248341599945e46150d6ff7002f043d7a5375..bc4c813665babcc7c3abfa26d9ffed47f653c64d 100644
--- a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_newton.xml
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_newton.xml
@@ -1,6 +1,14 @@
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProjectDiff base_file="flow_fully_saturated_gas.prj">
<remove msel="/*/*/process/jacobian_assembler"/>
- <replace sel="/*/time_loop/output/prefix/text()">flow_fully_saturated_gas_newton</replace>
+ <replace msel="/*/time_loop/outputs/output[1]/prefix/text()">flow_fully_saturated_gas_newton</replace>
+ <replace msel="/*/time_loop/outputs/output[2]/prefix/text()">flow_fully_saturated_gas_{:meshname}_newton</replace>
+ <replace msel="/*/time_loop/outputs/output[3]/prefix/text()">empty_variables_flow_fully_saturated_gas_{:meshname}_newton</replace>
+
<replace msel="/*/test_definition/vtkdiff/regex/text()">flow_fully_saturated_gas_newton_ts_.*.vtu</replace>
+
+ <replace msel="/*/test_definition/vtkdiff[last()-3]/regex/text()">flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_.*.vtu</replace>
+ <replace msel="/*/test_definition/vtkdiff[last()-2]/regex/text()">flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_.*.vtu</replace>
+ <replace msel="/*/test_definition/vtkdiff[last()-1]/regex/text()">empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_.*.vtu</replace>
+ <replace msel="/*/test_definition/vtkdiff[last()]/regex/text()">empty_variables_flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_.*.vtu</replace>
</OpenGeoSysProjectDiff>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b3f385d6d2581d5c5b9e72b58ae2dcd5dab6f5a3
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_0_t_0.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="152" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="212" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="272" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="0" offset="336" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="396" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="520" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="584" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="0" offset="676" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="736" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="796" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="888" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="948" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="1040" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1100" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15" offset="1164" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15" offset="1252" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1320" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1384" />
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="1448" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1628" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1752" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1844" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAHwAAAAAAAAA=eF5LSwOCRUUOaVCaAQ2kockTUkeufnT16HwA5KYkhw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAEgAAAAAAAAA=eF5LSwOCRUUOaYOcBgAG90W7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8635fea0de8fc1a39bb50a1d5f3ee29f3ad8d89d
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_1_t_1.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5.000000031e-05" RangeMax="5.000000031e-05" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.03056133107" RangeMax="0.028248973988" offset="212" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="352" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190561" RangeMax="0.050476190468" offset="468" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-5.9428570479e-10" RangeMax="5.942856941e-10" offset="840" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-5.9428570479e-10" RangeMax="5.942856941e-10" offset="1012" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142719e-12" offset="1260" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="2.7527042325e-26" RangeMax="3.7142857143e-11" offset="1484" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2208" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2300" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2468" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2724" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2816" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2992" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3084" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="5.4782669632e-16" RangeMax="0.58466805513" offset="3144" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15583048" offset="3836" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15583048" offset="3996" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.00099999999703" RangeMax="0.0010000000045" offset="4232" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4604" />
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="4976" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5156" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5280" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5372" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5rt7/69pEZlz0DGPDtYcAKLuxuhKjbD+H/sMau7sVu7OIPdqPKv8GhDg4aULkXdgMAKW8UaA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAASAAAAAAAAAA=eF5LNt3Q6/N2jj0DEDR80rBlwAa0mG27Xt0zrfCevx+sTlYLu7odHHuxil/Ug6g/828PmHYWw64fBj5csEHhh3DbAgA88RREAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANgAAAAAAAAA=eF4zdpZ5/ciMy54BDPj2MGAFF3brQtTth/B/WGNX92I3dvEHu1Hl3+BQBwcNqNwbuwH3kRJaAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9gAAAAAAAAA=eF7jEQSBxv1Gay9Krb3YZc8ABk17lha98eR4PQ3Cb6jdo59Tupbn32KovMUe82pD111ci6D8CpuEp4GVySyT7fkg5tlLXnGUveKYC5WP2PPo04yHW79PgZpXtedLduRLLY85UPlmm1ubDonGf1wKU2/zh2P+/b72JijfwKbnp3y5zZQJUH6cTZLqpGfMN1ZC+UV7pMH2ToTyA2xY4rM3GOe1Qflle9rmBm2KujcDan+zzeNlH+6kiy+Aykft4TRh9XhYBXWfQ4ONfId83cYny6HyDTaZcUKNN0Jg4ROw59C9fRJtiTPg7o2rXBBYa91nDwDqkWKHAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYICAl4oVBSnZLPsY0EDtZJuKlOwve2H8Jo/tj7/fZcRQx/Jn2nKxhf9s0cXPuWQdqW/N2jNjuW8W0Hy7ucZ7VgPV7WXAAQKbkm4BzbeD8bWLOOOA9tsCAFsPJC4=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAlwAAAAAAAAA=eF5jYICAl4oVBSnZLPscJl2VXyX8Z298cH1rSvaTvbWTbSpSsr/shSpjaPLY/vj7XcZ9adbpnEB1tmZbxbXWqTDvY/kzbbnYwn+2EPqd7VzjPauB9N5zLllH6luz9sxY7psFNN8OIv5vr+emYwFA823PTuK5DzTP7o6sNtjcwKakWyD+F5EIJaC5dtpFnHFA+20BIP9Kmw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhQAAAAAAAAA=eF5jYEAFO9K3f+JXnrUXxldnyA17O2URnG98oitUznwJnL/r9v2v3rMWw/noYFOfYTmT+kK4vBqaeVOLIzaILZwH52fzuX4sTO+F8939pn1YxFkP58ei8WNy8kpTshHu2Y3m/nw089DdrydYs/XD5SlwviYafyOa+yeiuXcnmv8BPjVObw==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA/gEAAAAAAAA=eF7bvuNDaUr2EVux33W5cq927mKAAusE+U0FLces/AxkP99r2GE7adcXO8avK+Hyf3MXV3203WyVznT+072GFbYL/1z3MCmbaAWTL5yrdpDz1jKrBh/l7y9DJ9j27GM3k3jVAJfPi9m/Su920a7dYPkFtgkH93f5W2yGmy/m8HjaW6Nlu65A3Lf35pynOg/zY+HyXdf+zDhzp2TXy2hboP4Ntmuensh7t/Q03PwJP50vX7ONsvKHuG+vw9yCjJ7W2XD9nWD7m3bJuOkCzd9i27REQj94+XS4/k9vMqxag7utjGcuBZp/YK+Pc03rsi0FcHn21EfPTtZN28WgedpsecKhvfl8fHeyi6Ph8juvH2p9kZNpxQSRt20P19V/V70Rbn8jT8r+3NedYH7D9ElWf2L09h5ZOQuu/4m6pvaO0pZdupDw31sRorfyu8EEuH6mrNC8b7ELrTQg7rON2skgof9uAVz/rE96/llT91nthYT/3jO7xepOHkqG678RPE3dNWy9laehSWHktD17rzcmngzkr4TrT9mRXeXAkLjLBiJvayrwIcLbtwcu7/e3I+ZsQeUulhhw+O/dt3ragWeKbnDzuU/YsT09NHGXGiR89249dOpTiEUdXP7oHQ8JW7O5uxh8wfG/9238mlTH6/lwedGLJwt6xCZbAQDvKgHBAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXQAAAAAAAAA=eF5jYACBD/bsb9xnzpz50p4BDWTFcBkbGz+Gi7e/k0tLS3uGoY5ldfPsmTN3YogzgskH9vf9H8yaOfOk/YKu/jNnzrzBUAcDH5UOAe27DJfXdn8NdNdNewC2pCLTAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAngAAAAAAAAA=eF5jYACBD/bsb9xnzpz50r5i6ZQzZ848sTfTawXyH9pnxXAZGxs/tmeAgvZ3cmlpac/s+bZeA9LX7N1zN/////+5Pcvq5tkzZ+6E0ivtA466A+Xf2TOCdT2wv+//YNbMmSftF3T1A81/Yz9t8Vcgfcd+g1CgibHxYfuQ3a+NjI1f239UOgS077J9gepMoL4L9trur4HuuGkPAPHLS0A=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYQAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL6Dqd0LplVD5d/aMYPkH9jPA4ieh5r+xPwem79gbgc0/DLXnNZR/2f7ffxA4DzXvpj0AbEpLzA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA5gEAAAAAAAA=eF47+R8E7tt3bcu9vS33tP1iKC1aG84tpHXFRuZGv5Gx8WV7f6aX/uylB+yjobR7he4RQ+uDNp+mBgHlD9tPbtgCFN9g3wWlxbJl7HaWbbHhehw0Y+bMlfa/Q26KxH+cav8KSi9P5v104FHDHm+w/E777WDxpfZ/oPSz2G31wg837/kDcd/+mxB3wWn/SsY3m9Y17VHeVQ3Uf9L+S9gfoL6t9hNg9GL+4uv16TaPIO7b3w1x1/46KL0HbP+EPSpbz/77//+8/Rzd9FvbcnfbT4TSPsHn6hsrZth4icrPnDnz5v6k6A7R+I9H4fRVjYv5TKpL9szr6j9z5syd/dKc1T+61p/YLwWljSWtQ/cdKbOZBZG3Z4KI2wtC6QiXa9+4zKbvYQCChkPbbUD0AfPdYJqhbJ1NvctWlpSgSXs4IOG/PwQS7vsDoPTH1ZO27ZVab2MFcZ99AsRd9pFQuuDfvwNOyadt0iDhv98qFBzu+02g9N4lV8u/8e+xWVZ2PTUt7dr+G+t2gMMJRnvxHZv2vChvz1yIvP01iLj9JSj9IcmY0fpE2x41SPjvz4KEO5zm6+Y+3HR+/h59SPjuh4bv/ulQ+pJz5r/jGqv2nIfE//7JYHct3d8PpcPXP5r2SnyhDQBdKF+tAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVgAAAAAAAAA=eF5LSwOCRUUOh7Zz9lUAaQY0cHytpEcdUDwNqk716be9NVjUcU9qLcnBIj5T0XdmLVDcWf6aRhGQvnlmk0shFnUw8Nz5u0Y5kvz1G/skqoF8AEpEKWo=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAkAAAAAAAAAA=eF5LSwOCRUUOh7Zz9lUA6TJPjuZaIL0hffMHEH18raRHHZBOg6pTffptbw2Qjue7rVAJpEv8GJZVAWnuSa0lOUDa4k9pbiaQXjCz52o2kJ6p6DsTZI6z/DWNIiB988wml0KQur+skSBzlt9Y65cPpP2c4l6WAOnnzt81yoH01I3PVoD412/sk6gG0gCeiU+CAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9gAAAAAAAAA=eF475Kd9xTcpYL+Xd+8n0wcKe3jPcl8G8Q9e5Jrych3jntnXQy6B+JP/5lissrKwKZmoAOa3PX1++k+n2p7OLVZg/teyq76Fz79a74GaxwAFu2Q7wfKLjnoa3E9R2zMDah7v3W+bbUtZ9lxJ3gvmz3v05nD6DdU91YcOge1/bF74v8X2p/UuNPP2Q/mrIeptqqHuYWC5d8evr82aG+r+/H1eXl9KXu5uhpqnfuXd94k/rG0aoOp/rfwIpA1sLMyiwPIRLx5dnmj9fbcTlL8Jav5uqPsLTObxLXz9xfo21L3rwPJfd/dB/R/3cEXax3QBGwCDyKo7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9AAAAAAAAAA=eF5jXvXxkm9SwP7vk+0+mz5Q2PN3JYTvxZD8fVcg054fUL6/9o2lnn/Nbb5A+YZyj40uzVHb8xPKdznWkXxB8os1TD8Dw4F1Vo6Cu2Dyq5zuxa5VUIebt7nsnjyzI/MemHmlrO/ZVnqp7fkN5e/TjK8zXvfT+g+U/3jp7CMLtHusGaDuXbKi/nbGZlWbXzD73C5fe7pmJVz9r8qWG65Zr3azQNVP+ZAgOPmytc1/qPzSkDmC9asMbGDu0Us886vP4ftuJqj6fUUZG/nmq9rA3JOxaCnfurCP1jD+1lWrVpk++L4bZt+erZe2ZzQK2QAAOJO7aA==AQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d113ea3be156797067ad6248bf7f668043e27e8d
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_newton_ts_2_t_2.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="26" format="appended" RangeMin="45" RangeMax="121" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="92" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.029894096593" RangeMax="0.029723375654" offset="208" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="348" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190476" RangeMax="0.050476190476" offset="464" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-3.3646614802e-18" RangeMax="3.9588653941e-18" offset="824" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-3.3646614802e-18" RangeMax="3.9588653941e-18" offset="996" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142857e-12" offset="1256" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="7.505834107e-26" RangeMax="3.7142857143e-11" offset="1480" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2188" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2280" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2412" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2584" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2676" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2848" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="2940" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.0356994338e-15" RangeMax="0.58466805513" offset="3000" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15053296" offset="3656" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15053296" offset="3808" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4032" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4372" />
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="4740" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="4920" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5044" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5136" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9E1NLSw1E1PSkpMTLG0MNdLySwqqQQAVuIHew==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANAAAAAAAAAA=eF6TcJZ5/ciMy54BAvYwYAUfrPUg6vZDBRqwq7thjV38wW4oowGNxgHg6qHghzUAw9gPAA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAARgAAAAAAAAA=eF5bHFP882z+PHsGEGCUt2XABuRlbfnkb97pmzNvP4jboKe4F6u6w7rYxU+JQcQTeCC0FxN2dTAQ8sEGhX9Oci8AxoET8A==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANAAAAAAAAAA=eF6TcJZ5/ciMy54BDP7tZsAKvljrQdTthwo0YFd3wxq7+AOYuQ1oNA7wAM0d36wBXh4P/w==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA7QAAAAAAAAA=eF7jEASBxv16ay9Krb3YZc8ABgk2rdZyXhyvp0H5VTbuy6av5fm3GC5f8NbWdRfXIig/Zc+5k1mVySyT7fkg5tmLXXGUveKYC5V32LOq6cjDrd+nQPgNtXv2fRd5qeUxBypfvEc3/qNo/MelUD5Dw/LokPt97U1w/sNFv8pspkyA8mNstl5Y8oz5xkqY/TYSYHsnQvkZNo5n2TcY57XB+WpfyjdF3ZsB5SfvWW7y5E66+AK4+TUi69wfVkHdx5Cz57G4d93GJ8uh/AabB3a7G26EwMIndY9EWyIQwcxrtKnZkBdYa91nDwAz8GGoAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYAAAAAAAAAA=eF5jYICAAxJtj69cd7VhQAMOF24sPesYBBfnUJnNLWFgiqFu2k5xy1X7vDDEpV7OFZzK4b8nRPfd9dbdZjbnJCW43q10wlAHAz/FSwMXywTC5bUmepv/53DcAwCnGx6BAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAogAAAAAAAAA=eF5jYICAAxJtj69cd7Vxmrl/xeVed5sZd/1eSLwRtXG4cGPpWccgG6gyBg6V2dwSBqY2J+bNy3Y3MLRJurqy5IO5g820neKWq/Z5QWkrm3OSElzvVhrZSL2cKziVw39PiO676627zaDiTjZKe0tUTnF47FHIaJ9kxuds4y55vuyhg4vNT/HSwMUygTZ+iwuOt5x2s9Ga6G3+n8NxDwDO00FrAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhwAAAAAAAAA=eF5jYEAFr5bIf+ZXnrUXxneuvxH2dsoiOH9CJVeYnPkSOH/z7vBv3rMWw/nowO5ITzmT+kK4vAOaeR6iLzeILZwH57fObP5YmN4L51+ZvPXDIs56OP8tGv/Wjg+lKdkI97xBc38/mnn9aO7PNt+79cPlKXB+BhrfBs39bmju3YjmfwAOJVQ2AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA8QEAAAAAAAA=eF7bvuNDaUr2EduoT2stV83YtIsBCqbn1VfMcTthtWHV28/3GnbYnlY7Pill2SK4vMnhrg97WbdbdYLlV9gq/CiudYpcaAWTlyp9f1eKfYKV9PHP31+GTrD9V+zU9doPIf9HdGf/zjnTdhmB5RfYimWo/bGxWAY333KK8oPCRUutrkPct3eRgPWvWReT4PLNpksYm/+k7OIC699gu8ypdopH6Cm4+TXOPwMEJRutlkDct5fpfnOQfkw/XH+07RFVN+cpu+6Czd9ia6J4X/dO3AS4/ic6S5j/qU6zkgSbf2Dvis0riw58jYHLT1lzckGL1NRdn+w/mS1POLT3Rnonn5dhOFy+5xKfgcnOBquvEHlbN6Fk7mTmNXD7D7QLNjZxLAfzDxxabPVt6+P10gyL4foZvfx+qx7cZNUHCf+9DesPr97fjLD/r2ejtif/UisRiPtsszQf3OPvWASX/yfjpqhiv8uKDRL+e2cZ2B6P/lsMlz8WmBPw4PlCq+g5rwojp+3Zuy71uxa7bS9c/qN1fdq8BWlWgRB5203MIs5MMp1weS2WkHUn73RZKUDCf+9lhfuylxfUw+VzRTznFF9bb/UYEr57zyq8eMzzuwGh32d+x7O9LVaskPjfe6Dm6TJ7Vxe4/JHtphZOmvOtAOXxBgk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAQQAAAAAAAAA=eF5jYACBD/YzweClPQMaMAaDx3Dx9DQQeIahbgZY/04McUYw+QBq/kn7s2dA4A2GOhiA2HcZLg/Rd9MeAPoQIv0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXgAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL+xngNXvhNIrofLv7BnB8g+g5p+Emv/G/hyYvmNvBDb/MNSe11D6MtTcC1B9N+0BGpVF0Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIgAAAAAAAAA=eF5jYACBD/YMUPrffxB4D6cZ0OQJqSNXPy5zYHwAZHtUjA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXgAAAAAAAAA=eF5jYACBD/YzweCl/dkzIPDEfhaY/9DeGAwe2zNAQXoaCDyzB1Np16DiL+xngNXvhNIrofLv7BnB8g+g5p+Emv/G/hyYvmNvBDb/MNSe11D6MtTcC1B9N+0BGpVF0Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF6bORMEdtrPgtLT0GiY+Gwi1eESR6fR1eEyB8YHAKJbX5Q=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAAyQEAAAAAAAA=eF678B8E7tvP2pZ7e1vuafs1UHpxsqei8d6rNs7GIHDZPoi9NIC99IB9FJRmW+vsdWfuERs9sPxh+2Sw+Ab7QCid880zY92WOTaHZoLASvsJ8R9F4z9OtQ+B0k2/vJwv2i/ZcxAsv9P+XxxIfKm9ajyEtue48WJD/GYbVgYQeLD/AcRdcNrhruIvidNFezaA9Z+03wTWt9U+FEobLFzXGzm/z8YW4r79KRB37Y+B0pzFrNXmqov2cILNv2D/BWzubvs3UNonecHP2quLbZaBzb+5XwFs7lE43VeY/vj700V7Tp8BgTv7f3StB6ITcHrzbk/XFZ97bc5B5O2fQ8Ttv0LpAx07zoct2LIHbP2lozYgKmEKhGb4vtNmUjbPloLTB2y0IeG/3xMS7vs9oHSceN+LF7c32SyCuM9eE+IueyUorVC17Q1j43GbVkj47/8JDt+p+xkg4b+//9jhmPscG2xS00Dg2n5oPO43gNLTXJzf+zqX2MRB5O31ofGsBaVj9yib6fHNsIGG/341SLjv14DSZ7oZn+5v2mPDDQnf/a8g4br/LZQWCVlRwD55kk0HJP73s0HifT8HlL6+aNuy5c9X2wAAhqxO9g==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAUQAAAAAAAAA=eF5LSwOCRUUOc3+vOpoOpBnQQKnV/tYMoHgaVN2L2Rn2GVjU/Xf5zIlN/21P+6kg9e/f75oPkpfeYKeETR0M1NQfkUI2f8rzhfkgPgD1QypeAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAhQAAAAAAAAA=eF5LSwOCRUUOc3+vOpoOpPVeT07KANIaTvW9ILrUan8riE6DqnsxO8MexO8wNnYF0YdOdzKB6P8unzlB+jdd1dkBUufQHXQERN/2tJ8Kkn//ftd8kLz0BjslEL1jRkETSLxc6noIiH/n6pdiEF1Tf0QKJL5tu/1dEH/K84X5ID4AFNpQaQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA3QAAAAAAAAA=eF5jXPXxkm9SwH5+RyuGyZ/09nxfCeHP4Lf7mpyjZvMXyg+RfrGl0c/a5guUv/HRm8PpN0T3/IHyM66vqyg8ZmTzG8pnmLkhMGrRlt0/oPzL9q0v7uY423yF8g95fG91/8e75z2UvxJing3MvG/3lQzOn3y4+xeUDxH/sPsflM8KdreBDUw+SV9vj/JdNbj9IT8mFLLG/LBmgPpvGfeKpZoFDntg+g+B7WOFu//T2xuPd897sRvGh/hPdc9PKF+kXyLqxWZOGxifAQw+7oCpz+HidpmwWc0GAOpfv6c=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF5jXPXxkm9SwH4BRyuGyZ/09nxfCeFbMvQIiexXs/kL5UtnBCmuvmdl8wXK72d0VC34KLrnD5RvZzXNpPOzoc1vKJ9j5obAqEVbdv+A8o1u/y6rTnO2+Qrl+wmbGCpU8+x5D+WfcBCV7hMUs4GZd3L/PsNctee7f0H5HxdGBJfu/Lj7H5Rv4qJw6Fi5gQ1M/oYX9+PXmRpw+8+8mabNa/3HmgHqv7LCym0bohz2wPQ7qV1INtzBAHd/6ssFWsx5z3fD+EzbXlyTbVHd8xPKv/Du5r1XazhsYPwlNt0WZwTWWsPUm/va+D7epm4DAAtitE0=AQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..11539ae284f22bc4fd11244ab5d2480e3a054314
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_0_t_0.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="144" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="0" RangeMax="0" offset="204" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="264" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="328" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="0" offset="448" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="508" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="568" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="632" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="696" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="0" offset="792" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="852" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="912" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="1008" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="1068" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="1160" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1220" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15" offset="1284" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15" offset="1372" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1440" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="1504" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="1568" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="1660" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="1720" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1900" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2024" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2116" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw0E03TE42MDe2NAQALXoEiQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAADwAAAAAAAAA=eF5jYBgFo2DwAgACoAABAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAHwAAAAAAAAA=eF5LSwOCRUUOaVCaAQ2kockTUkeufnT16HwA5KYkhw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAEgAAAAAAAAA=eF5LSwOCRUUOaYOcBgAG90W7AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAADQAAAAAAAAA=eF5jYBgF1AQAAVAAAQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7b6d273dea266394083594e95dfe2e8645a6cf86
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_1_t_1.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5.000000031e-05" RangeMax="5.000000031e-05" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.03056133107" RangeMax="0.028248973988" offset="204" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="344" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190561" RangeMax="0.050476190468" offset="464" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="836" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="-5.9428564254e-10" RangeMax="5.9428574849e-10" offset="956" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="-5.9428564254e-10" RangeMax="5.9428574849e-10" offset="1128" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142719e-12" offset="1384" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="1.0423184196e-25" RangeMax="3.7142857143e-11" offset="1608" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2328" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2424" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2592" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2848" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2944" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="3116" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3208" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.1047392924e-15" RangeMax="0.58466805513" offset="3268" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15583048" offset="3968" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15583048" offset="4128" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.00099999999703" RangeMax="0.0010000000045" offset="4368" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4736" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="5104" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="5196" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="5256" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5436" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5560" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5652" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw0E03TE42MDe2NAQALXoEiQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOQAAAAAAAAA=eF5rt7/69pEZlz0DGLDtYcAKPlg3QtTth/B/WGNXd2E3dvEPUPELUH1PcOiHgQNo8hd2AwC86xNTAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAARwAAAAAAAAA=eF4LD93Q6/N2jj0DGNy2YcAGjOX3Ttt1z7TCe/5+MJ+H3Rarui2f9mAVVxYFq3f4a74XzFf5gd0eGIhUQDX/gNFeAN8AE70=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOAAAAAAAAAA=eF4zcZZ5/ciMy54BDFj3MGAFH6y1Ier2Q/i/rLGru7Abu/gHmHgDhHqBQz8MHEAz58FuAHSjEUY=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA9QAAAAAAAAA=eF7jEASBxv16ay9Krb3YZc8ABil7Fhe98eR4PQ3Cb6jdo5dTupbn32KofMEes2pD111ci6D8GpuEp4GVySyT7QUg5tnLXHGUveKYC5XPsHn0acbDrd+nQM3r3/M5O/KllsccqLyBzY1Nh0TjPy6F8hP2/OWYf7+vvQnKD7Dp+SlfbjNlApQfsydeddIz5hsrofw6GwmwvRNh9u1hi8/eYJzXBtO/p2Fu0KaoezOg/OI9D5Z9uJMuvgDKz7VhN2H1eFgFdR9Dso10h3zdxifLofxmm5w4ocYbIbDwybc5cG+fRFsizLwWm+TKBYG11n32ANpRYfk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAYAAAAAAAAAA=eF5jYIAAmetTfFKyWfYxoIFZS3hZUrK/7IXxL14o2PL9LiOGOustO8+JLfxniy7eKCNy62PJbJsXS1d0As23u3ujvw2obi8DDpD+7hDHj7uMdjD+Hdv0M0D7bQF9zyZFAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAnQAAAAAAAAA=eF5jYIAAmetTfFKyWfbJLS7ftlL4z16ju5bTkrOf7J21hJclJfvLXqgyhosXCrZ8v8u4L7T3W94q4T+23y43MaxTYd5nvWXnObGF/2wh9Dvbuzf624D03kYZkVsfS2bbvFi6ohNovh1E/N9eF1Zrq9TsJ7bMB5p+AM2zO5r47jjI3PR3hzh+APkap0w9geba3bFNPwO03xYAf8hPgw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhgAAAAAAAAA=eF5jYEAFO9O3f+JXnrUXxtdkyA17O2URnG92oitUznwJnL/n9v2v3rMWw/noYEufYTmT+kK4vDqaedOKIzaILZwH5+fyuX4sTO+F8338pn1YxFkP5yeh8eNy8kpTshHu2Yvm/hI080zQ3K8vWLP1w+UpcL42Gn8TmvunoLl3F5r/AVkVTpk=AQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA+QEAAAAAAAA=eF47uONDaUr2Edtd7hdvJF4+sIsBClaU+25sPHncKsZA9vO9hh22004HagZlLofLrxI5FcHTu9cqi+n8p3sNK2x/JbrLp5ROsILJ/8u5f6bsR4NVvY/y95ehE2y/bc7wdfs+ES4PARXmm8DyC2yZrhjUK+1fDzff4deLzISflVa3Ie7byym1fcMmzmlw/c0vmjsKLPN2ccbYAvVvsF2+fakxZ8dpuPwDN88pYqlJVm4Q9+3tqs1pfnOqH27+nV//7tZqx+1SctMFmr/FtmjH+0YO9ilw/U+953c829tiZThzKdD8A3vLLdYJuIq1wuWLrWKtFhbO3sWsedpsecKhvds3u4bdOd4Pl59XEVeiLtFtxQ6Rt80Ns1K8sHcb3P6Vm5beihddBOY7fNm563hGcejMrQvh+me52qsc/rPBygUS/nu9b/gxC1wrhsvP9J3h7/dymZUxxH22fm1skSLls+Dy7NE2ng08u60eQcJ/b/mTQqM9J7IR4S+2s3/nnGVWHoYmhZHT9uz97zntyWmXqXB5zw+fdlyZUWZlB5G3fbVEceXLAz1w+X7BV4GHtvRafY8Gh//eeSe951i/6IXLL/97PNQ/cNIuWUj47t06O3dJdQIi/CYzCN+W2Vayawsk/vfO41JREM8xgMuLvfp4pqarYRcAKDL5yw==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXQAAAAAAAAA=eF5jYACBD/bsb9xnzpz50p4BDWTGcBkbGz+Gi7e/k0tLS3uGoY55dfPsmTN3YogzgskH9g/8H8yaOfOk/fyu/jNnzrzBUAcDH5UOAe27DJfXdn8NdNdNewC1xCLRAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAngAAAAAAAAA=eF5jYACBD/bsb9xnzpz50r5i6ZQzZ848sTfVawXyH9pnxnAZGxs/tmeAgvZ3cmlpac/s+bZeA9LX7N1zN/////+5PfPq5tkzZ+6E0ivtA466A+Xf2TOCdT2wf+D/YNbMmSft53f1A81/Yz9t8Vcgfcd+o1CgibHxYfvg3a+NjI1f239UOgS077J9gepMoL4L9trur4HuuGkPAPALSzw=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYACBD/YzZoLAS/szYPAEyn9ob2gMAo/tGaAgDQye2aeA6Wv2/8Hguf00sPqdUHqlPUTdO6i+B/bTweInoea/gdJ3oOYftjcC06+h/Mv2f8Dmnoe646Y9AMFVUag=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA6gEAAAAAAAA=eF67+R8E7tv3bMu9vS33tP1qKL1poqaQtu9VG/Ub/UbGxpfto5he+rOXHrBPhdJh+hN+tJSdsvkyNQgof9h+UsMWoPgG+24o/ebZk2dfUvpshB8HzZg5c6W9ROhNkfiPU+2/h0BoBjBosXQCy++0PwAWX2r/E0p7vNNo7HBqt2EBq3uw/w7EXfsfQmmrDYwlTwWq9wTtqgbqP2kvE/4HqG+r/e4wCJ0hbcfJUFlgcxPivv1VEHftL4bSsWzrX6azZO/R2nr23///5+3n66bf2pa7234alL4evKKAffIkG1dR+ZkzZ97cHxPdIRr/8SicdnhuzVExecWexV39Z86cubNfnLP6R9f6E/ulofRT9p+3wp1n2CyEyNszQsTtRaH0oYuJh+cabtgD8p1CykUwzVA9G0K/3LZnaVTCvOU5+2zEIOG/Hxr++6Oh9JfVIbGpLVtsPCHus8+BuMs+DUqbyq/JfmFxwmYeJPz3R0HCf388lO745eV80X6LzbKy66lpadf2X1i3AxxOV6D0ug1B00NsWmzmQeTtr0LE7S9D6W95LAtEZs2y4YeE/34fSLjvD4bS+/bnP/ZMX7BHDhK++5sg4bq/DUo7HdguwdfTtCcEEv/71cDuWrpfBUrPcf/y6c/83j0ABH5TbQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVgAAAAAAAAA=eF5LSwOCRUUOh7dz9lUAaQY0cHStpEcdUDwNqk716be9NVjUcU9qLcnBIj5D0XdmLVDcWf6aRhGQvnlmk0shFnUw8Nz5u0Y5kvy1G/skqoF8AEmEKWc=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAkQAAAAAAAAA=eF5LSwOCRUUOh7dz9lUA6VJPjuZaIL0uffMHEH10raRHHZBOg6pTffptbw2QjuO7rVAJpEv8GJZVAWnuSa0lOUDa4k9pbiaQXjCz52o2kJ6h6DsTZI6z/DWNIiB988wml0Igbf6XNRJkzvIba/3ygbSfU9zLEiD93Pm7RjmQnrrx2QoQ/9qNfRLVQBoAm3FPeg==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF4746d9xTcpYP8EV7t7Lwy89wic5b4M4nMJ/Jmi2yGxZ+H1kEsg/oZHbw6n32DdUzRRAcxngIK+LVZgfpmA+Os1ri579kDNg8nvl+0Ey3s/UveW3aNvMwtq3ofIhP8LjI5ZX0/eC+avAZsvalN36BDYfruOe2/cHzDZ7EIz7xiUX3SdYYmwkPCeOqh7nBM6ZoivtrHhg7r/zdHrD6UTFG3aoOZ5HHQN6lSzsmmAqveSfrGl0U/bxsosCix/7/7H3ZfUWPfYQfkbwe5RtdkHdT+D7sSoLoWzu+9A3QuR/7p7EtT/+6TDnOdOUNwDABS+mjw=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8wAAAAAAAAA=eF7jW/Xxkm9SwP7tmfn3Xhh472GC8qfsv/vyHYvknv8rIfzGWw7pdwQZ93yA8s+2LuHMnfHQ+g+Uv5zhxJS8N857fkH5DAwMf8z0VHcxQM1juvftvEmsgc13qLz8j/XsVlM+W/+A8ntUZ+QHzxCz+Qfls3z/657JyGoDk3/2e+VHLs6D1hxQ82Iur7o6467gHpj9Wspn9ioK2drA3FtV5FU6O0bFhhOq/oVQ6e5jB6zg8l4fEirUNuvawNx7JbDrTKvVn90w964xv3SjJFgN7h6HN6dXf5Xbtxum/ponD9NGrkdw9TEspgIXtKT3AABXfLS2AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..324ac934602eb081446f79f093487b13ed62bd8e
--- /dev/null
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas_square_1x1_quad8_1e2_bottom_ts_2_t_2.000000.vtu
@@ -0,0 +1,48 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="GasMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="84" />
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0.029894096594" RangeMax="0.029723375655" offset="204" />
+ <DataArray type="Float64" Name="LiquidMassFlowRate" format="appended" RangeMin="-5e-05" RangeMax="5e-05" offset="352" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.0090476190476" RangeMax="0.050476190476" offset="468" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="139" offset="836" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="2.6164164567e-17" offset="956" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="2.6164164567e-17" offset="1128" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="9.2857142857e-12" offset="1388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="5.1677569392e-26" RangeMax="3.7142857143e-11" offset="1612" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="2308" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="1" offset="2404" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2532" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="1" RangeMax="1" offset="2692" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0" RangeMax="1" offset="2788" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2948" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3040" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="3.9998130556e-16" RangeMax="0.58466805513" offset="3100" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="298.15053296" offset="3740" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="298.15" RangeMax="298.15053296" offset="3892" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4116" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0.001" RangeMax="0.001" offset="4436" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="4792" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="4884" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="4944" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5124" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5248" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5340" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9E1NLSw0E03TE42MDe2NAQALXoEiQ==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAOQAAAAAAAAA=eF4zcJZ5/ciMy54BDDj2MGAFH3ZrQdTth/Af7MaujqEBu/APawj9AqrvgzV2dTBwAF2+AQCdchA3AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAATQAAAAAAAAA=eF5zry77eTZ/nj0DEDQYhO1lwAZWBtim9t260zdn3n4Qd0GAJ1Z1B74722ITL3iQAhafYZ4Kpi+0xGK3Bwoc9piiyL+4m2ELABr8Gsk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANQAAAAAAAAA=eF7Td5Z5/ciMy54BDNj2MGAFH3ZrQ9Tth/Gxq7tgjV38D1T8BVTfJxzqYADdnAvWANhAEXE=AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA8gAAAAAAAAA=eF7jEASBxv1Gay9Krb3YZc8ABiU2bdZyXhyvp0H5rXvclk1fy/NvMZSfZJP/1tZ1F9ciKL/O5tzJrMpklsn2bBDz7AWuOMpeccyFyjfYrGg68nDr9ylQfrnNvu8iL7U85sDs26MV/1E0/uNSKD9hz6rokPt97U1Qfr7Ng0W/ymymTIDy0222XVjyjPnGSii/eY8E2N6JUH7GHoez7BuM89pg5tlofinfFHVvBoSvMMlmhcmTO+niC+D+qRFZ5/6wCua+mj2PxL3rNj5ZDjfvpt3uhhshsPDJt5FpS5RoS4Sax9C5p25DXmCtdZ89AItoYvI=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAANwAAAAAAAAA=eF5VxTEKgCAAAEBLS+3HERESEtGnG7qlW24Jn+Lu03P4yz58efDqx5NHVzdvvh29O/kFt9cFTA==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAXwAAAAAAAAA=eF5jYICANQdKQ78dybZhQANKzHd8WsJL4OIb4h7fqHb0wFBn94ZNcl1OJYb4mQeyF24x1NpobZH+82NCrQ3T14dbOrgyMdTBgI5SU6KPVh1cfsox/k6rA1U2AGjeILk=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAogAAAAAAAAA=eF5jYICANQdKQ78dybb5kffj+Ob4dJvvH2v69NdU2Cgx3/FpCS+xgSpj2BD3+Ea1o4dNwhePUufSGhsdn+tLV71KsLF7wya5LqcSSmfaMH19uKWDK9LmzAPZC7cYam20tkj/+TGhFiqeabPh8rQ1XQ+qbUwCJnZdrqu2uX5Lu3V/epaNjlJToo9WnY3263PrFt2ttZlyjL/T6kCVDQAMNEn4AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAhwAAAAAAAAA=eF5jYEAFb5bIf+ZXnrUXxneuvxH2dsoiOH9iJVeYnPkSOH/L7vBv3rMWw/nowP5ITzmT+kK4vCOaed6iLzeILZwH53fNbP5YmN4L59+evPXDIs56OP8TGv/Wjg+lKdkI97xFc/9ENPMmoLk/z3zv1g+Xp8D5mWh8GzT3u6K5dyOa/wEllVRXAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAA6QEAAAAAAAA=eF47vONDaUr2Eduqi0ZtMh8O7GKAAiufBmlb2eNWa1e9/XyvYYctS4yi0221tXB5/hStGUKle63mgeVX2B5ZYb8inWmBFUze1md+x7O9U6xeHfv8/WXoBNtLtRzfViQg5J3B8i1WPMdB8gtsHaZEHVyVsRpuPovify7JJRus7kDctzf98i67c5oT4PrjQsr/LJYqsfoLNn+D7a6tyv5CP07C5XWqb4WGPi60mgFx3163vmmS2++Wws3XkTu5IHb11F03wOZvsT3c861r9ukmuH4waJi3XQLsvgN7l4hEFah/74fLn5qwPanhV+uur/afzJYnHNqbqBIXFNHdDZd/nXP/TNmPgF0/IPK2HpVGPgtDd8Htj1AROCEYuADMd/i8aFfaw6BPyhxL4PpfrFugdMd0gdUaSPjvVREPt2vRb4Hrn3HO8NzX5MVWQhD32eoFvH0mapcCl68I27ZwZf9+Kw2w/IS9n1SMxYP6euHmx/z6d7dWu84qbc6rwshpe/au2Flh2vMWIQ8BBeYhEHlbt/YiTV3ZVrj5fJZfn/tFV1tBwmfD3vnyrzmi+ifC9T/5bTTt9P8lVvch4bt31uUQw8W+9XD5E4E5AQ+eL7Rih8T/XvvHHH83/euAy98r9H+9tr9hFwDCfglUAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAPQAAAAAAAAA=eF5jYACBD/YzweClPQMaMAaDx3DxNDB4hqFuFlj/TgxxCHhgD5E/aX8GDN7gUAez7zJcHuKum/YA+hgi/Q==AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVwAAAAAAAAA=eF5jYACBD/YzweCl/RkweGI/C8x/aG8MBo/tGaAgDQye2UPoa1DxF1D1O6H0Sqj8O6j8A6j4Saj5b6D0HXsTsPmH7Y3A9GuofZeh+i5A3XXTHgAahUXRAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAJgAAAAAAAAA=eF5jZQCBD/YsUPrvfxB4D6dh4kxEqsMljk6jq8NlDowPAHETVK0=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAVwAAAAAAAAA=eF5jYACBD/YzweCl/RkweGI/C8x/aG8MBo/tGaAgDQye2UPoa1DxF/YzwOp3QumVUPl3UPkHUPNPQs1/A6Xv2BuBzT8MpV9D7bsM1XcBqu+mPQAYhUXKAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAIwAAAAAAAAA=eF5bOBMEdtrPg9LT0WiY+Cwi1eESR6fR1eEyB8YHALXLX8k=AQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAADAAAAAAAAAA=eF5jYBgaAAAAqAABAQAAAAAAAAAAgAAAAAAAAKACAAAAAAAAvwEAAAAAAAA=eF67/x8E7ttP2pZ7e1vuafvNUFrS0OXxs01XbOyMQeCyvRt7aQB76QH7YCgd+dp4Jk/iKRtbsPxh+2yw+Ab7OCitFbKigH3yIpvSmSCw0t4z/qNo/Mep9nxQWgcsP8mmCSy/035XHEh8qf0bKH1yl47Z7bN7bVgYQODB/kcQd+1/AqXfTtkfeIWhyWYeWP9J+/lgc7faO0Hp4oLyl2dZGmxMIO7bHwJx1/5AKB3R9Jd5XfriPYxg8y/YPwWbu9v+AZSGWLtm53Kw+Tf3y4DNPbpfDkrXzag60b1u8p4LZ0Dgzv5vXet/dK0/sf8HlD737MmzLylxe85B5O1vQcTtP0BpXq0H79Qs1+4B2/N9PYQW2mwDpitybbYtkcnN5lhr4wIJ//1JkHDfHw+l9+/bImYbu9FmFsR99mIQd9mLQ2nvn/Lzop3P2GyHhP9+Tki471eG0o9Z179MZ+m2yUoDgWv7DaHhZAylGSAOsUyAyNurQ+NZG0over7velV9h81iSPjvF4KE+35pKK2+g2UH57lNNhyQ8N3/GBKu+59B6UnHDsfc59gAi//9XyDxvv8vlD7LY5rhubJ3DwBNXUclAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAUQAAAAAAAAA=eF5LSwOCRUUOc36vOpoOpBnQQJHV/tYMoHgaVN2z2Rn2GVjU/XX5zIlN/01P+6kg9e/e75oPkpfeYKeETR0MVNcfkUI2f9LzhfkgPgDwYypQAQAAAAAAAAAAgAAAAAAAAKgAAAAAAAAAhQAAAAAAAAA=eF5LSwOCRUUOc36vOpoOpHVeT07KANJqTvW9ILrIan8riE6Dqns2O8MexG8zNnYF0QdPdzKB6L8unzlB+jdd1dkBUufQHXQERN/0tJ8Kkn/3ftd8kLz0BjslEL1tRkETSLxM6noIiH/n6pdiEF1df0QKJL5lu/1dEH/S84X5ID4AC7pQTw==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAAzQAAAAAAAAA=eF7jW/Xxkm9SwH4V6565H/fy2fxdCeFHZM37rrpBcs9PKP/b2xuPd8/jsIHxGcH6FPb8g/I9IyUKk9+LwPkMUPAHylfJyEqVXv5t93co/4TH91b3f7x7vkD5yx69OZx+w9SGAeqeX8vPbWYP/Gz9B808jlXI/B87YPaZrYiO8BI22fMfym9yW9qaxypiwwZVv6vg97qrz533/ILKrwLbx7oHZj4rWB0DnL8FLC8K50v1S0S92MwJ9z/E/g0WMPvXvuk/vmO+sg0AmeGyAg==AQAAAAAAAAAAgAAAAAAAAFABAAAAAAAA6gAAAAAAAAA=eF7jW/Xxkm9SwH5V6565H/fy2fxdCeFzLJoYW10svec3lP8mcdqlDyoCNj+h/CzmdTkaXkp7GKD6p1lK93JkSuz5A5VnAIMGCxh/RcEHfoE3X62/Q/kur9+7HKrn2PMFyi8Nulk0jcXc5j+Uf0Fm1RW7Z8+tf8HUpz1TXVh+djcH1D4tL2Yvs+/ndv+Dynd4iq/eLWO6B6af55bp0lfFAjZsUPV+5f3SR1Rd9sDMK1VOjmYuYoW7917N/qPepxH8azPMNM9sFIfzJQ6alc6fwQ33/6XZTNOdfl6B2/+ez2fh2aeqNgDdrLMOAQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIgAAAAAAAAA=eF4txbcBACAIADAsiP7/sAPJkog2PL28nT4uXz9/BXgALg==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAADAAAAAAAAAA=eF5jYKAuAAAAUAABAQAAAAAAAAAAgAAAAAAAAPgBAAAAAAAAZgAAAAAAAAA=eF5jYACBD/YMaGDWTBB4iSF+Bgye4FD/EEPcGAweY4jjAmlg8AxDPUT8Gg5zXuBwz04c4itxmP8Oh/kPcJhzEkP8LDh83uAQv4MhbgIOn8M4xF/jCM/LONx5AYc7b2KIAwBUVUXXAQAAAAAAAAAAgAAAAAAAAPAAAAAAAAAAOwAAAAAAAAA=eF5dxzcWgCAAwFALYgEbcv+zMhiWZPkvy/C346E/8dEX/PQVV33GWx9w1kec9CNu+gtffS9hA6UoATQ=AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF5jZoAANijNCaV5oDQ/lBaC0qJQWgJKS0NpOSgNABTwAKY=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TFYUBAASNANM=
+ </AppendedData>
+</VTKFile>