From 65a8c10929aafd6ede4eed07c454e729660df6c5 Mon Sep 17 00:00:00 2001
From: Norbert Grunwald <Norbert.Grunwald@ufz.de>
Date: Sat, 10 Apr 2021 00:06:45 +0200
Subject: [PATCH] <Tests/TH2M> add 1D heat transport Dirichlet BC
---
ProcessLib/TH2M/Tests.cmake | 15 +
.../TH2M/T/T_1d_dirichlet/T_1d_dirichlet.prj | 532 ++++++++++++++++++
.../T/T_1d_dirichlet/boundary8_bottom.vtu | 24 +
.../TH2M/T/T_1d_dirichlet/boundary8_left.vtu | 24 +
.../TH2M/T/T_1d_dirichlet/boundary8_right.vtu | 24 +
.../TH2M/T/T_1d_dirichlet/boundary8_top.vtu | 24 +
Tests/Data/TH2M/T/T_1d_dirichlet/domain.vtu | 24 +
.../result_TH2M_T_dirichlet.pvd | 10 +
...esult_TH2M_T_dirichlet_ts_0_t_0.000000.vtu | 44 ++
...ult_TH2M_T_dirichlet_ts_10_t_10.000000.vtu | 44 ++
...lt_TH2M_T_dirichlet_ts_20_t_200.000000.vtu | 44 ++
...t_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu | 44 ++
...t_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu | 44 ++
Tests/Data/TH2M/T/T_1d_dirichlet/square.vtu | 22 +
.../TH2M/T/T_1d_dirichlet/square_quad8.vtu | 22 +
15 files changed, 941 insertions(+)
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/T_1d_dirichlet.prj
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_bottom.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_left.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_right.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_top.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/domain.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet.pvd
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_0_t_0.000000.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_10_t_10.000000.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_20_t_200.000000.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/square.vtu
create mode 100644 Tests/Data/TH2M/T/T_1d_dirichlet/square_quad8.vtu
diff --git a/ProcessLib/TH2M/Tests.cmake b/ProcessLib/TH2M/Tests.cmake
index a743dc07926..68065f1c200 100644
--- a/ProcessLib/TH2M/Tests.cmake
+++ b/ProcessLib/TH2M/Tests.cmake
@@ -1,3 +1,18 @@
+# TH2M 1d heat diffusion w/ Dirichlet-BC
+AddTest(
+ NAME TH2M_T_1d_dirichlet
+ PATH TH2M/T/T_1d_dirichlet
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS T_1d_dirichlet.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+
+ result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu
+ temperature_interpolated temperature_interpolated 1e-8 1e-8
+)
+
# TH2M Heatpipe w/ static gas phase in radial domain
AddTest(
NAME TH2M_TH_unsaturated_heatpipe_radial
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/T_1d_dirichlet.prj b/Tests/Data/TH2M/T/T_1d_dirichlet/T_1d_dirichlet.prj
new file mode 100644
index 00000000000..8794bd52e03
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/T_1d_dirichlet.prj
@@ -0,0 +1,532 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh axially_symmetric="false">square_quad8.vtu</mesh>
+ <mesh axially_symmetric="false">domain.vtu</mesh>
+ <mesh axially_symmetric="false">boundary8_left.vtu</mesh>
+ <mesh axially_symmetric="false">boundary8_right.vtu</mesh>
+ <mesh axially_symmetric="false">boundary8_top.vtu</mesh>
+ <mesh axially_symmetric="false">boundary8_bottom.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>TH2M</name>
+ <type>TH2M</type>
+ <integration_order>4</integration_order>
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+
+ <jacobian_assembler>
+ <type>CentralDifferences</type>
+ <component_magnitudes>
+ 1e5 1e5 1e5 1e5
+ 1e5 1e5 1e5 1e5
+ 1e2 1e2 1e2 1e2
+ 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2
+ 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2 1e-2
+ </component_magnitudes>
+ <relative_epsilons>
+ 1e-08 1e-08 1e-08 1e-08
+ 1e-08 1e-08 1e-08 1e-08
+ 1e-08 1e-08 1e-08 1e-08
+ 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08
+ 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08</relative_epsilons>
+ </jacobian_assembler>
+
+ <reference_temperature>T0</reference_temperature>
+ <process_variables>
+ <gas_pressure>gas_pressure</gas_pressure>
+ <capillary_pressure>capillary_pressure</capillary_pressure>
+ <temperature>temperature</temperature>
+ <displacement>displacement</displacement>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="velocity_gas" output_name="velocity_gas"/>
+ <secondary_variable internal_name="velocity_liquid" output_name="velocity_liquid"/>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="liquid_pressure" output_name="liquid_pressure"/>
+ <secondary_variable internal_name="liquid_density" output_name="liquid_density"/>
+ <secondary_variable internal_name="gas_density" output_name="gas_density"/>
+ <secondary_variable internal_name="porosity" output_name="porosity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ </secondary_variables>
+ <specific_body_force>0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>Gas</type>
+ <properties>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>100</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>10</value>
+ </property>
+ <property>
+ <name>molar_mass</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1.0e-5</value>
+ </property>
+ <property>
+ <name>compressibility</name>
+ <type>Constant</type>
+ <value>1.0e-5</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>200.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>20</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>20</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>2.0e-6</value>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>2.e-3</value>
+ </property>
+ <property>
+ <name>compressibility</name>
+ <type>Constant</type>
+ <value>2.0e-9</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>25.0</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>30</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>3000</value>
+ </property>
+
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>3.0e-8</value>
+ </property>
+ <property>
+ <name>compressibility</name>
+ <type>Constant</type>
+ <value>3.e-10</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>2.5e-12</value>
+ </property>
+
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+
+ <property>
+ <name>saturation</name>
+ <type>SaturationLiakopoulos</type>
+ </property>
+
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+
+ <property>
+ <name>relative_permeability_nonwetting_phase</name>
+ <type>Constant</type>
+ <value>1.0</value>
+ </property>
+
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.3</value>
+ </property>
+
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="TH2M">
+ <nonlinear_solver>basic_newton</nonlinear_solver>
+ <compensate_non_equilibrium_initial_residuum>false</compensate_non_equilibrium_initial_residuum>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <!-- <abstols>1e-09 1e-09 1e-10 1e-6 1e-6</abstols> -->
+ <reltols>1e-5 1e-5 1e-6 1e-8 1e-8</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>4000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>10</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>9</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ <pair>
+ <repeat>9</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ <pair>
+ <repeat>9</repeat>
+ <delta_t>500</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>result_TH2M_T_dirichlet</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>10</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_pressure</variable>
+ <variable>liquid_density</variable>
+ <variable>gas_density</variable>
+ <variable>porosity</variable>
+ <variable>saturation</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <parameters>
+
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>21.e8</value>
+ </parameter>
+
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>0.3</value>
+ </parameter>
+ <parameter>
+ <name>T0</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+
+ <parameter>
+ <name>pGR_ic</name>
+ <type>Constant</type>
+ <value>100000</value>
+ </parameter>
+
+ <parameter>
+ <name>pGR_bc_bottom</name>
+ <type>Constant</type>
+ <value>100100</value>
+ </parameter>
+
+
+ <parameter>
+ <name>pCap_ic</name>
+ <type>Constant</type>
+ <value>15000</value>
+ </parameter>
+
+ <parameter>
+ <name>T_ic</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </parameter>
+
+ <parameter>
+ <name>T_bc</name>
+ <type>Constant</type>
+ <value>333.15</value>
+ </parameter>
+
+ <parameter>
+ <name>ux_bc_right</name>
+ <type>Constant</type>
+ <value>0.4</value>
+ </parameter>
+
+ <parameter>
+ <name>sigma_xx_right</name>
+ <type>Constant</type>
+ <value>100.0</value>
+ </parameter>
+
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+
+ <parameter>
+ <name>neumann</name>
+ <type>Constant</type>
+ <value>-100100</value>
+ </parameter>
+
+ <parameter>
+ <name>zero</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+
+ <parameter>
+ <name>displacementRight</name>
+ <type>Constant</type>
+ <value>-0.05</value>
+ </parameter>
+
+ <parameter>
+ <name>displacementRamp</name>
+ <type>CurveScaled</type>
+ <curve>timeRamp</curve>
+ <parameter>displacementRight</parameter>
+ </parameter>
+
+ </parameters>
+
+ <curves>
+
+ <curve>
+ <name>timeRamp_temperature</name>
+ <coords>0 86400 1.0e99 </coords>
+ <values>1. 1. 1. </values>
+ </curve>
+
+ <curve>
+ <name>timeRamp_pressure</name>
+ <coords>0 1 3 5</coords>
+ <values>1. 1. 1.05 1.</values>
+ </curve>
+
+ <curve>
+ <name>timeRamp</name>
+ <coords>0 1000 10000</coords>
+ <values>0 1 1</values>
+ </curve>
+
+ </curves>
+
+ <!-- ## Displacement equation ##############################################-->
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+
+ <boundary_condition>
+ <mesh>domain</mesh>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+
+ <boundary_condition>
+ <mesh>domain</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+
+ </boundary_conditions>
+ </process_variable>
+
+ <!--- ## Gas pressure equation #############################################-->
+ <process_variable>
+ <name>gas_pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pGR_ic</initial_condition>
+ <boundary_conditions>
+
+ <boundary_condition>
+ <mesh>domain</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>pGR_ic</parameter>
+ </boundary_condition>
+
+ </boundary_conditions>
+ </process_variable>
+
+ <!--- ## Capillary pressure equation #######################################-->
+ <process_variable>
+ <name>capillary_pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>pCap_ic</initial_condition>
+ <boundary_conditions>
+
+ <boundary_condition>
+ <mesh>domain</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>pCap_ic</parameter>
+ </boundary_condition>
+
+ </boundary_conditions>
+ </process_variable>
+
+ <!-- ## Temperature equation ##############################################-->
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_ic</initial_condition>
+ <boundary_conditions>
+
+ <boundary_condition>
+ <mesh>boundary8_left</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>T_bc</parameter>
+ </boundary_condition>
+
+ <boundary_condition>
+ <mesh>boundary8_right</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>T_ic</parameter>
+ </boundary_condition>
+
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+
+
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_newton</name>
+ <type>Newton</type>
+ <max_iter>50</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i bicgstab -p ilu -tol 1e-16 -maxiter 10000</lis>
+ <!-- <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-25</error_tolerance>
+ </eigen>
+-->
+
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ <!--
+ <solver_type>PardisoLU</solver_type>
+ -->
+ </eigen>
+
+ <!--
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <scaling>true</scaling>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-26</error_tolerance>
+ </eigen>
+-->
+ <petsc>
+ <prefix>sd</prefix>
+ <parameters>-sd_ksp_type cg -sd_pc_type bjacobi -sd_ksp_rtol 1e-16 -sd_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_bottom.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_bottom.vtu
new file mode 100644
index 00000000000..618529cff1c
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_bottom.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="50" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="232" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="344" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1024" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1352" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1464" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _qAAAADIAAAAAAAAACgAAAAAAAAAsAAAAAAAAAAgAAAAAAAAAKQAAAAAAAAAHAAAAAAAAACYAAAAAAAAAIwAAAAAAAAAAAAAAAAAAAC8AAAAAAAAAAQAAAAAAAAAGAAAAAAAAAB0AAAAAAAAAFwAAAAAAAAACAAAAAAAAAAMAAAAAAAAAGgAAAAAAAAAEAAAAAAAAAAkAAAAAAAAAIAAAAAAAAAAFAAAAAAAAAA==UAAAAAAAAAAAAAAAAQAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAA+AEAAGZmZmZmZu4/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOg/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAAAAAAAAAAAAAM7MzMzMzOQ/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmeE/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzNw/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADQzMzMzM+s/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAAAAAAAAANA/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAAAAAAAAAAAAADQzMzMzM9M/AAAAAAAAAAAAAAAAAAAAADQzMzMzM8M/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZtY/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAAA==8AAAAAgAAAAAAAAACgAAAAAAAAANAAAAAAAAAAoAAAAAAAAADgAAAAAAAAAQAAAAAAAAAA4AAAAAAAAADwAAAAAAAAAMAAAAAAAAAA8AAAAAAAAAEQAAAAAAAAATAAAAAAAAABEAAAAAAAAAFAAAAAAAAAAHAAAAAAAAABQAAAAAAAAACwAAAAAAAAAGAAAAAAAAAAsAAAAAAAAABQAAAAAAAAAEAAAAAAAAAAUAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAMAAAAAAAAAEgAAAAAAAAAJAAAAAAAAABIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==UAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAACgAAABUVFRUVFRUVFRU=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_left.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_left.vtu
new file mode 100644
index 00000000000..93ef579c517
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_left.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="3" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="22" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="40" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="56" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="160" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="200" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="216" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _GAAAABYAAAAAAAAAAAAAAAAAAAALAAAAAAAAAA==CAAAAAAAAAAAAAAASAAAAAAAAAAAAAAAAAAAAAAA4D8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAA==GAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==CAAAAAMAAAAAAAAAAQAAABU=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_right.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_right.vtu
new file mode 100644
index 00000000000..457ab170ac8
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_right.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="3" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="10" RangeMax="52" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="9" RangeMax="9" offset="40" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1.4142135624" offset="56" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="160" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="200" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="216" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _GAAAADQAAAAAAAAAFQAAAAAAAAAKAAAAAAAAAA==CAAAAAkAAAAAAAAASAAAAAAAAAAAAPA/AAAAAAAA4D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAA==GAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==CAAAAAMAAAAAAAAAAQAAABU=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_top.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_top.vtu
new file mode 100644
index 00000000000..be942cffb0c
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/boundary8_top.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="21" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="11" RangeMax="51" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="232" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1.4142135624" offset="344" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1024" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1352" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1464" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _qAAAABgAAAAAAAAACwAAAAAAAAAbAAAAAAAAAAwAAAAAAAAAHgAAAAAAAAAhAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAkAAAAAAAAAA8AAAAAAAAAFQAAAAAAAAAUAAAAAAAAADMAAAAAAAAAKgAAAAAAAAATAAAAAAAAADAAAAAAAAAAEgAAAAAAAAAtAAAAAAAAABEAAAAAAAAAJwAAAAAAAAAQAAAAAAAAAA==UAAAAAkAAAAAAAAACAAAAAAAAAAHAAAAAAAAAAYAAAAAAAAABQAAAAAAAAAEAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAABAAAAAAAAAAAAAAAAAAAA+AEAAJqZmZmZmak/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAADQzMzMzM8M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D8AAAAAAAAAAAAAAAAAANA/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZtY/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzNw/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA8D8AAAAAAAAAAGZmZmZmZu4/AAAAAAAA8D8AAAAAAAAAAM7MzMzMzOQ/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+s/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOg/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmeE/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA8D8AAAAAAAAAAA==8AAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAAsAAAAAAAAADgAAAAAAAAAPAAAAAAAAAA4AAAAAAAAAEAAAAAAAAAARAAAAAAAAABAAAAAAAAAAEgAAAAAAAAANAAAAAAAAABIAAAAAAAAAFAAAAAAAAAATAAAAAAAAABQAAAAAAAAACQAAAAAAAAAIAAAAAAAAAAkAAAAAAAAABgAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAEAAAAAAAAAAcAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAMAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==UAAAAAMAAAAAAAAABgAAAAAAAAAJAAAAAAAAAAwAAAAAAAAADwAAAAAAAAASAAAAAAAAABUAAAAAAAAAGAAAAAAAAAAbAAAAAAAAAB4AAAAAAAAACgAAABUVFRUVFRUVFRU=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/domain.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/domain.vtu
new file mode 100644
index 00000000000..765deaf04ba
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/domain.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="52" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="572" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="684" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2388" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="3248" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="3360" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _qAEAAAAAAAAAAAAAAQAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAACgAAAAAAAAALAAAAAAAAAAwAAAAAAAAADQAAAAAAAAAOAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAARAAAAAAAAABIAAAAAAAAAEwAAAAAAAAAUAAAAAAAAABUAAAAAAAAAFgAAAAAAAAAXAAAAAAAAABgAAAAAAAAAGQAAAAAAAAAaAAAAAAAAABsAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAB4AAAAAAAAAHwAAAAAAAAAgAAAAAAAAACEAAAAAAAAAIgAAAAAAAAAjAAAAAAAAACQAAAAAAAAAJQAAAAAAAAAmAAAAAAAAACcAAAAAAAAAKAAAAAAAAAApAAAAAAAAACoAAAAAAAAAKwAAAAAAAAAsAAAAAAAAAC0AAAAAAAAALgAAAAAAAAAvAAAAAAAAADAAAAAAAAAAMQAAAAAAAAAyAAAAAAAAADMAAAAAAAAANAAAAAAAAAA=UAAAAAAAAAAAAAAAAQAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAA+AQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAAAAAAAAAAAAADQzMzMzM9M/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA4D8AAAAAAAAAAJqZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA4D8AAAAAAAAAADQzMzMzM8M/AAAAAAAAAAAAAAAAAAAAADQzMzMzM8M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA4D8AAAAAAAAAAAAAAAAAANA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAANA/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA4D8AAAAAAAAAAGdmZmZmZtY/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZtY/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA4D8AAAAAAAAAAM3MzMzMzNw/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzNw/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA4D8AAAAAAAAAAJqZmZmZmeE/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmeE/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA4D8AAAAAAAAAAM7MzMzMzOQ/AAAAAAAAAAAAAAAAAAAAAM7MzMzMzOQ/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA4D8AAAAAAAAAAAAAAAAAAOg/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOg/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA4D8AAAAAAAAAADQzMzMzM+s/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+s/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA4D8AAAAAAAAAAGZmZmZmZu4/AAAAAAAAAAAAAAAAAAAAAGZmZmZmZu4/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA4D8AAAAAAAAAAA==gAIAAAAAAAAAAAAAAQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAAFwAAAAAAAAAZAAAAAAAAABgAAAAAAAAAFgAAAAAAAAABAAAAAAAAAAIAAAAAAAAADQAAAAAAAAAMAAAAAAAAABoAAAAAAAAAHAAAAAAAAAAbAAAAAAAAABkAAAAAAAAAAgAAAAAAAAADAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAdAAAAAAAAAB8AAAAAAAAAHgAAAAAAAAAcAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAAIAAAAAAAAAAiAAAAAAAAACEAAAAAAAAAHwAAAAAAAAAEAAAAAAAAAAUAAAAAAAAAEAAAAAAAAAAPAAAAAAAAACMAAAAAAAAAJQAAAAAAAAAkAAAAAAAAACIAAAAAAAAABQAAAAAAAAAGAAAAAAAAABEAAAAAAAAAEAAAAAAAAAAmAAAAAAAAACgAAAAAAAAAJwAAAAAAAAAlAAAAAAAAAAYAAAAAAAAABwAAAAAAAAASAAAAAAAAABEAAAAAAAAAKQAAAAAAAAArAAAAAAAAACoAAAAAAAAAKAAAAAAAAAAHAAAAAAAAAAgAAAAAAAAAEwAAAAAAAAASAAAAAAAAACwAAAAAAAAALgAAAAAAAAAtAAAAAAAAACsAAAAAAAAACAAAAAAAAAAJAAAAAAAAABQAAAAAAAAAEwAAAAAAAAAvAAAAAAAAADEAAAAAAAAAMAAAAAAAAAAuAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAVAAAAAAAAABQAAAAAAAAAMgAAAAAAAAA0AAAAAAAAADMAAAAAAAAAMQAAAAAAAAA=UAAAAAgAAAAAAAAAEAAAAAAAAAAYAAAAAAAAACAAAAAAAAAAKAAAAAAAAAAwAAAAAAAAADgAAAAAAAAAQAAAAAAAAABIAAAAAAAAAFAAAAAAAAAACgAAABcXFxcXFxcXFxc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet.pvd b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet.pvd
new file mode 100644
index 00000000000..4e9d0973785
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet.pvd
@@ -0,0 +1,10 @@
+<?xml version="1.0"?>
+<VTKFile type="Collection" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
+ <Collection>
+ <DataSet timestep="0" group="" part="0" file="result_TH2M_T_dirichlet_ts_0_t_0.000000.vtu"/>
+ <DataSet timestep="10" group="" part="0" file="result_TH2M_T_dirichlet_ts_10_t_10.000000.vtu"/>
+ <DataSet timestep="200" group="" part="0" file="result_TH2M_T_dirichlet_ts_20_t_200.000000.vtu"/>
+ <DataSet timestep="2000" group="" part="0" file="result_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu"/>
+ <DataSet timestep="4000" group="" part="0" file="result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu"/>
+ </Collection>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_0_t_0.000000.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_0_t_0.000000.vtu
new file mode 100644
index 00000000000..ca3a18c9263
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_0_t_0.000000.vtu
@@ -0,0 +1,44 @@
+<?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="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="0" RangeMax="0" offset="92" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="156" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="15000" offset="224" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="15000" RangeMax="15000" offset="296" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="368" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="436" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="512" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="100000" offset="620" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="100000" RangeMax="100000" offset="692" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="20" RangeMax="20" offset="764" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="85000" RangeMax="85000" offset="868" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.3" RangeMax="0.3" offset="940" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="1036" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="1140" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="293.15" offset="1216" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15" RangeMax="293.15" offset="1288" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="1.1118626977e-18" RangeMax="1.1986176893e-17" offset="1360" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="5.6231417867e-21" RangeMax="1.5765461726e-19" offset="2540" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="3700" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="3772" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="4096" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="4328" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="4420" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9A1MjXQTU8zSjQ1TDJNM9VLySwqqQQAU3wHOg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAADgAAAAAAAAA=eF5jYBgFQwEAAAGoAAE=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYAACn7MODEOJHoEAALN+Hac=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYAACn7MODKP0oKYBnvRHbg==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALwAAAAAAAAA=eF77/R8E3tszMoDAB4ppVij9HWouE5RPKc0JpZnRaAYoDbN/uIn/g4YjAGQNUO0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACCrB8ODEOJHoEAAOXeI+0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYACCrB8ODKP0oKYBC7dWiw==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALQAAAAAAAAA=eF77+x8EjB1YGEDAhGKaC0r/hprLBOVTSsPMZUajGaE0wzAV/wsNRwDY8ipdAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5jYGBgaDjwxYFhlB7UNAD+poIFAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAJgAAAAAAAAA=eF4zMgaBy/ZmVKKtobQelDanEu2AxreA0iZo9HATN4LSAOEteKs=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAKwAAAAAAAAA=eF7L8f40x8XvuX0xlegqKJ1FpHpi6RooXQKlS6F0ARo93MRzoDQAsPTUqg==AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5LSwOCoCKHtCFCM4xQAAANBUItAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5LSwOCoCKHtFF6UNMAeGqfaw==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAUwMAAAAAAAA=eF49knlM0wccxQdqU1HCXYI6N+qBx9wsk9LBl5+ztBAsVfpri6SUQmlpGY6BOHBs4BAQ2RQVpFgOKx54xA7LJYjfOfmp4EVgpiMOhEIsTmHLQLmsCWyJxPff++993nszYicubU4JxyOK7t3t3o6RCw5IevtVQPZcrOy0hkKzS4Y594kGQu5zzc63uFA/2MeOBC0cZhI3Hax8cL4juySf/AgspaRY/kIGcqr3d56FDepF/xZOJAcD7WXGqTPffAU2Kmd85pwU0jJ6jPpkFXxi+7XLw3crcOMYYz/FxoOD9egd3T4pDFb1fr12KhaiFN6sFvs4SNZzcs8+I2GZbWXW4t0imHyXF3mmDkO4KQKvBf3FSp5UIvnkD3rlrBAVbhpiRbkKmf8kWpxNfJSvCHowYlCgV0m7/GFpJKavfKyjFsZgWRI192cFgVuurFMUTGiw/UhT/E6RHbJ/kAwUrIrHRx43rk+zJSDb0FesD4hBTuX5UgZ3IzpecOkdKpRjFF3ZFUV+jk4/WkraCAV+W+1myiFYYBd/dVXgczX6HKSX5TH48MH/aq/78gvb5jNFw3wSZ+f7buseqjlQ5o1v53myytty2E9paI3o1HYM++HHCxl7hM22wOv1rzS1mAAjEquxgynATJ/fvM2jSszwy9fmjgjA8e1nlFs5G/UJ7vSLh+yASn2zbdeEFsx1TYFjs64omFsarQ5JwIdlDErvI4Cqaf90brgfRlFe5J69HlgbcXOjNF8D2p6/xwPMJPoGxA0srSPxGaclu1EtBdHlZYY3n67GmlRjjXMfC1POc5Rz+WrkORTl7TV6ocbu7FowqrBxh6OW/lgM1Pixk6/s10CFrmtn640QcJBe1XVDEhT36xuqO7bi4nceq07IM9fVO8H+MG2QaXoLtlaM1fk3cGBJNlGTVp0ILVlNwdFEd2Af97Xwe3cJrn/08zHXTaHwUtSvY4bxMGX58xecTj+gLE8Ts9PlUDA6pQwtYuJ04VBS9GYVGtxFvpzTIhgUl+bf9eLhoTxuyul+4j3vLpfDnpc9PnzP+52MRisQiKHWOBOyr4KFhDazWHh7G6TP/y+tcNb1Xhcf7ef3/SU13HV/qxB8Fg0bpl57gn7UKXjJehn8B7OqkpY=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAQwMAAAAAAAA=eF41kntQCwAcxxniQqeWhGudHpujPK5S/EqiGOWtiIREWdfyajcymaK288hW5r1W5pG6jGrpl2Wpptp6LG21XkNemd3KHeZR3Jnvf9//vp/v90szF3GvVCjBZZrQ0ThGjRpWR6vdIjVk280iW/kgktpIaW9uvYMnSx3s1hdIYMjz2O215SZYJrO1yd6uBv5x+Wib+BHIvFQ3LonTD7MD3/is+22Gfts7DUslNSBSqDPEWR+gJuSX+76+B0DjnbIVfm8EE6EuSeQbAze6+SdTJsphnYCqoAjegskwQCZZiVHNdswoDtPCNmnHSLN4EN1cjeLJ+jbY/i8v5qy6pKTEdqCV91GzE1+JLDdidXWaGKdvmnOQW6zHA9miqMBHTai80UldOPwchPXN9jX+FUj3SmTmuprAhhcTrHeWQVbbnbm7Q1thwOv9w0xqNpxYY6RsZavhl2RLeyG7FPIKWK7tyS3QPF6wqy2hAUXE99fe+mphdl5v3GJtAxTvdDeXpGqgICqXrgvQwKOhGXudM7rwDI/sQvRshVF/lX63zE9R150vD9FhqqXvsMs6qjRPACkWHj+Ll+YUVZ1JL8Pb2sXMFeRCMA9HRnfaaOAoX2eI6C1DOX1JfH6CEq/2rP5JYd/DsWk9fnClHD8tNwmjdzyG2M2Zcz1nDsKraU39jMC7cG7CSiuP4AHcynBkzk+VQhIvuqQ0uRe5zA6K6mkRcibQnOjhJtjA9lmje6FAVkJkeN6sIZhU+TlVQbgA5wn5XQLja0wLrv9G+VGK4YxJXyICxvtHDEYdDhLJQfV1tG9AqAFarMM+LVNVw/5yVfpNSRUy6/suTm5UgGRlmUPFSDt4ZHESC5NrsfifR7bG0MkgNKDskPYlqasH5zkZWzZd64SfPu4E6ZQvINjY1E1NVGB34movb/ZnmD71Hk1n/QzxZqjef78eHRQ/+O1KFaALN4Xk/gHJKTGc66dlmHFW4Cxz0CHT3hDEtJZC/IKPubaH9FilInD3vMuB/7y8Y/fTV1SWwH9ejziGW+OtXBg36voR+0odLjEP04M4Koi0/I8YS5u6OVyN9y37sh7E1TouKkGv8kqiUPwEMrdciOTs7IM/P/u41Q==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAEwAAAAAAAAA=eF4r8P40x8XvuX0BlWgAWQooCw==AQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA0AAAAAAAAAA=eF51kkEKwkAQBPMnvfub+Cb9QZ7g2VNAJSAIAQ8bIpKgSJ4QNoMMPZ3uY2Wp3fRMUazneMg57dZ5TXy7ybmL8y1xSxKenvi+zHkL/0j8esn5insnwfk79hB5TRx7iOdb4pYkPD1x7CH6R+LYQ7x3EvyfFPyV6L8iD/bmHvuvM3mQq57jvjTi3Q15cC7usT4f5EGu5uge6/lJHuRq7rHnTvTckQf3xD235d4XeZCrvYo9D6LngTy4h3HuH/IgV3vrnuWZ5Y88yNWeu2cG/AjT6w==AQAAAAAAAAAAgAAAAAAAAIACAAAAAAAAiwAAAAAAAAA=eF41zMtCgQEAhNGfUFSkkNxF7iG8/7O1cGY2ZzVfUdxX4guf+ckv9tljfmW+Mp0BRxwyvfwe2GQ6Y045YXr5VdhiOjN+c8708qvyjeks+MMl08uvxjbTWXHDNdPL75HvTGfLPXdML78nfjCdA4/8ZXr51dlhOif+8cz08muwy3QuvPHK9P4B1RkHbQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF7jYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TF4cBAAT7AOc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_10_t_10.000000.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_10_t_10.000000.vtu
new file mode 100644
index 00000000000..2436955558a
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_10_t_10.000000.vtu
@@ -0,0 +1,44 @@
+<?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="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-0" RangeMax="-0" offset="92" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149" RangeMax="-nan" offset="164" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="15000" offset="244" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="15000" RangeMax="15000" offset="316" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="532" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="100000" offset="640" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="100000" RangeMax="100000" offset="712" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="20" RangeMax="20" offset="784" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="85000" RangeMax="85000" offset="888" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.29999748" RangeMax="0.30000001729" offset="960" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="1284" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.3978998243e-19" RangeMax="10911.920088" offset="1388" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="333.15" offset="2292" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="292.87554191" RangeMax="333.15" offset="2480" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="1.1118626977e-18" RangeMax="1.1986176893e-17" offset="2812" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="5.6231417867e-21" RangeMax="1.5765461726e-19" offset="3992" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="5152" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5224" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5548" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5780" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5872" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9A1MjXQTU8zSjQ1TDJNM9VLySwqqQQAU3wHOg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACDBoYhRf+oH2k0AOdeOGo=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAGQAAAAAAAAA=eF5jYACDBgj1o55hlE8D/ig9UmgA7P2GDQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYAACn7MODEOJHoEAALN+Hac=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYAACn7MODKP0oKYBnvRHbg==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALwAAAAAAAAA=eF77/R8E3tszMoDAB4ppVij9HWouE5RPKc0JpZnRaAYoDbN/uIn/g4YjAGQNUO0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACCrB8ODEOJHoEAAOXeI+0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYACCrB8ODKP0oKYBC7dWiw==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALQAAAAAAAAA=eF77+x8EjB1YGEDAhGKaC0r/hprLBOVTSsPMZUajGaE0wzAV/wsNRwDY8ipdAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5jYGBgaDjwxYFhlB7UNAD+poIFAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA0AAAAAAAAAA=eF67cz9uoYbxZfuO6x3/9YH0UaELf4yAtJvwEVdjIC1z5KgRiD7dvA1Mnzrnbwyi1x02ANOXLhmB6QNJxmDa1hhC34Ka245mrjuauWeg5p5GM/cKmrn2UHMfQc3ddrs4QAdIb4HSbVB7VEt16gyBtAqUPgS1l13YSxGknw1Ku0Ddsem1iw2I3gClpaDu+qK8D0x/gNKnoO701mQGu8MDSh+HutvW0wGFXgP1h8NpQzBtB6UvQP11Qgrin+NQeh/Un7+9IPRPKG0E9TcA9XusEg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAKwAAAAAAAAA=eF7L8f40x8XvuX0xlegqKJ1FpHpi6RooXQKlS6F0ARo93MRzoDQAsPTUqg==AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAhQIAAAAAAAA=eF5d1F1IU2Ecx/Hj2UEs7MKEoDeCXDUpYpW20MXZhXmRFivtjS6K0F1EoNHFLiojs4TeLyp0L25FSBClpZXl9IyCKNNlFwYSBC6tJUVUY6bz7Gk75/cf7nluPvz3ZTf/c86jsuRx9wRVyKCAEyleqbz44A2Sk5D6/pamq1X3DgXt8ACkfttXV5Tv3WbzQ5qp5y50hZ2jvxQJxvJ0qTddWqyGjX8UXuqGwZKCcPuwLEIJUld93RfyO24pcTgLqTPWWPF9wTFlFiYgdVGIOrbcPyhnQwlSt4c2f7bcSJRed3tfVbVM9dqGHBNLpZ+91OPY5zQndS/22cZJfdil7zMEaab+0KXvswPSTP1ls0fbZzfsg9TPvl6i7bORk3rinb5PFdJMPYZ9RjmpS8I5bZ8iFNLSaahJ7dMAhbT6Gdiw82vbwHrrssKCmROTi6zrvnQc/7HCaBXmHg/2mTRr7pw8/eXj+yw5ncE+GIDUuzA/5qT+DM/jKezhno/jzPx/8uC1IC91GTMv9SOt+vN0wNrWzOfb8KYy+na0LHgSnoLp/w/p82HoGMrsJq/+PhRyUt9rmah+sLzYtoeT+hWTPl/mpL57h1t7n+xwF6S+dfsm5/uc3wov9fLzRdpMlkHqCZP+PvKme+fNSHXsk0IKjyAOW92izWQcUhe5+0Pk7o+EwTcy5uyXSUGCOCp+J+OGzJ7r178Hcp4/8/sQhFWdw+aLisRJNc6M2sw4qUdwX43DCHd/jbE8s6u+UvkGI5C6ypo3puYZOA2pi8Lz2tT3yEudsb/J+7BCztLuxQpZhNRn2Oma1MygCqmHjnbdbR9ZE/DESkqz602BJ+apuo931gao/wc2CqplAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAbAAAAAAAAAA=eF5LSwOCSyUOxxrNDrmKFDtM8tCe5hxf5KB418aCya/IoebzVo+GoCIHazn9J81AeoqL3qJ4IM069+HPNCC9X8G4FkTPUNaIBtEg4+A00NwjUHMnQ81VgppbBzXXloC5M9HMZRihAABqI0VYAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA1gAAAAAAAAA=eF5LSwOCSyUOxxrNDrmKFDtM8tCe5hxf5KB418aCya/IoebzVo+GoCIHazn9J81AeoqL3qJ4IM069+HPNCC9X8G4FkTPUNaIBtEg4+A00NwjUHMnQ81VgppbBzXXloC5M3GYK1biJxL2udhBFEofhdqj89RgzZGdRXAa5p9bQibPF4UhaJg7/DIqHRwDihwCoHQt1F13OorEmoD0PShtA3Vnxka9w4VINMzdZz6w+yUj0TB/PIrr2o1Mw/yls0g3B5mGhd9/l+MJIJrBFULD/A0Az7Gr0A==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAUwMAAAAAAAA=eF49knlM0wccxQdqU1HCXYI6N+qBx9wsk9LBl5+ztBAsVfpri6SUQmlpGY6BOHBs4BAQ2RQVpFgOKx54xA7LJYjfOfmp4EVgpiMOhEIsTmHLQLmsCWyJxPff++993nszYicubU4JxyOK7t3t3o6RCw5IevtVQPZcrOy0hkKzS4Y594kGQu5zzc63uFA/2MeOBC0cZhI3Hax8cL4juySf/AgspaRY/kIGcqr3d56FDepF/xZOJAcD7WXGqTPffAU2Kmd85pwU0jJ6jPpkFXxi+7XLw3crcOMYYz/FxoOD9egd3T4pDFb1fr12KhaiFN6sFvs4SNZzcs8+I2GZbWXW4t0imHyXF3mmDkO4KQKvBf3FSp5UIvnkD3rlrBAVbhpiRbkKmf8kWpxNfJSvCHowYlCgV0m7/GFpJKavfKyjFsZgWRI192cFgVuurFMUTGiw/UhT/E6RHbJ/kAwUrIrHRx43rk+zJSDb0FesD4hBTuX5UgZ3IzpecOkdKpRjFF3ZFUV+jk4/WkraCAV+W+1myiFYYBd/dVXgczX6HKSX5TH48MH/aq/78gvb5jNFw3wSZ+f7buseqjlQ5o1v53myytty2E9paI3o1HYM++HHCxl7hM22wOv1rzS1mAAjEquxgynATJ/fvM2jSszwy9fmjgjA8e1nlFs5G/UJ7vSLh+yASn2zbdeEFsx1TYFjs64omFsarQ5JwIdlDErvI4Cqaf90brgfRlFe5J69HlgbcXOjNF8D2p6/xwPMJPoGxA0srSPxGaclu1EtBdHlZYY3n67GmlRjjXMfC1POc5Rz+WrkORTl7TV6ocbu7FowqrBxh6OW/lgM1Pixk6/s10CFrmtn640QcJBe1XVDEhT36xuqO7bi4nceq07IM9fVO8H+MG2QaXoLtlaM1fk3cGBJNlGTVp0ILVlNwdFEd2Af97Xwe3cJrn/08zHXTaHwUtSvY4bxMGX58xecTj+gLE8Ts9PlUDA6pQwtYuJ04VBS9GYVGtxFvpzTIhgUl+bf9eLhoTxuyul+4j3vLpfDnpc9PnzP+52MRisQiKHWOBOyr4KFhDazWHh7G6TP/y+tcNb1Xhcf7ef3/SU13HV/qxB8Fg0bpl57gn7UKXjJehn8B7OqkpY=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAQwMAAAAAAAA=eF41kntQCwAcxxniQqeWhGudHpujPK5S/EqiGOWtiIREWdfyajcymaK288hW5r1W5pG6jGrpl2Wpptp6LG21XkNemd3KHeZR3Jnvf9//vp/v90szF3GvVCjBZZrQ0ThGjRpWR6vdIjVk280iW/kgktpIaW9uvYMnSx3s1hdIYMjz2O215SZYJrO1yd6uBv5x+Wib+BHIvFQ3LonTD7MD3/is+22Gfts7DUslNSBSqDPEWR+gJuSX+76+B0DjnbIVfm8EE6EuSeQbAze6+SdTJsphnYCqoAjegskwQCZZiVHNdswoDtPCNmnHSLN4EN1cjeLJ+jbY/i8v5qy6pKTEdqCV91GzE1+JLDdidXWaGKdvmnOQW6zHA9miqMBHTai80UldOPwchPXN9jX+FUj3SmTmuprAhhcTrHeWQVbbnbm7Q1thwOv9w0xqNpxYY6RsZavhl2RLeyG7FPIKWK7tyS3QPF6wqy2hAUXE99fe+mphdl5v3GJtAxTvdDeXpGqgICqXrgvQwKOhGXudM7rwDI/sQvRshVF/lX63zE9R150vD9FhqqXvsMs6qjRPACkWHj+Ll+YUVZ1JL8Pb2sXMFeRCMA9HRnfaaOAoX2eI6C1DOX1JfH6CEq/2rP5JYd/DsWk9fnClHD8tNwmjdzyG2M2Zcz1nDsKraU39jMC7cG7CSiuP4AHcynBkzk+VQhIvuqQ0uRe5zA6K6mkRcibQnOjhJtjA9lmje6FAVkJkeN6sIZhU+TlVQbgA5wn5XQLja0wLrv9G+VGK4YxJXyICxvtHDEYdDhLJQfV1tG9AqAFarMM+LVNVw/5yVfpNSRUy6/suTm5UgGRlmUPFSDt4ZHESC5NrsfifR7bG0MkgNKDskPYlqasH5zkZWzZd64SfPu4E6ZQvINjY1E1NVGB34movb/ZnmD71Hk1n/QzxZqjef78eHRQ/+O1KFaALN4Xk/gHJKTGc66dlmHFW4Cxz0CHT3hDEtJZC/IKPubaH9FilInD3vMuB/7y8Y/fTV1SWwH9ejziGW+OtXBg36voR+0odLjEP04M4Koi0/I8YS5u6OVyN9y37sh7E1TouKkGv8kqiUPwEMrdciOTs7IM/P/u41Q==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAEwAAAAAAAAA=eF4r8P40x8XvuX0BlWgAWQooCw==AQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA0AAAAAAAAAA=eF51kkEKwkAQBPMnvfub+Cb9QZ7g2VNAJSAIAQ8bIpKgSJ4QNoMMPZ3uY2Wp3fRMUazneMg57dZ5TXy7ybmL8y1xSxKenvi+zHkL/0j8esn5insnwfk79hB5TRx7iOdb4pYkPD1x7CH6R+LYQ7x3EvyfFPyV6L8iD/bmHvuvM3mQq57jvjTi3Q15cC7usT4f5EGu5uge6/lJHuRq7rHnTvTckQf3xD235d4XeZCrvYo9D6LngTy4h3HuH/IgV3vrnuWZ5Y88yNWeu2cG/AjT6w==AQAAAAAAAAAAgAAAAAAAAIACAAAAAAAAiwAAAAAAAAA=eF41zMtCgQEAhNGfUFSkkNxF7iG8/7O1cGY2ZzVfUdxX4guf+ckv9tljfmW+Mp0BRxwyvfwe2GQ6Y045YXr5VdhiOjN+c8708qvyjeks+MMl08uvxjbTWXHDNdPL75HvTGfLPXdML78nfjCdA4/8ZXr51dlhOif+8cz08muwy3QuvPHK9P4B1RkHbQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF7jYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TF4cBAAT7AOc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_20_t_200.000000.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_20_t_200.000000.vtu
new file mode 100644
index 00000000000..dc651e7596e
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_20_t_200.000000.vtu
@@ -0,0 +1,44 @@
+<?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="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-0" RangeMax="-0" offset="92" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149" RangeMax="-nan" offset="164" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="15000" offset="244" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="15000" RangeMax="15000" offset="316" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="532" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="100000" offset="640" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="100000" RangeMax="100000" offset="712" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="20" RangeMax="20" offset="784" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="85000" RangeMax="85000" offset="888" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.29999748" RangeMax="0.3" offset="960" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="1312" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.5499675644e-14" RangeMax="10911.920088" offset="1416" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="333.15" offset="2180" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15" RangeMax="333.15" offset="2392" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="1.1118626977e-18" RangeMax="1.1986176893e-17" offset="2788" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="5.6231417867e-21" RangeMax="1.5765461726e-19" offset="3968" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="5128" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5200" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5524" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5756" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5848" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9A1MjXQTU8zSjQ1TDJNM9VLySwqqQQAU3wHOg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACDBoYhRf+oH2k0AOdeOGo=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAGQAAAAAAAAA=eF5jYACDBgj1o55hlE8D/ig9UmgA7P2GDQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYAACn7MODEOJHoEAALN+Hac=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYAACn7MODKP0oKYBnvRHbg==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALwAAAAAAAAA=eF77/R8E3tszMoDAB4ppVij9HWouE5RPKc0JpZnRaAYoDbN/uIn/g4YjAGQNUO0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACCrB8ODEOJHoEAAOXeI+0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYACCrB8ODKP0oKYBC7dWiw==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALQAAAAAAAAA=eF77+x8EjB1YGEDAhGKaC0r/hprLBOVTSsPMZUajGaE0wzAV/wsNRwDY8ipdAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5jYGBgaDjwxYFhlB7UNAD+poIFAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA5QAAAAAAAAA=eF67dz9uoYbxZXu1LwprtYD0DMXH03SAdAOfmosekO7lqF+mD6QZDFj3GADpI00lPYZAWtSXRdUISAvVHJgCohPX+r8C0TbGIHDZ/i4Bc3ug5jJBzT2MZq4A1NwkqLmOUHOfQM1lfmW/XBNIM0JpFag997mYFmmD7IfSU6H28sxoydUF0pxQug7qjrbuoK8guhVKd0Pd5SrjZAhylxuU/qcPcWfyTRsVkDsToPRBqLvfpO+5AaJfQmkhqD+mpDyKAdEToDQf1F87RDr2g+htUDoe6k/fgoN8IH/6QGkjqL8BAgqpow==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAKwAAAAAAAAA=eF7L8f40x8XvuX0xlegqKJ1FpHpi6RooXQKlS6F0ARo93MRzoDQAsPTUqg==AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAGwIAAAAAAAA=eF5l1E1I0wEYx3EZWjClNyizQ+yikGGhWehoXTp46NDLSEQ7RPQ/GB7Mi6CZMKOhUATVQBpNTExcbb6gDPFvDxLDiELUwERSiRKGL5c12bIt6Xl+Ix53+fDjix6e/39LpHY+L0KUEH+LGfKJj2R+cc8PE0yI6PfsVXsbpgLUIt4X0Zc7fro74r2kRbd6wlea3V0Es0X0nOkbk1djHoLZIrpxZNN/fvIxadELW+58f9TXRieU6GWXrKu9+xrJLmKjb72eneo+ZxCMiejT9VkXflisJrvHtBVFOm1lG+PocbnnlhhT952Qe2rRnXJPeE3d96HcU4s+9JzvqUU/LvfUoj/I43tq03/v4ntq093J99SiTwT5nlr0Y+FIyuYsNtkSs3a93N7UV2Kip/6/pzdEFhE9o8YVuWgfpV3KZ6ma9zcRG/2rPA+4oJ6Pt7aq7ubLQdKi35WtRb8tz1OLPnDLV5Fc7qeAEt2QrUUPyvsQUKI3Pqlcc9Er0qLnytaiL8j7pEU/lIwW91i8BA+I6ON/eGvRD8r7qEX3FRXmhE4+JS163SneWvT1XH6f10Rs9Pcf3lQYZ9tJi/7xE28telcrfx98Inb6/3dvVq6cbiUteo+ftxY9/zJ/nwpEbPTO/WP523P1pEVvPsy7SYk+4+fvoxY9+u/38Dr9ErHRR4O8R0Rs9LF3nqTx7IyDLXWEF6NHP78tdaD/Bdlm2KI=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAfQAAAAAAAAA=eF5LSwOCSyUOYjtfFreElTionjMyXfuo2GFKhoW0am2xQ8ScbIlQzWIHHueacxueFzk8qM99nbytyKHKis3Ib1KRwz4Xh94H5UUOnCKN9xYnFzmAjEsLgtJAcwWh5ipDze2DmhuMZu4rqLldUHNPQ80VQDOXYYQCAMdxSkc=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAABgEAAAAAAAA=eF5LSwOCSyUOYjtfFreElTionjMyXfuo2GFKhoW0am2xQ8ScbIlQzWIHHueacxueFzk8qM99nbytyKHKis3Ib1KRwz4Xh94H5UUOnCKN9xYnFzmAjEsLgtJAcwWh5ipDze2DmhuMZu4rqLldUHNPQ80VwGHuPn6N3NIpJQ57oLQI1J55zn1XZsiUOMyC0ipQe2Nmm2q8XA/0D5SeCHVHGdMlSdvgYocCKB0Gddem98WfmziKHTZAaZg7yx5+ieE6X+RQDaWfQt2te3cnX+SSIgcrKN0E9ccc+8Xx4q1FDqug9BGov1LWPDA9lFvkkAeleaD+tNhbrMQaXeRgDaVh/gYAL8O1KQ==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAUwMAAAAAAAA=eF49knlM0wccxQdqU1HCXYI6N+qBx9wsk9LBl5+ztBAsVfpri6SUQmlpGY6BOHBs4BAQ2RQVpFgOKx54xA7LJYjfOfmp4EVgpiMOhEIsTmHLQLmsCWyJxPff++993nszYicubU4JxyOK7t3t3o6RCw5IevtVQPZcrOy0hkKzS4Y594kGQu5zzc63uFA/2MeOBC0cZhI3Hax8cL4juySf/AgspaRY/kIGcqr3d56FDepF/xZOJAcD7WXGqTPffAU2Kmd85pwU0jJ6jPpkFXxi+7XLw3crcOMYYz/FxoOD9egd3T4pDFb1fr12KhaiFN6sFvs4SNZzcs8+I2GZbWXW4t0imHyXF3mmDkO4KQKvBf3FSp5UIvnkD3rlrBAVbhpiRbkKmf8kWpxNfJSvCHowYlCgV0m7/GFpJKavfKyjFsZgWRI192cFgVuurFMUTGiw/UhT/E6RHbJ/kAwUrIrHRx43rk+zJSDb0FesD4hBTuX5UgZ3IzpecOkdKpRjFF3ZFUV+jk4/WkraCAV+W+1myiFYYBd/dVXgczX6HKSX5TH48MH/aq/78gvb5jNFw3wSZ+f7buseqjlQ5o1v53myytty2E9paI3o1HYM++HHCxl7hM22wOv1rzS1mAAjEquxgynATJ/fvM2jSszwy9fmjgjA8e1nlFs5G/UJ7vSLh+yASn2zbdeEFsx1TYFjs64omFsarQ5JwIdlDErvI4Cqaf90brgfRlFe5J69HlgbcXOjNF8D2p6/xwPMJPoGxA0srSPxGaclu1EtBdHlZYY3n67GmlRjjXMfC1POc5Rz+WrkORTl7TV6ocbu7FowqrBxh6OW/lgM1Pixk6/s10CFrmtn640QcJBe1XVDEhT36xuqO7bi4nceq07IM9fVO8H+MG2QaXoLtlaM1fk3cGBJNlGTVp0ILVlNwdFEd2Af97Xwe3cJrn/08zHXTaHwUtSvY4bxMGX58xecTj+gLE8Ts9PlUDA6pQwtYuJ04VBS9GYVGtxFvpzTIhgUl+bf9eLhoTxuyul+4j3vLpfDnpc9PnzP+52MRisQiKHWOBOyr4KFhDazWHh7G6TP/y+tcNb1Xhcf7ef3/SU13HV/qxB8Fg0bpl57gn7UKXjJehn8B7OqkpY=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAQwMAAAAAAAA=eF41kntQCwAcxxniQqeWhGudHpujPK5S/EqiGOWtiIREWdfyajcymaK288hW5r1W5pG6jGrpl2Wpptp6LG21XkNemd3KHeZR3Jnvf9//vp/v90szF3GvVCjBZZrQ0ThGjRpWR6vdIjVk280iW/kgktpIaW9uvYMnSx3s1hdIYMjz2O215SZYJrO1yd6uBv5x+Wib+BHIvFQ3LonTD7MD3/is+22Gfts7DUslNSBSqDPEWR+gJuSX+76+B0DjnbIVfm8EE6EuSeQbAze6+SdTJsphnYCqoAjegskwQCZZiVHNdswoDtPCNmnHSLN4EN1cjeLJ+jbY/i8v5qy6pKTEdqCV91GzE1+JLDdidXWaGKdvmnOQW6zHA9miqMBHTai80UldOPwchPXN9jX+FUj3SmTmuprAhhcTrHeWQVbbnbm7Q1thwOv9w0xqNpxYY6RsZavhl2RLeyG7FPIKWK7tyS3QPF6wqy2hAUXE99fe+mphdl5v3GJtAxTvdDeXpGqgICqXrgvQwKOhGXudM7rwDI/sQvRshVF/lX63zE9R150vD9FhqqXvsMs6qjRPACkWHj+Ll+YUVZ1JL8Pb2sXMFeRCMA9HRnfaaOAoX2eI6C1DOX1JfH6CEq/2rP5JYd/DsWk9fnClHD8tNwmjdzyG2M2Zcz1nDsKraU39jMC7cG7CSiuP4AHcynBkzk+VQhIvuqQ0uRe5zA6K6mkRcibQnOjhJtjA9lmje6FAVkJkeN6sIZhU+TlVQbgA5wn5XQLja0wLrv9G+VGK4YxJXyICxvtHDEYdDhLJQfV1tG9AqAFarMM+LVNVw/5yVfpNSRUy6/suTm5UgGRlmUPFSDt4ZHESC5NrsfifR7bG0MkgNKDskPYlqasH5zkZWzZd64SfPu4E6ZQvINjY1E1NVGB34movb/ZnmD71Hk1n/QzxZqjef78eHRQ/+O1KFaALN4Xk/gHJKTGc66dlmHFW4Cxz0CHT3hDEtJZC/IKPubaH9FilInD3vMuB/7y8Y/fTV1SWwH9ejziGW+OtXBg36voR+0odLjEP04M4Koi0/I8YS5u6OVyN9y37sh7E1TouKkGv8kqiUPwEMrdciOTs7IM/P/u41Q==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAEwAAAAAAAAA=eF4r8P40x8XvuX0BlWgAWQooCw==AQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA0AAAAAAAAAA=eF51kkEKwkAQBPMnvfub+Cb9QZ7g2VNAJSAIAQ8bIpKgSJ4QNoMMPZ3uY2Wp3fRMUazneMg57dZ5TXy7ybmL8y1xSxKenvi+zHkL/0j8esn5insnwfk79hB5TRx7iOdb4pYkPD1x7CH6R+LYQ7x3EvyfFPyV6L8iD/bmHvuvM3mQq57jvjTi3Q15cC7usT4f5EGu5uge6/lJHuRq7rHnTvTckQf3xD235d4XeZCrvYo9D6LngTy4h3HuH/IgV3vrnuWZ5Y88yNWeu2cG/AjT6w==AQAAAAAAAAAAgAAAAAAAAIACAAAAAAAAiwAAAAAAAAA=eF41zMtCgQEAhNGfUFSkkNxF7iG8/7O1cGY2ZzVfUdxX4guf+ckv9tljfmW+Mp0BRxwyvfwe2GQ6Y045YXr5VdhiOjN+c8708qvyjeks+MMl08uvxjbTWXHDNdPL75HvTGfLPXdML78nfjCdA4/8ZXr51dlhOif+8cz08muwy3QuvPHK9P4B1RkHbQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF7jYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TF4cBAAT7AOc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu
new file mode 100644
index 00000000000..e141eddf91b
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_30_t_2000.000000.vtu
@@ -0,0 +1,44 @@
+<?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="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-0" RangeMax="-0" offset="92" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149" RangeMax="-nan" offset="164" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="15000" offset="244" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="15000" RangeMax="15000" offset="316" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="532" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="100000" offset="640" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="100000" RangeMax="100000" offset="712" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="20" RangeMax="20" offset="784" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="85000" RangeMax="85000" offset="888" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.29999748" RangeMax="0.3" offset="960" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="1308" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="4.9595887094e-14" RangeMax="10911.920088" offset="1412" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="333.15" offset="2172" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15" RangeMax="333.15" offset="2380" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="1.1118626977e-18" RangeMax="1.1986176893e-17" offset="2768" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="5.6231417867e-21" RangeMax="1.5765461726e-19" offset="3948" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="5108" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5180" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5504" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5736" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5828" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9A1MjXQTU8zSjQ1TDJNM9VLySwqqQQAU3wHOg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACDBoYhRf+oH2k0AOdeOGo=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAGQAAAAAAAAA=eF5jYACDBgj1o55hlE8D/ig9UmgA7P2GDQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYAACn7MODEOJHoEAALN+Hac=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYAACn7MODKP0oKYBnvRHbg==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALwAAAAAAAAA=eF77/R8E3tszMoDAB4ppVij9HWouE5RPKc0JpZnRaAYoDbN/uIn/g4YjAGQNUO0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACCrB8ODEOJHoEAAOXeI+0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYACCrB8ODKP0oKYBC7dWiw==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALQAAAAAAAAA=eF77+x8EjB1YGEDAhGKaC0r/hprLBOVTSsPMZUajGaE0wzAV/wsNRwDY8ipdAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5jYGBgaDjwxYFhlB7UNAD+poIFAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA4gAAAAAAAAA=eF67cz9uoYbxZXtt670bNIG0wHP2A1pAeoe87nltIJ14Ru6eDpBe/ODyG10g/c3O87cekJ7E1sRpAKR/ubSLGwLpL0LpqkZA2t4YBC7b34Kaq4PD3CQc5k6GmvsbzVwHqLkPoeZ28PJpgsxth9IaUHtkJiZZgOyRhtI8UHtTmme5g+xNgNLboO7g/rY0DOQODiidAHVXQ1hlKshd9VB6EdSdJ/r5SkDuPA6lP0HdfXhRapM+kN4PpSdA/XFiacsEZPoH1F/ftMvngejPUPoT1J8iyr5rQLQwlDaC+hsAy9ev8A==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAKwAAAAAAAAA=eF7L8f40x8XvuX0xlegqKJ1FpHpi6RooXQKlS6F0ARo93MRzoDQAsPTUqg==AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAGQIAAAAAAAA=eF511E9I03EYx3GDGCwI7DCWAxFhBblDf9ymIEIReFBEvUQRTKkZYZTSpkIURVFQ0hb4B6MifyuGuCL6Q5uwnz4Hgy4dokMeLFSoxBAFqU2RLPw+H8EP9Lu8eHgfBs/3+11u/d/3ICM5Na8W6HenwM74vG/lNolu5bM1N5+9FphQ0ZeKxh3BwEth0R1VxQP25efCos8O7ardd29YZtRpFb3t886U/fSJsOi+ynRoMvpIYJmKfvjM45Ar0SssetDTXhSybkiFihn9wNErKxFfoT1oja2Nurx2Ybb829RQqY2+qvtc+c9+23SfLPot3SeLPqf7ZNGdus8dqpP2+1v3yaJ7dZ8sej5o9smix8NmnzAW3rrfEbfZ5zCJXv5+z/2J9kbb2GRfX4x+Kn7XuLnfje9hRraRSHXXQlOulrSw6AGd/WqAepeeR6capfMZWMtePT73RmC/it6pM4se0/OEcTrfOudY8kvHK6kl0Q/pDA9Sn9f7wKJ3pz987Ft+ISz6WZ1Z9O16n1j03Hl3R8SREvhLRf+jM4v+U+8ji743G78U8SSFRXfqzKLX6H1m0cdH/Qt9+y1h0et1ZtFjFeY9sOiVDamvXT2Dwm7+fr2ZWfRwq3lPp9VTrVvf14nlY+cmLt4VeFJFb14yM4v+Y7d5j99J9MDG/2FE/CR6wm1mS8WMnlw90nOhtKEatnjWm0tGmqrR/wJ1l6VkAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAewAAAAAAAAA=eF5LSwOCSyUOIqUGV8wmlTjIXHj7kjuoxKFE1vf2S8EShy+R9/9evljscMG2q/TMxGKHV4ve/7scWOzgUs5696VgsUPJsaWvuS8VOYj+tb9qNqnIAWRcWhCUBporBDVXCmpuEQ5zn0HNtYaamwk1lwfNXIYRCgDZ1FXKAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAAgEAAAAAAAA=eF5LSwOCSyUOIqUGV8wmlTjIXHj7kjuoxKFE1vf2S8EShy+R9/9evljscMG2q/TMxGKHV4ve/7scWOzgUs5696VgsUPJsaWvuS8VOYj+tb9qNqnIAWRcWhCUBporBDVXCmpuEQ5zn0HNtYaamwk1lweHuXtzvef6bSpx2AOlhaH2zFDij1csKnGYBqWlofaeKJv76JdhicMxKF0MdccW62lv7n0sdtgMpWHuenR6i+WFjQga5s7YDzZWFwqLHaKh9Auou8V7f729Z1jsIAClHaD+uDMv9Mmvj0UOl6B0HtRfRx4WJShuKnLYBaUFoP7cZ3h5rl9RkcNOKA3zNwA3k9PsAQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAUwMAAAAAAAA=eF49knlM0wccxQdqU1HCXYI6N+qBx9wsk9LBl5+ztBAsVfpri6SUQmlpGY6BOHBs4BAQ2RQVpFgOKx54xA7LJYjfOfmp4EVgpiMOhEIsTmHLQLmsCWyJxPff++993nszYicubU4JxyOK7t3t3o6RCw5IevtVQPZcrOy0hkKzS4Y594kGQu5zzc63uFA/2MeOBC0cZhI3Hax8cL4juySf/AgspaRY/kIGcqr3d56FDepF/xZOJAcD7WXGqTPffAU2Kmd85pwU0jJ6jPpkFXxi+7XLw3crcOMYYz/FxoOD9egd3T4pDFb1fr12KhaiFN6sFvs4SNZzcs8+I2GZbWXW4t0imHyXF3mmDkO4KQKvBf3FSp5UIvnkD3rlrBAVbhpiRbkKmf8kWpxNfJSvCHowYlCgV0m7/GFpJKavfKyjFsZgWRI192cFgVuurFMUTGiw/UhT/E6RHbJ/kAwUrIrHRx43rk+zJSDb0FesD4hBTuX5UgZ3IzpecOkdKpRjFF3ZFUV+jk4/WkraCAV+W+1myiFYYBd/dVXgczX6HKSX5TH48MH/aq/78gvb5jNFw3wSZ+f7buseqjlQ5o1v53myytty2E9paI3o1HYM++HHCxl7hM22wOv1rzS1mAAjEquxgynATJ/fvM2jSszwy9fmjgjA8e1nlFs5G/UJ7vSLh+yASn2zbdeEFsx1TYFjs64omFsarQ5JwIdlDErvI4Cqaf90brgfRlFe5J69HlgbcXOjNF8D2p6/xwPMJPoGxA0srSPxGaclu1EtBdHlZYY3n67GmlRjjXMfC1POc5Rz+WrkORTl7TV6ocbu7FowqrBxh6OW/lgM1Pixk6/s10CFrmtn640QcJBe1XVDEhT36xuqO7bi4nceq07IM9fVO8H+MG2QaXoLtlaM1fk3cGBJNlGTVp0ILVlNwdFEd2Af97Xwe3cJrn/08zHXTaHwUtSvY4bxMGX58xecTj+gLE8Ts9PlUDA6pQwtYuJ04VBS9GYVGtxFvpzTIhgUl+bf9eLhoTxuyul+4j3vLpfDnpc9PnzP+52MRisQiKHWOBOyr4KFhDazWHh7G6TP/y+tcNb1Xhcf7ef3/SU13HV/qxB8Fg0bpl57gn7UKXjJehn8B7OqkpY=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAQwMAAAAAAAA=eF41kntQCwAcxxniQqeWhGudHpujPK5S/EqiGOWtiIREWdfyajcymaK288hW5r1W5pG6jGrpl2Wpptp6LG21XkNemd3KHeZR3Jnvf9//vp/v90szF3GvVCjBZZrQ0ThGjRpWR6vdIjVk280iW/kgktpIaW9uvYMnSx3s1hdIYMjz2O215SZYJrO1yd6uBv5x+Wib+BHIvFQ3LonTD7MD3/is+22Gfts7DUslNSBSqDPEWR+gJuSX+76+B0DjnbIVfm8EE6EuSeQbAze6+SdTJsphnYCqoAjegskwQCZZiVHNdswoDtPCNmnHSLN4EN1cjeLJ+jbY/i8v5qy6pKTEdqCV91GzE1+JLDdidXWaGKdvmnOQW6zHA9miqMBHTai80UldOPwchPXN9jX+FUj3SmTmuprAhhcTrHeWQVbbnbm7Q1thwOv9w0xqNpxYY6RsZavhl2RLeyG7FPIKWK7tyS3QPF6wqy2hAUXE99fe+mphdl5v3GJtAxTvdDeXpGqgICqXrgvQwKOhGXudM7rwDI/sQvRshVF/lX63zE9R150vD9FhqqXvsMs6qjRPACkWHj+Ll+YUVZ1JL8Pb2sXMFeRCMA9HRnfaaOAoX2eI6C1DOX1JfH6CEq/2rP5JYd/DsWk9fnClHD8tNwmjdzyG2M2Zcz1nDsKraU39jMC7cG7CSiuP4AHcynBkzk+VQhIvuqQ0uRe5zA6K6mkRcibQnOjhJtjA9lmje6FAVkJkeN6sIZhU+TlVQbgA5wn5XQLja0wLrv9G+VGK4YxJXyICxvtHDEYdDhLJQfV1tG9AqAFarMM+LVNVw/5yVfpNSRUy6/suTm5UgGRlmUPFSDt4ZHESC5NrsfifR7bG0MkgNKDskPYlqasH5zkZWzZd64SfPu4E6ZQvINjY1E1NVGB34movb/ZnmD71Hk1n/QzxZqjef78eHRQ/+O1KFaALN4Xk/gHJKTGc66dlmHFW4Cxz0CHT3hDEtJZC/IKPubaH9FilInD3vMuB/7y8Y/fTV1SWwH9ejziGW+OtXBg36voR+0odLjEP04M4Koi0/I8YS5u6OVyN9y37sh7E1TouKkGv8kqiUPwEMrdciOTs7IM/P/u41Q==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAEwAAAAAAAAA=eF4r8P40x8XvuX0BlWgAWQooCw==AQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA0AAAAAAAAAA=eF51kkEKwkAQBPMnvfub+Cb9QZ7g2VNAJSAIAQ8bIpKgSJ4QNoMMPZ3uY2Wp3fRMUazneMg57dZ5TXy7ybmL8y1xSxKenvi+zHkL/0j8esn5insnwfk79hB5TRx7iOdb4pYkPD1x7CH6R+LYQ7x3EvyfFPyV6L8iD/bmHvuvM3mQq57jvjTi3Q15cC7usT4f5EGu5uge6/lJHuRq7rHnTvTckQf3xD235d4XeZCrvYo9D6LngTy4h3HuH/IgV3vrnuWZ5Y88yNWeu2cG/AjT6w==AQAAAAAAAAAAgAAAAAAAAIACAAAAAAAAiwAAAAAAAAA=eF41zMtCgQEAhNGfUFSkkNxF7iG8/7O1cGY2ZzVfUdxX4guf+ckv9tljfmW+Mp0BRxwyvfwe2GQ6Y045YXr5VdhiOjN+c8708qvyjeks+MMl08uvxjbTWXHDNdPL75HvTGfLPXdML78nfjCdA4/8ZXr51dlhOif+8cz08muwy3QuvPHK9P4B1RkHbQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF7jYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TF4cBAAT7AOc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu
new file mode 100644
index 00000000000..effdaf4db0e
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/result_TH2M_T_dirichlet_ts_34_t_4000.000000.vtu
@@ -0,0 +1,44 @@
+<?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="53" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="-0" RangeMax="-0" offset="92" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149" RangeMax="-nan" offset="164" />
+ <DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="15000" offset="244" />
+ <DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="15000" RangeMax="15000" offset="316" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="388" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="456" />
+ <DataArray type="Float64" Name="gas_density" format="appended" RangeMin="1" RangeMax="1" offset="532" />
+ <DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="100000" offset="640" />
+ <DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="100000" RangeMax="100000" offset="712" />
+ <DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="20" RangeMax="20" offset="784" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="85000" RangeMax="85000" offset="888" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.29999748" RangeMax="0.3" offset="960" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="1308" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="1.5745564387e-13" RangeMax="10911.920088" offset="1412" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="333.15" offset="2180" />
+ <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15" RangeMax="333.15" offset="2388" />
+ <DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="1.1118626977e-18" RangeMax="1.1986176893e-17" offset="2764" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="5.6231417867e-21" RangeMax="1.5765461726e-19" offset="3944" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.72830420164" RangeMax="0.72830420164" offset="5104" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5176" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5500" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="5732" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="5824" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABoAAAAAAAAAIgAAAAAAAAA=eF4z0zPRM9A1MjXQTU8zSjQ1TDJNM9VLySwqqQQAU3wHOg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACDBoYhRf+oH2k0AOdeOGo=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAGQAAAAAAAAA=eF5jYACDBgj1o55hlE8D/ig9UmgA7P2GDQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYAACn7MODEOJHoEAALN+Hac=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYAACn7MODKP0oKYBnvRHbg==AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAEQAAAAAAAAA=eF5jYBgFo2AUUAsAAANQAAE=AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAFgAAAAAAAAA=eF5jYBgFo2AUjIJRMAoGHwAABqAAAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALwAAAAAAAAA=eF77/R8E3tszMoDAB4ppVij9HWouE5RPKc0JpZnRaAYoDbN/uIn/g4YjAGQNUO0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFAAAAAAAAAA=eF5jYACCrB8ODEOJHoEAAOXeI+0=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5jYACCrB8ODKP0oKYBC7dWiw==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAALQAAAAAAAAA=eF77+x8EjB1YGEDAhGKaC0r/hprLBOVTSsPMZUajGaE0wzAV/wsNRwDY8ipdAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAFQAAAAAAAAA=eF5jYGBgaDjwxYFhlB7UNAD+poIFAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA5AAAAAAAAAA=eF67ez9uoYbxZXulxq/rNYH0ZIaefVpA+oSt0lltIP1NbsdtHSA9d7LvK10gzT7t0Q89IP1Cu5zdAKQulVvMEEinx85XMQLSdsYgcNn+Npq5U3CYOw9qLgfU3Fc4zHWCmvsYam6zwSoNkLkNUFoBas/dAw7mIHtuQekJUHt15K+7guzVgtJHoO6I08sNAbkjGkp/gbrL6yZTMshdHlB6NtSdF0VmFILceR5Ks0Hd/f2BToM+yF9Q+hnUH9c9DvaB6GtQ+ijUX9MSQ+eC6ElQOhXqT1+Pl6tBtDeUNoT6GwA0ZrLQAQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAKwAAAAAAAAA=eF7L8f40x8XvuX0xlegqKJ1FpHpi6RooXQKlS6F0ARo93MRzoDQAsPTUqg==AQAAAAAAAAAAgAAAAAAAAKAGAAAAAAAAHQIAAAAAAAA=eF511F1IU3EYx/EMekHIitidLwOVCrowRoHKeiEKEVRYaHaVSUiE5nwptCtDcBAIQlTiDBJBhIpD9ibUoQe66UKwl4PNbIbKJvhyZ5trma3+zy/YDzoX+/DwPTd7zvmf+GbqCo5LTF1Xt+jV6rjf/Cx8ISz6xOe1iYT1VFj0LFfmtXjxE9mt7lHRa2OPosuBx8Ki7xwY+7oUHJXtJPpoOHW/NSws+vnv+51w4L7Ukuh3Ih89i9ZtYdEHQtsSEatbgipm9P7E6cm34SO2sdR2erpOVQ+X2OhJ3ecPNUn7vaH77FQ7aL9Tuk8W/ajuk0Vv0X1CP+33hO6TRc+dNfuEebPp+y2Lm32y6DWLZp8s+oEZs08WvaHBNb1Q5LONZ+3nvZf2eXt8//ab8edncFy2kujHr/Rd+FX3UuAxFb1UZ1hC3a/Pw/+f9z9nKvNMcumZZJPo+TrDAuqT+jzhe3q+zdc3H6y3jwm8qqLf1Bl2UXfr+8Cin+yru7y6YQmL3qQzi16v7xOL7nlnXVze+1DgYRX9nM6whrpb30cW/dOr1P2FI+KomNEPvTYziz6v3wsW/WBj6v8WDwmLPqczi14ZM+cBVsXSz0dn00b3fLBfYIeKXtFsZha9PWrOE2yLpp+vbx+qi74EegXOqegRx8wsemjanEcW/d7f72GbwLsqet6MmXNVzOgVt0KRHYkqL/TUl+9acfm86L8B7VodAQ==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAeQAAAAAAAAA=eF5LSwOCSyUOVkWnBNMmlTjsNdA9mhpU4nCn07ojVajEIeAmT2LqpWKHcq+u4NRJxQ58N4D8oGKH4naQfLHDMl2g+ktFDsvzQfqLHEDGpQVBaaC51lBz90PNvQ81NxBqbiXUXBGoufVQczdBzV2PZi7DCAUAH0VIRg==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAA+AAAAAAAAAA=eF5LSwOCSyUOVkWnBNMmlTjsNdA9mhpU4nCn07ojVajEIeAmT2LqpWKHcq+u4NRJxQ58N4D8oGKH4naQfLHDMl2g+ktFDsvzQfqLHEDGpQVBaaC51lBz90PNvQ81NxBqbiXUXBGoufVQczdBzV2Pw9yAHAmbtE0IGub+movVr1OLShxqofQ+qL0+sVuWpRqVOPhD6XtQd0wzXFKS+qnYYQaUhvn3yUTvqNRNxQ5PoXQF1J1OIH5RsYMblBaEuvuAPlC/UbHDKShdCfVHTxTQvk9FDjOg9Bqov5adA7pvU5HDBii9GupPtmygf4qKHLigNMzfAL1JtB0=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAUwMAAAAAAAA=eF49knlM0wccxQdqU1HCXYI6N+qBx9wsk9LBl5+ztBAsVfpri6SUQmlpGY6BOHBs4BAQ2RQVpFgOKx54xA7LJYjfOfmp4EVgpiMOhEIsTmHLQLmsCWyJxPff++993nszYicubU4JxyOK7t3t3o6RCw5IevtVQPZcrOy0hkKzS4Y594kGQu5zzc63uFA/2MeOBC0cZhI3Hax8cL4juySf/AgspaRY/kIGcqr3d56FDepF/xZOJAcD7WXGqTPffAU2Kmd85pwU0jJ6jPpkFXxi+7XLw3crcOMYYz/FxoOD9egd3T4pDFb1fr12KhaiFN6sFvs4SNZzcs8+I2GZbWXW4t0imHyXF3mmDkO4KQKvBf3FSp5UIvnkD3rlrBAVbhpiRbkKmf8kWpxNfJSvCHowYlCgV0m7/GFpJKavfKyjFsZgWRI192cFgVuurFMUTGiw/UhT/E6RHbJ/kAwUrIrHRx43rk+zJSDb0FesD4hBTuX5UgZ3IzpecOkdKpRjFF3ZFUV+jk4/WkraCAV+W+1myiFYYBd/dVXgczX6HKSX5TH48MH/aq/78gvb5jNFw3wSZ+f7buseqjlQ5o1v53myytty2E9paI3o1HYM++HHCxl7hM22wOv1rzS1mAAjEquxgynATJ/fvM2jSszwy9fmjgjA8e1nlFs5G/UJ7vSLh+yASn2zbdeEFsx1TYFjs64omFsarQ5JwIdlDErvI4Cqaf90brgfRlFe5J69HlgbcXOjNF8D2p6/xwPMJPoGxA0srSPxGaclu1EtBdHlZYY3n67GmlRjjXMfC1POc5Rz+WrkORTl7TV6ocbu7FowqrBxh6OW/lgM1Pixk6/s10CFrmtn640QcJBe1XVDEhT36xuqO7bi4nceq07IM9fVO8H+MG2QaXoLtlaM1fk3cGBJNlGTVp0ILVlNwdFEd2Af97Xwe3cJrn/08zHXTaHwUtSvY4bxMGX58xecTj+gLE8Ts9PlUDA6pQwtYuJ04VBS9GYVGtxFvpzTIhgUl+bf9eLhoTxuyul+4j3vLpfDnpc9PnzP+52MRisQiKHWOBOyr4KFhDazWHh7G6TP/y+tcNb1Xhcf7ef3/SU13HV/qxB8Fg0bpl57gn7UKXjJehn8B7OqkpY=AQAAAAAAAAAAgAAAAAAAAFADAAAAAAAAQwMAAAAAAAA=eF41kntQCwAcxxniQqeWhGudHpujPK5S/EqiGOWtiIREWdfyajcymaK288hW5r1W5pG6jGrpl2Wpptp6LG21XkNemd3KHeZR3Jnvf9//vp/v90szF3GvVCjBZZrQ0ThGjRpWR6vdIjVk280iW/kgktpIaW9uvYMnSx3s1hdIYMjz2O215SZYJrO1yd6uBv5x+Wib+BHIvFQ3LonTD7MD3/is+22Gfts7DUslNSBSqDPEWR+gJuSX+76+B0DjnbIVfm8EE6EuSeQbAze6+SdTJsphnYCqoAjegskwQCZZiVHNdswoDtPCNmnHSLN4EN1cjeLJ+jbY/i8v5qy6pKTEdqCV91GzE1+JLDdidXWaGKdvmnOQW6zHA9miqMBHTai80UldOPwchPXN9jX+FUj3SmTmuprAhhcTrHeWQVbbnbm7Q1thwOv9w0xqNpxYY6RsZavhl2RLeyG7FPIKWK7tyS3QPF6wqy2hAUXE99fe+mphdl5v3GJtAxTvdDeXpGqgICqXrgvQwKOhGXudM7rwDI/sQvRshVF/lX63zE9R150vD9FhqqXvsMs6qjRPACkWHj+Ll+YUVZ1JL8Pb2sXMFeRCMA9HRnfaaOAoX2eI6C1DOX1JfH6CEq/2rP5JYd/DsWk9fnClHD8tNwmjdzyG2M2Zcz1nDsKraU39jMC7cG7CSiuP4AHcynBkzk+VQhIvuqQ0uRe5zA6K6mkRcibQnOjhJtjA9lmje6FAVkJkeN6sIZhU+TlVQbgA5wn5XQLja0wLrv9G+VGK4YxJXyICxvtHDEYdDhLJQfV1tG9AqAFarMM+LVNVw/5yVfpNSRUy6/suTm5UgGRlmUPFSDt4ZHESC5NrsfifR7bG0MkgNKDskPYlqasH5zkZWzZd64SfPu4E6ZQvINjY1E1NVGB34movb/ZnmD71Hk1n/QzxZqjef78eHRQ/+O1KFaALN4Xk/gHJKTGc66dlmHFW4Cxz0CHT3hDEtJZC/IKPubaH9FilInD3vMuB/7y8Y/fTV1SWwH9ejziGW+OtXBg36voR+0odLjEP04M4Koi0/I8YS5u6OVyN9y37sh7E1TouKkGv8kqiUPwEMrdciOTs7IM/P/u41Q==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAEwAAAAAAAAA=eF4r8P40x8XvuX0BlWgAWQooCw==AQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA0AAAAAAAAAA=eF51kkEKwkAQBPMnvfub+Cb9QZ7g2VNAJSAIAQ8bIpKgSJ4QNoMMPZ3uY2Wp3fRMUazneMg57dZ5TXy7ybmL8y1xSxKenvi+zHkL/0j8esn5insnwfk79hB5TRx7iOdb4pYkPD1x7CH6R+LYQ7x3EvyfFPyV6L8iD/bmHvuvM3mQq57jvjTi3Q15cC7usT4f5EGu5uge6/lJHuRq7rHnTvTckQf3xD235d4XeZCrvYo9D6LngTy4h3HuH/IgV3vrnuWZ5Y88yNWeu2cG/AjT6w==AQAAAAAAAAAAgAAAAAAAAIACAAAAAAAAiwAAAAAAAAA=eF41zMtCgQEAhNGfUFSkkNxF7iG8/7O1cGY2ZzVfUdxX4guf+ckv9tljfmW+Mp0BRxwyvfwe2GQ6Y045YXr5VdhiOjN+c8708qvyjeks+MMl08uvxjbTWXHDNdPL75HvTGfLPXdML78nfjCdA4/8ZXr51dlhOif+8cz08muwy3QuvPHK9P4B1RkHbQ==AQAAAAAAAAAAgAAAAAAAAFAAAAAAAAAAIwAAAAAAAAA=eF7jYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=AQAAAAAAAAAAgAAAAAAAAAoAAAAAAAAACwAAAAAAAAA=eF4TF4cBAAT7AOc=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/square.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/square.vtu
new file mode 100644
index 00000000000..6c9a931359a
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/square.vtu
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="22" NumberOfCells="10" >
+ <PointData>
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="0" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="712" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1144" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1256" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _EAIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAAAAAAAAAAAAADQzMzMzM9M/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAAAAAA==QAEAAAAAAAAAAAAAAQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAAAQAAAAAAAAACAAAAAAAAAA0AAAAAAAAADAAAAAAAAAACAAAAAAAAAAMAAAAAAAAADgAAAAAAAAANAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAABAAAAAAAAAAFAAAAAAAAABAAAAAAAAAADwAAAAAAAAAFAAAAAAAAAAYAAAAAAAAAEQAAAAAAAAAQAAAAAAAAAAYAAAAAAAAABwAAAAAAAAASAAAAAAAAABEAAAAAAAAABwAAAAAAAAAIAAAAAAAAABMAAAAAAAAAEgAAAAAAAAAIAAAAAAAAAAkAAAAAAAAAFAAAAAAAAAATAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAVAAAAAAAAABQAAAAAAAAAUAAAAAQAAAAAAAAACAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAACgAAAAAAAAACgAAAAkJCQkJCQkJCQk=
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/T/T_1d_dirichlet/square_quad8.vtu b/Tests/Data/TH2M/T/T_1d_dirichlet/square_quad8.vtu
new file mode 100644
index 00000000000..1122f60f752
--- /dev/null
+++ b/Tests/Data/TH2M/T/T_1d_dirichlet/square_quad8.vtu
@@ -0,0 +1,22 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="53" NumberOfCells="10" >
+ <PointData>
+ </PointData>
+ <CellData>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="0" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1704" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2564" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2676" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _+AQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAAAAAAAAAAAAADQzMzMzM9M/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAAAAAAAAAAAA4D8AAAAAAAAAAJqZmZmZmak/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmbk/AAAAAAAA4D8AAAAAAAAAADQzMzMzM8M/AAAAAAAAAAAAAAAAAAAAADQzMzMzM8M/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmck/AAAAAAAA4D8AAAAAAAAAAAAAAAAAANA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAANA/AAAAAAAA8D8AAAAAAAAAADQzMzMzM9M/AAAAAAAA4D8AAAAAAAAAAGdmZmZmZtY/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZtY/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmdk/AAAAAAAA4D8AAAAAAAAAAM3MzMzMzNw/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzNw/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAOA/AAAAAAAA4D8AAAAAAAAAAJqZmZmZmeE/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmeE/AAAAAAAA8D8AAAAAAAAAADQzMzMzM+M/AAAAAAAA4D8AAAAAAAAAAM7MzMzMzOQ/AAAAAAAAAAAAAAAAAAAAAM7MzMzMzOQ/AAAAAAAA8D8AAAAAAAAAAGdmZmZmZuY/AAAAAAAA4D8AAAAAAAAAAAAAAAAAAOg/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOg/AAAAAAAA8D8AAAAAAAAAAJqZmZmZmek/AAAAAAAA4D8AAAAAAAAAADQzMzMzM+s/AAAAAAAAAAAAAAAAAAAAADQzMzMzM+s/AAAAAAAA8D8AAAAAAAAAAM3MzMzMzOw/AAAAAAAA4D8AAAAAAAAAAGZmZmZmZu4/AAAAAAAAAAAAAAAAAAAAAGZmZmZmZu4/AAAAAAAA8D8AAAAAAAAAAAAAAAAAAPA/AAAAAAAA4D8AAAAAAAAAAA==gAIAAAAAAAAAAAAAAQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAAFwAAAAAAAAAZAAAAAAAAABgAAAAAAAAAFgAAAAAAAAABAAAAAAAAAAIAAAAAAAAADQAAAAAAAAAMAAAAAAAAABoAAAAAAAAAHAAAAAAAAAAbAAAAAAAAABkAAAAAAAAAAgAAAAAAAAADAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAdAAAAAAAAAB8AAAAAAAAAHgAAAAAAAAAcAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAAIAAAAAAAAAAiAAAAAAAAACEAAAAAAAAAHwAAAAAAAAAEAAAAAAAAAAUAAAAAAAAAEAAAAAAAAAAPAAAAAAAAACMAAAAAAAAAJQAAAAAAAAAkAAAAAAAAACIAAAAAAAAABQAAAAAAAAAGAAAAAAAAABEAAAAAAAAAEAAAAAAAAAAmAAAAAAAAACgAAAAAAAAAJwAAAAAAAAAlAAAAAAAAAAYAAAAAAAAABwAAAAAAAAASAAAAAAAAABEAAAAAAAAAKQAAAAAAAAArAAAAAAAAACoAAAAAAAAAKAAAAAAAAAAHAAAAAAAAAAgAAAAAAAAAEwAAAAAAAAASAAAAAAAAACwAAAAAAAAALgAAAAAAAAAtAAAAAAAAACsAAAAAAAAACAAAAAAAAAAJAAAAAAAAABQAAAAAAAAAEwAAAAAAAAAvAAAAAAAAADEAAAAAAAAAMAAAAAAAAAAuAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAVAAAAAAAAABQAAAAAAAAAMgAAAAAAAAA0AAAAAAAAADMAAAAAAAAAMQAAAAAAAAA=UAAAAAgAAAAAAAAAEAAAAAAAAAAYAAAAAAAAACAAAAAAAAAAKAAAAAAAAAAwAAAAAAAAADgAAAAAAAAAQAAAAAAAAABIAAAAAAAAAFAAAAAAAAAACgAAABcXFxcXFxcXFxc=
+ </AppendedData>
+</VTKFile>
--
GitLab