From d17ad6fc7b638abe04be7e9638f78b12d94afe7e Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Wed, 6 Apr 2022 17:54:44 +0200
Subject: [PATCH] [T/TH2M] Swelling test from RM
The test's data is the same but few variable names were updated.
Temperature and gas pressure are held fast.
---
ProcessLib/TH2M/Tests.cmake | 1 +
.../TH2M/H2M/OrthotropicSwelling/square.prj | 506 ++++++++++++++++++
.../H2M/OrthotropicSwelling/square_1x1.gml | 31 ++
.../square_1x1_quad8_1e0.vtu | 25 +
.../square_ts_104_t_3.000000.vtu | 47 ++
.../square_ts_25_t_0.500000.vtu | 47 ++
.../square_ts_50_t_1.000000.vtu | 47 ++
.../square_ts_54_t_2.000000.vtu | 47 ++
.../square_ts_5_t_0.100000.vtu | 47 ++
.../square_ts_79_t_2.500000.vtu | 47 ++
.../square_ts_99_t_2.900000.vtu | 47 ++
11 files changed, 892 insertions(+)
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1.gml
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1_quad8_1e0.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_104_t_3.000000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_25_t_0.500000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_50_t_1.000000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_54_t_2.000000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_5_t_0.100000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_79_t_2.500000.vtu
create mode 100644 Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_99_t_2.900000.vtu
diff --git a/ProcessLib/TH2M/Tests.cmake b/ProcessLib/TH2M/Tests.cmake
index ecf9700f004..a7bce2f1964 100644
--- a/ProcessLib/TH2M/Tests.cmake
+++ b/ProcessLib/TH2M/Tests.cmake
@@ -10,6 +10,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE TH2M/TH/idealGasLaw/compression_gas.prj RUNTIME 1)
OgsTest(PROJECTFILE TH2M/H2M/Liakopoulos/liakopoulos_TH2M.prj RUNTIME 15)
OgsTest(PROJECTFILE TH2M/H2M/Liakopoulos/liakopoulos_newton.xml RUNTIME 5)
+ OgsTest(PROJECTFILE TH2M/H2M/OrthotropicSwelling/square.prj RUNTIME 1)
OgsTest(PROJECTFILE TH2M/H2/mcWorther/mcWorther_h2.prj RUNTIME 55)
OgsTest(PROJECTFILE TH2M/H2/mcWorther/mcWorther_h2_newton.xml RUNTIME 20)
OgsTest(PROJECTFILE TH2M/TH2/unicube/unicube.prj RUNTIME 25)
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
new file mode 100644
index 00000000000..513444bbb8e
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
@@ -0,0 +1,506 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+ <mesh>square_1x1_quad8_1e0.vtu</mesh>
+ <geometry>square_1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>process</name>
+ <type>TH2M</type>
+ <integration_order>3</integration_order>
+ <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
+ 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
+ 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08
+ </relative_epsilons>
+ </jacobian_assembler>
+ <!--
+ -->
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <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="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="gas_density" output_name="gas_density"/>
+ <secondary_variable internal_name="liquid_density" output_name="liquid_density"/>
+ <secondary_variable internal_name="mass_fraction_gas" output_name="xmCG"/>
+ <secondary_variable internal_name="mass_fraction_liquid" output_name="xmWL"/>
+ <secondary_variable internal_name="mole_fraction_gas" output_name="xnCG"/>
+ <secondary_variable internal_name="porosity" output_name="porosity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="swelling_stress" output_name="swelling_stress"/>
+ <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_variables>
+ <reference_temperature>T_0</reference_temperature>
+ <specific_body_force>0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Linear</type>
+ <reference_value>1</reference_value>
+ <independent_variable>
+ <variable_name>phase_pressure</variable_name>
+ <reference_condition>0</reference_condition>
+ <slope>0.1</slope>
+ </independent_variable>
+ </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>1</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>swelling_stress_rate</name>
+ <type>SaturationDependentSwelling</type>
+ <swelling_pressures>.1 .2 .3</swelling_pressures>
+ <exponents>1 1 1</exponents>
+ <lower_saturation_limit>0.3</lower_saturation_limit>
+ <upper_saturation_limit>0.8</upper_saturation_limit>
+ </property>
+ <property>
+ <name>thermal_expansivity</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>SaturationVanGenuchten</type>
+ <residual_liquid_saturation>0.1</residual_liquid_saturation>
+ <residual_gas_saturation>0.05</residual_gas_saturation>
+ <exponent>0.8</exponent>
+ <p_b>0.5</p_b>
+ </property>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>relative_permeability_nonwetting_phase</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>EffectiveThermalConductivityPorosityMixing</type>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsPowerLaw</type>
+ <exponent>1</exponent>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="process">
+ <nonlinear_solver>nonlinear_solver</nonlinear_solver>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-15 1e-15 1e-15 1e-15 1e-15</abstols>
+ <reltols>1e-15 1e-15 1e-15 1e-15 1e-15</reltols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>3</t_end>
+ <timesteps>
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+ <pair>
+ <repeat>4</repeat>
+ <delta_t>.25</delta_t>
+ </pair>
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>square</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>capillary_pressure</variable>
+ <variable>displacement</variable>
+ <variable>epsilon</variable>
+ <variable>gas_density</variable>
+ <variable>gas_pressure</variable>
+ <variable>liquid_density</variable>
+ <variable>porosity</variable>
+ <variable>saturation</variable>
+ <variable>sigma</variable>
+ <variable>swelling_stress</variable>
+ <variable>temperature</variable>
+ <variable>vapour_pressure</variable>
+ <variable>velocity_gas</variable>
+ <variable>velocity_liquid</variable>
+ <variable>xmCG</variable>
+ <variable>xmWL</variable>
+ <variable>xnCG</variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>1 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>0 1</values>
+ </parameter>
+ <!-- Mechanics -->
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>T_0</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0</values>
+ </parameter>
+ <parameter>
+ <name>gas_pressure</name>
+ <type>Constant</type>
+ <value>1e-16</value>
+ </parameter>
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>-0.1</value>
+ </parameter>
+ <parameter>
+ <name>pressure_bc</name>
+ <type>CurveScaled</type>
+ <curve>pressure_ramp</curve>
+ <parameter>p0</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>pressure_ramp</name>
+ <coords>0 1 2 3</coords>
+ <values>1 -10 -10 1</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>2</components>
+ <order>2</order>
+ <compensate_non_equilibrium_initial_residuum>
+ true</compensate_non_equilibrium_initial_residuum>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>capillary_pressure</name>
+ <components>1</components>
+ <compensate_non_equilibrium_initial_residuum>
+ true</compensate_non_equilibrium_initial_residuum>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>gas_pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>gas_pressure</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>gas_pressure</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>square_1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>gas_pressure</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T_0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>square_1x1_quad8_1e0</mesh>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>T_0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ <source_terms/>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>nonlinear_solver</name>
+ <type>Newton</type>
+ <max_iter>100</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>liquid_pressure_interpolated</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>saturation</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>1e-15</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>displacement</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>sigma</field>
+ <absolute_tolerance>2e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>swelling_stress</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>epsilon</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>velocity_liquid</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>square_ts_.*.vtu</regex>
+ <field>porosity</field>
+ <absolute_tolerance>1e-14</absolute_tolerance>
+ <relative_tolerance>1e-15</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1.gml b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1.gml
new file mode 100644
index 00000000000..d959c7c34c6
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1.gml
@@ -0,0 +1,31 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<?xml-stylesheet type="text/xsl" href="OpenGeoSysGLI.xsl"?>
+
+<OpenGeoSysGLI xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ogs="http://www.opengeosys.org">
+ <name>square_1x1_geometry</name>
+ <points>
+ <point id="0" x="0" y="0" z="0" name="origin"/>
+ <point id="1" x="0" y="1" z="0"/>
+ <point id="2" x="1" y="0" z="0"/>
+ <point id="3" x="1" y="1" z="0"/>
+ </points>
+
+ <polylines>
+ <polyline id="0" name="left">
+ <pnt>0</pnt>
+ <pnt>1</pnt>
+ </polyline>
+ <polyline id="1" name="right">
+ <pnt>2</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ <polyline id="2" name="bottom">
+ <pnt>0</pnt>
+ <pnt>2</pnt>
+ </polyline>
+ <polyline id="3" name="top">
+ <pnt>1</pnt>
+ <pnt>3</pnt>
+ </polyline>
+ </polylines>
+</OpenGeoSysGLI>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1_quad8_1e0.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1_quad8_1e0.vtu
new file mode 100644
index 00000000000..bd86778b4a8
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_1x1_quad8_1e0.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="7" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="96" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="112" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="136" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="404" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="500" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="524" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _QAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAABAAAAAAAAAAAAAAACAAAAAAAAAAAAAAAAAAAAA==wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAOA/AAAAAAAAAAAAAAAAAADgPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADgPwAAAAAAAPA/AAAAAAAAAAAAAAAAAADwPwAAAAAAAOA/AAAAAAAAAAA=QAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAwAAAAAAAAACAAAAAAAAAAUAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAQAAAAAAAAACAAAAAAAAAAIAAAAAAAAAA==AQAAAAAAAAAX
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_104_t_3.000000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_104_t_3.000000.vtu
new file mode 100644
index 00000000000..6eee1996226
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_104_t_3.000000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.00018393210664" RangeMax="0.00018393210664" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="580" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.95" RangeMax="0.95" offset="648" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.00026900080031" RangeMax="0.00026900080031" offset="716" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.00047928991978" RangeMax="0.00047928991978" offset="1000" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1152" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="0.0782365375" RangeMax="0.0782365375" offset="1220" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1308" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.00018393210664" offset="1464" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.00018393210664" RangeMax="0.00018393210664" offset="1604" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1848" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="0" RangeMax="0.1" offset="1936" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="0.1" RangeMax="0.1" offset="2008" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.95" RangeMax="0.95" offset="2076" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.00026900080031" RangeMax="0.00026900080031" offset="2156" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.00047928991978" RangeMax="0.00047928991978" offset="2528" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2716" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2776" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2836" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.95" RangeMax="0.95" offset="2900" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.00026900080031" RangeMax="0.00026900080031" offset="2964" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="3060" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="3148" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="3232" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="3292" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAApAAAAAAAAAA=eF77+DLWX6rso50Zl6JuCpOGPQMUKKx+KL2T2dlmw3GIfAsHqnzM65v/1miU2WxdBpE/x4YqL/DDfOvp1dU26/ZA5D8IoMqzeSr6c4T47VGBmr+SG1V+BU/lF7l/BTaPL0HkE9HMZ90YeX9XRrlNVj9E/oQ4mnymX//fUxk2M49B5J/xocrPjTzQMU200Yb/P0SeGc385ZenzWw+02IDACqFUxY=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5LSwOBd/ZpFNIAjO0gGg==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAtAAAAAAAAAA=eF5jYAABBxswxVC1Z0PipcWiCw33y0xsYj7nYAsVz9kDoWfZOCRA5DvEuY7M4C6EykdA5Vv2PEiGyE/n1b4l5FkKkW+YBVXXbfOhHCL/wltXMzTEHaYPKh9go5AGkX9Xmp1Sp5IF1T8Pqg4IkiDyYr6Gl08xFkHkHRZA5Btm7zGAus94c/ITH6FkqP7pUPM79hikQ+T//k+6NMOjCireBjU/Zc8DqH71x4VJ/f/rbQBvU1UhAQAAAAAAAAAAgAAAAAAAACABAAAAAAAAUQAAAAAAAAA=eF5bUNX9V+eUwv4FYNpgf4NmYI2zvuV+BihYUAmVr4TKa6DKJ1RA5CG0wX4HNPkFUPkFUPkGI9LMJyhPyP3o9uuguR+qPwGq/wCafgC5XGnFAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF47pvfK0599i/0xKH0SSh+H0gxg8KMeFw0A824ToQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAUgAAAAAAAAA=eF4zupDJPNeswV4fSptC6P22UL49lDaAisPkjaG0BZheYM/29P23KHv7PQ1NF29satazMYeKLwmyd5QS8rSxgvD3W0LpvSoPbVQ8wmwAXLE1hg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAARgAAAAAAAAA=eF5jYEAFD3bG+kuVfbRDE2aYw66om8KkYd99DSJ/mx/Ch8lzcID4Evb190HyD+H6F2yE8FdyQ9SrHIfodxaDqAcAmJkanQ==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAlgAAAAAAAAA=eF6TYI3zlyr7aFfKoaibwqRhzwADCiHSIZzTdq/MiwXLb5dBle++7FgcEJZhs60HIn+YDVV+Ifte5adf620SFCDmczKjyisc0vIxDuq0uXkcon8+O6p8y9mlvIIPq20e7YTI1wihym/Z/5r9JUfmnp5rEPkbLKjySWHsU7lmltkUHIPI3+ZHlb+0pyRWdmOrDQCAjkUxAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF6bNRMEdtrPwkEzEAAAeIUSPQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF6bNRMEdtrPIpMGAKksJHk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGwAAAAAAAAA=eF5LSgOBd/b5UDoDjU5B46dC6TQoDQCa3RyRAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAA9AAAAAAAAAA=eF5bqX57ydvvE/bs60w5NnH1oj1V4ZcWiy403M/A4HSs+vOE3Q3ncrtPJ+zaY5d4wSysZOmej74Q+bhtuct7diTZbOT7cZdj3wKbdbbxk4LDZ+3ZkgeR/+Ag/eXw1AqbmilNLhwn+/ccU5sk6a7Xb7MsCSIv/UThhl9Vs03xY2UjVrbuPWlyWxamuk+22Qq1v/GogRxvUamN5w0L2cNbp9jclJ67Mz4xw+ZQJUT+zq47Ie8ck/c0Hfz42Kxg1Z5THwVnMh3u2ZMfC5EvcLqnzJFUaMPinv4uPGKqTTHI/ONz96xIgMgLS7vlL1jfYAMA9TV/Sg==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAawAAAAAAAAA=eF6rr+7+q3NKYX89mDbYf1EtsMZZ33I/AxTcr4HIQ2iD/ceUUeW9yiDyENpgv6IuqvyvKog8hDbY36WDKh8DlY+BykuoosoLV0LkITTQfGNU+Wio/miofjM09wtA9QtA9TtoocoDAKm0WtY=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5LTwOBd/YAD1EDkw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJgAAAAAAAAA=eF6TOV4oc7wwYQ8DGCTsUUi9tFh0oeH+VUtr+0845NkAAL3kDI4=AQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_25_t_0.500000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_25_t_0.500000.vtu
new file mode 100644
index 00000000000..c52291a0bea
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_25_t_0.500000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.17481899845" RangeMax="0.17481899845" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="504" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="572" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.29488098968" RangeMax="0.29488098968" offset="652" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.085010243865" RangeMax="0.085010243865" offset="868" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="984" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0.30583802835" RangeMax="-0.30583802835" offset="1052" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1140" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.17481899845" offset="1268" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.17481899845" RangeMax="0.17481899845" offset="1388" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1608" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-0.45" RangeMax="0" offset="1696" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-0.45" RangeMax="-0.45" offset="1768" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="1836" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.29488098968" RangeMax="0.29488098968" offset="1920" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.085010243865" RangeMax="0.085010243865" offset="2156" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2312" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2372" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2432" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="2496" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.29488098968" RangeMax="0.29488098968" offset="2564" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2656" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2744" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2828" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2888" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAbAAAAAAAAAA=eF6r5X7+Pmfx7v11OleOdc88uJ8BCs5Y16fNW5BmUwWVb0CT388lz3DMI8amDCpfj1V/254aHPJM4tsl99yrgpuPLn8z5/6Zsh8FNtVEyqO7Lzhk/cyuM5NtcNkPkS+Gy9dgdX+ZDQBeX3wyAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAGgAAAAAAAAA=eF7bHS2w6++3p/Z70OjdaDQucRgNAFYALLM=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAgAAAAAAAAAA=eF4rOruHQ+fmyf1FMPozw9ZDZzftn8t2w+9OeqJNAVS8GEqXQ+V9OPqSLh8MtclH018MlZ8N1t+4pxBNvgQqf7N9wdzQ1SVw/cVo+l9zupSdks6ygekvgdJlOORh+guh8hO5f3UeLO3D0A+zvw8snwf3Xy6UzkPxf6ENAIVIofU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAANAAAAAAAAAA=eF771r1w2lPn6fbfwPRy+58aPw3eFG20Z4CCD1D5D1D5LyTKfyNgPiF5WptPSB4ApbJ+gg==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF779jp91c1Jl/f/gNJvoPRPKM0ABj/qcdEA0UAdlA==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPQAAAAAAAAA=eF4zupDJPNeswV4JSptB6P0wvgkqvd8FSjtBaRswvcCeAQzqbSC0wR4TFPEKGyh/vwWUhohPsAEAuLsqxg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOQAAAAAAAAA=eF5jYEAFNdzP3+cs3r0fTZihXufKse6ZB/dXQeXroHyYPIS/cX8tWH41XLwCykfXXw9VDwB9iTDDAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAggAAAAAAAAA=eF6r5X7+Pmfx7v3VOleOdc88uJ8BCrbkFsyanh9u0wKVb0OTF1p8l+/GvIU2ZVD5ejR51v77WydenbWnHipfgybPd/DSjslHqm2qoPKNaPJcCz7N3rjVzqYOh/lXZD8vvbOqDq6/Dk3ewTb4teLjQptaHPoN1DLK70YW2QAA/uJu5g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF47fwYE7uw/j4NmIAAA2ggZmQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF47fwYE7uw/TyYNAICyMzE=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHQAAAAAAAAA=eF7bFS2w6++3p/aHoPReKL0PBw2T3wWlAQU7J8s=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAjgAAAAAAAAA=eF4rPLuHQ+fmyf2FULrgM8PWQ2c37b+2jV1DiD3IphoqXgOly6Dyn0zZeg9Jz7YpgIoXw2io/JZGz80eS6bsKYKKZ0HpHKi8x+MPd5RiSm1g9pZA6VKY/RKy70rmm9uUosmXQOXTK2UfyU8vt4GZX4Qmn3Ll3I6LLdlweZj+Qqh8hLCWVI11jg0AdMON1g==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAUwAAAAAAAAA=eF772r1w2lPn6fZfwfRy+x8aPw3eFG20Z4CCn1D5n1D5v2jyH6HyH6Hy39Dkf0Hlf0Hl/6DJf4LKf4LKfyFRP7r96Pp/QOV/QOV/o8kDAIyucJs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF7bHS2w6++3p/YAFT4E+A==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF4rOLuHQ+fmyf1FULrkM8PWQ2c37X9Y/WzKhOlZNgAy/xLeAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_50_t_1.000000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_50_t_1.000000.vtu
new file mode 100644
index 00000000000..3f41b5c6c98
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_50_t_1.000000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.22090710401" RangeMax="0.22090710401" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="504" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="572" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="644" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.37368644876" RangeMax="0.37368644876" offset="864" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="980" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0.027773126425" RangeMax="-0.027773126425" offset="1048" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1136" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.22090710401" offset="1264" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.22090710401" RangeMax="0.22090710401" offset="1388" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1616" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-1" RangeMax="0" offset="1704" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-1" RangeMax="-1" offset="1772" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="1840" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="1924" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.37368644876" RangeMax="0.37368644876" offset="2168" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2312" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2372" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2432" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="2496" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="2564" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2656" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2744" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2828" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2888" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAagAAAAAAAAA=eF77Z3e2+eCBzftPOb1+fKju1H4GKDhjXZ82b4Hbnj9Q+XNo8paf3p25KRNv8w0qfwarfL/NHyLNP4smfzPn/pmyHwU2P3DYzyS+XXLPvS4bXO6D6G/A6X6I/WF7fuCwHyI/zQYAMyaLxA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEwAAAAAAAAA=eF5bVr3vk0v+YftlFNIA3Bgotw==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAhAAAAAAAAAA=eF67fN13q9WM/ftPQ+m9cwUV/scetJ/LdsPvTrrjnotQcRi9Byqvy5n55ty7cJuzUPGzaPolwPLtNheg4keg9HY088/gMP8lp0vZKeksm8tQcZg5+6Dyt9sXzA1d3QKXv4Rm/0ew/qo9MPGraPYvAtsfuAfm7lNQ+gCK+ybaAADOMp/GAQAAAAAAAAAAgAAAAAAAACABAAAAAAAANQAAAAAAAAA=eF67e1mq2WTiTvu7YPqk/Qu5kKTnOpftGaDgLgH5e1D5e0TKPyPRfErl0e1/TqJ+APdrbDE=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF6T33LhW1/BnP0GUJoXSktCaQYw+FGPiwYACpQZng==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPgAAAAAAAAA=eF5TuZDJPNeswV4BSttA6P06UL4FlHaFijtBaS8oDZFfYM8ABv02ELrVBlWcocEMwt9vCKUhwh02ALbJKoU=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPAAAAAAAAAA=eF5jYEAFv+3ONh88sHk/mjDDWafXjw/Vndr/DSoP48Pkz4D5u/b/Acsvhot/hfJh6v+g6N+1HwDUyjdhAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAiQAAAAAAAAA=eF77Z3e2+eCBzfuPOb1+fKju1H4GKNgfcuBB04fCPYz2EPlbaPKLNn+QiS7t3fMFqv8UmnygNPMZwQ9TbX5A5S+gyfMzvJaaabHI5i9U/hKa/Dv27af49DJs/kHlz6LJt22/y//NLs7mG1T+PJp8Qar2hKuqU+Dmn0OT3/PHrmUzS+geACM4dqY=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEgAAAAAAAAA=eF5jYACBD/sZ8NK4AQAzGwa9AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEAAAAAAAAAA=eF5jYACBD/sZyKQBjmcNeQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHAAAAAAAAAA=eF5bVr3vk0v+YfvVUHo5lF6Bg78MTRwAppQkPQ==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAlQAAAAAAAAA=eF67ed13q9WM/ftPQ+l9cwUV/scetLfjmfc4bWb2nkdQ8YdQ+jhU3uXoRCbulW17zkLFT0LpXVD5hW/nb+B9MMHmHFT8EJQ+CpU/KczO5Mcxx+YSVPwalN4Plb8iWWk5dWaSzRU0/Xug8oWV+XvlXcLh+mH2HITKx3VM1FH83weXvwGld0Llta6F2zFoB+wBAMbZies=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAASQAAAAAAAAA=eF67e1mq2WTiTvu7YPqk/TO5kKTnOpftGaDgMVT+MVT+A5r8Q6j8Q6j8SzT5B1D5B1D5V2jy96Dy93DIo+t/QUA/qfIAEq5gWg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF5bXr3vk0v+YXsAFe0EiA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF67cN13q9WM/fvPQek9cwUV/scetP9xxCJO6VS+DQA5qBHAAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_54_t_2.000000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_54_t_2.000000.vtu
new file mode 100644
index 00000000000..3015a4b0fe9
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_54_t_2.000000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.22090710401" RangeMax="0.22090710401" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="500" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="568" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="644" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.37368644876" RangeMax="0.37368644876" offset="876" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="992" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-6.7010057245e-17" RangeMax="4.4609450118e-17" offset="1060" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1168" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.22090710401" offset="1304" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.22090710401" RangeMax="0.22090710401" offset="1428" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1648" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-1" RangeMax="0" offset="1736" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-1" RangeMax="-1" offset="1808" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="1876" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="1956" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.37368644876" RangeMax="0.37368644876" offset="2200" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2344" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2404" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2464" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.15183317044" RangeMax="0.15183317044" offset="2528" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.22107160338" RangeMax="0.22107160338" offset="2596" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2688" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2776" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2860" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2920" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAZwAAAAAAAAA=eF77YXe2+eCBzfvPOr1+fKju1H4GKDhjXZ82b4GbzQ+o/Hk0+ZKo6mWt4sU2v6HyZ9DkYeAHDvmbOffPlP0I2PMTKn8Oq/1lNr9w2A/R32CD3/yCPb9wyEPMb4Prx+7/tj0A21iHSQ==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAGAAAAAAAAAA=eF5bVb3vk0v+YfuVBOhVaDS6PADg0CjVAQAAAAAAAAAAgAAAAAAAACABAAAAAAAAjQAAAAAAAAA=eF47dd13q9WM/ftPQ+l9cwUV/scetJ/LdsPvTrqjzRmo+HkovRcqb9zNc9Q5PtfmEpr8Hqj8Gev6tHkL3CzPQsVPQOltUPnXnC5lp6S99sDMvwald0Pl54DtL7S5ARW/jGb/O7D+Krj7YO7fD5X/BpbP2nMKKg5zx04U8xttYOIwdQdQ/N+4BwCpDJzwAQAAAAAAAAAAgAAAAAAAACABAAAAAAAANQAAAAAAAAA=eF67e1mq2WTiTvu7YPqk/Qu5kKTnOpftGaDgLgH5e1D5e0TKPyPRfErl0e1/TqJ+APdrbDE=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAMAAAAAAAAAA=eF67dzcuwXB7p43AxAv+6SqT9wjk/ZrDH99p8/1Ek2FFwOQ9DGDwox4XDQBC5xVtAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAARAAAAAAAAAA=eF6zuJDJPNeswd4FSltC6P0BUL4JlFaCimtBaQUobQamF9gzgIBD3x4wzZCyxwpZnMHBxg7C3w8V3w8Rn7wHAPPXLCs=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOwAAAAAAAAA=eF5jYEAF3+3ONh88sHk/mjDDeafXjw/Vndr/Cyp/BspHld+1H6J/MVz8F5R/Dqr+J1T/aah6ANTiN2E=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAggAAAAAAAAA=eF77bne2+eCBzftPOL1+fKju1H4GKIj7fLKybF6FzW+o/FU0+Vv7m16yhrXv+Y9DP7eVS+6phpI9P6DyF9DkV7l7bOc7vmDPL6j8RTT529qfjgvxBcHlz6LJLyv9ZHBsUssemPvOo8l3CenkffVOgLv/DJq8xxbxT4bOtTYAk3F4pg==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF778x8E3u//g4NmIAAAz4selQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF778x8E3u//QyYNAMwpPSk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGwAAAAAAAAA=eF5bWb3vk0v+YfuNUHollF5NgIapAwCqpCRZAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAlgAAAAAAAAA=eF47ed13q9WM/fuPQumDcwUV/scetP/BuT6W50ORzQ2o+AUofRoqP2NyClP658Y9V6HiJ6H0Tqi8W+QELoHu/D1noeKHoPQRqHx/1NN7x8/P3HMCKn4Opg4qr3NW9eLDY957zkPFT6CZf4y59OCuZfV77kLFr6G5f9OkSU9ShSJtYPrOQOk9UPkLDIkGpU1lNgCSeY81AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAASQAAAAAAAAA=eF67e1mq2WTiTvu7YPqk/TO5kKTnOpftGaDgMVT+MVT+A5r8Q6j8Q6j8SzT5B1D5B1D5V2jy96Dy93DIo+t/QUA/qfIAEq5gWg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF5bWb3vk0v+YXsAFf0Eig==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF47c913q9WM/fvPQek9cwUV/scetC9i33zj+pRgGwA3CxHRAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_5_t_0.100000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_5_t_0.100000.vtu
new file mode 100644
index 00000000000..acdc6e75f58
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_5_t_0.100000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.038424182482" RangeMax="0.038424182482" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="484" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="552" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.080267973049" RangeMax="0.080267973049" offset="624" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0" RangeMax="0" offset="816" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="880" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-0.077526734592" RangeMax="-0.077526734592" offset="948" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1032" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.038424182482" offset="1160" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.038424182482" RangeMax="0.038424182482" offset="1272" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1476" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-0.01" RangeMax="0" offset="1564" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-0.01" RangeMax="-0.01" offset="1640" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="1712" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.080267973049" RangeMax="0.080267973049" offset="1796" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2028" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2088" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2148" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2208" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="2272" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.080267973049" RangeMax="0.080267973049" offset="2340" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2432" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2520" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2604" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2664" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAXAAAAAAAAAA=eF4zubDwavGl2fstoDQDFJyxrk+bt8DNxgQqboYmHxyyfmbXGWe4vClW+WAbc7zmq+E0H2a/MQ75mzn3z5T9cNhjjkMeYn8yTv1uyv7Gm0o8bQxwyMMAAGL8eig=AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEwAAAAAAAAA=eF4zDbVOTUt7Z29KIQ0A9egcPA==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAbQAAAAAAAAA=eF5buSpMqPnAqv2roLSRcXqhNO+C/XPZbvjdSXe0WQEVX4kmP5n7V+fBUjub5VDxFWjyE8HyfjYwc1dDaRMU8xVtVpJp/xtOl7JT0lZ70M1HdV8cTv2F7JtvXJ/ibLMMzf0GUHkGKAAAwy6FtA==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAADQAAAAAAAAA=eF5jYBgF+AAAASAAAQ==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHQAAAAAAAAA=eF6bI7r46Kkbm/fPQqNnQGkGMPhRj4sGADHUGok=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF7TvZDJPNeswd4ISptC6P1Q2l4fShugiu+Hiu+3BtML7BkgoAFMHkizsUIWd0ixgarbbwmlIcpdbADGFSrgAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYEAFZhcWXi2+NHs/mjCDOVTcCEqbo6mzAPO795tCaZi4CZQPM9cESptBxQHoGzLKAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAeAAAAAAAAAA=eF4zurDwavGl2ftNoTQDFEy2KNruH+dp4wgVt0WTP2opoBe+M9fGHCpuiSbf1Xsl6SZLrI0hDvl7E4+enJljuMcIh/2zZgZlLdAIt7GGiluhyT8L/RJhlO9oA9NvhiZ/jHdlnWif+R5c8oZbo3kzzkfZAABzl24jAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAFgAAAAAAAAA=eF5rEFnn/rCqZX8DDpqBAAAAHi0QnQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF5rEFnn/rCqZX8DmTQAKHwhOQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHQAAAAAAAAA=eF4zDbVOTUt7Z28DpS3QaBMobY7GN4PSACN9GSc=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAiwAAAAAAAAA=eF5btipMqPnAqv1LobSBcXqhNO+C/V5VBifmNzvbbIKKb4HSplD5Xyd/HAi7m2azAiq+CkqbQOWdpty7Wh8UZrMcTd4QKs+xs66w7JzWHlzyHXOrUqsvBNqshYqvQ7Nf0HP+9ZJLVjYroeIw2ggqv87VeEObnNGeFWjyMPNLjn2X7Z4fYgMA49978g==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAADAAAAAAAAAA=eF5jYBjZAAABAAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF4zDbVOTUt7Zw8ADEMDJA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF5bsSpMqPnAqv0robShcXqhNO+C/cummBqIvXGwAQAE9Q6uAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_79_t_2.500000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_79_t_2.500000.vtu
new file mode 100644
index 00000000000..0b5d8bd31c5
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_79_t_2.500000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.17466467123" RangeMax="0.17466467123" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="516" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="584" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.29494811038" RangeMax="0.29494811038" offset="660" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.084530953945" RangeMax="0.084530953945" offset="884" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1032" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="0.31047806824" RangeMax="0.31047806824" offset="1100" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1184" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.17466467123" offset="1316" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.17466467123" RangeMax="0.17466467123" offset="1440" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1664" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-0.45" RangeMax="0" offset="1752" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-0.45" RangeMax="-0.45" offset="1824" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="1892" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.29494811038" RangeMax="0.29494811038" offset="1980" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.084530953945" RangeMax="0.084530953945" offset="2220" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2380" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2440" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2500" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.68640028325" RangeMax="0.68640028325" offset="2564" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.29494811038" RangeMax="0.29494811038" offset="2632" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2724" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2812" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2896" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="2956" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAdQAAAAAAAAA=eF5bIurUF7po9361Jyd0uycf3M8ABWes69PmLQizWQKV10CT388lz3DMo8ZmPlReE02+5l/bs/ObeuH6VbCan7ZnMVReHat8mc0CqLwWVvt74Oaju/9mzv0zZT8K9izC4f53jt7FM1+32szD4X6Y/QADAHFjAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAFgAAAAAAAAA=eF47Gy2w6++3p/ZnyaTPQ2kAbsgtVA==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAhgAAAAAAAAA=eF4rO7uHQ+fmyf0lUFri2JJTyfc37Z/LdsPvTnqgTTlUvBhKi0LlvTn6ki4fLLUpQtMvBpXf2FGtHnijDS5fAKU1oPLzweYn7ilBMx/V/kKbCqh4GZQWgcq7gO1vtamEilegyf/kdCk7JZ21pxhNXgYqn+KWmDP3TgPcfBgth2Y/AOQloOY=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAATAAAAAAAAAA=eF7rL1lmE6g03b4fTC+3r7he83ZP5EZ7BiiYAJWfAJWvQpOfDJWfjEN+ElR+ElS+gETzCcn3E3A/uv1laPJNUPkmqHw+mjwAC0pkZw==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5zyPBLv3/vsr0DlHaC0h5QmgEMftTjogFzARaXAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQAAAAAAAAAA=eF5zuJDJPNeswd4IShtA6P1WqHx7Pag4TB6qfr8ZmF5gzwAGCXsgdMUeExTxhj0Q8xZA9S/YDxEvsAEA6e0r0Q==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYEAFi0Wd+kIX7d6PJsyg8eSEbvfkg/vnQuVhfJi8Gpi/cT9E/2q4+FwoXx2qHmY+TD0APcMqqA==AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAhwAAAAAAAAA=eF6bL+rUF7po9361Jyd0uycf3M8ABVPOCp9KFq+wWQuV10GTT92wa+3TLS175kLltdDktTjfGk4rn2yzACqvjSZ/oqojbWFm2p7FUHlNNHnhGzvXTd1RZrMUKq+KJj/hEp+Byc6GPfNwmD+5Iq5EXaLbZiFUXgNN/tn3rdffp7XaAAAoM2WHAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF7bcwYE7uzfg4NmIAAAypgZTQ==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEQAAAAAAAAA=eF7bcwYE7uzfQyYNAGtSMpk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAIQAAAAAAAAA=eF47GS2w6++3p/ZXoPRFKH0JSp+F0ueh9Dk0GgAY8yhZAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAkwAAAAAAAAA=eF4rPbuHQ+fmyf3lUJr/2JJTyfc37Z/EEJeg2lFk0wUVb4HSXFB5vbip3y1v1+/JgYqXQGlFqPz36YYlfJ97beqg4jVQWgsqf+fA9hv9tol7KqDipVBaCCp/0rn6oHdVoQ3M3CIorQOV3zYpdU35rso9VWjuh+lfo7j6fMP1Frh+mP2SUPnYS6IyE8sabADbipGLAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAVQAAAAAAAAA=eF7rK1lmE6g03b4PTC+3r7pe83ZP5EZ7BiiYBJWfBJVvRJOfBZWfBZWvQJOvgMpXQOVT0eT70OyvQZOfDJWfDJUvIKAf3f39UPl+qHw5mjwA7pNZLg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF47Fy2w6++3p/YAFdYFCw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF4rObuHQ+fmyf2lUFry2JJTyfc37RcQVf46manIBgAkLRCmAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_99_t_2.900000.vtu b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_99_t_2.900000.vtu
new file mode 100644
index 00000000000..001b76f3793
--- /dev/null
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square_ts_99_t_2.900000.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="465" format="appended" RangeMin="34" RangeMax="125" offset="0" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="232" />
+ <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.038283063011" RangeMax="0.038283063011" offset="316" />
+ <DataArray type="Float64" Name="porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="520" />
+ <DataArray type="Float64" Name="saturation_ip" NumberOfTuples="9" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="588" />
+ <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.080373417303" RangeMax="0.080373417303" offset="660" />
+ <DataArray type="Float64" Name="swelling_stress_ip" NumberOfComponents="4" NumberOfTuples="9" format="appended" RangeMin="0.00047928991978" RangeMax="0.00047928991978" offset="896" />
+ <DataArray type="Float64" Name="transport_porosity_ip" NumberOfTuples="9" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1048" />
+ </FieldData>
+ <Piece NumberOfPoints="8" NumberOfCells="1" >
+ <PointData>
+ <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="0.077531194322" RangeMax="0.077531194322" offset="1116" />
+ <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0.011195857369" RangeMax="0.031666666667" offset="1200" />
+ <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0.038283063011" offset="1328" />
+ <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.038283063011" RangeMax="0.038283063011" offset="1452" />
+ <DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.5" RangeMax="0.5" offset="1676" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-0.01" RangeMax="0" offset="1764" />
+ <DataArray type="Float64" Name="liquid_pressure_interpolated" format="appended" RangeMin="-0.01" RangeMax="-0.01" offset="1840" />
+ <DataArray type="Float64" Name="saturation" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="1912" />
+ <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0.080373417303" RangeMax="0.080373417303" offset="1996" />
+ <DataArray type="Float64" Name="swelling_stress" NumberOfComponents="4" format="appended" RangeMin="0.00047928991978" RangeMax="0.00047928991978" offset="2240" />
+ <DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2428" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="2488" />
+ <DataArray type="Float64" Name="porosity_avg" format="appended" RangeMin="0.5" RangeMax="0.5" offset="2548" />
+ <DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0.94999999782" RangeMax="0.94999999782" offset="2612" />
+ <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="0.080373417303" RangeMax="0.080373417303" offset="2680" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2772" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2860" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2944" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="3004" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAANEBAAAAAAAAiwAAAAAAAAA=eF6t0LEKwlAMBdB/ydxF7NRfEQlRYwn0JY8kRYr0332Di0s71PHeC2e4bxBNHp1STLFaS0jutAQMl/fPaP5gh+HcgVJhGCBkLIRSoTVzubGjPfFupZqyZgP6tdsjKOfvtuGcdp1qbiG5HFPSSaNRif/x4sXTJDpipHPEsau4hkzbP/Xrdf0A+2CrFQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPWM9Y1MjbWTTc3TktKNTc3tjQCADOEBPU=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAeAAAAAAAAAA=eF6L/cHrL3Zu9v4Tve68+Ytm72eAgps598+U/QiwSYHKn8IqX2BTgEP/O0fv4pmvS22yofK7sOp32JMBlT+Mw/x0qPw5NPntacIHqtcW2lRA5VegyRv6Pe7pbyyB278PTf5v9LpMl1dNNjFQ+QtY7W+wAQD4knG2AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEwAAAAAAAAA=eF4zDbVOTUt7Z29KIQ0A9egcPA==AQAAAAAAAAAAgAAAAAAAACABAAAAAAAAkAAAAAAAAAA=eF4LWBUm1Hxg1f4QKH3w+ro5sQYL9r/mdCk7Je1lEwAVd4HSe6HyP8HyWXD5YCh9BCqf5ZaYM/dOgQ1Mnw+UfgeVZ+AC6bfaAxMPhNKnUOzPskmCigdB6UNQ+aAuQ2sPo2y4fByU3geV/9j93M++J9/GHSoeCaXPQuXZtC8GL5GstYlA8/8BFPurbAAHwpEfAQAAAAAAAAAAgAAAAAAAACABAAAAAAAAUQAAAAAAAAA=eF5bUNX9V+eUwv4FYNpgf4NmYI2zvuV+BihYUAmVr4TKa6DKJ1RA5CG0wX4HNPkFUPkFUPkGI9LMJyhPyP3o9uuguR+qPwGq/wCafgC5XGnFAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAEAAAAAAAAAA=eF5jYACBB/YMFNIAVS4KGA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHAAAAAAAAAA=eF5TyuGYIXpzs708lJZDoxnA4Ec9LhoAlsISCg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPgAAAAAAAAA=eF7TvpDJPNeswV4XSptC6P1WUL4DmrgZlDZG4S+wZwCBhiQbCJ1qY44szsBgYwnh74eK74cIp9gAANBTKsY=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOwAAAAAAAAA=eF5jYEAF2T94/cXOzd6PJsxwqtedN3/R7P2pUPldUD6qfPf+RLB8N1w8F8o/DFWfAdW/BqoeAMpqJ9o=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAhQAAAAAAAAA=eF4L/sHrL3Zu9v7Tve68+Ytm72eAgs4/Nx6fE4q06YDKL0GT32UXab3QvcymCip/HE1+pczia342VTYhUPk7aPLHdOpzWqc12iThsF/XdIFgin6ZTR5UfieavM+Xx/e/rUrekwGVv4Amz9Q3m+tLcb5NLlR+F5p8yAZGo019HTYACfJeCA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHwAAAAAAAAA=eF77+x8E7tuzMoDAA3tGKM0EpRnQxGE0TBwAfAsO+g==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAFgAAAAAAAAA=eF6bJbjO/WFVy/5ZOGgGAgAAMOkQ+Q==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEwAAAAAAAAA=eF6bJbjO/WFVy/5ZZNIAQnQh8Q==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHQAAAAAAAAA=eF4zDbVOTUt7Z28DpS3QaBMobY7GN4PSACN9GSc=AQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAlgAAAAAAAAA=eF6LWBUm1Hxg1f54KL3m+ro5sQYL9mseMijp9gy26YOKZ0PpuVD5jZM4QzOUC20coeLRUPo2VH6v5kT2FL0SmxiouD+UPgKVv8T5etkXvmqbUKi4B5r9nzQcFA8LFNp4QcVDoPRnqLx2SlGEpFPcnjyoeCSU3gGVP8I1P+u8YYYNzNxEKH0AKu83+23H68dNNgDMboLIAQAAAAAAAAAAgAAAAAAAAAABAAAAAAAAawAAAAAAAAA=eF6rr+7+q3NKYX89mDbYf1EtsMZZ33I/AxTcr4HIQ2iD/ceUUeW9yiDyENpgv6IuqvyvKog8hDbY36WDKh8DlY+BykuoosoLV0LkITTQfGNU+Wio/miofjM09wtA9QtA9TtoocoDAKm0WtY=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBB/YAAgcBIA==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF4zDbVOTUt7Zw8ADEMDJA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF4LWhUm1Hxg1f4gKH3i+ro5sQYL9hfrvY0S3ppjAwAQKhBOAQAAAAAAAAAAgAAAAAAAAMAAAAAAAAAAIQAAAAAAAAA=eF5jYMAHPtjjlcaQh/Fx6XuAJo7ORxfHZT5CHwB6OwuXAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHgAAAAAAAAA=eF5jYIAARijNDKWZoDQrlGaH0mxQmgVKAwADIAAdAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF7jYIAAAABIAAk=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF4TBwAAGAAY
+ </AppendedData>
+</VTKFile>
--
GitLab