diff --git a/ProcessLib/TH2M/Tests.cmake b/ProcessLib/TH2M/Tests.cmake
index ea7d68981fe3a102d2993751c6a9be87ce3f1dd4..f71fe0f127fde1a5b12972c5da407a5ca6ccb300 100644
--- a/ProcessLib/TH2M/Tests.cmake
+++ b/ProcessLib/TH2M/Tests.cmake
@@ -8,6 +8,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE TH2M/THM/Confined_Compression/THM_confined_compression_gas.prj RUNTIME 55)
OgsTest(PROJECTFILE TH2M/THM/Confined_Compression/THM_confined_compression_liquid.prj RUNTIME 55)
OgsTest(PROJECTFILE TH2M/TH/idealGasLaw/compression_gas.prj RUNTIME 1)
+ OgsTest(PROJECTFILE TH2M/H2M/Liakopoulos/liakopoulos_TH2M.prj RUNTIME 15)
endif()
# TH2M 1d heat diffusion w/ Dirichlet-BC
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/boundary_bottom.vtu b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_bottom.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c340e6425f717b5fbae401835f2313b2cfb901f1
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_bottom.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="2" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="1" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="32" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.1" offset="56" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="132" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="164" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="188" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _EAAAAAAAAAABAAAAAAAAAAAAAAAAAAAACAAAAAAAAAAAAAAAAAAAAA==MAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA=EAAAAAAAAAABAAAAAAAAAAAAAAAAAAAACAAAAAAAAAACAAAAAAAAAA==AQAAAAAAAAAD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/boundary_left.vtu b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_left.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..56d44b8d67bc9f69b19e189c6479487ac46b7392
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_left.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="200" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="1088" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="2168" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5412" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7556" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="8636" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _KAMAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAAFQAAAAAAAAAVgAAAAAAAABYAAAAAAAAAFoAAAAAAAAAXAAAAAAAAABeAAAAAAAAAGAAAAAAAAAAYgAAAAAAAABkAAAAAAAAAGYAAAAAAAAAaAAAAAAAAABqAAAAAAAAAGwAAAAAAAAAbgAAAAAAAAAAAAAAAAAAAHAAAAAAAAAAcgAAAAAAAAB0AAAAAAAAAHYAAAAAAAAAeAAAAAAAAAB6AAAAAAAAAHwAAAAAAAAAfgAAAAAAAAAQAAAAAAAAAIAAAAAAAAAAEgAAAAAAAACCAAAAAAAAABQAAAAAAAAAhAAAAAAAAAAWAAAAAAAAAIYAAAAAAAAAGAAAAAAAAACIAAAAAAAAABoAAAAAAAAAigAAAAAAAAAcAAAAAAAAAIwAAAAAAAAAHgAAAAAAAACOAAAAAAAAACAAAAAAAAAAkAAAAAAAAAAiAAAAAAAAAJIAAAAAAAAAJAAAAAAAAACUAAAAAAAAACYAAAAAAAAAlgAAAAAAAAAoAAAAAAAAAJgAAAAAAAAAKgAAAAAAAACaAAAAAAAAACwAAAAAAAAAnAAAAAAAAAAuAAAAAAAAAJ4AAAAAAAAAMAAAAAAAAACgAAAAAAAAADIAAAAAAAAAogAAAAAAAAA6AAAAAAAAAKoAAAAAAAAAPAAAAAAAAACsAAAAAAAAAD4AAAAAAAAArgAAAAAAAABAAAAAAAAAALAAAAAAAAAAQgAAAAAAAACyAAAAAAAAAEQAAAAAAAAAtAAAAAAAAABGAAAAAAAAALYAAAAAAAAASAAAAAAAAAC4AAAAAAAAAEoAAAAAAAAAugAAAAAAAABMAAAAAAAAALwAAAAAAAAATgAAAAAAAAC+AAAAAAAAAFAAAAAAAAAAwAAAAAAAAABSAAAAAAAAAMIAAAAAAAAANAAAAAAAAACkAAAAAAAAAMQAAAAAAAAANgAAAAAAAACmAAAAAAAAAMYAAAAAAAAAOAAAAAAAAACoAAAAAAAAAMgAAAAAAAAAIAMAAAAAAABjAAAAAAAAAGIAAAAAAAAAYQAAAAAAAABgAAAAAAAAAF8AAAAAAAAAXgAAAAAAAABdAAAAAAAAAFwAAAAAAAAAWwAAAAAAAABaAAAAAAAAAFkAAAAAAAAAWAAAAAAAAABXAAAAAAAAAFYAAAAAAAAAVQAAAAAAAABUAAAAAAAAAFMAAAAAAAAAUgAAAAAAAABRAAAAAAAAAFAAAAAAAAAATwAAAAAAAABOAAAAAAAAAE0AAAAAAAAATAAAAAAAAABLAAAAAAAAAEoAAAAAAAAASQAAAAAAAABIAAAAAAAAAEcAAAAAAAAARgAAAAAAAABFAAAAAAAAAEQAAAAAAAAAQwAAAAAAAABCAAAAAAAAAEEAAAAAAAAAQAAAAAAAAAA/AAAAAAAAAD4AAAAAAAAAPQAAAAAAAAA8AAAAAAAAADsAAAAAAAAAOgAAAAAAAAA5AAAAAAAAADgAAAAAAAAANwAAAAAAAAA2AAAAAAAAADUAAAAAAAAANAAAAAAAAAAzAAAAAAAAADIAAAAAAAAAMQAAAAAAAAAwAAAAAAAAAC8AAAAAAAAALgAAAAAAAAAtAAAAAAAAACwAAAAAAAAAKwAAAAAAAAAqAAAAAAAAACkAAAAAAAAAKAAAAAAAAAAnAAAAAAAAACYAAAAAAAAAJQAAAAAAAAAkAAAAAAAAACMAAAAAAAAAIgAAAAAAAAAhAAAAAAAAACAAAAAAAAAAHwAAAAAAAAAeAAAAAAAAAB0AAAAAAAAAHAAAAAAAAAAbAAAAAAAAABoAAAAAAAAAGQAAAAAAAAAYAAAAAAAAABcAAAAAAAAAFgAAAAAAAAAVAAAAAAAAABQAAAAAAAAAEwAAAAAAAAASAAAAAAAAABEAAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA4AAAAAAAAADQAAAAAAAAAMAAAAAAAAAAsAAAAAAAAACgAAAAAAAAAJAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAUAAAAAAAAABAAAAAAAAAADAAAAAAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==eAkAAAAAAAAAAAAAAAAAAHsUrkfheoQ/AAAAAAAAAAAAAAAAAAAAAHsUrkfhepQ/AAAAAAAAAAAAAAAAAAAAALgehetRuJ4/AAAAAAAAAAAAAAAAAAAAAHsUrkfheqQ/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAAAAAAAAAAAAALgehetRuK4/AAAAAAAAAAAAAAAAAAAAAOxRuB6F67E/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4do/AAAAAAAAAAAAAAAAAAAAAIXrUbgehds/AAAAAAAAAAAAAAAAAAAAAClcj8L1KNw/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzNw/AAAAAAAAAAAAAAAAAAAAAHE9CtejcN0/AAAAAAAAAAAAAAAAAAAAABWuR+F6FN4/AAAAAAAAAAAAAAAAAAAAALgehetRuN4/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXN8/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUeA/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo+A/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9eA/AAAAAAAAAAAAAAAAAAAAAEjhehSuR+E/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmeE/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOxRuB6F6+E/AAAAAAAAAAAAAAAAAAAAAD4K16NwPeI/AAAAAAAAAAAAAAAAAAAAAI/C9Shcj+I/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4eI/AAAAAAAAAAAAAAAAAAAAADMzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAIXrUbgeheM/AAAAAAAAAAAAAAAAAAAAANejcD0K1+M/AAAAAAAAAAAAAAAAAAAAAClcj8L1KOQ/AAAAAAAAAAAAAAAAAAAAAHsUrkfherQ/AAAAAAAAAAAAAAAAAAAAAHsUrkfheuQ/AAAAAAAAAAAAAAAAAAAAAArXo3A9Crc/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOQ/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAB+F61G4HuU/AAAAAAAAAAAAAAAAAAAAAClcj8L1KLw/AAAAAAAAAAAAAAAAAAAAAHE9CtejcOU/AAAAAAAAAAAAAAAAAAAAALgehetRuL4/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwuU/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo8A/AAAAAAAAAAAAAAAAAAAAABWuR+F6FOY/AAAAAAAAAAAAAAAAAAAAAOxRuB6F68E/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAAAAAAAAAAAAADMzMzMzM8M/AAAAAAAAAAAAAAAAAAAAALgehetRuOY/AAAAAAAAAAAAAAAAAAAAAHsUrkfhesQ/AAAAAAAAAAAAAAAAAAAAAArXo3A9Cuc/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwsU/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXOc/AAAAAAAAAAAAAAAAAAAAAArXo3A9Csc/AAAAAAAAAAAAAAAAAAAAAK5H4XoUruc/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUcg/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOg/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUeg/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4co/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo+g/AAAAAAAAAAAAAAAAAAAAAClcj8L1KMw/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9eg/AAAAAAAAAAAAAAAAAAAAAHE9CtejcM0/AAAAAAAAAAAAAAAAAAAAAEjhehSuR+k/AAAAAAAAAAAAAAAAAAAAALgehetRuM4/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAAAAAAAAAAAAAAAAAAAAANA/AAAAAAAAAAAAAAAAAAAAAOxRuB6F6+k/AAAAAAAAAAAAAAAAAAAAAI/C9Shcj9I/AAAAAAAAAAAAAAAAAAAAADMzMzMzM+s/AAAAAAAAAAAAAAAAAAAAADMzMzMzM9M/AAAAAAAAAAAAAAAAAAAAAIXrUbgehes/AAAAAAAAAAAAAAAAAAAAANejcD0K19M/AAAAAAAAAAAAAAAAAAAAANejcD0K1+s/AAAAAAAAAAAAAAAAAAAAAHsUrkfhetQ/AAAAAAAAAAAAAAAAAAAAAClcj8L1KOw/AAAAAAAAAAAAAAAAAAAAAB+F61G4HtU/AAAAAAAAAAAAAAAAAAAAAHsUrkfheuw/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwtU/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZtY/AAAAAAAAAAAAAAAAAAAAAB+F61G4Hu0/AAAAAAAAAAAAAAAAAAAAAArXo3A9Ctc/AAAAAAAAAAAAAAAAAAAAAHE9CtejcO0/AAAAAAAAAAAAAAAAAAAAAK5H4XoUrtc/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwu0/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUdg/AAAAAAAAAAAAAAAAAAAAABWuR+F6FO4/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9dg/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZu4/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAAAAAAAAAAAAALgehetRuO4/AAAAAAAAAAAAAAAAAAAAAD4K16NwPdo/AAAAAAAAAAAAAAAAAAAAAArXo3A9Cu8/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo9A/AAAAAAAAAAAAAAAAAAAAAD4K16NwPeo/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXO8/AAAAAAAAAAAAAAAAAAAAAEjhehSuR9E/AAAAAAAAAAAAAAAAAAAAAJDC9Shcj+o/AAAAAAAAAAAAAAAAAAAAAK5H4XoUru8/AAAAAAAAAAAAAAAAAAAAAOxRuB6F69E/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4eo/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAA=QAYAAAAAAABkAAAAAAAAAGEAAAAAAAAAYQAAAAAAAABeAAAAAAAAAF4AAAAAAAAAWwAAAAAAAABbAAAAAAAAAFkAAAAAAAAAWQAAAAAAAABXAAAAAAAAAFcAAAAAAAAAVQAAAAAAAABVAAAAAAAAAFMAAAAAAAAAUwAAAAAAAABRAAAAAAAAAFEAAAAAAAAATwAAAAAAAABPAAAAAAAAAE0AAAAAAAAATQAAAAAAAABLAAAAAAAAAEsAAAAAAAAASQAAAAAAAABJAAAAAAAAAEcAAAAAAAAARwAAAAAAAABFAAAAAAAAAEUAAAAAAAAAQwAAAAAAAABDAAAAAAAAAGMAAAAAAAAAYwAAAAAAAABgAAAAAAAAAGAAAAAAAAAAXQAAAAAAAABdAAAAAAAAAEEAAAAAAAAAQQAAAAAAAAA/AAAAAAAAAD8AAAAAAAAAPQAAAAAAAAA9AAAAAAAAADsAAAAAAAAAOwAAAAAAAAA5AAAAAAAAADkAAAAAAAAANwAAAAAAAAA3AAAAAAAAADUAAAAAAAAANQAAAAAAAAAzAAAAAAAAADMAAAAAAAAAMQAAAAAAAAAxAAAAAAAAAC8AAAAAAAAALwAAAAAAAAAtAAAAAAAAAC0AAAAAAAAAKwAAAAAAAAArAAAAAAAAACkAAAAAAAAAKQAAAAAAAAAnAAAAAAAAACcAAAAAAAAAJQAAAAAAAAAlAAAAAAAAACMAAAAAAAAAIwAAAAAAAAAhAAAAAAAAACEAAAAAAAAAHwAAAAAAAAAfAAAAAAAAAB0AAAAAAAAAHQAAAAAAAAAcAAAAAAAAABwAAAAAAAAAGwAAAAAAAAAbAAAAAAAAABoAAAAAAAAAGgAAAAAAAAAZAAAAAAAAABkAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABcAAAAAAAAAFwAAAAAAAAAWAAAAAAAAABYAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABMAAAAAAAAAEwAAAAAAAAASAAAAAAAAABIAAAAAAAAAEQAAAAAAAAARAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA8AAAAAAAAADgAAAAAAAAAOAAAAAAAAAA0AAAAAAAAADQAAAAAAAAAMAAAAAAAAAAwAAAAAAAAACwAAAAAAAAALAAAAAAAAAAoAAAAAAAAACgAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAcAAAAAAAAABwAAAAAAAABaAAAAAAAAAFoAAAAAAAAAWAAAAAAAAABYAAAAAAAAAFYAAAAAAAAAVgAAAAAAAABUAAAAAAAAAFQAAAAAAAAAUgAAAAAAAABSAAAAAAAAAFAAAAAAAAAAUAAAAAAAAABOAAAAAAAAAE4AAAAAAAAATAAAAAAAAABMAAAAAAAAAEoAAAAAAAAASgAAAAAAAABIAAAAAAAAAEgAAAAAAAAARgAAAAAAAABGAAAAAAAAAEQAAAAAAAAARAAAAAAAAABCAAAAAAAAAEIAAAAAAAAAYgAAAAAAAABiAAAAAAAAAF8AAAAAAAAAXwAAAAAAAABcAAAAAAAAAFwAAAAAAAAAQAAAAAAAAABAAAAAAAAAAD4AAAAAAAAAPgAAAAAAAAA8AAAAAAAAADwAAAAAAAAAOgAAAAAAAAA6AAAAAAAAADgAAAAAAAAAOAAAAAAAAAA2AAAAAAAAADYAAAAAAAAANAAAAAAAAAA0AAAAAAAAADIAAAAAAAAAMgAAAAAAAAAwAAAAAAAAADAAAAAAAAAALgAAAAAAAAAuAAAAAAAAACwAAAAAAAAALAAAAAAAAAAqAAAAAAAAACoAAAAAAAAAKAAAAAAAAAAoAAAAAAAAACYAAAAAAAAAJgAAAAAAAAAkAAAAAAAAACQAAAAAAAAAIgAAAAAAAAAiAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAeAAAAAAAAAB4AAAAAAAAABgAAAAAAAAAGAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAEAAAAAAAAAAQAAAAAAAAAAwAAAAAAAAADAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABUAAAAAAAAAIAMAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAABYAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABwAAAAAAAAAHgAAAAAAAAAgAAAAAAAAACIAAAAAAAAAJAAAAAAAAAAmAAAAAAAAACgAAAAAAAAAKgAAAAAAAAAsAAAAAAAAAC4AAAAAAAAAMAAAAAAAAAAyAAAAAAAAADQAAAAAAAAANgAAAAAAAAA4AAAAAAAAADoAAAAAAAAAPAAAAAAAAAA+AAAAAAAAAEAAAAAAAAAAQgAAAAAAAABEAAAAAAAAAEYAAAAAAAAASAAAAAAAAABKAAAAAAAAAEwAAAAAAAAATgAAAAAAAABQAAAAAAAAAFIAAAAAAAAAVAAAAAAAAABWAAAAAAAAAFgAAAAAAAAAWgAAAAAAAABcAAAAAAAAAF4AAAAAAAAAYAAAAAAAAABiAAAAAAAAAGQAAAAAAAAAZgAAAAAAAABoAAAAAAAAAGoAAAAAAAAAbAAAAAAAAABuAAAAAAAAAHAAAAAAAAAAcgAAAAAAAAB0AAAAAAAAAHYAAAAAAAAAeAAAAAAAAAB6AAAAAAAAAHwAAAAAAAAAfgAAAAAAAACAAAAAAAAAAIIAAAAAAAAAhAAAAAAAAACGAAAAAAAAAIgAAAAAAAAAigAAAAAAAACMAAAAAAAAAI4AAAAAAAAAkAAAAAAAAACSAAAAAAAAAJQAAAAAAAAAlgAAAAAAAACYAAAAAAAAAJoAAAAAAAAAnAAAAAAAAACeAAAAAAAAAKAAAAAAAAAAogAAAAAAAACkAAAAAAAAAKYAAAAAAAAAqAAAAAAAAACqAAAAAAAAAKwAAAAAAAAArgAAAAAAAACwAAAAAAAAALIAAAAAAAAAtAAAAAAAAAC2AAAAAAAAALgAAAAAAAAAugAAAAAAAAC8AAAAAAAAAL4AAAAAAAAAwAAAAAAAAADCAAAAAAAAAMQAAAAAAAAAxgAAAAAAAADIAAAAAAAAAA==ZAAAAAAAAAADAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/boundary_right.vtu b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_right.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b07c3e541b937d654cd512ca723d8e2f4eec4672
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_right.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="1" RangeMax="201" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="1088" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0.1" RangeMax="1.0049875621" offset="2168" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="5412" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7556" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="8636" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _KAMAAAAAAAB9AAAAAAAAAHsAAAAAAAAAeQAAAAAAAAB3AAAAAAAAAHUAAAAAAAAAcwAAAAAAAABxAAAAAAAAAG8AAAAAAAAAbQAAAAAAAABrAAAAAAAAAGkAAAAAAAAAZwAAAAAAAABlAAAAAAAAAGMAAAAAAAAAYQAAAAAAAABfAAAAAAAAAF0AAAAAAAAAWwAAAAAAAABZAAAAAAAAAFcAAAAAAAAAwwAAAAAAAABTAAAAAAAAAMEAAAAAAAAAUQAAAAAAAAC/AAAAAAAAAE8AAAAAAAAAvQAAAAAAAABNAAAAAAAAALsAAAAAAAAASwAAAAAAAAC5AAAAAAAAAEkAAAAAAAAAtwAAAAAAAABHAAAAAAAAALUAAAAAAAAARQAAAAAAAACzAAAAAAAAAEMAAAAAAAAAsQAAAAAAAABBAAAAAAAAAK8AAAAAAAAAPwAAAAAAAACtAAAAAAAAAD0AAAAAAAAAqwAAAAAAAAA7AAAAAAAAAMkAAAAAAAAAqQAAAAAAAAA5AAAAAAAAAMcAAAAAAAAApwAAAAAAAAA3AAAAAAAAAMUAAAAAAAAApQAAAAAAAAA1AAAAAAAAAKMAAAAAAAAAMwAAAAAAAAChAAAAAAAAADEAAAAAAAAAnwAAAAAAAAAvAAAAAAAAAJ0AAAAAAAAALQAAAAAAAACbAAAAAAAAACsAAAAAAAAAlwAAAAAAAAAnAAAAAAAAAJUAAAAAAAAAJQAAAAAAAACTAAAAAAAAACMAAAAAAAAAkQAAAAAAAAAhAAAAAAAAAI8AAAAAAAAAHwAAAAAAAACNAAAAAAAAAB0AAAAAAAAAiwAAAAAAAAAbAAAAAAAAAIkAAAAAAAAAGQAAAAAAAACHAAAAAAAAABcAAAAAAAAAhQAAAAAAAAAVAAAAAAAAAIMAAAAAAAAAEwAAAAAAAACBAAAAAAAAABEAAAAAAAAAfwAAAAAAAAAPAAAAAAAAAA0AAAAAAAAACwAAAAAAAAAJAAAAAAAAAAcAAAAAAAAABQAAAAAAAABVAAAAAAAAAAMAAAAAAAAAmQAAAAAAAAApAAAAAAAAAAEAAAAAAAAAIAMAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAA==eAkAAAAAAACamZmZmZm5P9ejcD0K1+M/AAAAAAAAAACamZmZmZm5P4XrUbgeheM/AAAAAAAAAACamZmZmZm5PzMzMzMzM+M/AAAAAAAAAACamZmZmZm5P+F6FK5H4eI/AAAAAAAAAACamZmZmZm5P4/C9Shcj+I/AAAAAAAAAACamZmZmZm5Pz4K16NwPeI/AAAAAAAAAACamZmZmZm5P+xRuB6F6+E/AAAAAAAAAACamZmZmZm5P5qZmZmZmeE/AAAAAAAAAACamZmZmZm5P0jhehSuR+E/AAAAAAAAAACamZmZmZm5P/YoXI/C9eA/AAAAAAAAAACamZmZmZm5P6RwPQrXo+A/AAAAAAAAAACamZmZmZm5P1K4HoXrUeA/AAAAAAAAAACamZmZmZm5PwAAAAAAAOA/AAAAAAAAAACamZmZmZm5P1yPwvUoXN8/AAAAAAAAAACamZmZmZm5P7gehetRuN4/AAAAAAAAAACamZmZmZm5PxWuR+F6FN4/AAAAAAAAAACamZmZmZm5P3E9CtejcN0/AAAAAAAAAACamZmZmZm5P83MzMzMzNw/AAAAAAAAAACamZmZmZm5Pylcj8L1KNw/AAAAAAAAAACamZmZmZm5P4XrUbgehds/AAAAAAAAAACamZmZmZm5PwrXo3A9Cu8/AAAAAAAAAACamZmZmZm5Pz4K16NwPdo/AAAAAAAAAACamZmZmZm5P7gehetRuO4/AAAAAAAAAACamZmZmZm5P5qZmZmZmdk/AAAAAAAAAACamZmZmZm5P2dmZmZmZu4/AAAAAAAAAACamZmZmZm5P/YoXI/C9dg/AAAAAAAAAACamZmZmZm5PxWuR+F6FO4/AAAAAAAAAACamZmZmZm5P1K4HoXrUdg/AAAAAAAAAACamZmZmZm5P8P1KFyPwu0/AAAAAAAAAACamZmZmZm5P65H4XoUrtc/AAAAAAAAAACamZmZmZm5P3E9CtejcO0/AAAAAAAAAACamZmZmZm5PwrXo3A9Ctc/AAAAAAAAAACamZmZmZm5Px+F61G4Hu0/AAAAAAAAAACamZmZmZm5P2dmZmZmZtY/AAAAAAAAAACamZmZmZm5P83MzMzMzOw/AAAAAAAAAACamZmZmZm5P8P1KFyPwtU/AAAAAAAAAACamZmZmZm5P3sUrkfheuw/AAAAAAAAAACamZmZmZm5Px+F61G4HtU/AAAAAAAAAACamZmZmZm5Pylcj8L1KOw/AAAAAAAAAACamZmZmZm5P3sUrkfhetQ/AAAAAAAAAACamZmZmZm5P9ejcD0K1+s/AAAAAAAAAACamZmZmZm5P9ejcD0K19M/AAAAAAAAAACamZmZmZm5P4XrUbgehes/AAAAAAAAAACamZmZmZm5PzMzMzMzM9M/AAAAAAAAAACamZmZmZm5PzMzMzMzM+s/AAAAAAAAAACamZmZmZm5P4/C9Shcj9I/AAAAAAAAAACamZmZmZm5PwAAAAAAAPA/AAAAAAAAAACamZmZmZm5P+F6FK5H4eo/AAAAAAAAAACamZmZmZm5P+xRuB6F69E/AAAAAAAAAACamZmZmZm5P65H4XoUru8/AAAAAAAAAACamZmZmZm5P5DC9Shcj+o/AAAAAAAAAACamZmZmZm5P0jhehSuR9E/AAAAAAAAAACamZmZmZm5P1yPwvUoXO8/AAAAAAAAAACamZmZmZm5Pz4K16NwPeo/AAAAAAAAAACamZmZmZm5P6RwPQrXo9A/AAAAAAAAAACamZmZmZm5P+xRuB6F6+k/AAAAAAAAAACamZmZmZm5PwAAAAAAANA/AAAAAAAAAACamZmZmZm5P5qZmZmZmek/AAAAAAAAAACamZmZmZm5P7gehetRuM4/AAAAAAAAAACamZmZmZm5P0jhehSuR+k/AAAAAAAAAACamZmZmZm5P3E9CtejcM0/AAAAAAAAAACamZmZmZm5P/YoXI/C9eg/AAAAAAAAAACamZmZmZm5Pylcj8L1KMw/AAAAAAAAAACamZmZmZm5P6RwPQrXo+g/AAAAAAAAAACamZmZmZm5P+F6FK5H4co/AAAAAAAAAACamZmZmZm5PwAAAAAAAOg/AAAAAAAAAACamZmZmZm5P1K4HoXrUcg/AAAAAAAAAACamZmZmZm5P65H4XoUruc/AAAAAAAAAACamZmZmZm5PwrXo3A9Csc/AAAAAAAAAACamZmZmZm5P1yPwvUoXOc/AAAAAAAAAACamZmZmZm5P8P1KFyPwsU/AAAAAAAAAACamZmZmZm5PwrXo3A9Cuc/AAAAAAAAAACamZmZmZm5P3sUrkfhesQ/AAAAAAAAAACamZmZmZm5P7gehetRuOY/AAAAAAAAAACamZmZmZm5PzMzMzMzM8M/AAAAAAAAAACamZmZmZm5P2dmZmZmZuY/AAAAAAAAAACamZmZmZm5P+xRuB6F68E/AAAAAAAAAACamZmZmZm5PxWuR+F6FOY/AAAAAAAAAACamZmZmZm5P6RwPQrXo8A/AAAAAAAAAACamZmZmZm5P8P1KFyPwuU/AAAAAAAAAACamZmZmZm5P7gehetRuL4/AAAAAAAAAACamZmZmZm5P3E9CtejcOU/AAAAAAAAAACamZmZmZm5Pylcj8L1KLw/AAAAAAAAAACamZmZmZm5Px+F61G4HuU/AAAAAAAAAACamZmZmZm5P5qZmZmZmbk/AAAAAAAAAACamZmZmZm5P83MzMzMzOQ/AAAAAAAAAACamZmZmZm5PwrXo3A9Crc/AAAAAAAAAACamZmZmZm5P3sUrkfheuQ/AAAAAAAAAACamZmZmZm5P3sUrkfherQ/AAAAAAAAAACamZmZmZm5Pylcj8L1KOQ/AAAAAAAAAACamZmZmZm5P+xRuB6F67E/AAAAAAAAAACamZmZmZm5P7gehetRuK4/AAAAAAAAAACamZmZmZm5P5qZmZmZmak/AAAAAAAAAACamZmZmZm5P3sUrkfheqQ/AAAAAAAAAACamZmZmZm5P7gehetRuJ4/AAAAAAAAAACamZmZmZm5P3sUrkfhepQ/AAAAAAAAAACamZmZmZm5P+F6FK5H4do/AAAAAAAAAACamZmZmZm5P3sUrkfheoQ/AAAAAAAAAACamZmZmZm5P1K4HoXrUeg/AAAAAAAAAACamZmZmZm5P5qZmZmZmck/AAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAAA=QAYAAAAAAABkAAAAAAAAAGEAAAAAAAAAYQAAAAAAAABfAAAAAAAAAF8AAAAAAAAAXgAAAAAAAABeAAAAAAAAAF0AAAAAAAAAXQAAAAAAAABcAAAAAAAAAFwAAAAAAAAAWwAAAAAAAABbAAAAAAAAAFoAAAAAAAAAWgAAAAAAAABYAAAAAAAAAFgAAAAAAAAAVgAAAAAAAABWAAAAAAAAAFQAAAAAAAAAVAAAAAAAAABSAAAAAAAAAFIAAAAAAAAAUAAAAAAAAABQAAAAAAAAAE4AAAAAAAAATgAAAAAAAABMAAAAAAAAAEwAAAAAAAAASgAAAAAAAABKAAAAAAAAAEgAAAAAAAAASAAAAAAAAABGAAAAAAAAAEYAAAAAAAAARAAAAAAAAABEAAAAAAAAAEIAAAAAAAAAQgAAAAAAAABjAAAAAAAAAGMAAAAAAAAAQAAAAAAAAABAAAAAAAAAAD4AAAAAAAAAPgAAAAAAAAA8AAAAAAAAADwAAAAAAAAAOgAAAAAAAAA6AAAAAAAAADgAAAAAAAAAOAAAAAAAAAA2AAAAAAAAADYAAAAAAAAAMwAAAAAAAAAzAAAAAAAAADAAAAAAAAAAMAAAAAAAAAAtAAAAAAAAAC0AAAAAAAAAKwAAAAAAAAArAAAAAAAAACkAAAAAAAAAKQAAAAAAAAAnAAAAAAAAACcAAAAAAAAAJQAAAAAAAAAlAAAAAAAAACMAAAAAAAAAIwAAAAAAAAAhAAAAAAAAACEAAAAAAAAAHwAAAAAAAAAfAAAAAAAAAB0AAAAAAAAAHQAAAAAAAAAbAAAAAAAAABsAAAAAAAAAGQAAAAAAAAAZAAAAAAAAABcAAAAAAAAAFwAAAAAAAAAVAAAAAAAAABUAAAAAAAAAYAAAAAAAAABgAAAAAAAAABMAAAAAAAAAEwAAAAAAAAASAAAAAAAAABIAAAAAAAAAEQAAAAAAAAARAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA8AAAAAAAAADgAAAAAAAAAOAAAAAAAAAA0AAAAAAAAADQAAAAAAAAAMAAAAAAAAAAwAAAAAAAAACwAAAAAAAAALAAAAAAAAAAoAAAAAAAAACgAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAcAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABQAAAAAAAAAFAAAAAAAAAAQAAAAAAAAABAAAAAAAAAADAAAAAAAAAAMAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAEAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAWQAAAAAAAABZAAAAAAAAAFcAAAAAAAAAVwAAAAAAAABVAAAAAAAAAFUAAAAAAAAAUwAAAAAAAABTAAAAAAAAAFEAAAAAAAAAUQAAAAAAAABPAAAAAAAAAE8AAAAAAAAATQAAAAAAAABNAAAAAAAAAEsAAAAAAAAASwAAAAAAAABJAAAAAAAAAEkAAAAAAAAARwAAAAAAAABHAAAAAAAAAEUAAAAAAAAARQAAAAAAAABDAAAAAAAAAEMAAAAAAAAAQQAAAAAAAABBAAAAAAAAAGIAAAAAAAAAYgAAAAAAAAA/AAAAAAAAAD8AAAAAAAAAPQAAAAAAAAA9AAAAAAAAADsAAAAAAAAAOwAAAAAAAAA5AAAAAAAAADkAAAAAAAAANwAAAAAAAAA3AAAAAAAAADUAAAAAAAAANQAAAAAAAAAyAAAAAAAAADIAAAAAAAAALwAAAAAAAAAvAAAAAAAAACwAAAAAAAAALAAAAAAAAAAqAAAAAAAAACoAAAAAAAAAKAAAAAAAAAAoAAAAAAAAACYAAAAAAAAAJgAAAAAAAAAkAAAAAAAAACQAAAAAAAAAIgAAAAAAAAAiAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAeAAAAAAAAAB4AAAAAAAAAHAAAAAAAAAAcAAAAAAAAABoAAAAAAAAAGgAAAAAAAAAYAAAAAAAAABgAAAAAAAAAFgAAAAAAAAAWAAAAAAAAABQAAAAAAAAAFAAAAAAAAAA0AAAAAAAAADQAAAAAAAAAMQAAAAAAAAAxAAAAAAAAAC4AAAAAAAAAIAMAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAABYAAAAAAAAAGAAAAAAAAAAaAAAAAAAAABwAAAAAAAAAHgAAAAAAAAAgAAAAAAAAACIAAAAAAAAAJAAAAAAAAAAmAAAAAAAAACgAAAAAAAAAKgAAAAAAAAAsAAAAAAAAAC4AAAAAAAAAMAAAAAAAAAAyAAAAAAAAADQAAAAAAAAANgAAAAAAAAA4AAAAAAAAADoAAAAAAAAAPAAAAAAAAAA+AAAAAAAAAEAAAAAAAAAAQgAAAAAAAABEAAAAAAAAAEYAAAAAAAAASAAAAAAAAABKAAAAAAAAAEwAAAAAAAAATgAAAAAAAABQAAAAAAAAAFIAAAAAAAAAVAAAAAAAAABWAAAAAAAAAFgAAAAAAAAAWgAAAAAAAABcAAAAAAAAAF4AAAAAAAAAYAAAAAAAAABiAAAAAAAAAGQAAAAAAAAAZgAAAAAAAABoAAAAAAAAAGoAAAAAAAAAbAAAAAAAAABuAAAAAAAAAHAAAAAAAAAAcgAAAAAAAAB0AAAAAAAAAHYAAAAAAAAAeAAAAAAAAAB6AAAAAAAAAHwAAAAAAAAAfgAAAAAAAACAAAAAAAAAAIIAAAAAAAAAhAAAAAAAAACGAAAAAAAAAIgAAAAAAAAAigAAAAAAAACMAAAAAAAAAI4AAAAAAAAAkAAAAAAAAACSAAAAAAAAAJQAAAAAAAAAlgAAAAAAAACYAAAAAAAAAJoAAAAAAAAAnAAAAAAAAACeAAAAAAAAAKAAAAAAAAAAogAAAAAAAACkAAAAAAAAAKYAAAAAAAAAqAAAAAAAAACqAAAAAAAAAKwAAAAAAAAArgAAAAAAAACwAAAAAAAAALIAAAAAAAAAtAAAAAAAAAC2AAAAAAAAALgAAAAAAAAAugAAAAAAAAC8AAAAAAAAAL4AAAAAAAAAwAAAAAAAAADCAAAAAAAAAMQAAAAAAAAAxgAAAAAAAADIAAAAAAAAAA==ZAAAAAAAAAADAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/boundary_top.vtu b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_top.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..13e747d6f55b590cb14fd287f0ceb234b8bfd519
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/boundary_top.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="2" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="200" RangeMax="201" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="99" RangeMax="99" offset="32" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1.0049875621" offset="56" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="132" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="164" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="188" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _EAAAAAAAAADIAAAAAAAAAMkAAAAAAAAACAAAAAAAAABjAAAAAAAAAA==MAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAACamZmZmZm5PwAAAAAAAPA/AAAAAAAAAAA=EAAAAAAAAAABAAAAAAAAAAAAAAAAAAAACAAAAAAAAAACAAAAAAAAAA==AQAAAAAAAAAD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/domain.vtu b/Tests/Data/TH2M/H2M/Liakopoulos/domain.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..68770c8bcb178632f77eb91f404f424ae8f3f05f
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/domain.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="202" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="201" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="2168" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.0049875621" offset="3248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="9724" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="14004" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="15084" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _UAYAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAGQAAAAAAAAAZQAAAAAAAABmAAAAAAAAAGcAAAAAAAAAaAAAAAAAAABpAAAAAAAAAGoAAAAAAAAAawAAAAAAAABsAAAAAAAAAG0AAAAAAAAAbgAAAAAAAABvAAAAAAAAAHAAAAAAAAAAcQAAAAAAAAByAAAAAAAAAHMAAAAAAAAAdAAAAAAAAAB1AAAAAAAAAHYAAAAAAAAAdwAAAAAAAAB4AAAAAAAAAHkAAAAAAAAAegAAAAAAAAB7AAAAAAAAAHwAAAAAAAAAfQAAAAAAAAB+AAAAAAAAAH8AAAAAAAAAgAAAAAAAAACBAAAAAAAAAIIAAAAAAAAAgwAAAAAAAACEAAAAAAAAAIUAAAAAAAAAhgAAAAAAAACHAAAAAAAAAIgAAAAAAAAAiQAAAAAAAACKAAAAAAAAAIsAAAAAAAAAjAAAAAAAAACNAAAAAAAAAI4AAAAAAAAAjwAAAAAAAACQAAAAAAAAAJEAAAAAAAAAkgAAAAAAAACTAAAAAAAAAJQAAAAAAAAAlQAAAAAAAACWAAAAAAAAAJcAAAAAAAAAmAAAAAAAAACZAAAAAAAAAJoAAAAAAAAAmwAAAAAAAACcAAAAAAAAAJ0AAAAAAAAAngAAAAAAAACfAAAAAAAAAKAAAAAAAAAAoQAAAAAAAACiAAAAAAAAAKMAAAAAAAAApAAAAAAAAAClAAAAAAAAAKYAAAAAAAAApwAAAAAAAACoAAAAAAAAAKkAAAAAAAAAqgAAAAAAAACrAAAAAAAAAKwAAAAAAAAArQAAAAAAAACuAAAAAAAAAK8AAAAAAAAAsAAAAAAAAACxAAAAAAAAALIAAAAAAAAAswAAAAAAAAC0AAAAAAAAALUAAAAAAAAAtgAAAAAAAAC3AAAAAAAAALgAAAAAAAAAuQAAAAAAAAC6AAAAAAAAALsAAAAAAAAAvAAAAAAAAAC9AAAAAAAAAL4AAAAAAAAAvwAAAAAAAADAAAAAAAAAAMEAAAAAAAAAwgAAAAAAAADDAAAAAAAAAMQAAAAAAAAAxQAAAAAAAADGAAAAAAAAAMcAAAAAAAAAyAAAAAAAAADJAAAAAAAAAA==IAMAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAA==8BIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACamZmZmZm5PwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAHsUrkfheoQ/AAAAAAAAAACamZmZmZm5P3sUrkfheoQ/AAAAAAAAAAAAAAAAAAAAAHsUrkfhepQ/AAAAAAAAAACamZmZmZm5P3sUrkfhepQ/AAAAAAAAAAAAAAAAAAAAALgehetRuJ4/AAAAAAAAAACamZmZmZm5P7gehetRuJ4/AAAAAAAAAAAAAAAAAAAAAHsUrkfheqQ/AAAAAAAAAACamZmZmZm5P3sUrkfheqQ/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmak/AAAAAAAAAACamZmZmZm5P5qZmZmZmak/AAAAAAAAAAAAAAAAAAAAALgehetRuK4/AAAAAAAAAACamZmZmZm5P7gehetRuK4/AAAAAAAAAAAAAAAAAAAAAOxRuB6F67E/AAAAAAAAAACamZmZmZm5P+xRuB6F67E/AAAAAAAAAAAAAAAAAAAAAHsUrkfherQ/AAAAAAAAAACamZmZmZm5P3sUrkfherQ/AAAAAAAAAAAAAAAAAAAAAArXo3A9Crc/AAAAAAAAAACamZmZmZm5PwrXo3A9Crc/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmbk/AAAAAAAAAACamZmZmZm5P5qZmZmZmbk/AAAAAAAAAAAAAAAAAAAAAClcj8L1KLw/AAAAAAAAAACamZmZmZm5Pylcj8L1KLw/AAAAAAAAAAAAAAAAAAAAALgehetRuL4/AAAAAAAAAACamZmZmZm5P7gehetRuL4/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo8A/AAAAAAAAAACamZmZmZm5P6RwPQrXo8A/AAAAAAAAAAAAAAAAAAAAAOxRuB6F68E/AAAAAAAAAACamZmZmZm5P+xRuB6F68E/AAAAAAAAAAAAAAAAAAAAADMzMzMzM8M/AAAAAAAAAACamZmZmZm5PzMzMzMzM8M/AAAAAAAAAAAAAAAAAAAAAHsUrkfhesQ/AAAAAAAAAACamZmZmZm5P3sUrkfhesQ/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwsU/AAAAAAAAAACamZmZmZm5P8P1KFyPwsU/AAAAAAAAAAAAAAAAAAAAAArXo3A9Csc/AAAAAAAAAACamZmZmZm5PwrXo3A9Csc/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUcg/AAAAAAAAAACamZmZmZm5P1K4HoXrUcg/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmck/AAAAAAAAAACamZmZmZm5P5qZmZmZmck/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4co/AAAAAAAAAACamZmZmZm5P+F6FK5H4co/AAAAAAAAAAAAAAAAAAAAAClcj8L1KMw/AAAAAAAAAACamZmZmZm5Pylcj8L1KMw/AAAAAAAAAAAAAAAAAAAAAHE9CtejcM0/AAAAAAAAAACamZmZmZm5P3E9CtejcM0/AAAAAAAAAAAAAAAAAAAAALgehetRuM4/AAAAAAAAAACamZmZmZm5P7gehetRuM4/AAAAAAAAAAAAAAAAAAAAAAAAAAAAANA/AAAAAAAAAACamZmZmZm5PwAAAAAAANA/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo9A/AAAAAAAAAACamZmZmZm5P6RwPQrXo9A/AAAAAAAAAAAAAAAAAAAAAEjhehSuR9E/AAAAAAAAAACamZmZmZm5P0jhehSuR9E/AAAAAAAAAAAAAAAAAAAAAOxRuB6F69E/AAAAAAAAAACamZmZmZm5P+xRuB6F69E/AAAAAAAAAAAAAAAAAAAAAI/C9Shcj9I/AAAAAAAAAACamZmZmZm5P4/C9Shcj9I/AAAAAAAAAAAAAAAAAAAAADMzMzMzM9M/AAAAAAAAAACamZmZmZm5PzMzMzMzM9M/AAAAAAAAAAAAAAAAAAAAANejcD0K19M/AAAAAAAAAACamZmZmZm5P9ejcD0K19M/AAAAAAAAAAAAAAAAAAAAAHsUrkfhetQ/AAAAAAAAAACamZmZmZm5P3sUrkfhetQ/AAAAAAAAAAAAAAAAAAAAAB+F61G4HtU/AAAAAAAAAACamZmZmZm5Px+F61G4HtU/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwtU/AAAAAAAAAACamZmZmZm5P8P1KFyPwtU/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZtY/AAAAAAAAAACamZmZmZm5P2dmZmZmZtY/AAAAAAAAAAAAAAAAAAAAAArXo3A9Ctc/AAAAAAAAAACamZmZmZm5PwrXo3A9Ctc/AAAAAAAAAAAAAAAAAAAAAK5H4XoUrtc/AAAAAAAAAACamZmZmZm5P65H4XoUrtc/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUdg/AAAAAAAAAACamZmZmZm5P1K4HoXrUdg/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9dg/AAAAAAAAAACamZmZmZm5P/YoXI/C9dg/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmdk/AAAAAAAAAACamZmZmZm5P5qZmZmZmdk/AAAAAAAAAAAAAAAAAAAAAD4K16NwPdo/AAAAAAAAAACamZmZmZm5Pz4K16NwPdo/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4do/AAAAAAAAAACamZmZmZm5P+F6FK5H4do/AAAAAAAAAAAAAAAAAAAAAIXrUbgehds/AAAAAAAAAACamZmZmZm5P4XrUbgehds/AAAAAAAAAAAAAAAAAAAAAClcj8L1KNw/AAAAAAAAAACamZmZmZm5Pylcj8L1KNw/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzNw/AAAAAAAAAACamZmZmZm5P83MzMzMzNw/AAAAAAAAAAAAAAAAAAAAAHE9CtejcN0/AAAAAAAAAACamZmZmZm5P3E9CtejcN0/AAAAAAAAAAAAAAAAAAAAABWuR+F6FN4/AAAAAAAAAACamZmZmZm5PxWuR+F6FN4/AAAAAAAAAAAAAAAAAAAAALgehetRuN4/AAAAAAAAAACamZmZmZm5P7gehetRuN4/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXN8/AAAAAAAAAACamZmZmZm5P1yPwvUoXN8/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAACamZmZmZm5PwAAAAAAAOA/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUeA/AAAAAAAAAACamZmZmZm5P1K4HoXrUeA/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo+A/AAAAAAAAAACamZmZmZm5P6RwPQrXo+A/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9eA/AAAAAAAAAACamZmZmZm5P/YoXI/C9eA/AAAAAAAAAAAAAAAAAAAAAEjhehSuR+E/AAAAAAAAAACamZmZmZm5P0jhehSuR+E/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmeE/AAAAAAAAAACamZmZmZm5P5qZmZmZmeE/AAAAAAAAAAAAAAAAAAAAAOxRuB6F6+E/AAAAAAAAAACamZmZmZm5P+xRuB6F6+E/AAAAAAAAAAAAAAAAAAAAAD4K16NwPeI/AAAAAAAAAACamZmZmZm5Pz4K16NwPeI/AAAAAAAAAAAAAAAAAAAAAI/C9Shcj+I/AAAAAAAAAACamZmZmZm5P4/C9Shcj+I/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4eI/AAAAAAAAAACamZmZmZm5P+F6FK5H4eI/AAAAAAAAAAAAAAAAAAAAADMzMzMzM+M/AAAAAAAAAACamZmZmZm5PzMzMzMzM+M/AAAAAAAAAAAAAAAAAAAAAIXrUbgeheM/AAAAAAAAAACamZmZmZm5P4XrUbgeheM/AAAAAAAAAAAAAAAAAAAAANejcD0K1+M/AAAAAAAAAACamZmZmZm5P9ejcD0K1+M/AAAAAAAAAAAAAAAAAAAAAClcj8L1KOQ/AAAAAAAAAACamZmZmZm5Pylcj8L1KOQ/AAAAAAAAAAAAAAAAAAAAAHsUrkfheuQ/AAAAAAAAAACamZmZmZm5P3sUrkfheuQ/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOQ/AAAAAAAAAACamZmZmZm5P83MzMzMzOQ/AAAAAAAAAAAAAAAAAAAAAB+F61G4HuU/AAAAAAAAAACamZmZmZm5Px+F61G4HuU/AAAAAAAAAAAAAAAAAAAAAHE9CtejcOU/AAAAAAAAAACamZmZmZm5P3E9CtejcOU/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwuU/AAAAAAAAAACamZmZmZm5P8P1KFyPwuU/AAAAAAAAAAAAAAAAAAAAABWuR+F6FOY/AAAAAAAAAACamZmZmZm5PxWuR+F6FOY/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZuY/AAAAAAAAAACamZmZmZm5P2dmZmZmZuY/AAAAAAAAAAAAAAAAAAAAALgehetRuOY/AAAAAAAAAACamZmZmZm5P7gehetRuOY/AAAAAAAAAAAAAAAAAAAAAArXo3A9Cuc/AAAAAAAAAACamZmZmZm5PwrXo3A9Cuc/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXOc/AAAAAAAAAACamZmZmZm5P1yPwvUoXOc/AAAAAAAAAAAAAAAAAAAAAK5H4XoUruc/AAAAAAAAAACamZmZmZm5P65H4XoUruc/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOg/AAAAAAAAAACamZmZmZm5PwAAAAAAAOg/AAAAAAAAAAAAAAAAAAAAAFK4HoXrUeg/AAAAAAAAAACamZmZmZm5P1K4HoXrUeg/AAAAAAAAAAAAAAAAAAAAAKRwPQrXo+g/AAAAAAAAAACamZmZmZm5P6RwPQrXo+g/AAAAAAAAAAAAAAAAAAAAAPYoXI/C9eg/AAAAAAAAAACamZmZmZm5P/YoXI/C9eg/AAAAAAAAAAAAAAAAAAAAAEjhehSuR+k/AAAAAAAAAACamZmZmZm5P0jhehSuR+k/AAAAAAAAAAAAAAAAAAAAAJqZmZmZmek/AAAAAAAAAACamZmZmZm5P5qZmZmZmek/AAAAAAAAAAAAAAAAAAAAAOxRuB6F6+k/AAAAAAAAAACamZmZmZm5P+xRuB6F6+k/AAAAAAAAAAAAAAAAAAAAAD4K16NwPeo/AAAAAAAAAACamZmZmZm5Pz4K16NwPeo/AAAAAAAAAAAAAAAAAAAAAJDC9Shcj+o/AAAAAAAAAACamZmZmZm5P5DC9Shcj+o/AAAAAAAAAAAAAAAAAAAAAOF6FK5H4eo/AAAAAAAAAACamZmZmZm5P+F6FK5H4eo/AAAAAAAAAAAAAAAAAAAAADMzMzMzM+s/AAAAAAAAAACamZmZmZm5PzMzMzMzM+s/AAAAAAAAAAAAAAAAAAAAAIXrUbgehes/AAAAAAAAAACamZmZmZm5P4XrUbgehes/AAAAAAAAAAAAAAAAAAAAANejcD0K1+s/AAAAAAAAAACamZmZmZm5P9ejcD0K1+s/AAAAAAAAAAAAAAAAAAAAAClcj8L1KOw/AAAAAAAAAACamZmZmZm5Pylcj8L1KOw/AAAAAAAAAAAAAAAAAAAAAHsUrkfheuw/AAAAAAAAAACamZmZmZm5P3sUrkfheuw/AAAAAAAAAAAAAAAAAAAAAM3MzMzMzOw/AAAAAAAAAACamZmZmZm5P83MzMzMzOw/AAAAAAAAAAAAAAAAAAAAAB+F61G4Hu0/AAAAAAAAAACamZmZmZm5Px+F61G4Hu0/AAAAAAAAAAAAAAAAAAAAAHE9CtejcO0/AAAAAAAAAACamZmZmZm5P3E9CtejcO0/AAAAAAAAAAAAAAAAAAAAAMP1KFyPwu0/AAAAAAAAAACamZmZmZm5P8P1KFyPwu0/AAAAAAAAAAAAAAAAAAAAABWuR+F6FO4/AAAAAAAAAACamZmZmZm5PxWuR+F6FO4/AAAAAAAAAAAAAAAAAAAAAGdmZmZmZu4/AAAAAAAAAACamZmZmZm5P2dmZmZmZu4/AAAAAAAAAAAAAAAAAAAAALgehetRuO4/AAAAAAAAAACamZmZmZm5P7gehetRuO4/AAAAAAAAAAAAAAAAAAAAAArXo3A9Cu8/AAAAAAAAAACamZmZmZm5PwrXo3A9Cu8/AAAAAAAAAAAAAAAAAAAAAFyPwvUoXO8/AAAAAAAAAACamZmZmZm5P1yPwvUoXO8/AAAAAAAAAAAAAAAAAAAAAK5H4XoUru8/AAAAAAAAAACamZmZmZm5P65H4XoUru8/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAACamZmZmZm5PwAAAAAAAPA/AAAAAAAAAAA=gAwAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAwAAAAAAAAAFAAAAAAAAAAQAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAcAAAAAAAAABgAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACQAAAAAAAAAIAAAAAAAAAAgAAAAAAAAACQAAAAAAAAALAAAAAAAAAAoAAAAAAAAACgAAAAAAAAALAAAAAAAAAA0AAAAAAAAADAAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADwAAAAAAAAAOAAAAAAAAAA4AAAAAAAAADwAAAAAAAAARAAAAAAAAABAAAAAAAAAAEAAAAAAAAAARAAAAAAAAABMAAAAAAAAAEgAAAAAAAAASAAAAAAAAABMAAAAAAAAAFQAAAAAAAAAUAAAAAAAAABQAAAAAAAAAFQAAAAAAAAAXAAAAAAAAABYAAAAAAAAAFgAAAAAAAAAXAAAAAAAAABkAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGwAAAAAAAAAaAAAAAAAAABoAAAAAAAAAGwAAAAAAAAAdAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAB8AAAAAAAAAHgAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIQAAAAAAAAAgAAAAAAAAACAAAAAAAAAAIQAAAAAAAAAjAAAAAAAAACIAAAAAAAAAIgAAAAAAAAAjAAAAAAAAACUAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJwAAAAAAAAAmAAAAAAAAACYAAAAAAAAAJwAAAAAAAAApAAAAAAAAACgAAAAAAAAAKAAAAAAAAAApAAAAAAAAACsAAAAAAAAAKgAAAAAAAAAqAAAAAAAAACsAAAAAAAAALQAAAAAAAAAsAAAAAAAAACwAAAAAAAAALQAAAAAAAAAvAAAAAAAAAC4AAAAAAAAALgAAAAAAAAAvAAAAAAAAADEAAAAAAAAAMAAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMwAAAAAAAAAyAAAAAAAAADIAAAAAAAAAMwAAAAAAAAA1AAAAAAAAADQAAAAAAAAANAAAAAAAAAA1AAAAAAAAADcAAAAAAAAANgAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOQAAAAAAAAA4AAAAAAAAADgAAAAAAAAAOQAAAAAAAAA7AAAAAAAAADoAAAAAAAAAOgAAAAAAAAA7AAAAAAAAAD0AAAAAAAAAPAAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPwAAAAAAAAA+AAAAAAAAAD4AAAAAAAAAPwAAAAAAAABBAAAAAAAAAEAAAAAAAAAAQAAAAAAAAABBAAAAAAAAAEMAAAAAAAAAQgAAAAAAAABCAAAAAAAAAEMAAAAAAAAARQAAAAAAAABEAAAAAAAAAEQAAAAAAAAARQAAAAAAAABHAAAAAAAAAEYAAAAAAAAARgAAAAAAAABHAAAAAAAAAEkAAAAAAAAASAAAAAAAAABIAAAAAAAAAEkAAAAAAAAASwAAAAAAAABKAAAAAAAAAEoAAAAAAAAASwAAAAAAAABNAAAAAAAAAEwAAAAAAAAATAAAAAAAAABNAAAAAAAAAE8AAAAAAAAATgAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUQAAAAAAAABQAAAAAAAAAFAAAAAAAAAAUQAAAAAAAABTAAAAAAAAAFIAAAAAAAAAUgAAAAAAAABTAAAAAAAAAFUAAAAAAAAAVAAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVwAAAAAAAABWAAAAAAAAAFYAAAAAAAAAVwAAAAAAAABZAAAAAAAAAFgAAAAAAAAAWAAAAAAAAABZAAAAAAAAAFsAAAAAAAAAWgAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXQAAAAAAAABcAAAAAAAAAFwAAAAAAAAAXQAAAAAAAABfAAAAAAAAAF4AAAAAAAAAXgAAAAAAAABfAAAAAAAAAGEAAAAAAAAAYAAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYwAAAAAAAABiAAAAAAAAAGIAAAAAAAAAYwAAAAAAAABlAAAAAAAAAGQAAAAAAAAAZAAAAAAAAABlAAAAAAAAAGcAAAAAAAAAZgAAAAAAAABmAAAAAAAAAGcAAAAAAAAAaQAAAAAAAABoAAAAAAAAAGgAAAAAAAAAaQAAAAAAAABrAAAAAAAAAGoAAAAAAAAAagAAAAAAAABrAAAAAAAAAG0AAAAAAAAAbAAAAAAAAABsAAAAAAAAAG0AAAAAAAAAbwAAAAAAAABuAAAAAAAAAG4AAAAAAAAAbwAAAAAAAABxAAAAAAAAAHAAAAAAAAAAcAAAAAAAAABxAAAAAAAAAHMAAAAAAAAAcgAAAAAAAAByAAAAAAAAAHMAAAAAAAAAdQAAAAAAAAB0AAAAAAAAAHQAAAAAAAAAdQAAAAAAAAB3AAAAAAAAAHYAAAAAAAAAdgAAAAAAAAB3AAAAAAAAAHkAAAAAAAAAeAAAAAAAAAB4AAAAAAAAAHkAAAAAAAAAewAAAAAAAAB6AAAAAAAAAHoAAAAAAAAAewAAAAAAAAB9AAAAAAAAAHwAAAAAAAAAfAAAAAAAAAB9AAAAAAAAAH8AAAAAAAAAfgAAAAAAAAB+AAAAAAAAAH8AAAAAAAAAgQAAAAAAAACAAAAAAAAAAIAAAAAAAAAAgQAAAAAAAACDAAAAAAAAAIIAAAAAAAAAggAAAAAAAACDAAAAAAAAAIUAAAAAAAAAhAAAAAAAAACEAAAAAAAAAIUAAAAAAAAAhwAAAAAAAACGAAAAAAAAAIYAAAAAAAAAhwAAAAAAAACJAAAAAAAAAIgAAAAAAAAAiAAAAAAAAACJAAAAAAAAAIsAAAAAAAAAigAAAAAAAACKAAAAAAAAAIsAAAAAAAAAjQAAAAAAAACMAAAAAAAAAIwAAAAAAAAAjQAAAAAAAACPAAAAAAAAAI4AAAAAAAAAjgAAAAAAAACPAAAAAAAAAJEAAAAAAAAAkAAAAAAAAACQAAAAAAAAAJEAAAAAAAAAkwAAAAAAAACSAAAAAAAAAJIAAAAAAAAAkwAAAAAAAACVAAAAAAAAAJQAAAAAAAAAlAAAAAAAAACVAAAAAAAAAJcAAAAAAAAAlgAAAAAAAACWAAAAAAAAAJcAAAAAAAAAmQAAAAAAAACYAAAAAAAAAJgAAAAAAAAAmQAAAAAAAACbAAAAAAAAAJoAAAAAAAAAmgAAAAAAAACbAAAAAAAAAJ0AAAAAAAAAnAAAAAAAAACcAAAAAAAAAJ0AAAAAAAAAnwAAAAAAAACeAAAAAAAAAJ4AAAAAAAAAnwAAAAAAAAChAAAAAAAAAKAAAAAAAAAAoAAAAAAAAAChAAAAAAAAAKMAAAAAAAAAogAAAAAAAACiAAAAAAAAAKMAAAAAAAAApQAAAAAAAACkAAAAAAAAAKQAAAAAAAAApQAAAAAAAACnAAAAAAAAAKYAAAAAAAAApgAAAAAAAACnAAAAAAAAAKkAAAAAAAAAqAAAAAAAAACoAAAAAAAAAKkAAAAAAAAAqwAAAAAAAACqAAAAAAAAAKoAAAAAAAAAqwAAAAAAAACtAAAAAAAAAKwAAAAAAAAArAAAAAAAAACtAAAAAAAAAK8AAAAAAAAArgAAAAAAAACuAAAAAAAAAK8AAAAAAAAAsQAAAAAAAACwAAAAAAAAALAAAAAAAAAAsQAAAAAAAACzAAAAAAAAALIAAAAAAAAAsgAAAAAAAACzAAAAAAAAALUAAAAAAAAAtAAAAAAAAAC0AAAAAAAAALUAAAAAAAAAtwAAAAAAAAC2AAAAAAAAALYAAAAAAAAAtwAAAAAAAAC5AAAAAAAAALgAAAAAAAAAuAAAAAAAAAC5AAAAAAAAALsAAAAAAAAAugAAAAAAAAC6AAAAAAAAALsAAAAAAAAAvQAAAAAAAAC8AAAAAAAAALwAAAAAAAAAvQAAAAAAAAC/AAAAAAAAAL4AAAAAAAAAvgAAAAAAAAC/AAAAAAAAAMEAAAAAAAAAwAAAAAAAAADAAAAAAAAAAMEAAAAAAAAAwwAAAAAAAADCAAAAAAAAAMIAAAAAAAAAwwAAAAAAAADFAAAAAAAAAMQAAAAAAAAAxAAAAAAAAADFAAAAAAAAAMcAAAAAAAAAxgAAAAAAAADGAAAAAAAAAMcAAAAAAAAAyQAAAAAAAADIAAAAAAAAAA==IAMAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAAAAQAAAAAAAAQBAAAAAAAACAEAAAAAAAAMAQAAAAAAABABAAAAAAAAFAEAAAAAAAAYAQAAAAAAABwBAAAAAAAAIAEAAAAAAAAkAQAAAAAAACgBAAAAAAAALAEAAAAAAAAwAQAAAAAAADQBAAAAAAAAOAEAAAAAAAA8AQAAAAAAAEABAAAAAAAARAEAAAAAAABIAQAAAAAAAEwBAAAAAAAAUAEAAAAAAABUAQAAAAAAAFgBAAAAAAAAXAEAAAAAAABgAQAAAAAAAGQBAAAAAAAAaAEAAAAAAABsAQAAAAAAAHABAAAAAAAAdAEAAAAAAAB4AQAAAAAAAHwBAAAAAAAAgAEAAAAAAACEAQAAAAAAAIgBAAAAAAAAjAEAAAAAAACQAQAAAAAAAA==ZAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/Liakopoulos/liakopoulos_TH2M.prj b/Tests/Data/TH2M/H2M/Liakopoulos/liakopoulos_TH2M.prj
new file mode 100644
index 0000000000000000000000000000000000000000..3def3a94de52eb770b819a77928d16e7e0d891f0
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/Liakopoulos/liakopoulos_TH2M.prj
@@ -0,0 +1,507 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh axially_symmetric="true">domain.vtu</mesh>
+ <mesh axially_symmetric="true">boundary_left.vtu</mesh>
+ <mesh axially_symmetric="true">boundary_right.vtu</mesh>
+ <mesh axially_symmetric="true">boundary_top.vtu</mesh>
+ <mesh axially_symmetric="true">boundary_bottom.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>TH2M</name>
+ <type>TH2M</type>
+ <integration_order>2</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>
+ 1.0 1.0 1.0 1.0
+ 1.0 1.0 1.0 1.0
+ 1.0 1.0 1.0 1.0
+ 1.0 1.0 1.0 1.0
+ 1.0 1.0 1.0 1.0
+ </component_magnitudes>
+ <relative_epsilons>
+ 1e-03 1e-03 1e-03 1e-03
+ 1e-03 1e-03 1e-03 1e-03
+ 1e-06 1e-06 1e-06 1e-06
+ 1e-08 1e-08 1e-08 1e-08
+ 1e-08 1e-08 1e-08 1e-08
+ </relative_epsilons>
+ </jacobian_assembler>
+ <reference_temperature>T_0</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="vapour_pressure" output_name="vapour_pressure"/>
+ <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_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_variable internal_name="mole_fraction_gas" output_name="xnCG"/>
+ <secondary_variable internal_name="mass_fraction_gas" output_name="xmCG"/>
+ <secondary_variable internal_name="mass_fraction_liquid" output_name="xmWL"/>
+ </secondary_variables>
+ <initial_stress>sigma0</initial_stress>
+ <specific_body_force>0 -9.81</specific_body_force>
+ <mass_lumping>false</mass_lumping>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1e3</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>4280.0</value>
+ </property>
+ <property>
+ <name>molar_mass</name>
+ <type>Constant</type>
+ <value>.018</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>0.6</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Gas</type>
+ <properties>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>800.0</value>
+ </property>
+ <property>
+ <name>molar_mass</name>
+ <type>Constant</type>
+ <value>.028949</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>0.6</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>IdealGasLaw</type>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1.8e-5</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2e3</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.838</value>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>917.654</value>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>4.5e-13</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>RelPermLiakopoulos</type>
+ </property>
+ <property>
+ <name>relative_permeability_nonwetting_phase</name>
+ <type>RelPermBrooksCoreyNonwettingPhase</type>
+ <residual_liquid_saturation>0.2</residual_liquid_saturation>
+ <residual_gas_saturation>0.0</residual_gas_saturation>
+ <min_relative_permeability>1.0e-4</min_relative_permeability>
+ <lambda>3.0</lambda>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.2975</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-15 1e-15 1e-05 1e-05 1.e-05</abstols>
+ <reltols>1e-12 1e-12 1e-15 1e-10 1.e-10</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0.0 </t_initial>
+ <t_end>7200</t_end>
+ <timesteps>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>5</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>50</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ <pair>
+ <repeat>5</repeat>
+ <delta_t>200</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>result_liakopoulos</prefix>
+ <suffix>_t_{:gtime}</suffix>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1000</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>gas_pressure</variable>
+ <variable>gas_pressure_interpolated</variable>
+ <variable>capillary_pressure</variable>
+ <variable>capillary_pressure_interpolated</variable>
+ <variable>liquid_pressure_interpolated</variable>
+ <variable>temperature</variable>
+ <variable>temperature_interpolated</variable>
+ <variable>displacement</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>vapour_pressure</variable>
+ <variable>velocity_gas</variable>
+ <variable>velocity_liquid</variable>
+ <variable>liquid_density</variable>
+ <variable>gas_density</variable>
+ <variable>porosity</variable>
+ <variable>saturation</variable>
+ <variable>xnCG</variable>
+ <variable>xmCG</variable>
+ <variable>xmWL</variable>
+ </variables>
+ <fixed_output_times>
+ 120.
+ 300.0
+ 4800.0
+ 7200.0
+ </fixed_output_times>
+ </output>
+ </time_loop>
+ <parameters>
+ <!-- Mechanics -->
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1.3e6</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.4</value>
+ </parameter>
+ <parameter>
+ <name>sigma0</name>
+ <type>Function</type>
+ <expression>101325</expression>
+ <expression>101325</expression>
+ <expression>0</expression>
+ <expression>0</expression>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>zero</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>temperature_ic</name>
+ <type>Constant</type>
+ <value>300</value>
+ </parameter>
+ <parameter>
+ <name>T_0</name>
+ <type>Constant</type>
+ <value>300</value>
+ </parameter>
+ <parameter>
+ <name>gas_pressure_ic</name>
+ <type>Constant</type>
+ <value>101325</value>
+ </parameter>
+ <parameter>
+ <name>capillary_pressure_ic</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ <parameter>
+ <name>capillary_pressure_bc</name>
+ <type>Constant</type>
+ <value>100</value>
+ </parameter>
+ </parameters>
+ <!-- ## Displacement equation #########################################-->
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>1</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>boundary_left</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>boundary_right</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>boundary_bottom</mesh>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>zero</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <!-- ## Gas pressure equation #############################################-->
+ <process_variable>
+ <name>gas_pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>gas_pressure_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>boundary_top</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>gas_pressure_ic</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>boundary_bottom</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>gas_pressure_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>capillary_pressure_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>boundary_bottom</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>capillary_pressure_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <!-- ## Temperature equation ##############################################-->
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>temperature_ic</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>domain</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>temperature_ic</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ <source_terms/>
+ </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>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </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>
+ <test_definition>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>gas_pressure_interpolated</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>capillary_pressure_interpolated</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>temperature_interpolated</field>
+ <absolute_tolerance>1e-6</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>displacement</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>sigma</field>
+ <absolute_tolerance>1e-07</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>velocity_gas</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>velocity_liquid</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>gas_density</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>liquid_density</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>porosity</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>result_liakopoulos_t_.*.vtu</regex>
+ <field>saturation</field>
+ <absolute_tolerance>1e-08</absolute_tolerance>
+ <relative_tolerance>0.</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>