diff --git a/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp b/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp
index 33a3ad57e1886a57ac80b19b2d96b1493a39b625..6e90fea693e62af4730849391d8947db08613684 100644
--- a/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp
+++ b/MaterialLib/MPL/Properties/SaturationDependentSwelling.cpp
@@ -83,6 +83,11 @@ PropertyDataType SaturationDependentSwelling::value(
double const delta_S_eff = S_eff - S_eff_prev;
+ if (delta_S_eff == 0.)
+ {
+ return delta_sigma_sw; // still being zero.
+ }
+
// \Delta\sigma_{sw} = - \sum_i k_i (\lambda p S_{eff}^{(\lambda_i - 1)}
// e_i \otimes e_i \Delta S_L / (S_{max} - S_{min}), where
// e_i \otimes e_i is a square matrix with e_i,0^2 e_i,0*e_i,1 etc.
diff --git a/ProcessLib/TH2M/IntegrationPointData.h b/ProcessLib/TH2M/IntegrationPointData.h
index 0ec28659c4162d5fed6c23fe5fa4a5d9a64958f3..87e933cd6cb98cd65990e8db8e69abcce574f825 100644
--- a/ProcessLib/TH2M/IntegrationPointData.h
+++ b/ProcessLib/TH2M/IntegrationPointData.h
@@ -40,6 +40,7 @@ struct IntegrationPointData final
static const int kelvin_vector_size =
MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
sigma_eff.setZero(kelvin_vector_size);
+ sigma_sw.setZero(kelvin_vector_size);
eps.setZero(kelvin_vector_size);
eps_m.setZero(kelvin_vector_size);
eps_m_prev.resize(kelvin_vector_size);
@@ -53,6 +54,7 @@ struct IntegrationPointData final
DisplacementDim, NPoints * DisplacementDim>
N_u_op;
typename BMatricesType::KelvinVectorType sigma_eff, sigma_eff_prev;
+ typename BMatricesType::KelvinVectorType sigma_sw, sigma_sw_prev;
typename BMatricesType::KelvinVectorType eps, eps_prev;
typename BMatricesType::KelvinVectorType eps_m, eps_m_prev;
@@ -183,6 +185,7 @@ struct IntegrationPointData final
eps_prev = eps;
eps_m_prev = eps_m;
sigma_eff_prev = sigma_eff;
+ sigma_sw_prev = sigma_sw;
s_L_prev = s_L;
rho_G_h_G_prev = rho_G_h_G;
@@ -199,6 +202,50 @@ struct IntegrationPointData final
material_state_variables->pushBackState();
}
+ MathLib::KelvinVector::KelvinMatrixType<DisplacementDim>
+ computeElasticTangentStiffness(
+ double const t, ParameterLib::SpatialPosition const& x_position,
+ double const dt, double const temperature_prev,
+ double const temperature)
+ {
+ namespace MPL = MaterialPropertyLib;
+
+ MPL::VariableArray variable_array;
+ MPL::VariableArray variable_array_prev;
+
+ auto const null_state = solid_material.createMaterialStateVariables();
+
+ using KV = MathLib::KelvinVector::KelvinVectorType<DisplacementDim>;
+
+ variable_array[static_cast<int>(MPL::Variable::stress)].emplace<KV>(
+ KV::Zero());
+ variable_array[static_cast<int>(MPL::Variable::mechanical_strain)]
+ .emplace<KV>(KV::Zero());
+ variable_array[static_cast<int>(MPL::Variable::temperature)]
+ .emplace<double>(temperature);
+
+ variable_array_prev[static_cast<int>(MPL::Variable::stress)]
+ .emplace<KV>(KV::Zero());
+ variable_array_prev[static_cast<int>(MPL::Variable::mechanical_strain)]
+ .emplace<KV>(KV::Zero());
+ variable_array_prev[static_cast<int>(MPL::Variable::temperature)]
+ .emplace<double>(temperature_prev);
+
+ auto&& solution =
+ solid_material.integrateStress(variable_array_prev, variable_array,
+ t, x_position, dt, *null_state);
+
+ if (!solution)
+ {
+ OGS_FATAL("Computation of elastic tangent stiffness failed.");
+ }
+
+ MathLib::KelvinVector::KelvinMatrixType<DisplacementDim> C =
+ std::move(std::get<2>(*solution));
+
+ return C;
+ }
+
typename BMatricesType::KelvinMatrixType updateConstitutiveRelation(
MaterialPropertyLib::VariableArray& variable_array,
double const t,
diff --git a/ProcessLib/TH2M/LocalAssemblerInterface.h b/ProcessLib/TH2M/LocalAssemblerInterface.h
index b66e9f7540f87905a41ad20071b46addfe72a518..53a5b52df1de21e3ad24cede912fae49afbd8299 100644
--- a/ProcessLib/TH2M/LocalAssemblerInterface.h
+++ b/ProcessLib/TH2M/LocalAssemblerInterface.h
@@ -32,6 +32,14 @@ struct LocalAssemblerInterface : public ProcessLib::LocalAssemblerInterface,
std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
std::vector<double>& cache) const = 0;
+ virtual std::vector<double> getSwellingStress() const = 0;
+
+ virtual std::vector<double> const& getIntPtSwellingStress(
+ const double t,
+ std::vector<GlobalVector*> const& x,
+ std::vector<NumLib::LocalToGlobalIndexMap const*> const& dof_table,
+ std::vector<double>& cache) const = 0;
+
virtual std::vector<double> getEpsilon() const = 0;
virtual std::vector<double> const& getIntPtEpsilon(
diff --git a/ProcessLib/TH2M/TH2MFEM-impl.h b/ProcessLib/TH2M/TH2MFEM-impl.h
index ce77e8a5eb9d008461367b5c47e936d8df608aab..9e9314ee5747ae52eb52629c7efe48b5f4ed6b45 100644
--- a/ProcessLib/TH2M/TH2MFEM-impl.h
+++ b/ProcessLib/TH2M/TH2MFEM-impl.h
@@ -160,6 +160,7 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
GlobalDimVectorType const gradT = gradNp * temperature;
MPL::VariableArray vars;
+ MPL::VariableArray vars_prev;
vars[static_cast<int>(MPL::Variable::temperature)] = T;
vars[static_cast<int>(MPL::Variable::phase_pressure)] = pGR;
vars[static_cast<int>(MPL::Variable::capillary_pressure)] = pCap;
@@ -177,6 +178,8 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
vars, pos, t, std::numeric_limits<double>::quiet_NaN());
vars[static_cast<int>(MPL::Variable::liquid_saturation)] = ip_data.s_L;
+ vars_prev[static_cast<int>(MPL::Variable::liquid_saturation)] =
+ ip_data.s_L_prev;
auto const chi = [&medium, pos, t, dt](double const s_L)
{
@@ -215,6 +218,7 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
MathLib::KelvinVector::kelvinVectorToSymmetricTensor(
sigma_total));
}
+ // Set volumetric strain rate for the general case without swelling.
vars[static_cast<int>(MPL::Variable::volumetric_strain)] =
Invariants::trace(eps);
vars[static_cast<int>(MPL::Variable::equivalent_plastic_strain)] =
@@ -251,6 +255,24 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
// solid phase compressibility
ip_data.beta_p_SR = (1. - ip_data.alpha_B) / K_S;
+ // If there is swelling stress rate, compute swelling stress.
+ if (solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate))
+ {
+ auto& sigma_sw = ip_data.sigma_sw;
+ auto const& sigma_sw_prev = ip_data.sigma_sw_prev;
+
+ sigma_sw = sigma_sw_prev;
+
+ using DimMatrix = Eigen::Matrix<double, 3, 3>;
+ auto const sigma_sw_dot =
+ MathLib::KelvinVector::tensorToKelvin<DisplacementDim>(
+ solid_phase
+ .property(MPL::PropertyType::swelling_stress_rate)
+ .template value<DimMatrix>(vars, vars_prev, pos, t,
+ dt));
+ sigma_sw += sigma_sw_dot * dt;
+ }
+
// solid phase linear thermal expansion coefficient
ip_data.alpha_T_SR = MathLib::KelvinVector::tensorToKelvin<
DisplacementDim>(MaterialPropertyLib::formEigenTensor<3>(
@@ -263,6 +285,8 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
ip_data.beta_T_SR = Invariants::trace(ip_data.alpha_T_SR);
double const T_dot = NT.dot(temperature_dot);
+ double const T_prev = T - T_dot * dt;
+
MathLib::KelvinVector::KelvinVectorType<DisplacementDim> const
dthermal_strain = ip_data.alpha_T_SR * T_dot * dt;
@@ -271,6 +295,15 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
auto& eps_m_prev = ip_data.eps_m_prev;
eps_m.noalias() = eps_m_prev + eps - eps_prev - dthermal_strain;
+
+ if (solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate))
+ {
+ auto const C_el =
+ ip_data.computeElasticTangentStiffness(t, pos, dt, T_prev, T);
+ eps_m.noalias() +=
+ C_el.inverse() * (ip_data.sigma_sw - ip_data.sigma_sw_prev);
+ }
+
vars[static_cast<int>(MaterialPropertyLib::Variable::mechanical_strain)]
.emplace<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
eps_m);
@@ -312,7 +345,6 @@ TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
auto const rhoSR = rho_ref_SR;
#endif // NON_CONSTANT_SOLID_PHASE_VOLUME_FRACTION
- auto const T_prev = T - T_dot * dt;
ip_cv.C = ip_data.updateConstitutiveRelation(vars, t, pos, dt, T_prev);
// constitutive model object as specified in process creation
@@ -814,6 +846,11 @@ std::size_t TH2MLocalAssembler<
return ProcessLib::setIntegrationPointScalarData(values, _ip_data,
&IpData::s_L);
}
+ if (name == "swelling_stress_ip")
+ {
+ return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
+ values, _ip_data, &IpData::sigma_sw);
+ }
if (name == "epsilon_ip")
{
return ProcessLib::setIntegrationPointKelvinVectorData<DisplacementDim>(
@@ -826,7 +863,7 @@ template <typename ShapeFunctionDisplacement, typename ShapeFunctionPressure,
typename IntegrationMethod, int DisplacementDim>
void TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
IntegrationMethod, DisplacementDim>::
- setInitialConditionsConcrete(std::vector<double> const& local_x,
+ setInitialConditionsConcrete(std::vector<double> const& local_x_data,
double const t,
bool const /*use_monolithic_scheme*/,
int const /*process_id*/)
@@ -835,14 +872,60 @@ void TH2MLocalAssembler<ShapeFunctionDisplacement, ShapeFunctionPressure,
gas_pressure_size + capillary_pressure_size + temperature_size +
displacement_size;
- assert(local_x.size() == matrix_size);
+ assert(local_x_data.size() == matrix_size);
- updateConstitutiveVariables(
- Eigen::Map<Eigen::VectorXd const>(local_x.data(), local_x.size()),
- Eigen::VectorXd::Zero(matrix_size), t, 0);
+ auto const local_x = Eigen::Map<Eigen::VectorXd const>(local_x_data.data(),
+ local_x_data.size());
+ auto const capillary_pressure =
+ local_x.template segment<capillary_pressure_size>(
+ capillary_pressure_index);
+
+ auto const temperature =
+ local_x.template segment<temperature_size>(temperature_index);
+
+ constexpr double dt = std::numeric_limits<double>::quiet_NaN();
+ auto const& medium = *_process_data.media_map->getMedium(_element.getID());
+ auto const& solid_phase = medium.phase("Solid");
unsigned const n_integration_points =
_integration_method.getNumberOfPoints();
+
+ ParameterLib::SpatialPosition pos;
+ pos.setElementID(_element.getID());
+ for (unsigned ip = 0; ip < n_integration_points; ip++)
+ {
+ pos.setIntegrationPoint(ip);
+ MPL::VariableArray vars;
+
+ auto& ip_data = _ip_data[ip];
+ auto const& Np = ip_data.N_p;
+ auto const& NT = Np;
+
+ double const pCap = Np.dot(capillary_pressure);
+ vars[static_cast<int>(MPL::Variable::capillary_pressure)] = pCap;
+
+ double const T = NT.dot(temperature);
+ vars[static_cast<int>(MPL::Variable::temperature)] = T;
+
+ ip_data.s_L_prev =
+ medium.property(MPL::PropertyType::saturation)
+ .template value<double>(
+ vars, pos, t, std::numeric_limits<double>::quiet_NaN());
+
+ // Set eps_m_prev from potentially non-zero eps and sigma_sw from
+ // restart.
+ auto const C_el =
+ ip_data.computeElasticTangentStiffness(t, pos, dt, T, T);
+ auto& eps = ip_data.eps;
+ auto& sigma_sw = ip_data.sigma_sw;
+ ip_data.eps_m_prev.noalias() =
+ solid_phase.hasProperty(MPL::PropertyType::swelling_stress_rate)
+ ? eps + C_el.inverse() * sigma_sw
+ : eps;
+ }
+ updateConstitutiveVariables(local_x, Eigen::VectorXd::Zero(matrix_size), t,
+ 0);
+
for (unsigned ip = 0; ip < n_integration_points; ip++)
{
auto& ip_data = _ip_data[ip];
diff --git a/ProcessLib/TH2M/TH2MFEM.h b/ProcessLib/TH2M/TH2MFEM.h
index ce55ca76ee0015a739ad34adc42d24b0d71ab1f4..1cb48cb87c1cde6c730c3ae362a69e2cf3523422 100644
--- a/ProcessLib/TH2M/TH2MFEM.h
+++ b/ProcessLib/TH2M/TH2MFEM.h
@@ -90,7 +90,7 @@ private:
double const* values,
int const integration_order) override;
- void setInitialConditionsConcrete(std::vector<double> const& local_x,
+ void setInitialConditionsConcrete(std::vector<double> const& local_x_data,
double const t,
bool const /*use_monolithic_scheme*/,
int const /*process_id*/) override;
@@ -215,6 +215,29 @@ private:
_ip_data, &IpData::sigma_eff, cache);
}
+ std::vector<double> getSwellingStress() const override
+ {
+ {
+ constexpr int kelvin_vector_size =
+ MathLib::KelvinVector::kelvin_vector_dimensions(
+ DisplacementDim);
+
+ return transposeInPlace<kelvin_vector_size>(
+ [this](std::vector<double>& values)
+ { return getIntPtSwellingStress(0, {}, {}, values); });
+ }
+ }
+
+ std::vector<double> const& getIntPtSwellingStress(
+ const double /*t*/,
+ std::vector<GlobalVector*> const& /*x*/,
+ std::vector<NumLib::LocalToGlobalIndexMap const*> const& /*dof_table*/,
+ std::vector<double>& cache) const override
+ {
+ return ProcessLib::getIntegrationPointKelvinVectorData<DisplacementDim>(
+ _ip_data, &IpData::sigma_sw, cache);
+ }
+
std::vector<double> getEpsilon() const override
{
constexpr int kelvin_vector_size =
diff --git a/ProcessLib/TH2M/TH2MProcess.cpp b/ProcessLib/TH2M/TH2MProcess.cpp
index 2084a2c522cbd1efdc8d4f313a336eb61595a5b6..1ac1744d9f118c977046e8cc3f68d89b9664452f 100644
--- a/ProcessLib/TH2M/TH2MProcess.cpp
+++ b/ProcessLib/TH2M/TH2MProcess.cpp
@@ -56,6 +56,13 @@ TH2MProcess<DisplacementDim>::TH2MProcess(
integration_order, _local_assemblers,
&LocalAssemblerInterface::getSigma));
+ _integration_point_writer.emplace_back(
+ std::make_unique<IntegrationPointWriter>(
+ "swelling_stress_ip",
+ static_cast<int>(mesh.getDimension() == 2 ? 4 : 6) /*n components*/,
+ integration_order, _local_assemblers,
+ &LocalAssemblerInterface::getSwellingStress));
+
_integration_point_writer.emplace_back(
std::make_unique<IntegrationPointWriter>(
"saturation_ip", 1 /*n components*/, integration_order,
@@ -177,6 +184,10 @@ void TH2MProcess<DisplacementDim>::initializeConcreteProcess(
MathLib::KelvinVector::KelvinVectorType<
DisplacementDim>::RowsAtCompileTime,
&LocalAssemblerInterface::getIntPtSigma);
+ add_secondary_variable("swelling_stress",
+ MathLib::KelvinVector::KelvinVectorType<
+ DisplacementDim>::RowsAtCompileTime,
+ &LocalAssemblerInterface::getIntPtSwellingStress);
add_secondary_variable("epsilon",
MathLib::KelvinVector::KelvinVectorType<
DisplacementDim>::RowsAtCompileTime,
diff --git a/ProcessLib/TH2M/Tests.cmake b/ProcessLib/TH2M/Tests.cmake
index ecf9700f0044c06d3c511f0385ef5b4d938bac5d..a7bce2f19647e12d9e0e31c4d75da38d04ff01f7 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 0000000000000000000000000000000000000000..513444bbb8ee16d86c352d4bf005737039c25d20
--- /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 0000000000000000000000000000000000000000..d959c7c34c6d9f3fd6e4b7e2c0bb1c774630c77e
--- /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 0000000000000000000000000000000000000000..bd86778b4a8b8f731417ab4dd03cdd3263fa15f5
--- /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 0000000000000000000000000000000000000000..6eee1996226dc110284bedbb1e9f985f42ccc4b6
--- /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 0000000000000000000000000000000000000000..c52291a0bea78f9719b672d940c4ec2205b0915d
--- /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 0000000000000000000000000000000000000000..3f41b5c6c986e5541b3b6d5436688152e3690c1c
--- /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 0000000000000000000000000000000000000000..3015a4b0fe961c5ee424d9ca8df9cdecf0a697ea
--- /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 0000000000000000000000000000000000000000..acdc6e75f585490e988b463b119abfe551c35aa1
--- /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 0000000000000000000000000000000000000000..0b5d8bd31c51dc4b8fd3f2f595e8974b5a9a7b04
--- /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 0000000000000000000000000000000000000000..001b76f3793455aa73108e3f4fabd82a50b4b27e
--- /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>