diff --git a/ChemistryLib/ChemicalSolverInterface.h b/ChemistryLib/ChemicalSolverInterface.h
index 50a123cf0a36635edf3a0d4fc847d07914fd7910..8596a8721387f54828966eefa919d663f03d2b90 100644
--- a/ChemistryLib/ChemicalSolverInterface.h
+++ b/ChemistryLib/ChemicalSolverInterface.h
@@ -61,7 +61,7 @@ public:
virtual void updateVolumeFractionPostReaction(
GlobalIndexType const& /*chemical_system_id*/,
- MaterialPropertyLib::Medium const* /*medium*/,
+ MaterialPropertyLib::Medium const& /*medium*/,
ParameterLib::SpatialPosition const& /*pos*/, double const /*porosity*/,
double const /*t*/, double const /*dt*/)
{
@@ -69,7 +69,7 @@ public:
virtual void updatePorosityPostReaction(
GlobalIndexType const& /*chemical_system_id*/,
- MaterialPropertyLib::Medium const* /*medium*/,
+ MaterialPropertyLib::Medium const& /*medium*/,
double& /*porosity*/)
{
}
diff --git a/ChemistryLib/PhreeqcIO.cpp b/ChemistryLib/PhreeqcIO.cpp
index a93dda604377772d05b0a7f273d14452be3d30de..b1b5333a92695c6c6df1875da2b43a829f3dc5dd 100644
--- a/ChemistryLib/PhreeqcIO.cpp
+++ b/ChemistryLib/PhreeqcIO.cpp
@@ -165,6 +165,64 @@ void setReactantMolality(Reactant& reactant,
(*reactant.molality)[chemical_system_id];
}
+template <typename Reactant>
+void updateReactantVolumeFraction(Reactant& reactant,
+ GlobalIndexType const& chemical_system_id,
+ MaterialPropertyLib::Medium const& medium,
+ ParameterLib::SpatialPosition const& pos,
+ double const porosity, double const t,
+ double const dt)
+{
+ auto const& solid_phase = medium.phase("Solid");
+ auto const& liquid_phase = medium.phase("AqueousLiquid");
+
+ MaterialPropertyLib::VariableArray vars;
+
+ auto const liquid_density =
+ liquid_phase[MaterialPropertyLib::PropertyType::density]
+ .template value<double>(vars, pos, t, dt);
+
+ auto const& solid_constituent =
+ solid_phase.component(reactant.name);
+
+ if (solid_constituent.hasProperty(
+ MaterialPropertyLib::PropertyType::molality))
+ {
+ return;
+ }
+
+ auto const molar_volume =
+ solid_constituent[MaterialPropertyLib::PropertyType::molar_volume]
+ .template value<double>(vars, pos, t, dt);
+
+ (*reactant.volume_fraction)[chemical_system_id] +=
+ ((*reactant.molality)[chemical_system_id] -
+ (*reactant.molality_prev)[chemical_system_id]) *
+ liquid_density * porosity * molar_volume;
+}
+
+template <typename Reactant>
+void setPorosityPostReaction(Reactant& reactant,
+ GlobalIndexType const& chemical_system_id,
+ MaterialPropertyLib::Medium const& medium,
+ double& porosity)
+{
+ auto const& solid_phase = medium.phase("Solid");
+
+ auto const& solid_constituent =
+ solid_phase.component(reactant.name);
+
+ if (solid_constituent.hasProperty(
+ MaterialPropertyLib::PropertyType::molality))
+ {
+ return;
+ }
+
+ porosity -=
+ ((*reactant.volume_fraction)[chemical_system_id] -
+ (*reactant.volume_fraction_prev)[chemical_system_id]);
+}
+
template <typename Reactant>
static double averageReactantMolality(
Reactant const& reactant,
@@ -670,61 +728,36 @@ std::vector<std::string> const PhreeqcIO::getComponentList() const
void PhreeqcIO::updateVolumeFractionPostReaction(
GlobalIndexType const& chemical_system_id,
- MaterialPropertyLib::Medium const* medium,
+ MaterialPropertyLib::Medium const& medium,
ParameterLib::SpatialPosition const& pos, double const porosity,
double const t, double const dt)
{
- auto const& solid_phase = medium->phase("Solid");
- auto const& liquid_phase = medium->phase("AqueousLiquid");
-
- MaterialPropertyLib::VariableArray vars;
-
- auto const liquid_density =
- liquid_phase[MaterialPropertyLib::PropertyType::density]
- .template value<double>(vars, pos, t, dt);
-
for (auto& kinetic_reactant : _chemical_system->kinetic_reactants)
{
- auto const& solid_constituent =
- solid_phase.component(kinetic_reactant.name);
-
- if (solid_constituent.hasProperty(
- MaterialPropertyLib::PropertyType::molality))
- {
- continue;
- }
-
- auto const molar_volume =
- solid_constituent[MaterialPropertyLib::PropertyType::molar_volume]
- .template value<double>(vars, pos, t, dt);
+ updateReactantVolumeFraction(kinetic_reactant, chemical_system_id, medium,
+ pos, porosity, t, dt);
+ }
- (*kinetic_reactant.volume_fraction)[chemical_system_id] +=
- ((*kinetic_reactant.molality)[chemical_system_id] -
- (*kinetic_reactant.molality_prev)[chemical_system_id]) *
- liquid_density * porosity * molar_volume;
+ for (auto& equilibrium_reactant : _chemical_system->equilibrium_reactants)
+ {
+ updateReactantVolumeFraction(equilibrium_reactant, chemical_system_id, medium,
+ pos, porosity, t, dt);
}
}
void PhreeqcIO::updatePorosityPostReaction(
GlobalIndexType const& chemical_system_id,
- MaterialPropertyLib::Medium const* medium,
+ MaterialPropertyLib::Medium const& medium,
double& porosity)
{
- auto const& solid_phase = medium->phase("Solid");
for (auto& kinetic_reactant : _chemical_system->kinetic_reactants)
{
- auto const& solid_constituent =
- solid_phase.component(kinetic_reactant.name);
-
- if (solid_constituent.hasProperty(
- MaterialPropertyLib::PropertyType::molality))
- {
- continue;
- }
+ setPorosityPostReaction(kinetic_reactant, chemical_system_id, medium, porosity);
+ }
- porosity -=
- ((*kinetic_reactant.volume_fraction)[chemical_system_id] -
- (*kinetic_reactant.volume_fraction_prev)[chemical_system_id]);
+ for (auto& equilibrium_reactant : _chemical_system->equilibrium_reactants)
+ {
+ setPorosityPostReaction(equilibrium_reactant, chemical_system_id, medium, porosity);
}
}
diff --git a/ChemistryLib/PhreeqcIO.h b/ChemistryLib/PhreeqcIO.h
index 2981106c3717bd7ca06f3784b7ca11686c311dca..22d31ff3f7f776d0c94127c52b258abe39c1f8f4 100644
--- a/ChemistryLib/PhreeqcIO.h
+++ b/ChemistryLib/PhreeqcIO.h
@@ -74,12 +74,12 @@ public:
void updateVolumeFractionPostReaction(
GlobalIndexType const& chemical_system_id,
- MaterialPropertyLib::Medium const* medium,
+ MaterialPropertyLib::Medium const& medium,
ParameterLib::SpatialPosition const& pos, double const porosity,
double const t, double const dt) override;
void updatePorosityPostReaction(GlobalIndexType const& chemical_system_id,
- MaterialPropertyLib::Medium const* medium,
+ MaterialPropertyLib::Medium const& medium,
double& porosity) override;
void computeSecondaryVariable(
diff --git a/ProcessLib/ComponentTransport/ComponentTransportFEM.h b/ProcessLib/ComponentTransport/ComponentTransportFEM.h
index f5740eb98a4d4dd5f74dfa68c9d0b8ec2cb98649..e077617de0d88760e7eb584314fe32105437f7f5 100644
--- a/ProcessLib/ComponentTransport/ComponentTransportFEM.h
+++ b/ProcessLib/ComponentTransport/ComponentTransportFEM.h
@@ -1164,7 +1164,8 @@ public:
{
if (_process_data.chemically_induced_porosity_change)
{
- auto const& medium = _process_data.media_map->getMedium(ele_id);
+ auto const& medium =
+ *_process_data.media_map->getMedium(ele_id);
ParameterLib::SpatialPosition pos;
pos.setElementID(ele_id);
diff --git a/ProcessLib/ComponentTransport/Tests.cmake b/ProcessLib/ComponentTransport/Tests.cmake
index 9b364a8e56916afa778d40d9a009910bff04241c..b9f3c2aa52e3a53f8251cbeb00cb6495cc9961dc 100644
--- a/ProcessLib/ComponentTransport/Tests.cmake
+++ b/ProcessLib/ComponentTransport/Tests.cmake
@@ -1026,6 +1026,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/ReactiveTransport/KineticReactant_AllAsComponents/KineticReactant2.prj RUNTIME 60)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/ReactiveTransport/SurfaceComplexation/RadionuclideSorption.prj RUNTIME 60)
OgsTest(PROJECTFILE Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calciteDissolvePrecipitateOnly.prj RUNTIME 25)
+ OgsTest(PROJECTFILE Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj RUNTIME 25)
endif()
AddTest(
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj
new file mode 100644
index 0000000000000000000000000000000000000000..72457ca07acb09c149007615cc1917c8cc4b551c
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj
@@ -0,0 +1,683 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>calcite.vtu</mesh>
+ <mesh>calcite_upstream.vtu</mesh>
+ <mesh>calcite_downstream.vtu</mesh>
+ <mesh>calcite_ReactiveDomain.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>hc</name>
+ <type>ComponentTransport</type>
+ <integration_order>2</integration_order>
+ <coupling_scheme>staggered</coupling_scheme>
+ <process_variables>
+ <concentration>C(4)</concentration>
+ <concentration>Ca</concentration>
+ <concentration>Mg</concentration>
+ <concentration>Cl</concentration>
+ <concentration>H</concentration>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <specific_body_force>0 0</specific_body_force>
+ <chemically_induced_porosity_change>true</chemically_induced_porosity_change>
+ <secondary_variables>
+ <secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <components>
+ <component>
+ <name>C(4)</name>
+ <properties>
+ <property>
+ <name>pore_diffusion</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>retardation_factor</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>decay_rate</name>
+ <type>Parameter</type>
+ <parameter_name>decay</parameter_name>
+ </property>
+ </properties>
+ </component>
+ <component>
+ <name>Ca</name>
+ <properties>
+ <property>
+ <name>pore_diffusion</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>retardation_factor</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>decay_rate</name>
+ <type>Parameter</type>
+ <parameter_name>decay</parameter_name>
+ </property>
+ </properties>
+ </component>
+ <component>
+ <name>Mg</name>
+ <properties>
+ <property>
+ <name>pore_diffusion</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>retardation_factor</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>decay_rate</name>
+ <type>Parameter</type>
+ <parameter_name>decay</parameter_name>
+ </property>
+ </properties>
+ </component>
+ <component>
+ <name>Cl</name>
+ <properties>
+ <property>
+ <name>pore_diffusion</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>retardation_factor</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>decay_rate</name>
+ <type>Parameter</type>
+ <parameter_name>decay</parameter_name>
+ </property>
+ </properties>
+ </component>
+ <component>
+ <name>H</name>
+ <properties>
+ <property>
+ <name>pore_diffusion</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>retardation_factor</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>decay_rate</name>
+ <type>Parameter</type>
+ <parameter_name>decay</parameter_name>
+ </property>
+ </properties>
+ </component>
+ </components>
+ <properties>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1e3</value>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <components>
+ <component>
+ <name>Calcite</name>
+ <properties>
+ <property>
+ <name>volume_fraction</name>
+ <type>Constant</type>
+ <value>2.446e-6</value>
+ </property>
+ <property>
+ <name>molar_volume</name>
+ <type>Constant</type>
+ <value>3.6934e-5</value>
+ </property>
+ </properties>
+ </component>
+ <component>
+ <name>Dolomite(dis)</name>
+ <properties>
+ <property>
+ <name>volume_fraction</name>
+ <type>Constant</type>
+ <value>0</value>
+ </property>
+ <property>
+ <name>molar_volume</name>
+ <type>Constant</type>
+ <value>6.439e-5</value>
+ </property>
+ </properties>
+ </component>
+ </components>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>permeability</name>
+ <type>Parameter</type>
+ <parameter_name>kappa</parameter_name>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Parameter</type>
+ <parameter_name>porosity</parameter_name>
+ </property>
+ <property>
+ <name>longitudinal_dispersivity</name>
+ <type>Constant</type>
+ <value>0.0067</value>
+ </property>
+ <property>
+ <name>transversal_dispersivity</name>
+ <type>Constant</type>
+ <value>0.1</value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <global_process_coupling>
+ <max_iter>6</max_iter>
+ <convergence_criteria>
+ <!-- convergence criterion for the first process (p) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process (C(4)) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process (Ca) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process (Mg) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process (Cl) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <!-- convergence criterion for the second process (H) -->
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ </convergence_criteria>
+ </global_process_coupling>
+ <processes>
+ <!-- convergence criterion for hydraulic equation -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ <!-- convergence criterion for component transport equation (C(4)) -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ <!-- convergence criterion for component transport equation (Ca) -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ <!-- convergence criterion for component transport equation (Mg) -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ <!-- convergence criterion for component transport equation (Cl) -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ <!-- convergence criterion for component transport equation (H) -->
+ <process ref="hc">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <reltol>1e-14</reltol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>21000</t_end>
+ <timesteps>
+ <pair>
+ <repeat>210</repeat>
+ <delta_t>100</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>calcitePorosityChange</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <timesteps>
+ <pair>
+ <repeat>5</repeat>
+ <each_steps>42</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>C(4)</variable>
+ <variable>Ca</variable>
+ <variable>Mg</variable>
+ <variable>Cl</variable>
+ <variable>H</variable>
+ <variable>pressure</variable>
+ <variable>darcy_velocity</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <chemical_system chemical_solver="Phreeqc">
+ <mesh>calcite_ReactiveDomain</mesh>
+ <database>PSINA_12_07_110615_DAV_s.dat</database>
+ <solution>
+ <temperature>25</temperature>
+ <pressure>1</pressure>
+ <pe>4</pe>
+ <components>
+ <component>C(4)</component>
+ <component>Ca</component>
+ <component>Mg</component>
+ <component>Cl</component>
+ </components>
+ <charge_balance>pH</charge_balance>
+ </solution>
+ <equilibrium_reactants>
+ <phase_component>
+ <name>Calcite</name>
+ <saturation_index>0.0</saturation_index>
+ </phase_component>
+ <phase_component>
+ <name>Dolomite(dis)</name>
+ <saturation_index>0.0</saturation_index>
+ </phase_component>
+ </equilibrium_reactants>
+ <knobs>
+ <max_iter>100</max_iter>
+ <relative_convergence_tolerance>1e-12</relative_convergence_tolerance>
+ <tolerance>1e-15</tolerance>
+ <step_size>100</step_size>
+ <scaling>0</scaling>
+ </knobs>
+ </chemical_system>
+ <parameters>
+ <parameter>
+ <name>kappa</name>
+ <type>Constant</type>
+ <values>1.157e-12</values>
+ </parameter>
+ <parameter>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.32</value>
+ </parameter>
+ <parameter>
+ <name>decay</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>p_upstream</name>
+ <type>Constant</type>
+ <value>1e5</value>
+ </parameter>
+ <parameter>
+ <name>p_downstream_Neumann</name>
+ <type>Constant</type>
+ <value>-2.9976852e-3</value>
+ </parameter>
+ <parameter>
+ <name>c0_C(4)</name>
+ <type>Constant</type>
+ <value>1.23e-4</value>
+ </parameter>
+ <parameter>
+ <name>c0_Ca</name>
+ <type>Constant</type>
+ <value>1.23e-4</value>
+ </parameter>
+ <parameter>
+ <name>c0_Mg</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </parameter>
+ <parameter>
+ <name>c0_Cl</name>
+ <type>Constant</type>
+ <value>1e-12</value>
+ </parameter>
+ <parameter>
+ <name>c0_H</name>
+ <type>Constant</type>
+ <!--pH=9.91-->
+ <value>1.2302687708123812e-10</value>
+ </parameter>
+ <parameter>
+ <name>c_C(4)</name>
+ <type>Constant</type>
+ <value>1e-10</value>
+ </parameter>
+ <parameter>
+ <name>c_Ca</name>
+ <type>Constant</type>
+ <value>1e-10</value>
+ </parameter>
+ <parameter>
+ <name>c_Mg</name>
+ <type>Constant</type>
+ <value>1e-3</value>
+ </parameter>
+ <parameter>
+ <name>c_Cl</name>
+ <type>Constant</type>
+ <value>2e-3</value>
+ </parameter>
+ <parameter>
+ <name>c_H</name>
+ <type>Constant</type>
+ <!--pH=7-->
+ <value>1e-7</value>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>p_upstream</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>calcite_downstream</mesh>
+ <type>Neumann</type>
+ <parameter>p_downstream_Neumann</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>C(4)</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0_C(4)</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_C(4)</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>Ca</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0_Ca</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_Ca</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>Mg</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0_Mg</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_Mg</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>Cl</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0_Cl</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_Cl</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>H</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>c0_H</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>calcite_upstream</mesh>
+ <type>Dirichlet</type>
+ <parameter>c_H</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>10</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i cg -p jacobi -tol 1e-16 -maxiter 20000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-14</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>hc</prefix>
+ <parameters>-hc_ksp_type bcgs -hc_pc_type bjacobi -hc_ksp_rtol 1e-8 -hc_ksp_max_it 20000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>pressure</field>
+ <absolute_tolerance>1e-6</absolute_tolerance>
+ <relative_tolerance>1e-10</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>Ca</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>Cl</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>H</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>Mg</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>Calcite_avg</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <regex>calcitePorosityChange_ts_[0-9]*_t_[0-9]*.000000.vtu</regex>
+ <field>Dolomite(dis)_avg</field>
+ <absolute_tolerance>1e-10</absolute_tolerance>
+ <relative_tolerance>1e-16</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9045e87a84f2c0bfb28513527691ba47919e532e
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_0_t_0.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0.00020719917662" RangeMax="0.00020719917662" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="0.00020695700439" RangeMax="0.00020695700439" offset="84" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0" offset="168" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0" offset="240" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="312" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="7.9018480282" RangeMax="7.9018554854" offset="392" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="2.4488622045e-06" RangeMax="2.4488622045e-06" offset="484" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="2.446e-06" RangeMax="2.446e-06" offset="568" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0" offset="652" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0" offset="724" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="0.00012275782751" RangeMax="0.00012275782751" offset="796" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="0.00012275782751" RangeMax="0.00012275782751" offset="872" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="9.9999999786e-13" RangeMax="9.9999999786e-13" offset="948" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.230841961e-10" RangeMax="1.230841961e-10" offset="1024" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="9.9999999786e-13" RangeMax="9.9999999786e-13" offset="1104" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="0" RangeMax="0" offset="1180" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1" RangeMax="1" offset="1248" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0.00020719917662" RangeMax="0.00020719917662" offset="1320" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0" offset="1396" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="1464" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="0" RangeMax="0" offset="1528" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="1596" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="2524" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="2824" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="3136" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAHQAAAAAAAAA=eF4T2fwyvlhD215klB6lR+lRepQepUmgAe1SZ78=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAHQAAAAAAAAA=eF7zyu1WDFbQtvcapUfpUXqUHqVHaRJoAIKYwhA=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMApIBwAGQAABAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMApIBwAGQAABAQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAIgAAAAAAAAA=eF4TYZdjbJgt75B0dA1P3Sg9So/So/QoPUoTSQMAhOOadA==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAHgAAAAAAAAA=eF4rOf7P+knXEbuSUXqUHqVH6VF6lCaBBgAmDNJOAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAHgAAAAAAAAA=eF47kdSyY1/LEbsTo/QoPUqP0qP0KE0CDQBFOqT+AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMApIBwAGQAABAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMApIBwAGQAABAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGAAAAAAAAAA=eF5bNatmq6C4gv2qUXqUHqVx0gDMdi18AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGAAAAAAAAAA=eF5bNatmq6C4gv2qUXqUHqVx0gDMdi18AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGAAAAAAAAAA=eF6bYtLTMHN6oe2UUXqUHqVx0gAQZnVKAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGwAAAAAAAAA=eF6bs1j+4eJXD2znQuk5o/QoPUpj0ACpQO/+AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGAAAAAAAAAA=eF6bYtLTMHN6oe2UUXqUHqVx0gAQZnVKAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAEAAAAAAAAAA=eF5jYBgFo2AUEAIAAygAAQ==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAFQAAAAAAAAA=eF5jYACBD/YMo/QoPUrjpAGtiHeMAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAGAAAAAAAAAA=eF4T2fwyvlhD215klB6lR2kMGgBlejPgAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAEAAAAAAAAAA=eF5jYBgFo2AU4AIAAyAAAQ==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAEAAAAAAAAAA=eF5jYBgFo2AU4AIAAyAAAQ==AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_126_t_12600.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_126_t_12600.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9c099a2745c0c64d0f4b749aeb2df06e28692fd8
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_126_t_12600.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.00020774045715" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-6.3649179011e-37" RangeMax="0.00020775286603" offset="868" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.00011424731226" offset="1824" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-6.2282857813e-120" RangeMax="0.00011418268655" offset="2128" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="2472" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="-5.5719283195" RangeMax="9.8236790107" offset="2552" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="-2.4163421816e-120" RangeMax="2.4552594854e-06" offset="4568" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-7.5226247636e-39" RangeMax="2.455406143e-06" offset="5592" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="-4.4094349429e-138" RangeMax="2.354043727e-06" offset="6612" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-1.2833356496e-121" RangeMax="2.3527121267e-06" offset="7008" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="9.900639174e-11" RangeMax="0.00014796621409" offset="7408" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="1.0000000001e-10" RangeMax="0.00022972306584" offset="8312" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="1.0000001482e-12" RangeMax="0.002" offset="9252" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.230841961e-10" RangeMax="1.0000614881e-07" offset="10388" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="1.0000000704e-12" RangeMax="0.0010238930009" offset="11324" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="12460" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="98704.543993" RangeMax="100000" offset="12964" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0" RangeMax="0.00020773251671" offset="14040" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0.00011404858608" offset="14520" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="14712" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="14776" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="15268" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="16196" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="16496" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="16808" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAaQIAAAAAAAA=eF7tz91L01EYwHG9KNSBxlqZLMWXSreTOqdWWu0IgS+ZzmnDtkVW1G9lpoTgbGL2Yqll+TKN1LXNZRQGpcOMmNOiYKm0gqALc6Zr5iAX5jDJi+o5Z/9A3XtuPnwPD4fz+Pisnv85xRK/948bt+DZ2oT1o4UB2OXcH70Sg/CQVbNpao6HJ0zmhuRUhMt1m0dCdyD8fXD5Kgea0TMLWtBP/pOl2onwL0HfKeMuhFns8iFzEsJRkT0XnMkIz2RKasVChNn596aliQibq3WOXAHC+eMVizEJCOMIpvp1HMzrTmwojUf4MlPK9McizFW42Hy4V7fMXgmEFilt4TKwss7jfrsd4YymGe4CaMh9NzcJXv/wau8UOBA3lbIN5gbUI+IgUKW/o2kDlePCgFLQFCY1cODdlMairK/Qv3nZt25DBwnv2hhw4PjKl0D4x8cnJ484oJsPBBkqoH3nc/l54EF72vAY2Ds/0fYITMzwLWDBPmFL8SoP9AOroCoVOrZbGcED+0NFMgWYlvcmT0z2Ngsyz4EXF13Zp8E1a39Ia8BD7lGsBktKonzrwZpvY45asHM5mH8T5HZZJTfAyZC6yGawWJ7OaQIZ5/OxFtBitycRy4Rn+lrBSvdwEfGpPJ+nAfsehvwh7c/ZpyS91fRCQvzcGXeN+HLjpSqi7ejuduJ97rpmIn/+mJZoPbvQQWRZTHpiB2rVUQ2J3UR0vt9A1NvstE01PdT63mgjcdoTSHW70qm6rj1UzacC6qBRTC3zl1GDZYXU2SwFdUkrp5a4vC2weM3JOUx1iryGF3t9pvDaULnqqv/uX4Au/do=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAArAIAAAAAAAA=eF7tz2tIU2EYwPHlskW5TAOdSEPDS3TyVpo35LUUS9KNjLwOMW/kBZUtM6zRzDDTeZukaTl1aolarouiQcIsc4TTSAvniImhZoUOWlppruc9ww8OBtFnny8//s95Oe85FAo5AspmN8awje0FBhrbGzu3MYZ7wzY2gv993hLzmeczIzTd6LC6PC1F6tLU9/Ld/gGKoFr5NctZbpnRtyR2l/kfVHSUO+Raan84IUehiSTKhoGCi56ekXoQqD50rX9KdQh9r1fcpPkT6FldqnSnL4GOnzPzewhqFKP7SvwIFFLo3Td3jEADurRFWx8C2Xcl+3l4EogVm/cm0otAlaKSwPfwvlGO/5+JIwTKoLJZC24EWh07vf7cnUD5OcO6XFcC0XlzPlawTx7sr813IZDFq8CWBdDbZMl38TCB5P3KTjdohXb+qgy6uG1cpAI144Nl06B59D02bo9MWbMXnONGzM/YgNbTlLhOUKU+JSsAK9STI65wX7hLtm4d+r7myicp9K49tN8C8EPrNQtn+A4mLZa9Cm3JXwkpheZrR+xSwdnta2dVoNj89uMXoFxitmgN/xHdvrtqG9j6Zbg2CAyzo5cdBbnFk5Lz4CfW2INIsCyTnpMLLgeoA7LAeEferwJw4Zs9lw/y37LEt8DOHUPqIpAlYDLLwXSKQ6EQvIgojypB67aUoQrcPVSGCOwOSrlTBTJ6TKZwXwg1rcGyo6xQNVjM/WiLZShjuNjspJoEbLKuuxTL+6m9jnWSz9Zi4+t6q7HKhNYGbDu1+y55ruhkEzb8RmIjtqHrUjN2MnEvqZizSip9PUpamOQpwaZx7EmpJ8JIx7KCSRvpkaSiuQhSmiaGdOVJNGl8UxxpB1MvMuWQTizruzdO36xwvcLLeg+kb+4tt/wX/wKaoB9AAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAwQAAAAAAAAA=eF5jYBjcYBnjZSG9vwftYPw1xn/7SpZJ2u91mmv0YJ+g/YaOd29W9sratzcfdWX7IGMvfOJ6Gu9DWfu//0J0vV7J2h9nv/3w8TlZe55CsYeMn2Xtwyt5T7sslLX38Dp7adkpWftA8ThZxwRZ+5apaVdvdsnal/+17rPglbV3Xnbjjr+NrL2orqVV4loZe916V51td2TsTb55TL0+77ed9VyVm7wLRO1h7mq3XXqWp+cK3J2jYBSMglEw3AEAW/FF7w==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA4AAAAAAAAAA=eF5jYEABDQzYAbHi6HxcoAGXeEnbO4fdgoWTYAK8M3kKGfne2vmxNd6wEzVY/iVNz1BZTMp+8uTujZ8VhO0vfN70KHyyrP1D3coSOSZZ+zmFi86135O1d+5/MEv6law9Y4nm5UUnZe2nLTWct/qdrP2BFWIOPLNk7UWvi08UPSJrX1G66LRhuKz9x82zvqc2ydr7LNjlfu6/jP0B7i6vPway9qZC2WkSC2Tst29ju9JwTsb+X//2c2/LP9p9kxHdpLeP3x7mTtH0vyZuSzfZYffWKBgFo2AUDD8AALxdT0o=AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAxwUAAAAAAAA=eF411Hk4VXkYB/AuhRoG95zfsUQdabVE2UL4FXEJGdkKLYYKdVNRmKi42SIlLbYpxSW50R6lXjUhuyaRYgoZuiRT2UrNuXee+evzfN/z/t73d57zPKfqU9Zw/BMFGBD4j5wJlIcRtZ4O02A2LA4xKjOPV4QcN1npibsEbJTd6T59ig0vLOdGjkWyIeavAOtygSI066sUy7si0Lq1cVVZvibmthgl6E6bh2tWODo/9GIDcbmFMolhQ6yd76OiNASCK8Ytf2QrQswjXlrAFhJcJ1e0aKykQE/R5fH6cKYvu/XCTDsSuicproMuG7qUJco8U0nIepZ8oKCCxnnKxwNTd9J4bLdDmsppNpQFG713FbIhjZhx1ruWxtK0Qs2sVzTuqSSOrlbVwNxDvEUBAzTWMDgrF/yExg/7hCtt3tA49UG/f/ggjW95zo6OSFQEbesbT2vbaNzkWtoQJqWBc/pmPLbtoHFCgozehzM03l2xS9DBzB96+ZDknaex0Ke1rekqCdmU56YNDTTe59RiUWRC48DtleqHztL4U9M1GfOPNFbWqWHt3UlAyQZCYzSOmb9eExdcIuBy9khi8xiNy/U4rzdH0tielP5Tm09jMi25yrqOACu1S4M6CTTOV7z982tm/+9vU/c3hRHQzqVDnaZo3LvACNvU01j7jUf4AhMS9JOyfPKiSCjYa1yrxSHgeP+7+bltBLR124/eoUkYE6oMCVoIMBRGc/TzCSi4JnO504CE/kJ1cxslErw9/jEuX0zCDc3V5FAEARNZp4ZLbpNw4ebHWgkvEmQLtzo07iRhyYe4OO2/SRgM6hrmVRDgMLBm7HMWCbuWkjP09pDw9jl/IQAJvpY+fSp8Eh6kmzn3TkOwW/uQ3KwGEraH+K12qidhNt8uVe0lY/TeZpnvJFTWc5Z+ySDhTd5hZZ48gpAwCo0KSPCsy0g5MklC8I4RmTrm3Cer125z2Qi8Hv9S7kogIKtzX4EGgoD31CZ5Zt4Ri0Xm5fMRSKnE5cYoIFhYzBrm2CII6zi2NWATAn7LuOoGTwTyTd3bDOwQ1FHxvX4uCNK6ii/2OiFYly6U0FqP4E5l+dXlvgiwkUJsymEEalKylpxQBO1Szt5oG4LOVn3X5DAE570LJWODmP1eawZjYhA0FJSYvCxBkDHLWih5gtl/sLOLz0OwbMmO730RCBLJb46lFxD4+rl3nElBcDO1/obSVQSmapy4bj6CqH2SJzqKmDy6QrCFyXwXtykfAQK7+yydJZnMe62ubFSvRsA7dZQz9zKC4uqLQqO7CCS/Puf8eg+BiaYKK6MBwbTirqrx+wheSKbm1FQhcL9xMyYQEKhMhC/SakFQP9T0Ra4ZgZVnSqxGL4KR+5/XXn+GoL+s++j2QQQna97HDvYgMNxyKyhlmPkeVrvSJrsRaEZmm8z5isBWKNsT3M/U++V10yYQuClwDXk/EEjkHD56fRzB85/dgr9NIjBQ28JdO50Co+Z01R/MnJi+dXvmyFGwvPu3zAoZCpQqQvJIRAHNmltlz+STplYJmaoUGHpY3nZk/HpFU12XcRfra4mxMgW2ul3xL2gKsF7QgemaFFBt6QYRWhRoS+qc5jCZ7W+eVLqMeZ7u3XVAlwJd5XL1KGMKWtsHHHkGFHg09lWZmVLgdq/If50JBXO4d8blmf+P+k8DErlmFChmVb4LwRQ8jbSIqrSkIMksYCbLloLNB0f0V1oz91Lx+XKPQ8FSqwa6w465f4FJfY49BRbtjacKGOXqlX+cZIyML5Q9xzjab1OXzNh5p3/yOOOmTHdClH25mjNFJutHLxHJqlmhJJIfKnBJYgyOvnj9GGNIxBrvRMZL8aqlIj9MalmI9Di8Q2xRoc4akdz9p6xEuqwdXCWujypaijw/T8lGZLVnJBafa/AX90k4J1uLzJ9SE5/3Mw4V1/NGKXE/nviv/n8+lyIttqLGSGzp5oVi/wUR9tY6AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA3gIAAAAAAAA=eF7t0OtLU3EYB/CJ5S3SmpfaUiQFE28TETWtE6RYGpIXCDe0vIw5pRBnambGDC3ZROd9ZipTM90QTDOct4meCeEbMY28QJQRzHKG2zSv9PzO+fmid/0B7s1n3+95nsPDCWFQP/EaVkrLKML5Ac4q7OG/84xknE9iOWb0n2ZacSzu/XDexnkSz2XQMAy4j8bZiOcVeO4F7q9ji3AfhecGcX90lwv2Bn5+dP8c1hnvzycWVuZrZ0o8cS+8syM2FSc5+18cmCASFHZZXxPX6uY0geIYVhnXdP4Nm71pm+ukIXTV7G7b9Q7iMNLNO1pKEt93apMF8inCkZu5FS8niR5upDC0niRK1b8PLkM+ZIb5KsC0uoCmaw0k0XrAt3ZrJIkrH3gdfnUkkcnJuGAH89pbSjG3hiSCmjYMrrUkYV7splmQkcRVVU7f42qScP4xcV9VRRK17NF33tAru5oTlythz0i2BUEv6el3T4Fsb9NXKwf1/dwoDpiwMO8eBU7fimgJBXt79G1h4MruNlMEsjY/jqE9vxZXqRZcl63xu8H9uQqBB7xXFZHKNQcFXnu95eBAvamKD7rHFg2sgp7ceItJ0NVY4xUAdy1LQuwdwUsc1kwOOL7wJ4QHJq+XdraDk6kVchm6f8R4exp86OTrrwaHZu/6fAEjT8cwP4Hto+VTG6BEH6bWgU7Cb/xd0KxBL90C00zFBgZ8l3NrrY2HkKeHlEsnIPOChXnm4OKrFZYleNYmJ8YCtHcY5lmBEYJ7OtS7NBRYWYON7qIZ1Ht7/ExHmZm1F450DtPWIM9wxM+Q4TJPJdL+l0M7smtWMogsWVW8RZ4SGYaRS3FsNTL4ad4YtR/SPIpcbJOPI3U+QZRJXJYGGcjQUzklOI7KK2UEZfqikLLUNpWyqVNEaZWbTfmkKp8y/2Ye5fPwR5TvJwoos0foPPqSVmNZSPl5n85SFzq/ZtKKPI899v/9CwrJgbk=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA3AIAAAAAAAA=eF7t0P8v1HEcB/Bjd6Sh7W5nEtaUuR+KW9Nw8vFtDS2GNnRd33w5EitfwnXJVZrma7lSHN2wvp1maczN17QPlroWmghNXb7F7ZavEdbr/fm8/eC3/gCfXx6f5/Pz+rw+731cGdQl+4XNo2VIcU7EuQa7sXWecRpnFtbJgL5R0MqCce+I8zLOb/FcLA1jHvfHcF7A85V4Lgf3Plgp7gPwXAPuN89lg/XDzzfP34e1xu/PSlfjzYbsWTzcn1NJchpSu5SuUm6p8rClXKLqtrpxlamWVTfNGZqLVdlevvvsznQQ+Waj+/MErwnRTAn3+z2S+BAaumchmCRyA9I0vx+RhN3gY/GlhyRRSuQOz4C23V+miqGvra1I7HtAEpN7OSEtJSSxxC751ikniaJ8zUz1fZIYt/JmrcO+22/0V9qKSSKR7TN95y5JcIr12Q7Qj6g7zS4XkYR2ws1zGpRHq7R1hSSxrAr5wYKs4Xs3+kFW1iy6SECXEYcxe/CCqVDhBp6wbjzrATaLb8ah3JnBa8oE6zZyW+NAixjTrgHQiykPagTLjMYOCWCvgJO+0wI0GWWXloPcopRTGWByeMD6Ivpu4IL2M9gQX/XOG85rqLN9ZQ/29ZSFXQerbN1nY8EemfTgS/CoO8dXAVqSEc4a0KvIgN2BFGh1WtBxvFw5jOYKPKfmwKGs3jQdGM4fjFoDu3brl/6AkSHMA4bwX/jXKiwZ4Mq8p5wFCtMjRUwwwl90yxg04l7sNwLNj/Af7ABjBM+dkU9TdA4mYOVc1gLKysHzqSgvJrWeRK5l8RRIrx51AVJTHVqL9NulfoZ8sV7YiBxui65HFgonm5ETtglNSJHJRCsy6pMr5dew923IseVMyvB6p3bkx1Uzyr8DQsq05OOUqz8TKNtXxJQtgamU/jNJlFMe6ZTcoDRKfXYGZb8drUcvbVgX7RO2hN5nTJvKoy2w2Zq33fZ//AdRdYmuAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAACAEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qPiDk+G/cizO8klD9c2tuXb+zY50u4ashhRnx6zFX2Rf889/ecAu9uDKuXvTttttc1svXlxz2C5JpWlOk8hhuzI5t668zYft7v+5sb9152G7cj+LxyXLD9v9ODjH0nX/YTs+R1fvhM7DdsdWvitZueSwndNBH9+mOKA+p8TAGVWH7Zh6Xp8/pH7YTtGqvfy822G7tZcflE58ccju8PLPh46xH7br1cro1jdbZOflpvogv2WPHdT5DPOe87uVS1TD+aNgFIyCUTDcAQAOa2VaAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAACQEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qPgB3+QDBU1ik6Sh+syZBMT9NSbbNZz2nlVtenWZ7zPbHCfjg3bLssSMV1nutLuW73LEsOGwXYjoh9sFCoftUr5l3j2y8bDd4oDdXOk7D9v9f3yJ5d/iw3bSp9q1M/ccttM/YF88t+Ww3Z7Q1y9DFhy2u6+4l6cz4rBd0OSCCaylh+0WvNqT8UrusN2Gpsb54faH7cSi62TO3z9kt+X6nrKMf4fsNHfnJ2hxTrcLX+07gf3dFjuo8xmWWxhI906LgPNHwSgYBaNguAMAAIBhDA==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAhQIAAAAAAAA=eF67Jtc644vZbdtVs2q2Coor2LOsW2TSvX2SXaz94usNQavtwtaqnbl0c6/dzb7beWIfz9hp3fEWuh1wy66rwnynL+9zu4REFRH39Z/t+E5ecFPxYbB/qKhSqLGf076hTdD1br+wPceES5JvUpTtZVQFXJRvKtpXzy+vMylTtG/02r/4HYOi/U3ORiGragV7B6NHPXevytsXvJlUm7BRzj7A/vvmW+2y9trR5dbbfGXs59WtM5rdIWNvlj/D/4+ErH1F4Prni1Jk7WelFZuy7Je192Pb524gLmd/8DPT4eJsOXvl4zn9h3bJ2cccSOP5zCxvP2nb6aefXeXtz6k/mfWrXt6+x+5H97SN8vZiaR6ZW2/J2x/4Vfh33m95+yo53lclXAr2q06eZNwqrGC/chOrJJ+Ugn2B/eft7HIK9gYrnO00FBTs7x7+5v0GSBupey06AKS/e2/f809ewb5L1/PLa6A6Y7u81btkFexv7RC60iOjYK8b2SZfLq1g/0nmsGsL0Ly0ay//b5dUsI/oXdsgCKRfPHkU1SuhYM8x68c9NSB95epFlpvA+Ag43v17HpB2mDGtvARIP1C6eTkUSE9dpCrmAKQ/qTI+NgDSvCLSmWpA2qnvs6c8kJ6wmL1ZCkh7u4RGigPpXV/+uYsCaQ1Xs0/CQPreFkMJEN259GifEJCOna5hCaJDqvJYQLTbkqAnoPTwcNKfCyBa/W35URC9x2fdPhB9z2LnLhCdsblgB4jWuaG0HUS/UE3bBqL/9fwBp6fK3wfB9MerC8H0h7gJYLpEqB1MKwk1genfOfVgOoC3Dkw3x9SCabdfkHRZcQpCu2+A0OwzIDQs3dKaBgBRmTZtAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAoAIAAAAAAAA=eF6l0v8v1AEcx3HNuGVD92XO5XZ97ug64jJ9k7rPu3FU3ETlNFpXHGbdGvLlmsuK61bSTiM5NA1dO80py4Wc0iSWXNEhGXO+3DorG5HN1vp8WX9BPz1+eL1/ez/jUme/zqi+iZprVO1UJgL2/R4ipesdNPiAOTO9vhE1rry9NSJuRyXaT5kj6X1or8vDJpZkBJVz7pXN3phBC87eL6yQOdD+he2tm4o19Hzy1MCJF04Q3h2uMwu2wsTx1wdzM5ig+lA/3FbEhNwv40XKaBaMVryb55h8wMimt3WtcqBBdihPUsuFX/rIcXGiHygzRWNSdwFER39MaF4JgpQcDjuMLYSwxpNLNGMgBC14WFhNuyHHwY2zawIga1/Ik8osf5D9GFNZ5QKIcNfUHknfBcJj0lj/Aj7kH97SUaLbCa6Pkparhv1Aywh2pjL9YHCVvky94gtpnwOUYOOBzH5zaDKVB9bO0OSxdS7E6pa7HTounNO62QwxXDDQW9jlblxwerCxSLMisED/o2gxIGCz2hwODQLe4Wz75UsIlN39HbWZgIDpDQ3NjcQclQbaQhHQU13qsvdgu9Lcmi1AoEfKS4v3RWCohc6s24HA0Qn3xjU2AuquULWXDwLTj+O5AywEGFLfvFlvBL77aPlizNKgcq9J7F/F789cqMB02mD2yjEtmm1JUZhiS4c5BNMWExDGx2wH/WkOZl/DXjkLs1fhn+iFyfA4dZWBWTa/nkjHrJlWN9PwO1F1BG6xM0LBTchYmcP7mFIbLbjPujh9uBJBlRl3mvKzE3fRsfgSN//VoAmXpxYSCvsnib5KKUZCRXwVoXjpNqExpYRQ5nmdkH+xiNwLrxHOeZLaJshOn/eQZulJ6QZS01NSSjXpv67/17+C5jkwAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/Iek8dJNYmA/Nx/qrJmXcT3UmxWrk2FgPxyUGw86YGA/VJWSletdYD+G2KbOM1pgP35niDJ5VGA/UwETvfdLYD96C1MGvT9gP/xN7cKmLmA/NbSlnGQXYD/yho98/PBfP45Udfi8oF8/8IVIlsI6Xz+WQOSQ17teP5ApDpIDIV4/XmDYU8NnXT9PSB8PP45cP2iLmmd9k1s/srGZyI13Wj9F2G+epTtZP0PKOcUs4lc/zKRITbZuVj/gG1HU5eVUPzyR/OJBTVM/prN2ffaqUT9XbUuijAVQPxpSF104x0w/+FvXbPyWST84gwJIC4ZGPz4lxkGGnkM/4vSIErDoQD9hx7tXatU8P+B4/QcvUTg/jwO1jndIND/n7tuYJrswP9tTLFVxSys/j4EH0sQDJj+rPd4ve44hP0L6vq/brxs/+5taZcSVFT/25EfaUaMQPwYA4YctXAk/Rm1yuggcAz9qsPxrqXn8PmTbgMFI+vQ+2iL+/RyQ7j56JK9THATmPrILYZ/rXd8+yl+QZRcZ1j6QZIAcVsvOPrD8EADbOMU+TLo17+DuvD6slO6HmIKzPrzfBGGsB6o+ukkTZ+ktoT6+Kn1QsG+WPn0wOE5+/ow+9CK55v2Jgj5a3ktib3Z3PiMX8JmTY20+lL2js/c3Yj4wnio1GVxWPm7Xzd7WKks+2HetbytXQD5KFZC9tHYzPhVksO5k9SY+thWZW0nSGj7nN79eDAoPPjq9DxSgywE+DUD7Nao59D1B48i1ZMzmPSAc4e7eg9k9JJvZQoZnzD3U2JRZL4+/PaQgoAd/mrE9P1Vsb3rsoz1lxAab30OXPYqPXmPws4w9ygpUn+g0gz3+v7X6HGN8Pdr2/pwmJHc9ICLvVaVrdD3UHENqjQVzPdYA+VkHT3I9DHhjQMnycT1bu9DOj8RxPbfSz0eXrXE9Do+pG0WicT02R1VtvJxxPXw5vYUNmnE9xp+QIMOYcT30WsN3JZhxPTh7GtPal3E94tJUwreXcT23NKFpp5dxPYlKqtmfl3E9EGjFYJyXcT3uZ7fLmpdxPYvpxhSal3E9KDyLxJmXcT0sTnikAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAnQIAAAAAAAA=eF6l0lswVAEcx/FsqVyKGdXMTojNZSw1NTKI/xkU06QwlKnEmq2YaCdFhYQyScuDOLMUm8y2EVltaVFbmGqHLugypRvn7OXsTi7tdGJ3GtQ5Z3vx3NPn5ff2+37YWqZKGclHxJJ1uOQ7Bq49+/wyt/MRw7sssrEhBukQnJLzPRCk2Z4juuS9CeEcPR12P4GDNK+9aGXXshoZdkKGfqK2yIm49LglbCskg3hcQEaZYDTaYRkvwAg1HkXD2ikNLCIdL8S362BmbgMrcVIHXzDhwKpgAvonxgrJKwRUemYH5pEEPJL2skP36mG5sPtckUIPaPXHdG+2AdDk9qbKNj3sKfD91uekh/67AbPDYgJkLWmLx10IUHIRfvZ1HQgHUHmUuw7cYx9s1DRpAR+KLHnqrwXHaa+GEJUGam//WsFL1UBuuTjIxNJASL5ReF6mhocjb8jNh9Tw8rCP0uymhsnpF47lehzOPp/LEyhwGHV7lvanAofekoSZ0gwc6vAeJDQGB4FLcFZCIA6+Wn6KtScOX4+7ckVrcNhyxtW83w6H1kLe604WDlIp+oo7i8Go3/z6CBMGZUu3sW6QGFjfOaYzGTFgBYkq7H9goCwaUrVNYlB8QLZTMYGBtILdYUOZOjjWXTWOwcF7gZJwyurwmHlbypy0ZJSgfiuNJ9BByjAuz/iEslAQ97mDsiB0R6aM8nfO+6QWyqnSrrdNlGb/mvBblIbyiF1SSm1XT+dNStG1vhTa2q4ZH9pih2lbWsXVCTPdR6rm8hRtfbXYQFtiE6mjtak6oqFFiCQ1rfOnAaYnSW4+o6oxljGxOYjRQ+7HaE7yYpw1cRiVu90ZnTPcGOUrLZ4ctPSpqrNY/6/XnFaLxqyFZrcu3P2vfwGzz8PyAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/M+r8dJNYlA/gFz4lZmXcT1UJ+L35mFQP0NC7X4wYVA/Kv3+4RBgUD/Wmj76eF5QPwoaq/t9XFA/tE+8fJBaUD9Eyuf/ullQP1LXcflTW1A/uOYNlLtfUD+X4E4nNmhQP8xgBL2ExlA/Bg93h5pgUD+qTK+hFfBPP4CqI1976k4/EB7DghqoTT8ildijJyJMPwX9z0HYU0o/buLNaTMlSD+GO6XQsUhFP0CBuqQx7EM/s5uVHFHmQj8cVK69c6lBP+RnzJVoY0A/hYmXuW4jPj94hUF3i3M7P+z2h7PKwjg/EmMiiWQcNj8BxtVlxoozP4SfN3E5FzE/RACzaBqTLT/pYioes08pP8vYcg+/bCU/C1grrRHvIT8AncaACrAdPwz3UQ+4Sxg/UDLdf4CoEz8k5DHqVnMPP3KE4kdB3wg/h9PCPhFyAz+2JbhLQg/+PtAP3vkV+PY+M7rK5gNa8T4S0VXdMuvpPi6n4wWYI+M+HPupWHLy2z64+Th85izUPssuScjyzcw+GN+In2hVxD5Iz12e9GO8PjIzLNkbmrM+qNmRQrzFqj5oS4oUnBWiPpMN2icJK5g+sxPBidXzjz64Qetnv+WEPhmiCjPHC3s+Ku2dyoVRcT5i+8vo1fJlPhG7aSuch1s+RFWk/KUWUT7oJT1C4P9EPr8azMPnizk+1mVD8OPELj6Su4Lv/1giPsba+RVrqxU+tPnB/5hZCT4JYP71vGH9Pd3NgvGX4PA9k5yAGD874z0AX1de5cLVPRys1e3egMg9nOB7KrSOuz3GESzT8ievPapZSIrU5qE96kiirY9GlT3LJDQUfsaKPcw8/tvKSYI9Ei1wMxCHez0yD84MOb92PVY/PEBnPnQ9jAZ8LdDxcj04RUABs0ZyPWqUYulr73E90Ns93EjDcT0eXI8TJ61xPTcCucIlonE9pDeFbbiccT0Pt2zlEJpxPZawSDnHmHE9chPckCiYcT1qNdbQ3JdxPbiZ6u64l3E9IsphEaiXcT2uV6szoJdxPTHlnI+cl3E91hd145qXcT2meEMgmpdxPRoREciZl3E9vqFeoZmXcT3XRnmqAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAVwEAAAAAAAA=eF4t0k8oA2AYx3Et4yBpDrIT+XPZgbWDg5ayRC3JQin/TkorkQNrrbVmB0Ii87/YMMNsQ5MyE0PhNErbYTVxUbNdlnKi9n1Pn56353mf31uvI5mOG6sfGzKB76yaD+r4Pfo8aHPj+gku7WNQnCt8WHWIDxvYtoaxZSyZR7+Yu9jDKTt6j9DqQrMDVbsYdeKKqEfFnhdxf8iKshHRP4DOXvHOITwdxh8tNivRWCPyqfGzXuypwIAEU7+prDlSak+lyFuGt4V4rsKuRlyQo/uL+YME+iMYfUbTK/rCqLrG4jO0eFF+g+Egqv1YcImLcYw9omsbbbNoHkO9BpV1KOvA2k7MaEWeVgzpRB4Fzv3xf1ok1JEisScXS6XYlKZP/4Tj4t8l3KgQ7u7g4CZ2r2LePuYfY3IP311ouMOeK+yzo2kG201YbsG3aVQbMDqB/ZOo28J/aQLFrw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAABgMAAAAAAAA=eF4N01tME2YUB3DiEjKE7FKhDKZWucigcldXQcu/AoLiZFRYAVvR4XR2QztdQAlqGAIiRsaUiBFEO53gDHuQMpEs4/LAstnKA9GtJkbHpSaLy6RR/8UCO9/LLyc5Od8558vn4yPHTCQEJvW8GUIc1ajuJOwjXlu0X+/6nKjsMsV9t5eYeVL1dGiPxKFtVvdnhEffb4wQjzQ6lQW7CQ57RmtLicPe9xp7PyVerdJkunYRFWWG+WDx5bXyvuydEj9qOXSkRPKVttgbOyQ/d8zlNEm9evcVf7FyQGFcZ5T7mKgs205UJeaNthcT3n2WU44i6c/alDFfSMw6u+fixWOLHLd3Goj5nGcHmz8hqk8ExA4VEG/8onZN5xN1LzZfCRf94szb87cRp/c0BNXqibc6Ou/Z8ojmByMNUx8TQe+40oPFC9m+c1m5xNLqyNuHtxLWvoyDXR8RUdOlK51biJsxNVMLxeRS6+XUHKLv4mDxl5uJtLHHge2biOEAn3v2bGJTpqphLotwHNWmx4v6XtNsyUbi/r9VP3+bSRij2r4azCDGS/rV0+mEudU5GSZOj3o6tm2QvfiFFJ/Qydw6TaANxPFKg2MyjVhwq/ykUqz7p2VDllbmirB5K9bLXMax3s51xNstbstfqcRZu0K9UFT6Jk2mpBCt2ryOL9YSiyssRW0a4vJPTYvsHxKRT7vts2uIH5c56uPE5KJnupLVxJ3mAG/TKkL3u7p3IJn4bUGO5XkSkZtqjgkT7x9qmNAnEqabnZdqEoipiZHCnnjiwBKXYjJO3rfA1x4kVp+JrN8YS/iPZOgqVhLn5ktfX1cTKk2N7c8Y4obFesBPXNM1GJ0STQw9eTxu/oDYGupz6WIU4dSrCu+uIPY2ahWzkYR72HQ3VjzurarbESF1V7ehKZw4X9Y/82sYEf6Ds+e/5UT3I8/+5WJKcEi0fhkxkqsZ/0ZF5J80tN9aSvw9UG6YWCL9elreDRK9ibY/MhcTp8xjteXvE8Hfu9OuhxJXHypmHsi/+B9SmpVCAQAAAAAAAAAAgAAAAAAAACADAAAAAAAARwEAAAAAAAA=eF5jYBicYPGmKyILlaTtg8Ley+d0aNkv08hZommhba+lrb87wkrbvsE1eBcPkC8+q1s2wFTbvsq3atMlI237tiVdr3YbaNtvjE/1lNTXtj+XsHrPW11t+5g/kZeMgXRC9yuW+zra9qITT066C6S3XDl3RAoonr/1qV8jkC4LvVXyF0iz6ZWsq9PTtnflNFj4A0h7yOinJgPN+/+kV3M3kNaYJ5H6H0hnmficMwTa98OMOTQUSL+8xuOVC6TXVd49VQOkv2l/L2wD0lEn39v2AGnlXWmz+4G00MPT9ROB9Bd+0Y5JQHoTo5nIZCAt/rgxBkRfuHG9AURb7Vs6CUT/0NGcDaKDSpTmg+hHPZ8WgGgN9YMLQbSJv/QiEO12yBFMq+wPBNPOhhFgOvdzFJg+dykaTEu4xYDpN0kQurNilMZGAwB+WP4bAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAbgAAAAAAAAA=eF5jYKAuWMZ4WUjv70Y7E+uzCas2idrLsEVJ3rCTte9IFju1+qmsfdaeioefHsjaX0h8e+7rViA/LKmyt1TW/m5gl+kRFVl7xTD7zdxHZexP3Crbd7idx77ddulZnp4jdoTsHQWjgBYAAMCUJ0k=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAATwEAAAAAAAA=eF4tzk0oA3AYx3EHWWs5KCvJaW0XJWlOK1qU1MqiuCi1klaUhiiNKLJhkyGvw+RtZTOvS63UlrJwkFg7LLSDmm2EwxyU2vd/+vR7np4XV/IjNqgKV8ejGDnCiT0ccePpBfYtYG4A9Vv4K+Zmx1Hqwf0RHHNi7AQdi2iew+cQqqziD3Fv2o4KP8oc+LSNdxb0mbHHgFVtKGnF134cbUdjBSrL8b4MbTrUlOJOMV4l0lnlUnIiB7vr0J6h/6MgD9SgWoupL/qGF/wUmoKofUT/Pe6e4eU1djhx6RKXxVzyHCfdOHWLK2ks2kH1Gj6aUWlEXwl69djbgu+N+Fcr/m7G4gaU5GHBbyrrkEzsLcT4D/UGOVkTJq/eYNMGdh3jmxMNbjTOYL0LQweo86BtHdWHGIhiZhODVqy0iTud+DCPkWF0mTBfZK8Fv8/wHxrTw/4=AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_168_t_16800.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_168_t_16800.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0bf6b692d3db5e7ecabaa340c208ba11263213c9
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_168_t_16800.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.00020791436765" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-5.3876319056e-54" RangeMax="0.00020790588733" offset="984" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.00011483497787" offset="2072" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-1.2098401028e-121" RangeMax="0.0001148156368" offset="2424" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="2840" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="-5.6381748341" RangeMax="10.19484368" offset="2920" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="-7.6396638713e-137" RangeMax="2.4573148974e-06" offset="5020" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-6.3675814519e-56" RangeMax="2.4572146695e-06" offset="6192" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="-4.2782287054e-139" RangeMax="2.3661523665e-06" offset="7360" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-2.492870412e-123" RangeMax="2.3657538475e-06" offset="7852" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="9.9950053039e-11" RangeMax="0.00013984891955" offset="8344" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="1e-10" RangeMax="0.0002563282497" offset="9364" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="1.0370224039e-12" RangeMax="0.002" offset="10420" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.230841961e-10" RangeMax="1.1234845032e-07" offset="11556" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="1.0176128318e-12" RangeMax="0.0010432249759" offset="12592" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="13728" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="98704.543993" RangeMax="100000" offset="14232" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0" RangeMax="0.00020789631148" offset="15308" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0.00011470893725" offset="15860" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="16076" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="16140" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="16632" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="17560" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="17860" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="18172" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvwIAAAAAAAA=eF7tz/8v1HEcwPGrMeRLRhErblTa3uN8X9y5t/IDE87XnC/l+9fLl67k0OZbly+Hc6IvVkSrZPnS3erW3LCFJb+IWKMZi3WhDSchqtf7c+tPaP3i/ctjz9fe7733otH2zv88964Pc30DK9h/O1DKM58yPIUn2jq7ut/R8Q5n+FWtO8KNTE1Es0RYvNb6nu+JsGlc7wMWRpg1qCknVve2i5fB6F1hq4qFsOzlOvubB8JpNfvkInivOYsbqpgIP2GHBLqeRtioiHOc4Ybwh8xzEi1X6L4Loyvgd1bxpL4zwuVXN7rmQF3xo7hQR4RnVwW2Bk5wXxps9tke4bPNjIpCB4RX6jn6vQyE5ZY+iUyYdwgKPBbsEPbkSFO9YF4yWdofCi1saPJuAHVqknw0QefcYFMEftW+ubZsi7C1KEpjB+TzAr1pMO8O++Sxa0v2euvjBJ0Ua/PLCnTRs8/OB634ieI40DH7Y+EIKMi5yOgEpXExJhbwb9VqpOI39Hl5/DYPutzkkAUHVGQZ0V+AQfT1lDvg5pv0K0vgUHxezwQYYLMrMIc9YrhdW9pgZk7xFBu0HmsJdwAP5z1VRYBHS54vcsB40UJmBji+mc1OBgO2A37mgwn7o2jXwGa33GM3QMmP+sEicCDV/kgluL7lv1gGHtDX+lINLu9whFXknWzOqxb0UyWkiUFFmYZUAha5cu2I83lFlXVgCX+XRWQOP9O6BXrwdKdJL80UBJDuQFwXYohhTw7RuH86hbjhpycizjFFpcQpafJt4tAor45oNht2n2gn9GokulfFNhPHBoKaiOYzyodE6WYs5UZ3ewvxjCxF7Y52KzGZr6R6fMSFavndk5QnlL6UBx97UfbJQikNQoIos3QiKBV14ZQZqkjK15fVXlJGUa5L1M6zoynbjNXS09X6+6utEOz5L/0DsnDZzg==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAADQMAAAAAAAA=eF7tz38s1HEcx/Eju3WTpd3klH6RnXzy84pL7j4hYjfKqlPOj6TSznYN1RS3EOoi52e6KM66VJN+7Eghq9EPhJguqhOVol9/XN1qZ+n9+d76I5u/Wv/1+eex5+v73fezL41GnXTan/4+03umPX2aM+0zvff7TN+n90wn/R89/+tdw6jN2O+zT172wSEiWJBr6h/w0S4hKr7m/GFz7UKbTsWYPi21ehkdmzXXDQiPvlLneG9jCp454u7Aen4jtse74zLVHmsR3jH88OA4QlgfyvmxYh3CfqNPat0xwnV329t7+QgPVbsePQF9ffZzdoYPwgf4ptZSHsKOiXkZLG+Eo71622zgO+e0Sl27F8JT8+nNLVyE58ns48+uRliZVeyd5Qn33Ot/V8tBeCnDyUq6CuHJbF/PKXeENQ9NgtUeCBvsb+eVuyEcdbo4zQP2g4GXnLNdEWZ0NRZOggKfcJ8mF4QHQpwimdBJIeY/uNCfTNjqZLBjYn2E3hnhSGnC+HzorgOyT++h7zs/ytERHZQjZrArVpQIfkInhDew+NBqyVVvBApvyexk4Ksl5eck4JhBa/YCTDEVjLWAg5Y2d13g3hEFK9AS5FdtLEgF/R98nhcDdrMrB++AO+WsXhVozav7ogcl5i8sRsF4/euL9vB/x7gqOhPcROe8DwBFG+o91oAr1SVOMSDWCPq2glbiFnYiuCxjTY4YFNCPnJSCC9rCaIdAVk3M02wwyNWyJxMUT7C9csFoOdNPBl5IeVyRD8pi6mNPgvI4W3UB+JlH18lBA0P1rRA8NOySSazoCrtRBHKVLyXE0s6yxcXgohyL76Ql2r0i0svPDAUQ72/RphFvhs9NJL4JGs8nlnYUHifuUSWfJva02pYSJ3nxZ4ma2CvlxF0BHyuJgi1zKDd/dVESL+9orSIW+T6m2lB2ijJJbVNNHBRNUa1p4lGtK3KnDKkMpZzoC6K03iWkFF7bTMnp306Zf34b5dtgEeWstxGUNcxISkWncQ8NMXaqm1E7sdEGkdG8lP/+S38BQJb2cA==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA5wAAAAAAAAA=eF5jYBjeoFCgpJLFU9J+Urklw9paAfue135Prs2WtX92+/H6Yy9l7O84L6r6xypn/yv7f538X1l7sflndNo55ex3XT/0tFBKzv7jpn1POb/L2lfV5hyu5pazlxORt5K7IWvf7nx17Zp3svYPZ62s2rUBqG9aZHTYSVn7V6rim31rgOKcjjEJU2XtmePkjhdby9rPUA9kSA6XtddVyz2x66mMfcas9nVKnLL2pnPs5F62yNinXdER4lotY68w6dHKCKtXdg5N3F3mx1ntYf4I9vC5/fNdvx0+v46CUTAKRgE9AQAPClNpAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQEAAAAAAAA=eF5jYEABDQzYAbHi6HxcoIFEcXR5OC3x4aWz7Nmv/DAFAduN/4l+Pmu3/NThqq5/l5d9iWr6s2iPpP2my42Z7q8E7VdH6+T7LJa1115122reXxn7uVs3dr9mlbP3Z0i4L/df1p4l+m/TMnY5+3Uq1+epSsrZf3rEeOzfZ1n7pqIvgQ845OzbTTbPOHJZ1j6GW7DwwitZ+8lOX8OOrpK1v1eRUvvjsKy9haXfs9QSWXueWLtvGn2y9uxqH30vG8naB3l1xs3xk7VfJXwvtu+2jP2bP8eXbf0vY5+a1rT3+i4x+3kiUmFzFsjYT9gwd+nVVY/tZF74vzg2h8Ee5q+qqMezTI/l2TGMglEwCkbBIAEAEYJt+w==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAABQYAAAAAAAA=eF4Nkn081ucex02mPN787t91lUi/8lxedx6GuMNv3M56mLFqoskdlSJs2uhFTBFuK846qVNqJvdI3DIrT4f0qfGKnXOkVltOD5slEQrFJHKuv76v7/v7uT6vz/e6rlK36F3Jcgn6d/aP2J41hkv9jRv3d5sg5W6b3KRHgqTJjlj3nabQXHulF3LLBH/8eiG7Qc7BZjpHtSjPFEJ4EzcXI0XZEq386zusRSMDPbX7AyIm3tzRObWFw4YQC9fJLA4BXkGuHc082j+fsTCtMEXpdxJuOIKHqedhHUs3CmGdU8GtKA4+UyeV/X48fGvHW1IsOPQWDixRl/IYPH7n4uOHy8RM26shH+ktE6tPO/UnSnio+tLj5KEEI9xs5pUeDlZdysSj902xbzKm9JJaEJuc9Rc8+UEQ/crim4oDBbH8UZJV+FNBnDhQ+3rvoCCGV2a3d54TxEML3+pEdgriSbrOyMNgmagd/tR21TNB3JJuKz1eKIjrvOtKNHcE0SLs/OOILEG09LH0Un/Lo2X9tF/8lCA6z7/inKDNoTXXP8b/N0FU5s5ZDUwLYhQMxtYOc8if3zgz0i6IN7VSqjr2SrHb2jNR/lIQvSyX6t+9JYjzvzLSOdgliIvOnDoSuVMKZ/W/HtVcFcTnzy01F1m+9fvLXo1HS/HfJLPYFrb/WTN17bwiQZw9dLBu2IqH2+Xk6917eLxXXHuhNIfDt84RbvodUuT05BX+LuWRZZlbad8nxe3oSy56zVIkf39k2+cWPCY7rUsub+chXCmqSHufh1xufKbHQIoBoyX+zqd57EkYmE0blsL7oHJFfRGP7aMaV69pHuFXJvIkLM+nud7aS1t59OoeaPw9lceAQb2Zz888anz9Wo/+k8eGxoZnp4wJZmvtXmxdSBBjXtzRVcZypgw9M+ji4UO/ia0e59FY4NGby/7HUq0/L334K4/QkOxzWawPKKn6rLKeR5HiRJb7v3k46ilpmB3BvL+CZGMrCewzg84kLiIYVTzsMdKw/Gt0AnUoQW5IzIBqM8GPTydWyF0J0vcN9ZdHEaQ5xmT8tZ1gQaNzf9EWgvFHwYqrYQSRkhcBTn4E3w/Xr3yrJPD2c1nk6EUQfPjiKZckArlO7NfclwS6r387lJZBILuV0GpZwHzH4mQmjFcFnGtZmc/O7wp/FdhI8DK0pHlaxfj0w+aILwg6PPNOjyYQTCka6lKOETicuH3b/iiB2QcZx3srmY+Miw6+SBCUrHW9rILAz0Xh+FRD8EfTxsDHbH598GX3B8UEx0x+kf/vJoFxuiLr3R8ITkx1R5/8iWB663ejHnUE1U26wWatBPd8viyMLCNYmq86cLKT4NJYTdtCppvV3BvxuEvwn/Z3PKqZj8Qsfq7vMYHNjw4Ruxgv/Kr8WOg4gcrqvrXZE4LLpYGq6DeMt92zqWG9uSgznBkgqBncNDQ6R+C7IKJo4xjBW32Z/IEhxcByzws3h9l9+s/waxdQxCbq2SYYUDzfKDvfzPSr2jZHzJNQrDpBi+KMKaK2Ojo81KdQGBPenPHTOudTG5ZT7Ev/JWUFR9GUcuQsb06xwS7YgUop1k/++Wb/MgoH3brKBiuKkYpYlT7TP9GSvShg+vS1SvXHMqbXWpV91Y7iPSctWYWcotKycajFk8LJ/lUgt5rNEZz+sSvFDLfE22YthZ3Xu8ldvhTt1ovf7vgbRZiiwP/MeorlGWvUuwMpQrfZu/euY31+0ArfDyn035TqbmaV35D73G8TxZ4ni/s0mynO/7w64vVGlu+Ln3yygykyHcorHT+lWM213ZgIoxiqML/WqKQIsnjjH7aNQq0afFa2h2L/ueZryYxP7Rx0j4+mKA5VRJZHUXwWsuN4eTxFxVRIahrTLX7Q90nRPorUaW2D3DiKMvdSQzGJ7W1o+DqT1XeMB8erUplvVUVnXDLFwgnfLslBCtu/N/2jJI29Q4hSwx2i2HRK59Eo42uqvcTuwxR5OXds9rKaYRz5slPF9tUOmkdyKIYbywPsjlLMubmvvPs1y/GN7eX0fHZvhk7h3Yz/H+GH3VU=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAATgMAAAAAAAA=eF7t0P9PzHEcB/Bzx3Wt5mu+RB11q/Vlq5TWFd4j4tSsE6VIMrnTWRRFKpxsRFpUSnVZOW3mlFM5LtWJt9HCfJm0ijbamrDyrYjh9f58Xv3gD/Cbzy+Pz/P5eb3fe+0jF3CPdgDN5RVkYk7CbEB//T0v2IR5Auo9jn8p59UqsffCPIL5Ns6peQSfsQ/F/AXnq3DuGPbBaCb2q3DuGvZjezmiK/H72P5PUQc838CjdcM+D/sw7BOx/4CO3Xcf5zSoN87jPoJWNLbJ9tadrKbcrrqtUSzvNK82eu68rI9eFnZRXS+e4XzoebbFRUaUraItav1XH41sbZ10wW2SPVtZ5qQ2k6pm04eSs5SkVA8rW45Q8uDNHVebCko+3VtzZJ6Oklk1Z4qZ4uqQ0mowPXiJRFVGyavGqx3O5ZREtsxyLYLzc5vbU6JKKUk1D4Q8LqZEHmDJOllCSW3us/PuZyjx89g1uQu03H9kbyikxGWyLDqsiJI6m0KvTQWU6Kb/vP4WXCW27Q45TYlMEuvdBJ7/Ujd/+ylKdpVvFtwDGz8m+9J8Sn57+SfZQJbIe/ujIPvYfDcdA0uN2e2zwamnaxMCwIgev60S0Klb0TYNTDPsE84ERw1qkx24Idq9ZDGoSzmwfD64IjAycQ84UTVcEg/23Rw0G0FjfE5eGdiR4us7CIa3Di7qBMVDvQ4esMfIxPUKO7AiLEMTB3roFw6HgvrjxTknwCwfyZwMcEnEC/taML279mUlaKjSnWsDC2LCVRbQQdbn+AqMGG042gFKC5Z2vQeD1jtp+sHR4bv2I+Br6bbIz+CkJ7HCX+A3z28ZP8CE7gU5QvhvgyflrwVg+IPCxAngaqM4bjyYGTLTzwrsfPg0TwxeKbBWSECPg8l9rK9Utb1j+dIUt6vM+ry0UGtQqr/mxqT+VlpmUaAoiflsyFTBtL8uL2T6Pu6oYS40uV9kigKTTEz5cU09c23kQCPzhmeMmdkfEdDMTDusaWLunhbfwtQpvnP9uiufuBwgusUZJwyyMC9sc+LU9m/kDN6u5NzYu4PTeYuK01Gxh3NSYzJnVv5ePqemcUqF6Zyiyn2cyTf5bJ3B2/mTV9bDm+u4n9Nixbvb7b//0j+9TF8QAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAASgMAAAAAAAA=eF7t0PtLU2EYB/B5WctmLOxiyWamXexmWZoj8sWStJaVEdVCpATT8laa6WyZW0Y2SjTNy7ygYrPIii6zULdEO6WEDQoJzTBQmSV5I010OHvecx5/6A/ot84vn/P9vs95eThSHvuoBtFbnDwl5njMtaj173leOGY+usWGeynlVIVi74V5EnMLzkVz8H5hL8M8jvNVOJeF/R5Uif1+nKvDfm4vCRqM53P7f0LF+L2eQ+WJfTb2B7A/h/0QOndfG87FoFtwHvfhNaPFisO3CysFOR0ea5VEt9p265vGF2ce2dw/3Zm3RpYRqZUdrF2kMVYShfdQrLdFoa+8bimfft9CGo629RhcjcRO33zJScuQNquf4lo2Qywjw44dZQxJi3u8azloqMgdNZUypFsuOh8G2fwi8YKohCEp2408E3jTlWfjX8yQ1lGxdhasSFq0T1nIkMSfSv2OIoaMuwqPdd1liEYxPpRVwJAzma2tEfkMsQ+xklnQLSW0RZzHELnbhnod2PVSahDdYYh6l2pADY7oH2zyz2WI1i9SWQSeG34Udy+HITqTZckA+JE8798NygNkRxLAdnOsSgQuvrF2lSc4Ie+R80FN0ObwhWBZUHO6C9jNj7u2FMw928nsBe8WRWX4gbH1PllXwRJhuyIGtF42SY3gzkz1hRp6HlnwwQJ65J1q6gcNPn1mX9grrDnqpBtY/mNd+1mQnzrVfQyMbqiU5oNtCxziM0G7958d68Dg4P7UB+DQ8uRxE1i6UX3iLZgrvCfqBd2Pi61fwOfV1ocj4KHomZpB8Mq7KuEUOL1QVjJB5y/1HJ4FZ5bF8WfAgsC9EXbw36r1AWYbUPjJ6RcfdD9gsKdKoi1ZArA2dHhiHigwb1XOB63zuj2oK14NODiAgYLezzRHPP4QRrOtqYZQB9abNGx2V6VRO1dIqqlE7VVC1fa1PKX2f/tdS5WETL+ivqy5WEf1Fdc3Ur/2xjRQZ22rjNSVFZMGqrNzwWtqyrMA1vQOxybqM2Ufm8fuB7H5SMM2Vs33CNaxOjnrxycJrPyQGNbVOcms5swkVhfPVFavwBTWG4EKVukk13t85XK1jrNJkMbdP8blbAmXdU6cyZ7//Zf+AedUcrY=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAATgEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qLgDVL00lH4NlTeCymdCuA2lUP6BuuPGmQu9eD5C1QeYNmf+O6hgrXCOI5c50WTxYf2D4cxrDtgtEd0mr7N/qx1XT9eN062H7XTKqt8kCxy2s/okwHvq9GE71uYlaz8cPWxnH+g+Qe7cYbv+o/KsbjcP23X+/rPu7cHDdr1TahZknD9sp7Sfob5r7WG7lfPv3M/cc9gupOD7yYmTD9s511jNPr34sN0lXrELDIWH7W5FWn3uaDxsx9hxMdLO7bCdWKpo48qIw3b8Lx6/38F72E5mlYiZoOphO6FP4hGPzh6ya+0JTWd7fsjOjff2OwXviXYHVq1tEtmxxg7qfYblqgbqxnymcP4oGAWjYBQMNAAAEC969w==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAATgEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qLgDVL00lH4NlTeCymdCuA2lUD6DE0fZy/1z+D9C1a8rjtsSvb/YzuFdW/+L2KplFuKCjvY/DtjtOazc/2v+druA699ebOs6bBfg8lG5UvawncBt84l7Tx+2e/fE7Lbj8cN2IpdfL3p45rBdA+fEd1k3DtvJ2Yp/f73/sJ1L38buDWcP2z3qMbFYs+qwXdPlg273dh62M0pZr7ek/7BdmfGbvY4LDtvtPTAz81T2YTubdzF3jtUctlv3NV+/z+Gw3f/D+3QfBAHlVxZtYWY/bPfaY7L2LbnDdlJTPvRG/dxrt0y7QkrtwSE7i/O8ApXvO+3yNC/rab9Yagf1PoP0vOzJbV4icP4oGAWjYBQMNAAAjVaEiQ==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAA3AIAAAAAAAA=eF4l030s1HEcB3BPK3fNXDiPv7v7HYsy3SQPG/p9JA/V2K5c0VykYzMz5NnC3aSGPJVCHmILI0IjXJ4iTMrT/VEZTgrlSnJJSOr7+/nr9c97n3/e789fuZm+MGn6RMxs8rODBjhULqyn/5yMIk525fscK80giNWDbn1qpQTfP9vBxv8xofk1V4OFdRCGi8OyOqtBIvRUsN9NAxnBoBcaRfnLiX0cu/Kxo8vExMKNlyO168TZ6z278w9UYKMkMF37FQ3SJPaM+kg9oI+M8Y3cTKFrUKiTIzCFxuy7te0GpqDUVHZqveFCQ+DHCe8oLvx2vs201eBCv7ltznQGDlqxtfZBKjjYx+jVpc6wIafje8BCNwuyywcujo1hgB/ZismINAHh6Lv6T1oYXGkpnE0WYeCOHeC1NmFQ5vdj+cwfDPixaSuWniyIFt9ZO1TAAtcSK0tfOQueuxa1SS3YoBhkVpyLYMN+dwbd7ikbTpuPVhWssmFmwKK89zAHhFNh83ZCDiQ0aKtVZXIgY1ttTdHMgYElWwFbxgHNPh2xZIUDHh4VL7pVcZCph6h70XGIFpilSBk44MyduEtMHAb4AUZiQxxWghI8k4xxmHSZC8szweGWQ1bkP+Sw2pBsB+m4WVjdjHTVEXk6IuXx34qaUV4lvEVkhpzjXZ0qM8JhDatZNEVim25BHehuniGDFoBU6W0DXaTxtdyHb1HPxTlyYQ0y9t7cfQly5bJ1nAip2STmeyNrwrdanJHnZ2Uh1kgXnluwBbJj1ucDjiwXcR0w5OsnIwpDZFZXWak+0pumVGUi81cVX3SRKVyaJak04FGdDtIr1UBAKlBUmpAy7Fs3yN1FNA7JSZ2KQsdJiy9UDlJ71I7tJY1PLOokTbSykZJG0qXtpL6/ZG1UrseTcub4OrXjifR+ygXjakpXZQGlkyKTMp+XTrk7LqEM/5xKqXyfQrnksaf79t5f/AcKaDg6AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAA9wIAAAAAAAA=eF4l038w03Ecx3GRo4uahQ1ZnxmjqUXkzJ3v+7L8qMPFya9Iym9adY22oyKUH1fyK47qLMrVEZWVqRMT6XKuEBMn9IsiaZjyq+/321+Pf15/vP94vp0qFYMfU4acE06fk+rREGxwSVuOzIvCKK5Z3DcnUjEBVlQx9yMfu/TVqHVToASj812G62LrsBeU6OoP/GfYGrd9o2S5A1uvqmW+XOvB4pry9syujWD2Xr3e99QmseDWSk5xpxJruHlLuzdsFeuB5Mrjj7SgVEVpPHZFHzqbfNViLxhA3DaP7SUCQ4iokSkDHtCAM2SlEgQYwaeXLbaibBOglpn15yYwgAILmTMmTLDKZh89eYMFAkefvTw5G6ZkcYZhTA48b9M7I+6ygZBoe2VO1C7o6FsuMA3ngmaa58G+mJ2gCC2ebju1A7IKM1s+i6yhjReknLrIgSO5Jd0x+dsh/ExgXt9dK9BonygOa7eEihxPrfGfbEgudb6qMGfDyrsleXiMBczFj89jTeawHPr+vsLYHFi9/rpuOSwYfXU4wV+bBUuFDcPlRWYgGLvGT+KaQXeu0J7xngny5Hpe3WUmcNMb3WxdmRAbrGEv0WFChMlmp8VhBL9HzibsliLw6Opx9CtCIFJOfheIEYzqH9lTdRxBjZ9ds5Uvggz538eu+xBM05LEQTwEpxp8UwdsEBQUmCa6cxBAfe1nvgUCrlA3U8TE929q9lcwEFRrREWmb0UwRp09NG+MoN9IpdNrhIBO/cNfh/uFy1OPpyMQ+G8MUMPVDIpzl+JdLGCeP1Jxp8o7pCG4M5KYcT5uCt1Q3Q5XSpmTWOIKlYlaCLd+1WnFGJeaRqml0Yi7qRkGuAyfJIY+7qqIytuCe3TNpJqK+3Bk8QDhN7GrHqFjXvw00aN6ts8AYU1s12tCB2FJK2HV4J3nhM3eD2WEOeG/nhLu0ph4Qu7FXqTDb5Vk10OWLaTWpbdJy8bySeXRWaRKVjpps0Mq6fXC86TJsv9/8Q8yQkXyAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/C6l8dJNYmA/wF0fblU+cj0Mm7QnMGJgPyKwn5L4YWA/r6INOphhYD+R63/u+WBgPxIiczsAYGA/+ROCJ4NeYD9bKMuvTVxgP4wJlhUbWWA//E73IZRUYD/+bpCQTE5gP7RgQd/ARWA/NVwJzFQ6YD8YFjDcUitgP8hbMEPtF2A/VLnFPX7+Xz+qG4MXpsBfPwKbdPxOdF8/GCh0x24XXz9RnnEuCKheP3bc45VBJF4/2r6D1H2KXT8mZj6YddlcPz4tE+JPEFw/pnlE67cuWz8cy6C67jRaP6G07kXXI1k/cqW9yvr8Vz8gu4tohsJWP6xgse4/d1U/v6XlD3MeVD9ETfyZ17tSP0j1pA1yU1E/RTX+sODSTz+yiq1wCQRNP6aT/9KBQko/DqCYY+eVRz+2eI5AFQVFP3PzsPzylUI/mKnw9VBNQD/v8I5epl08P/s2c3PVeTg/AD1LkbnxND+yyRGKq8UxP3wzO3oY6C0/3GCa+xjzKD+/G+CR+6IkP1W0lQBN7CA/jD0mqVOEGz84fCrxsy0WP4OnWJe3uBE/sTIeiWoTDD8lCtAZOgwGP3zUtRATKgE/7Otd2ZF++j5VHqIbjEX0Pqi0saX6wO4+RlFsAasg5z5SM9FhUD7hPheMv2klftk+UJepngav0j5VlMa6aifLPkDsXvOpkMM+Ox/8B9z0uz5cV6oJQs6zPmoIT1jw06s+Zgh4Pg5joz6IB0AU8cmaPkR9xd8QW5I+/GHAbcHyiD6KsbvA2tCAPl/LBxsOfHY+Ab1l+kvSbT5IuHqaC55jPotceUGCmlk+eofXilmTUD6OJR2w6EpFPupykoJFIzs+mg4nD8woMT6or+EuY4glPtVARqKf0Bo+hqDGfi+SED4wu6WKkFQEPuDycy+xw/g93PdkJPD17T2FCR2YuALiPUjCxtHvidU9cGIIlkKpyT0y15P0j4m+PQ5qA4KSOrI9+vHDkkT9pT3c6PFhPxmbPUhLniIbU5E93Ai8o+Rphz3a+rneuQCBPabmTHowvno9aMP8eqWodj3wZyCmMGF0PRuHQaReH3M9WvHN9BF2cj3CtXECAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAA5gIAAAAAAAA=eF4l0n0slAEcwPHbeQkbKRraEcnbXV5Clq3fgyPM2LyMluMq/zCSVnMt75XoaHNTU+2KmfcQzc7kWNS1jMasXKPU85y75zl31OHi1HV6nsdfnz++/37PcEzyGt8KJG73GNamQSGD2EY8IwSIPlipdcgoQIQ8XolTLh9pwmS/omeTEVF2koxvF43EzHYIRsLCEC1rtYorZCN+lvYu2gB3hCOewGLdHBHmouSOSmyNNDCsl4RqBsLC0/2QtB24Lq8cYjN1oN858RR3VQFXos4PKFBB2OATc1ajCoIk3Oav71QgqGorsmLiMN3vy3udhMO3dQ5h14XDkvFoc4wtAQkO4f7b5QQ8aHleGLJGQF5Xe64BUUNPlSZw7IMazr180bcaqgamR2vtzSkCgiZZjt/5BLh/liXzt3D4Y9MwcqAGh82G8PEWJxwGo1udU3pUYDarfpZ6VgXSsi+cj/NKyOnm1c3mK4Fl9HFIslKCaVw87NS7Aj6JTZnytBWIbBR5WpqtwJhngsQwooBR4UbcQrECzItLz4vDFeAZWdTmba6AeymlnNAFDDI/MeZz+zCIRVWxqUIMRIsXTIICDGo42VMpqRjohrx1PMBAc6i61f4kBu23tit0rmQ/Ep9VfRgDg/Oux0VrDJoUUo0bEwN2R05xoRGFvMfvbRgGFIZKXhma9CgkDzRMGjdQUAa4uFroULAwrWL1P1GYvz2nubGOgmhPadG5hoKR77RzkHSsxPiwTYtCvn7ubzopXuK15kaakRN0ykB+Ec+8FrhMOrG3eXeaVGDrPjJOyu7/bSEhlRUs+w2QBgvLI3pJh9e4Xd2kQVtjm52kzuhpeQcpNyrajZJlSuxuJ2W0crMpA5oJf0ppZZQ9pUgm/0f9uDhTtkmZJK3XUipC0gnKK8szSsrG/rgVyquXLisoDfJF+uMfVXW0T7qzaEf9uLT3bUJo32azabdqvWgfvTm+b7oHLSvfnXbJct///KyzpQ==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/KGr8dJNYlA/ZLGWt+vmcT3H1/tVRmJQP/2aYAk6YlA/MIsZtiliUD+yxdtjGmJQPzVbViIaYlA/GLClnEhiUD8cTkz45mJQP9BDwdN2ZFA/1ASM3u9nUD+MR8B37m5QP7Alhadte1A/qPxHBx6PUD+QR7F5urRQP3BB3EmaF1E/3rvZN9rxUD/A2LSTYLxQP0b24LrDeFA/ZzrfMUokUD+2qQ8z5ndPPzaAowLneE4/8KXpQ5pDTT9yIcmX6c1LPxQRdB9+F0o/WS3k8TlASD+csaAQozxFP3ogZf/wy0Q/epJKAiPxQz/qb7HaR/5CP3Sc6TB4+kE/7LTDPR/pQD+WUTg3Aps/PwbOAn/FVT0/AI5o12QJOz+VguoN7rw4P7s/RxJedzY/HOIWBmA/ND8ZesVGFBsyP1B/h/rgDzA/biHzQptELD/cqZ7K1asoP63SHDKOWiU/UNAQ9d1TIj+yAS0Q5jEfPxVEExp0Uho/lnWAhC0FFj8c+2/yS0MSPz6Ft5oRCA4/HF704ld6CD98xP/uJ8cDP3r57ezIrv8++NWO0Nkn+T6erzOknczzPo78tcU+5e4+MDHXKUDl5z6IPuLQKlLiPhFOjOql2ds+iAx5NRn81D40gYgvwFnPPsJskrCGN8c+cMFaZ8QLwT5yH33bB9G4PkDMpvYt6bE+jEeyCWaiqT4i4j9adDCiPjxPGhiPmJk+jhdyTcrbkT6cATUcb7aIPqwYaSX59IA+5qSSox4Udz6qIPClEidvPlrQMfyq2mQ+5LjMQKSxWz7uwgFOfj1SPoyw7Ju01Uc+7II0KbzlPj5OtGnzPd4zPpoj2DIAWik+AtsuB3QMID7UZjm0xikUPmATVhM9JAk+subbsUQe/z0IcpYwAx/zPbwrBeRyV+c97t4DHINS3D38ooA0QhnRPeC7/GEblcQ9UP4eUNrGuD1aOc+FdPqtPRYeoKRLYKI9A9bMX1Malz2UqXoSkESOPSDN/IUwCIU9O4PwKU9wfz0x6hwnHV95PQT3L0mp7nU9BkoA9S//cz16EUSGf+pyPdLLbuXHUXI9uHcT9GQBcj0HAmuoAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAVwEAAAAAAAA=eF4t0k8oA2AYx3Et4yBpDrIT+XPZgbWDg5ayRC3JQin/TkorkQNrrbVmB0Ii87/YMMNsQ5MyE0PhNErbYTVxUbNdlnKi9n1Pn56353mf31uvI5mOG6sfGzKB76yaD+r4Pfo8aHPj+gku7WNQnCt8WHWIDxvYtoaxZSyZR7+Yu9jDKTt6j9DqQrMDVbsYdeKKqEfFnhdxf8iKshHRP4DOXvHOITwdxh8tNivRWCPyqfGzXuypwIAEU7+prDlSak+lyFuGt4V4rsKuRlyQo/uL+YME+iMYfUbTK/rCqLrG4jO0eFF+g+Egqv1YcImLcYw9omsbbbNoHkO9BpV1KOvA2k7MaEWeVgzpRB4Fzv3xf1ok1JEisScXS6XYlKZP/4Tj4t8l3KgQ7u7g4CZ2r2LePuYfY3IP311ouMOeK+yzo2kG201YbsG3aVQbMDqB/ZOo28J/aQLFrw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAABgMAAAAAAAA=eF4N01tME2YUB3DiEjKE7FKhDKZWucigcldXQcu/AoLiZFRYAVvR4XR2QztdQAlqGAIiRsaUiBFEO53gDHuQMpEs4/LAstnKA9GtJkbHpSaLy6RR/8UCO9/LLyc5Od8558vn4yPHTCQEJvW8GUIc1ajuJOwjXlu0X+/6nKjsMsV9t5eYeVL1dGiPxKFtVvdnhEffb4wQjzQ6lQW7CQ57RmtLicPe9xp7PyVerdJkunYRFWWG+WDx5bXyvuydEj9qOXSkRPKVttgbOyQ/d8zlNEm9evcVf7FyQGFcZ5T7mKgs205UJeaNthcT3n2WU44i6c/alDFfSMw6u+fixWOLHLd3Goj5nGcHmz8hqk8ExA4VEG/8onZN5xN1LzZfCRf94szb87cRp/c0BNXqibc6Ou/Z8ojmByMNUx8TQe+40oPFC9m+c1m5xNLqyNuHtxLWvoyDXR8RUdOlK51biJsxNVMLxeRS6+XUHKLv4mDxl5uJtLHHge2biOEAn3v2bGJTpqphLotwHNWmx4v6XtNsyUbi/r9VP3+bSRij2r4azCDGS/rV0+mEudU5GSZOj3o6tm2QvfiFFJ/Qydw6TaANxPFKg2MyjVhwq/ykUqz7p2VDllbmirB5K9bLXMax3s51xNstbstfqcRZu0K9UFT6Jk2mpBCt2ryOL9YSiyssRW0a4vJPTYvsHxKRT7vts2uIH5c56uPE5KJnupLVxJ3mAG/TKkL3u7p3IJn4bUGO5XkSkZtqjgkT7x9qmNAnEqabnZdqEoipiZHCnnjiwBKXYjJO3rfA1x4kVp+JrN8YS/iPZOgqVhLn5ktfX1cTKk2N7c8Y4obFesBPXNM1GJ0STQw9eTxu/oDYGupz6WIU4dSrCu+uIPY2ahWzkYR72HQ3VjzurarbESF1V7ehKZw4X9Y/82sYEf6Ds+e/5UT3I8/+5WJKcEi0fhkxkqsZ/0ZF5J80tN9aSvw9UG6YWCL9elreDRK9ibY/MhcTp8xjteXvE8Hfu9OuhxJXHypmHsi/+B9SmpVCAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAewEAAAAAAAA=eF5jYBjaYGbNyQivgEq7Sbq8qUdzVexPuz4vOKGjbX955Ybt1+217Z84RiyeAqQ7IvcWFNhq25tvV1Djsda2j4kP3LXJQtv+87QrBjVm2vb7blgfLjfRti9z6HPfaKRtz63pGGhoqG0fcuv0jj/62vbfHosFiwBpFjNBllI9bft5THf0JYA026mC6i+62vYTE6Yd+gekK/krlLWB4jr2p/7nAumrRbOv7wHS8zVjX/IB9b9Ps+6JBdJsHswuS4B0vHDR9AdA+pp92E5BA237E4sFJS2A9Iy/mu9DgLSU0+prmUA64t+V2gqQuERbYROQjjG0Z+wE0r9btM16QfKXF13rB9Jy8XoFE4G08bd9yZOAtM2sqI8g2of1leNkkHjG6TwQnXPEux1EX3DMnAKil2zymQ2i2zRuzgPRaa+XLADRMkadC0F0yPYfYNrVSG8RiJ7l5g6m05iDwbSYXwSYnrwpCkwb7I0G026aMWB6WzSE7qwYWjQAGZvpIA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAfwAAAAAAAAA=eF5jYKAvYDpybU/EYxH7ncmSp845ydq/vn7hTx+jnH3GbvfOP4Jy9mYzpl5xAvIvaf8qTn0sa18zgzGbY6+s/TPxNezXO2Ttzx3anpztKWvvlVC1O++7jL3fDtMZP6fL2J9Qa9lyNvSPXbCHz+2f7+rtCLljFIwCYgAAeR4wGQ==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAATwEAAAAAAAA=eF4tzk0oA3AYx3EHWWs5KCvJaW0XJWlOK1qU1MqiuCi1klaUhiiNKLJhkyGvw+RtZTOvS63UlrJwkFg7LLSDmm2EwxyU2vd/+vR7np4XV/IjNqgKV8ejGDnCiT0ccePpBfYtYG4A9Vv4K+Zmx1Hqwf0RHHNi7AQdi2iew+cQqqziD3Fv2o4KP8oc+LSNdxb0mbHHgFVtKGnF134cbUdjBSrL8b4MbTrUlOJOMV4l0lnlUnIiB7vr0J6h/6MgD9SgWoupL/qGF/wUmoKofUT/Pe6e4eU1djhx6RKXxVzyHCfdOHWLK2ks2kH1Gj6aUWlEXwl69djbgu+N+Fcr/m7G4gaU5GHBbyrrkEzsLcT4D/UGOVkTJq/eYNMGdh3jmxMNbjTOYL0LQweo86BtHdWHGIhiZhODVqy0iTud+DCPkWF0mTBfZK8Fv8/wHxrTw/4=AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_210_t_21000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_210_t_21000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..7d992ec63379427792deb1ca53adaba35647a797
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_210_t_21000.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.00020798424948" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-6.8530538248e-70" RangeMax="0.00020798252291" offset="1084" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.0001149184777" offset="2276" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-1.2339193616e-136" RangeMax="0.00011493093047" offset="2696" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="3180" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="-5.6425412169" RangeMax="10.720415737" offset="3260" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="-6.7565904209e-39" RangeMax="2.4581408154e-06" offset="5348" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-8.0995469903e-72" RangeMax="2.4581204092e-06" offset="6680" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="-1.311157741e-154" RangeMax="2.3678728468e-06" offset="8008" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-2.5424856229e-138" RangeMax="2.3681294337e-06" offset="8584" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="9.9991896403e-11" RangeMax="0.0001420814591" offset="9164" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="1e-10" RangeMax="0.00025919448828" offset="10236" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="1.041488956e-10" RangeMax="0.002" offset="11340" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.2308419832e-10" RangeMax="1.1331948474e-07" offset="12476" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="4.9273677543e-11" RangeMax="0.0010527973507" offset="13572" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="14708" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="98704.543993" RangeMax="100000" offset="15212" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0" RangeMax="0.00020797256969" offset="16288" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0.00011487531132" offset="16908" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="17156" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="17220" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="17712" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="18640" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="18940" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="19252" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAACgMAAAAAAAA=eF7t0Osz1FEYwHGRoq1xyYtWkpJdOdZi3QoddBlW67KllchlSVmSy7ARkpKuNlsYk0iJWGOYpiH0IrchGowkY2xqKpqE3MqQnvPb6X9oms6bz3zPzHnmnKOk9H/9y0udvaau2kQX/+n6PI5OctVOfJDv+PSHBwMnx+eaTexB2Hhw1nSAhnDtfCj9pQvCk0z9GGVQqUWX1ueMcPSEU5EAWuO6pwfNCeG9X/mDneDJ4UOJenCee6G3egycPsOPZDjAvPfi03OgIE2+YL0bYR2tAzoq9ghXjjgY+tghHEG35OntQlhvl2bKDRuE9dN7mrm2CPcLd8x8tkLYu8a+JM8a4SA1N69IDsKNaW+YNNh/NF3YaGCJsFtLHUMGevcF+apbIJxr4i9NA5l8l2WmOcKLDWo2WeC+TCHrHBvu9XYysw3c/PCdhQoYW6Bdbwcqt3LD6swQXp4t/T4EapfF3b8DLr3ykZWBkh5eixSkt8tM80FXQbNjJcgZtdUuBVNLBfn94NB+15E2UFRhMqUGc4PZhnHz0Kpvnnm4QHc+Hk9igZ6cJx+TwZgBefUJsDefxawFQ+R9yvdAnuaHpvcgnXW6uAd0cRaK1sM7KsKUhpagHb8W9JtBb+w2+bkNFLjSD3DBbrlwrRMYNV0lDgQHltxyBCCWmwdEgzWV4fGnwLzAyzNJoP6opCwB/NjlYJYBHu4ONzwP1ocPj5J/+7V1pT0T/MRrb7oOzo/N8a+BEwsd7tmgqVWU/03QuLZxRgL6bdpylCjqoK/KAbWOc4pugUa3V8pJy6LNI4kpvs81pKBnVvkX0mdDMnikz8S52xJFiatjiS8SvoUS3cT2l4lbgsrSiLggQEp8dpGWTUxNtyggBuPWPKLRicVCYq/swV1iVOG6YmJWrFcRsa2TcZ8oUW+m9v1smqhGl7IoX3fRSoiRIVNU8+ysFG3EoBy+6kqpHOdM6dvApxSxPSmtxgWUxRuOUHoJ/ShTvY9SXss4RumkqRCr+ivmjymax1O01naFBhEKr4j/Tn8DVzaqYQ==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAWwMAAAAAAAA=eF7Nz2lQzGEcB/ClS70xYiJU1jDh6Vrb2tqNJ8c2k1q1bNqUGCJHKWFsrtlGIV1bLLuhjUmXDBlHasudmlYZa3JEB2O0thLJVTP8nmfHDDuzLxm/N5/5fp/n+T//h8GgI2P86a8xzeZ6mYnmenP7fo1pb5rNjewvrf+r3ty66ciy/BRRvxeVV6bpkr9lOwi8ztUbxqfb6rwtqkSGvEosay++ljShyCHx5kQL3hRcd4bNFD+o4bNcDG8v+M7CYdpHq7i3XLG+EAcx/RFuVlSH17qC/IJFmgUI90S3T1oM3k9QXs+Yj7C37aV4JuTPbVvk1zDCgS6ZWbvh3ERL9/I7cxH+IanoyJ2HsKrvfmwTH+Fkj/CNJ/wQNlTq7j33RXig+LJ9KQ/h7rqW+iEuws7dFcLbPginbvUY5s5BeLNE0TgAKorlfKU3wsElq4sEHIRjdLwJ09gIFy7v/XgDPKoPWt/FQnim1hFvnI3wvl3SNp0X/Pc2n8kB0DPu5j0d9kSYFXOhQwK90D1gTSRkXgjrwGkw5a0usMcD4bEOJQ1OkBlFI1LKIYcFR0Q3ggPP0nNyQLtQZUwBODjX2uko2Cb371SA0Z+0sovg4V5mUCm43JAV+BIMsfnu3Axaqk5dHAvf9R4aJWKA4v3d74Wg0tXqsg8Yxbm04yDoVd/qmwRqpicMasB9pRr7MjA2jR3RA46RKs49B6212/PHwTv2tJ6VW4GqHXYZHHBO6Ce/maBo4YfjIaA6cluHgPhkXvZaMEOmd44C4wQnNySBd1UvXseBr9Q3ZuwFfTiO7smgZWdxbSqoyWzqTwF5otDBdNAlUW17CGzo37wzC/xW3ueeCQ65cepzwC/PuH3Z4JtZLQtywcld06vkYPxXN3YeuC5N30b68/lODSSLtfw84lJxzdQj5J660SOJNr1vJMSm9gQBkbuzWkoUNIvjicq9LhnEgSU5qcSbDQoF8daMl7lEu4OJJ+j5ZStVxKqr5QXEd48OnCLqlZsKiflftWp6r3DFaeJDf2uqfexTatyiCmrasOMZYjzTgrp6N5/q+diTepIdTH0gCaD2HhNT08pE1MYWCfX84nBqfusKqk1/BJVVG0kVSoyOmRpFxVZGqyKN7vIyytxkNFP6f/oT2nvaLA==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAGgEAAAAAAAA=eF5jYBgFyGBpEJPzi80S9vME7jDe+ctr/7Y5mjW/Q9Y+sGDfVsaVMvZJfZ9qNnHJ2e/ZdfmS9GdZ++f3Y60cpeTsnZa8yeRRkLOf9HVycqyInH2i9/9DUTJy9n33S5cWscvZ/+KfP++FgJy9B5O78cGvsvacsxcf6/wva6/6Kbvi6C1Z+6ONdn9+PJO1F3i3596rnbL2Tfb7/qackLU/UXc8T7pH1v5u5wXdo3Nk7ft1fzYeC5K1t9OIFNubKWv/wqPN46iArP2MqJOmhzWAtIsGF+NeGXvxxP5T0+7J2D9Z7fD74Gxee+dc7nbVLhn7zCNLV6w9d9OurJops7Thix3M3z8fGjZtPrIczh8Fo2AUjAJCAAC5q26mAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAASQEAAAAAAAA=eF5jYEABDQzYAbHi6HxcoIFEcXR5XDQuAJe3WnnhrxzzOW4YP+bC45hq6UDrw7FXWVd93bZ4+75TKd3Zkvbype5SO+4I2J+0UuNPWShrX1/lyezxXMZ+c7e+thifnP3EmheL7f7I2pvXOCyRkZKzn+d9KvKRgpz9rbzyGf7CcvbSzJN9Z0nL2a97zdK3g1XOfu49h4RCfjl7/j3fXxl+lrV31hbX8/sra/+thLlH+bqsfVq7pOCBx7L2WvODOf9tkbU/Wjp1zofDsvaCXJqsR1tl7e+WR5q9mCZrP/+0E/NqL1n72wuS9zomydpXz5v2yJ1F1j6OefalGDlZe/2G60E31svY56T8+3r3goy95N8riRJs7PZFGv5/yhtl7Oef0fnzZuEFO1eJxxeXPH5jBwsHntYQz7enJsP5o2AUjIJRQAgAALNBd+Y=AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA+wUAAAAAAAA=eF4NkHs81mcYxpNonfD+3t/zlFnzIyRKB1aU8ksJeUtKDlsHomwOUWyyRKlMksiKlUaSTZmV5BTbpbypDbFMhA5aaqzXuQNe7fnr+7nu67nv574v8w1HVTSTNHApJP7fVn8NxFs97tumoonu2MuzP5RoYF9hZqiyUhOHOnN90xw10ds0qU11jxYG7qf7VszQgnJHWqzcVILXshbXRyGGIsmp9c6SG4ifjfb1OkVLYHap/dmi3RzmFU6xnb6Ng9I4w6A6h9X7HXp9rKSw2PDEM/UUj8y9EiuFNwc+vmNW1koeDes+/UOiy8HJlri1Z/F4ms+/mvpQX7T7vbh06lV90SbBLniOhIeq93lDj80Ej2deGzetlkPQBmnM4HQO6d08yW3WExNLr4WMq+SQ9PNnt7ebS1BRM3cgYqKeqJBPWuNGOfHG8aSbgesE0dN60DrygQSq+Q77/dsEUTLtsUFclxYCj6xuXZHCYShbYf9zkSDuyEsLPaAQxJvytIXkNA/l6xaPiCFBHAm/+mhHMweHJ8uagusF0UAh00lU0xO5wHpnYYUUOs3zJp91F8SQ3c96tGM5zHWqrFQw/47784D4fwXxly1eQ5kjgjjPz+SHLpkUK8tXNS75SxDfmprkZ7QLYpBmlGpTPQddy4QPvnp6onDMpFsRxOHDtimH8UKK2aPq1HIxjy9/2/j3G3cOQpZp69k8KeIC9H1ureNxtLQ8tS1RCqNXB+Md90txJVN58Lk6j5abQ8fsWN465ypkGTY81kWlaAhjHKKtnUcVATzSfB9qrm6UIufjm8XFyWxOh9bLSFWCMqdrLj+Ws39/km/dVczj4/SaBUtlPP5ben1F+m0enadG3b5l75sGNw/VSAne/aO7PWECwXfug7myMzwW3bt+90wnj5IpJooNjFnKfT4n8nj0SWTJeg94mLtNOOP4G4/flxcZJ97iMRirGL7DmOJn1Dlfn0D/XceJKSsJ5t33UsROI1AOT1L2s/7kL/w7jQSC0lqP+itfEaglOWq5uBKMu3ixac1BgoPFNy517iK47r2tamQTgdH37VUvmPZI+0Mvx5bALvzPclNPAnU/47oqBwKLnvwTvaEEkhGz4LGdBL+kfz0rPZzNkewptTlNcL5xOFf4huDJ1jCVsmMEo4Ffb9lVQOD4j1p/ShiBzYwFfSNsXvSas2fCfAmuRuYGaWawve5bqTskE5ysjqt5+ivBo1ulZT4VBGHLTr1LiCNwjUkocy0heO1reaU6h8DvtHPhEsZXtiZ1E5sJsptivn9TStBdvfxlSznBe6+Fmc5FBDqB9hP0qwjqohTrWy8S5JkZSHQaCFRS9yyKqCcI4Pi6E20Ehf6TvWrqCHbO0m5S6WV5el9eFfaI7VlnGfDFMAG/vu+TvNcEzwvcsjuGCO4F60Z4sn3GGu92FHexO5aMZZuoURzxdreMZv3Hc3ztnSZTKEJH59t0EzRcNgjYMpWiRB4UCClFs1rVXnUly7dI+1iBBsX4k/KKVFWKOfP3PU3SoqgbazU7xzje1LzZYiaFufz8IoFQfB77k72dhMJZWlrwdAaFZ22Jl7E2BT/e1s/MlGIkjFcv0KUoOBDxoJ/5Rp3R5/hlFKM9lUeGjFnfyPOXv66m+GiyWUuCCZs/1jhzwILC/4eet/sNKYrenHY5u4ZCWxYSHMpo2H/v2n47iquGZUVRMorFfQdcItZS3PBwCb/kSFF8q/bHuU4Us/eqjgtgtFwvzepyoFielG3xzJWiem1gSokHRYNOx4qQzeyuJbutDzMGGvbkh7lQtG43bSvcwnTfw7UxWykGHh4OivOnyCy/f87qc4q7L3UyqrwpbDYGLbywk+LtZs8FusyfuvRQWccOCi0d0SXzG5ZnYXuJ2lcUFg2Lu7wCKKA9g48Mp5DrtR/mYyhkc97vyWP15d0+w2ujKWzdOLPtERQLBrgvDzD/vWNNd2QUxdCFQuu7cRT36ucLb5luFEwcU45TuG1aam/A2G1829r2JMWL2f7LyhIpjCtmbXqTzHKWpn4kMv0/zv/1fQ==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAxQMAAAAAAAA=eF7N0ftTVGUYB/DDKrtBQxOo2QiieQGKRqhBFhTfEpKrxDpEOUyY6zAJ7MJq5coOtsuwIjo7YstlIVbCC3hrcCghII5c1mNry3iBZdYGKFKYKFZguAQOZfS85zz80H/Q+eVzvt/zvO9555wwhr/yx1CDIJOHORvz1+g//51n0jC7okEuwo1ZMF+G/RbM85gtOHdQgJnBPh7zLM6fx7ki7CPRPOzjcK4J+6VzrUVj8PnS+XtRH1zfKJAfgP1p7BOwz8R+HF3a7w7OZaFBOI/nYbrQapx/D/PbAkwE9kZcV4Y97sPgd2OW/kcFzmfiPnXYH2MjL1ft6fJK7dhv775ucVWOFdjmfhup/3WfIWRH2LvrdYN++inbDTJx6qHoSrqoxWkrSIpQ3iLTjtUX7EOdJCUpNnvEzJHd83v2vXUE9JSH9HzFkYzCgOxw0Lj3ZM9KUGk6M6ABZ1ad62fPcuTHIk82sZojl46nFEfD+k0dSskjsNBNet/2JUf0i77pB6s4smFXYHxyJUc4/WXVBOj83mQaNXFkZEEboa3giLxv3WRxOUeuP1Ws8YU+b2fXSlkZR+41SJ72gXGHPnwSXMqRhU6511VQ0bjeN7SEI+vs0pIaMK1t8c4BI0e214ulLNhtjk5u+YIjsdtLvV0gp7h3poVCnmu1VuWAHVn3+vvPcKSudOqDZZDvB3dFX4A84czIZMG/Lz5aKALfcas5Vw7GWGaNJ0C7Sj16Cnz/wfHkSnD8+c4rpaC1Zpu5GeyeDAqvB/PFw28OgV5rHoT2gJpbjgp3eE/MtMOxAPnTrWp/ep7Brck+fmDqct0bH4Eb5dI/k0Du9cBpPSjae4CoQZ9L5tyLoHOFx7cV4MtXn+hugibDrh9ugJYXJJG9oPXFkAIbKM1pin8Msrp52SCo1W6STdA5588P/wDH6oZXzIOv3R6yz4DBsxmJz8Dakt3MX2DrNk01/W5T1efFDOhvZU8uBz1EieploLo356gYVAzvt7iCcXe9WySgqdYvivr5qwP658ATv9wNp9pb+0VuYGXDQB/NL4WXyGgu8CoMobbvTNBRax8fy6GG/T5aRU2YdTFSz07OXKMGHw6ro3YrX2mkFlkHG6haQ0Ir1V+ubqZGKgJY6qEtm9v49XPNN6lRmwN5c8eK26lidRRv2zdMB1Xj+InPio0RfL72sT9v7HepvB6eMl79qizexeZ03jLbYV67t4o3t+kIL5v4Ga9P+VHeZ7fVvBtUubzjqwV39Aiq2gTbJRpe91EhG9YK+ZOA/6f/AqNZWpc=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAwwMAAAAAAAA=eF7N0ftPU3cYBvCq5WZApoLDCQmKESQLFaaNc/INKxFlXFZHmWYBxoSO2bUBMxFJbdNNXFVuZSKFlhbEbKFjiHaWIRUYeEywi6jMFcNgbUwY6hBhZTIo6faec15+2H+w/vI5z9PnXHLObg7zUz1Dy1k5cswyzG2o5797TjZmL5S3gj3Qs6qE2Mdgnsc8gLsCFo4L+3cwz+H+Eu7U2AtQOfbJuLNgv/xcYeh+/H/5+YfRUDz/OosqCvtK7FOwP4r9c3T5eoO4k6A83OPzcPpRA+5FmBNYOHuxr8HzarHH63DwvXGWv4cW99fCQ5ijr7E/m9H9w5QscsNSjm1DqHuLX9Vw5Ys1Ta9f63JI/LcpzZeTirc2Cyq6SKeXiNfok/lWQaF+Ntl+i5zJEPs+nu8nYiE5dbeRItGb58o+OE8Rl7Bk6baRIpoFW0shKL6y+pLDQJE/Pz81LoS8L87+hQL29YOT/WuhL1lnO7JKT5EF3URWE5jIVb4800ARq1abGqGjyGaj97tr6inygBfTYAL3q9o8l+soYt7K375HSxFlbKwg/SJFriYECsdASWLK4PpaipS5w111oK9xJGXxK4qcjDvIK7pAEdsW+S5vyFJB70spOOaoydhTQ5HOUXGfBjw37L6q1VAkwBqZNg7O6wzVYWDD6TLDe+DOGzvUd6opcii6cH4afBT5I6UFqarDYR3gE39F3mnwp2NafRXY2fzKL1+C4XMfRp8F252y9EbwnrMrpA5MHBXH9IJX+HqLGeSax0SToOZtafsjMDD6bt56uG/PYpBoFaj+Vvo7AW+XdttjwMi0Z9kFoLwyNeMwaH1Vwi8Hg8bfHFGAG3/ONJnA3npRshGsOH7rwgCo7PTZdwPU8SIO2EHPCvvzIbDjyMWPJsAc60S/A+Re/94wA2bnhxybAiUhsr6/wac+22/+BZoG1HIPeDQ9NXsJjA14KFsJ7zFY8XSBAw4PRQV7gQlD05u44KZA34feoNG/dQctL/cPlw/Ij/Arp100W02+YFZwdwEtPyX/NT/wTlPLLJ3Ph77IobPZsTuJ1n9kUk1b/TGvlFYpcTfRNiuytLTazIl2WkvRdCut7Y17FtpPWzxmWu6ve7tppwPaumhnBttv0rb+c99K26Hh9NJufP9cD+19p4XJ1rQTjE5BcB9zf6eLycW/JTH5u9pdjKqduYwr8w8xBo1KGXOlnzAK4z5jLHEXMT7OO8FYFV/MeHCqhPFJPGuP7iSjLJ119SRr/APWb9aVMs7MsrkijM3Ho/6f/gv6b1pTAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAjwEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qLgDVL00lH4NlTeCymdCuA2lUP5TKP8jVL09lLaG0gFQdWww/oEQwcMfjDpsIVyGBZ/VNP82zzJ3uHSueZXCn/kt8QIZmzwP2ElHL1cvPLbZTsMk35mh8rCd0ZT15uufHrIz/rr0xeNzh+2WrOkNCNt/2C7lxKq3qjcP25VzGt/9fuew3fr2c3ucrx62Kzn+NbLjFpAu+lSw/8xhO7Mjx2P/XALyjT5/ETl42O5bZkNEx/HDdgn72qO81h+2my364lro9sN2VobTri6Yedjugzdbmc7iw3adHi0xOdWH7UrrmX9ltB22q9RNVJ0fetju/cV6vcMph+30lz/4Jal12C4s4dOy15aH7QJsH+++9+mQnb478wpmzsN2N28tfL63dbNd4SO3oLCLh+y2+t+VU2ZptNvP3uKoxznbDhoMDA8zNmZ/CPeH80fBKBgFo4AQAADCgJ+4AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAkAEAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNJfoLQihGIohKozgfJPQ2l/qLgDVL00lH4NlTeCymdCuA2lUP5TKP8jVL09lLaG0gFQdWxQ/oNq8dXvThZz2zJAwZ8vmT7936wafq43/uUWt/j19aMZxw8dsHsg5WHSZLPNLtWx7llJ52E7jvcTru3jO2xnb5bDonDxsN3dx7I5MUcP281vFD3He/Ow3S+RmW9a7h62m5OuenPilcN2hl9/FXPeOmz3e7XXp9mnD9sd+7bxE8ulw3a5Xy4/Sdl/2M7yj3kZ/7HDdmp8dSZP1xy2aw1Xemu69bCdn8Pm+sRph+1eB6n02Cw4bOdRzjvPrvywXUL9d7/upsN24g85V7oFHLYL22nodSX+sN2FJ96H7yoetnuZ4G8kYHzYbknBzj7DV4fsOB+pl37/f8jOuvB15bSYdXaPtr7pnHv6kJ2a6wORr+9K7ZZp5DpE8Ey2gwVD/FXJn6ckreD8UTAKRsEoIAQASQKgLw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAAgMAAAAAAAA=eF4d03lIk3EYB3A7yLKlM9ec7XhfJ4qFttJOpd+TJ7WmZrHENLVEEZ3hsbwambfWZnmtxcRjlqHZSdgIJyWKB20YTleZoiKV2qllmpr93v76/PV8nwee5ykruzUaIBs+TN/horW2JSEiymr5QUwQkub/Eo9Ox6LEpT/FMkEmijim39KuKUaTgriZoWol0klknd6zDajja+9sHfEYBeruCuaNbSj+np9H70w3kifu2ndDM4CI/NJK0ZoxlFuR5JWzNIVyuVteF839RHJXX+Kp2gzC1H2LwZ2bINAFhbUbGJCT3kNmKO2h0bE9Ka2WDzTrehvFJT6sLbRw6Azgw7sC7cEzDD6MQnFNi94eppgx4dIMe3BNVo1bMuxB0tV+aLCOBKmwJ7iJJEGl4ZoidxGQd6cw0seSBwnam1ethjlworbM5PGWDXsK6l1CfNhQF64bmHjIhir/ylXxHBv66Ka4alcOBA6J+YsJHECnSmv6WjiQ5PhKxfnGgYk+vY+PGxd6G72Gai9yQSnPYFzTcmH+WWr66T9coAWlnD/pwYOREkWpMZ0HlmGWfugJD0zGynNrp3hwZFltaGYTIPwcTh85RgDLtyrzbCoBVcFRL+gqArbTmYalZwQsZzcavAYISAn90UJOExAiWhF9XiYgWhGUvn4DCb1F7qrZzSRoZFHbsugkOAeFmBttSJDsDTXsZ5LQOnxcbMR7NVux+viGRQLrfX9ZvB0JB4iD5k1YZRZz7D521TvTSYHV6PWbTmP7d1+7bYON25jwQo/r+mPd5CWUOZc9hViaJFxhjc2JKT80hvM/Pkrzb8XmeoKwAtvcY6dOw64UiZMjsYRTzV8RttLzgwxhdc2uB9yxEuOkYic1H/fEggP2sdnCOgKrKv9yko11u55yi4X1yopTM7Hmv7p2bKP67PQkGdhP46xoGyrnaN7IVqzF32/ZlE6D04gy/hNrK2VIaPIcdd8T2Q2jlCJpXT+lbYd3N6XgCu0lpZ3z7zbK1u+nnlNemG37/xf/AD07PqM=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGgMAAAAAAAA=eF4d0v0z03EcAHDGOSJXeyCPfTYTYytJcrLvm3Rqhh6Vhx41V01JIcpDOluXVRY92KkuiXI93bVj7eRk18OmB44T6+lYHlea0EQ36rN+ev0Dr54I7Yfe/E/cmMz0ZwudEZiYOyQvZiKJhqoJMcliJ3G77JitXn6EsFiRtNXwrpDo33SjfYHjeeLyZYnwvlhGhGTk8fZ01xDtlzLLjQmPiXzaqLzraBMxfe+VqMVRTVQIjJ/tL3QScqvgtWLFF+L12YJeZe0wodN6rt14doK4k/SLkp1nIpix13SW9Tbwp1LV0TdIhhE/dFwrpMLcMKOA/50KRUH6zI4rNLCRZt1n5DhBkcKyqu6mM3yLkpLPMV0gSLd9jkZ3A/cYbvYWpQcsfCnjXZcgyNJbVbMbGBAgCNLuEHtDgvTEeOBhFrToubxUEwecjIoPUlMA9P29OEkKCQCJFSnU6+oyUOYLd0VoloL7KG+2eJwDvk0ta97TOPBx7FTrx1Vs6BUza+WJ/nAyyV0TV+AHp/3DBs7XsMCxsHqI0+ELDqKeylZrXzAlpsa/DPcBW11MW+yZJdB0QFa/pdsb6ipTFJLl3kBsqH4qlDGB2mkwBdozYTwtt9RQ4gU/A0bGT1C9YJbSHVVyjwGUQ0ubH/EZMBjKstn7mw4D0ak5lg/pkMFu+6JLo0O4c5yat4IOQpU6No9Eh+EUSuITLQLrT66JEQ0INqe35KZUIGh0EPIOFSJ4xUgeajyIYJHdTakiAcGdjlvNXdEISquzG2WAwLBrt+JSMII4TW2JYBkCjaWe/ISFwIljlRzpjeCrsujcWzqCvtJ9bxYvRhBRm0mmeCCYv36/hcgNgZ/DjGqbK4IfX6fnZbsgsBM9Zr1bhH+FKVZuwpLi05Km8Leg7fyweqy1k22OCBuz2iMwBVs65tPKwwaX88Qh2Cn/KQEb+6BKNeaFdVXOUTyxWdF8B1esUZ7Gd8ZOG9wO07AvuI1+VOy6sqHNFKxn2e7nZGxYvFRgtm5y2xKzzT7qWfN77d3i/v//u2o6zY5oZtRmh9h1KrP/AHTaTYc=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/Kin8dJNYmA/ixistdOg3D3wh8i2SGJgP/MvQRY/YmA/UgbRPy5iYD9KNaFBEmJgP6iZzHTlYWA/8v646Z9hYD8kAnmxNmFgP/IEEwGbYGA/NAgWLblfYD/wEjZ+d15gP6C2vuK0XGA/M2GFhEdaYD+QMmVP+1ZgPwgyG3mQUmA/6lLKIbpMYD/4mNMyHUVgP0wMZ4RPO2A/anv9sNYuYD9zECn8KR9gP7YYcweyC2A/HSWvjZXnXz9WfueJjK1fP4wH55bgZ18/qEAW/yUVXz8NRC2/+LNePw9XgPsIQ14/Xp54hCjBXT8L1RzfWC1dP7BXnCzZhlw/MMeeQDPNWz/VGL0VRwBbPw3UDvxTIFo/evz4qf4tWT9+6fcrVCpYP1DwZqXIFlc/AwFCJzL1VT/whteQv8dUPygQscHrkFM/cbw4nW1TUj+6K3mPJRJRP0521M4SoE8/vJWh5R4gTT+/OVugM6pKP9jpPsjFQ0g/Hcu4k9fxRT8Pq/Dn2rhDP7BZCw2anEE/8ScdBU9APz8vsP5vqos7P8j9yw5mHjg/WkFN5Sf6ND/0UjOuMh8yP/0uh0EMGS8/WqrYeA2AKj+epPpDVW0mP1FYXBVS2SI/J9VcYah2Hz9qoCQE1BMaP3ga/chrdhU/4pLvsnWKET+FiW4NengMPwPtSeh98QY/Ls2+E/9bAj9kDJrfHi39PuhNgE84Bfc+djpj9jIJ8j5T4LdjORDsPtSwt/AWruU+/l0y366h4D4GuCU93lbZPrJXxJ4rK9M+mtthFbrMzD5Sh5GPInzFPgKpJFWR1b8+cCMI9fprtz7yQ/fpKB2xPhyd/DaW1qg+JnDMe7vmoT7SUIEbtqCZPm4brsRXOJI+naxvMma7iT6daQXAUwyCPhJc1LnPJXk+Hj3buyZncT7m3gJ6Ke1nPuKw2hOyVmA+nHOmPNAqVj7KYkiIseBNPtLrmOLqAEQ+fiYVR4mcOj66x4dIO5YxPrEm8QmwGCc+UWN1kbAjHj5NBRVCpooTPsRIMSM9Lwk+hGz+4SIiAD4X2Q3y1pL0PX9fwCHfR+o9whKUEsdv4T07EWJbAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAEwMAAAAAAAA=eF4V0udPE2EcAGAcUBSKKFSGBNQgMuoAREng99oqliKgBJWg1AaQoKEmLkwUFwnKSmRoRFQMoVKlOBhqwAVaK+BAgTiKJHJXenelMlpEioS0vvfp+Qce35nIk0LZefTiQAR9x0DA7AQHOe08i/xcng1ktZxC5+Itz/1qTqDOSc+V/K8yZNBcDFEa0xB/jEcr7JLQUl/ulKJxOwpwlC7+Hi5Ah6dLo+KoUFRziR49rvRH69USu1zGGxka+/MGbrqgYrlcw5EsQBVLbt1fPW6DlsY4ZIqkZsgSp0ubteMQLnmb2GiiYCZNN+/MZgroObLUyiQKfNXn1A+yKXB1+bpNfIuCJ++yelw/URC7d8Hcx/Y0VFZF37uYSENr3bXkpHoapqnYU48WMfDgtFODby4DKw5W6ov/MfBTScU4XNDDlMxdEqDXQ4JHRsU4dxj4fZrPqmo9SOosK0mRHsqCpsKHjAzwxEcVzhUM+CeTIX2bGJAKVLyAHzQUO3vnGbNp+NQj5loX0+CpDJKta6Dg1xHnfZJ4CqqUVaJ9ozoIPvuhwFSqA82dA4f6w3Rw3anf8nhwCBJNAYaFJUNQstXLMV44BKNUclz3Py0U9L5KVbVoYeOWtqLoM1rouPHXTxqFtWTZGJdogT6426uPJsG7h2PubCdhf4Z+Iv82CZ5XOJpnF0jQyhVofiYJgsv5atsEEiZdZickQEKr2w6Vcg0JhfKjj0qXk9DgsTq5hUdCbWj6ky+OJERGr+0tsSXhm8/LiGYrAdcHW3nWGQICc9KGhWYCGPQudsckAfkdTSkWEwF7hGGBIUYCVIL7bl/GCKg3b1jbPkrArqBjucYRAnISMsp3Y2dv9HYN/iagreP330KsKMqDFGP5bs9HPLBz74aqzfidvSbQfRDL55ZldmMLhfbdr7GByzfrWrDtRW3yJmytbL7oIbaOaRLUY/+U9/PrsBynq2/uYutP+nxUYLO7uKtY1QX001rsYZNVxurgHBzBqk4PWsbqHhdnx5pa836a/W+T4mVkLbNWG1hfpqUyrP8B1Qmnjg==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAMwMAAAAAAAA=eF4BKAPX/J2q8dJNYlA/jiKbU6cWyz3c3E1vTmJQP7E0sPJPYlA/gMYGpFNiUD+OTxcxXGJQP8IRlutuYlA/KMct9JViUD9cLFzs42JQPypf0x16Y1A/ihE2vJJkUD90BZDLkmZQPxd+t4sqalA/XiWkG4hwUD/wU/QUfntQP6gXNu9MjVA/tjyl8e6oUD+RrgocadxQP4Z9NIzAP1E/CNDTBhg8UT/GyE305h9RP9xsj4AN/FA/n1K+7ePOUD+WkfFXgZZQP9MkxC22UFA/ghFCBQf2Tz/4HnWqIyVPP8nqM0eMKE4/FOziXqn5TD/34p0byo1LP9IU92buxkk/RsOfTUHJRz8ILPcRaupFP7H2iN9OekU/2S6hneDGRD9sYzZj9/1DP5YRvsvmJ0M/2qB3K1ZFQj8/CmtyrFdBP2NFj4HHYEA/xN4LzZvFPj8tzdhyJ8A8PwEaInoRtjo/JEt6sU6sOD+LluEZxqc2PxBk5XIurTQ/GChTUu7AMj+Q3uJAAOcwP5rnlue3RS4/VB7rxcruKj/G+6hCxM0nP5Iv4wTq5SQ/Hm3Bhj85Ij/Meu+OKJEfPy0nLYQ2Jxs/ehavehMyFz9JwM26dq0TP4Nyj7OykxA/XsbCSwu8Cz9NnG3d0QkHPwCD+723AAM/wjjf4J4g/z5y2DVZbFD5PoIuzffucPQ+QRcfoLlj8D65uWlupRjqPhKBZ1YNoeQ+gXiQ6jEx4D6MxlBfYD3ZPjWF4TxmiNM+XMpOIxYFzj4uMjb6HujGPnuTr3BHW8E+3nJmo1Qeuj4a8BUIyIOzPlpYM49r9aw+Do76amJWpT6eOsjf/TmfPvRI6ckmsZY+4CB5M3dgkD5R3mbnAHqHPl7S1UphtoA+Eu6b1P2hdz6oz664qphwPkqlaOEjJ2c+h46+d58KYD7uSb0yjRRWPovAfKK5ME4+48GOA8qARD5Yd+QUSqo7PmAXGbHhijI+36gNatOxKD5kxyEjh1YgPqfzZ5QGexU+MHXRV2cQDD5yV+VnMjgCPqILkrXmg/c98o8yF2wu7j1pONdRlUnjPYRLPP5Uttg9G8287Y910D3eD23QAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAVwEAAAAAAAA=eF4t0k8oA2AYx3Et4yBpDrIT+XPZgbWDg5ayRC3JQin/TkorkQNrrbVmB0Ii87/YMMNsQ5MyE0PhNErbYTVxUbNdlnKi9n1Pn56353mf31uvI5mOG6sfGzKB76yaD+r4Pfo8aHPj+gku7WNQnCt8WHWIDxvYtoaxZSyZR7+Yu9jDKTt6j9DqQrMDVbsYdeKKqEfFnhdxf8iKshHRP4DOXvHOITwdxh8tNivRWCPyqfGzXuypwIAEU7+prDlSak+lyFuGt4V4rsKuRlyQo/uL+YME+iMYfUbTK/rCqLrG4jO0eFF+g+Egqv1YcImLcYw9omsbbbNoHkO9BpV1KOvA2k7MaEWeVgzpRB4Fzv3xf1ok1JEisScXS6XYlKZP/4Tj4t8l3KgQ7u7g4CZ2r2LePuYfY3IP311ouMOeK+yzo2kG201YbsG3aVQbMDqB/ZOo28J/aQLFrw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAABgMAAAAAAAA=eF4N01tME2YUB3DiEjKE7FKhDKZWucigcldXQcu/AoLiZFRYAVvR4XR2QztdQAlqGAIiRsaUiBFEO53gDHuQMpEs4/LAstnKA9GtJkbHpSaLy6RR/8UCO9/LLyc5Od8558vn4yPHTCQEJvW8GUIc1ajuJOwjXlu0X+/6nKjsMsV9t5eYeVL1dGiPxKFtVvdnhEffb4wQjzQ6lQW7CQ57RmtLicPe9xp7PyVerdJkunYRFWWG+WDx5bXyvuydEj9qOXSkRPKVttgbOyQ/d8zlNEm9evcVf7FyQGFcZ5T7mKgs205UJeaNthcT3n2WU44i6c/alDFfSMw6u+fixWOLHLd3Goj5nGcHmz8hqk8ExA4VEG/8onZN5xN1LzZfCRf94szb87cRp/c0BNXqibc6Ou/Z8ojmByMNUx8TQe+40oPFC9m+c1m5xNLqyNuHtxLWvoyDXR8RUdOlK51biJsxNVMLxeRS6+XUHKLv4mDxl5uJtLHHge2biOEAn3v2bGJTpqphLotwHNWmx4v6XtNsyUbi/r9VP3+bSRij2r4azCDGS/rV0+mEudU5GSZOj3o6tm2QvfiFFJ/Qydw6TaANxPFKg2MyjVhwq/ykUqz7p2VDllbmirB5K9bLXMax3s51xNstbstfqcRZu0K9UFT6Jk2mpBCt2ryOL9YSiyssRW0a4vJPTYvsHxKRT7vts2uIH5c56uPE5KJnupLVxJ3mAG/TKkL3u7p3IJn4bUGO5XkSkZtqjgkT7x9qmNAnEqabnZdqEoipiZHCnnjiwBKXYjJO3rfA1x4kVp+JrN8YS/iPZOgqVhLn5ktfX1cTKk2N7c8Y4obFesBPXNM1GJ0STQw9eTxu/oDYGupz6WIU4dSrCu+uIPY2ahWzkYR72HQ3VjzurarbESF1V7ehKZw4X9Y/82sYEf6Ds+e/5UT3I8/+5WJKcEi0fhkxkqsZ/0ZF5J80tN9aSvw9UG6YWCL9elreDRK9ibY/MhcTp8xjteXvE8Hfu9OuhxJXHypmHsi/+B9SmpVCAQAAAAAAAAAAgAAAAAAAACADAAAAAAAArgEAAAAAAAA=eF5jYBhZgFOfbcd6LS57plKLp5b/1O2Dlh1dm6iqbR/Q73g0zUnbfscN3asiQFrTJlA420HbPnL9pMlNdtr2PStT3/baaNubnRBqm2Olbc9VY75jh4W2vWxsgcIbM237SVwbTtmbatvzuQjs2W2sbf8janVhmpG2fUndvVIXQ217gfalkhEG2vZSGzSjFupr27vNr/wqC6QZYn+sPqmnbc+++8myeUD6inftj6lA2m7OX/UVQHrNV//dZ4H0szfLNv4H0hkOl56YA/U9YDnlVQikpWwYRFcA6c3zNNbfBNJllzcuYgHas+LOozQNIC00ZZ+aK5BW+l7dGg2kK1zPZOQA6SiFnrMVQFrG99KKRiB95/DbK+1AOm92HEcPkOaZJn2lD0g7Ps1ZNgFI75usdmEikH7BOKNrEpBOKzD+C6K/ZKU7TQbSYuWV2SDaOKewGUTX9JtMBNEFCb9ngOhqTv+5IPpX+9r5IFr71ZcFIPrQ3UULQfQC8X9g+o6e7iIQ/f6fC5hO+hoAprN+h4Hp9FuRYPqZZzSYvvYNQq83jwHTIYUQurOCOjQAO2Deog==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAmAAAAAAAAAA=eF5jYBgcoEr3h2L2TWF7+V89sTskZO37VR6rP/8nay/ipHpJTVbOvmqBp3CMhJz9kXKuRbo8cvYr/L7WLP0la//VMtTq/gNZ+5PTbN+pH5S1D2Lx8SuYImv/TEvzjGOcrP2Bi+r2R2Rk7W9cin51+qyM/ayHoU7Ma0TtJ7zRvqzw5ZXdz4eGTZuPTLcj5L5RMDIBANMoN0U=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAATwEAAAAAAAA=eF4tzk0oA3AYx3EHWWs5KCvJaW0XJWlOK1qU1MqiuCi1klaUhiiNKLJhkyGvw+RtZTOvS63UlrJwkFg7LLSDmm2EwxyU2vd/+vR7np4XV/IjNqgKV8ejGDnCiT0ccePpBfYtYG4A9Vv4K+Zmx1Hqwf0RHHNi7AQdi2iew+cQqqziD3Fv2o4KP8oc+LSNdxb0mbHHgFVtKGnF134cbUdjBSrL8b4MbTrUlOJOMV4l0lnlUnIiB7vr0J6h/6MgD9SgWoupL/qGF/wUmoKofUT/Pe6e4eU1djhx6RKXxVzyHCfdOHWLK2ks2kH1Gj6aUWlEXwl69djbgu+N+Fcr/m7G4gaU5GHBbyrrkEzsLcT4D/UGOVkTJq/eYNMGdh3jmxMNbjTOYL0LQweo86BtHdWHGIhiZhODVqy0iTud+DCPkWF0mTBfZK8Fv8/wHxrTw/4=AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_42_t_4200.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_42_t_4200.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6002741988bb42714e13f694238cc8a16ac1450b
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_42_t_4200.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.0002074413507" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-4.4215426928e-70" RangeMax="0.0002074413507" offset="548" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.0001083761108" offset="1168" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-9.9012113982e-38" RangeMax="0.00010809235643" offset="1340" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="1556" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="-5.2949227271" RangeMax="8.5090145293" offset="1636" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="-1.27584389e-38" RangeMax="2.4517244091e-06" offset="2864" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-5.2257718696e-72" RangeMax="2.4517244091e-06" offset="3512" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="-7.6868229821e-56" RangeMax="2.2330696344e-06" offset="4144" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-2.0401404186e-39" RangeMax="2.2272229284e-06" offset="4372" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="9.9935785911e-11" RangeMax="0.00013466972792" offset="4600" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="1.0000000022e-10" RangeMax="0.00019865155073" offset="5124" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="9.9999999786e-13" RangeMax="0.002" offset="5672" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.230841961e-10" RangeMax="1.0414881319e-08" offset="6440" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="9.9999999786e-13" RangeMax="0.001" offset="6976" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="7740" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="98704.543993" RangeMax="100000" offset="8244" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0" RangeMax="0.0002074413507" offset="9320" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0.00010784123157" offset="9640" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="9768" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="9832" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="10324" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="11252" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="11552" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="11864" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAegEAAAAAAAA=eF5jYKAvMKg7/4k5ScYext+ybEmX5Q0teyV2s0luP9Xtg2c/6ucX17aPNU9lXseubR9QZbIrUUHb/sdDl5u8ctr2D9rPGAcoa9s77f/EUKWkbV/OEJT0R1XbPk/C7Dc/kA76OEnrtIa2vcfz56eOqWvbL7HYNytaG6h/4prf2Vra9jzBNfHJutr2X+bOmLhMR9s+wb7Z8qyetv3ZjC43byB9movx+Bp9bft5j+ftjgXS78z4bTQMtO1P7dse/RfIvyn5emkqkN+Ut357AJBmn2iY3gikFz12l6sA0kncG5/1AOmGTC/5TiDNcOKP70Qgffdl38N+IO2+TubBJCAtzi+aAKK3Le5Nmwyk422ibEC0UPiTCSB63c1jzSCaY239XBDd+nDJTBD9hkVoIYiOW/5nPojOu8+xCER73tgHFj/r5AHm604yBdOJZ8LBdJl3MJjekRYNpncWRoHpJ3YxYHqFMITurIDQClmo/FF6lB6lR2l0GgDS6UjtAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAArgEAAAAAAAA=eF5jYACDBgZUGgbQ+RjiS3JVdxYopLHCxHvYJ/e6d92ZKSHCebDw3WX+VNcbJQcyOzYm9S0svDs9f8mU1f9v/JmmZP9yzdOE+n1brXMdZ8rJvdSyn7P9Wilfhaa9wexFbj0S2vZrjkuWcfBo2ytszF0/VUHbPtnALWmCnLZ9HFe57mxlbftvLgyfNyhp26+q2qJTpaZtfy7s5MpyVW37ubeuOMZqatsv39E3LUBD2z7vs/2qR9ra9nsUD/p91wKaw2rBfVZX217nfUSbNJB+ZPJMQFRf215ReMn9+Xra9q+PLrlyHcjXEz6eMRlIW3dPtbI30LY/+vX0FGkg/Sy/fmsBkNbb/u5nPJBO+Hjcrg1IN4WyS9YD6fSPiTf7gLSXvQF3D5CuXPeefRKQ5ixezjQRSF/XuqowGUjfnydyAiSe1LC5GsTP1U2JB9HdZzOmg2i5LY/7QHROWtp8EP193885IDrixoyFILpG9v0CEF1ZrLMIRM/PZQDTlmaBYHpOqRuYZv0QCaZflIaD6awX0WBaoRRCb4uOAdN+vhC6p2KUHqVH6VGaOBoAg0Rlrw==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAXgAAAAAAAAA=eF5jYKAM6K97dSP5IKc9jF+x5Oqp07ul7Nl/BlY7bxOz/9L9O/1/hoz9/JvNEqaKMvZL+Q/FMllw2S94kP1+WrQMXN/cyLd3k4Wf22G3ZRSMglEwCkbBYAMAgLkd/g==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAgAAAAAAAAAA=eF5jYEABDQxoQCY+mhWZP4FvSnB60+XFRv9kHhX+8up/+vHejAo3PvuXSWXfpt5w2NG/Kzp2ym8pe9VnwXs3rBC3X14g3P0lVMZ+1g0N6TOSMvYv3yqdCLnIYG+nYyPi4i1jDzO3bCvz3D/BV+3Q7R8Fo2AUjIJRMDgBAMjWKxo=AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAdwMAAAAAAAA=eF7t0u1Pk1cUAPCqVXSYafvcJxB17raIgPhIITIcBXaps5L4Ai2KYskQBAU3RSqwOWWMSQDtCL6QWl8oipiAQCWKojPUI1ZTFbYAAYSOVRCl4lAQmCww3X38tH9hiZ9+OSfn3HvuyT1r//5W6Vwx1Jn8FktqxZC+zZ44rZOB/gxDgrZEDNNP56fVHkNgPqQI6DAxcGtP4ZOrajFo1r2bcl4jgoGx8V3DcQhe53HFQXMxWRPj6fHQhEnM7ril8i9E8CsePdH7NSaj4Yq3tlAGZg8die1+gMlvL3+wZZVicjSqvfyOM4JFXHJBsVBCjnTtjVFqxXA423Vwx01MatbehMubMXGynRLsrsEkpeWzlFZvBnL6clszv8GkS3LdmDQghiWNS7nZuzBJWnjP7S7NC3ULt/5pw8TH6p6VSPtHtN3N2lZMKpLPV5YZMblgmX/fb+xTUnB2SYU0HpPxzevzilIw8YyXtS+LxaRbk/XCg85dGLDccboQk7ShTLlTESYhg2cCntC58gcuzhnMRxDoW90SX4VJfXf5DUMlJk8j+7qC6DmVgRMitzhMwv2TtKdUDFxb6eU/kkrfs2G9MV2JicTQr/xoDyb4TlPskEIEI/Unxy0Wes/PmU1CuifZXstjwR9iCCopUuhyMckLndkskzBgiZa7hkxjYH/Ms5a8Bkx8zZflfmEMdE401IdvRDD2VybjvpOBZPt+7hcdA5FRhpoz0XS/W195s+UMZOsqM2QcgraqR8JgKQJdz8G61CgEU/7ureOmInCet+PNd2sQbPKWOdT0vH86e3Y61iGYs0WwJXUFAnlFrVlkQiAafmZt3I6AnQxsupeDIDvD5vVVGYJcThHgX0L/i/65ZoLWPTQvcJY207pE25dtVxDMUPlOqh8jOMxNKqW/Izh+rj2huAdB/6wGu4oakSxE5DWdS68UlDho/kBwNTODhW97NW7GUQQH+vVlahcWBqodsxI/ZiFMlT4v3oOF+2Gh6qxPWFgb7PPoEseC+bbP54JlLPzkWnAyhcYVJq/lh6i1/rkFsdQwJ/2PCVQXY9qqKOpwXYJVQ3UXrN6upkqthsZIvm7ffJuKGrnP0Ma7Sm/o43Xmql7xmiJcXvK+m+75htcYLeviPTiz/n3f26n2dl7rinPv8zkhER28o6UuHf+NP/jB/6P/Ahm+r2E=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAxQEAAAAAAAA=eF6zYACDhldQuhtCM1RD+blQ/moo/Q9VPUNiw+e5a9o3MbJCxf0Z+2aUSBpOO+3r5/TumjaPwclwkc9pa9amlBfx371xROHFjYyJwhVb7T5u2HZ8VazgDpEFd8JM/Y7Y6eXb1noqHbJbzbfy3p7aI3aRs8PqpIqO2N32/nD0ZMsRu4OHCnnYm4/Y2azWymdsP2Lnfv1nZnzbEbslT5M/s3UesevsnasU0nHETjhs9WK17iN2HLFaGxO7jthtuLdPW7z3iN2CG+tfJ/YcsVsa6Bq8oO+IHfvJhGZBIJ26/x37pv4jdjuMNu7xB9Ip+xoVoicA9aekhskDaU6JdU07gPQHdmWRGUD6ZtyBb8+ANPeSw4rXgfT9QtOvv4G0e5Pv+y9A2nftuuUsE4/YOex4f5oRSG9ZGP+HHUgrHN7+gg1Ii1/r0uEE0mLmm09zAOkVP60ngfizrAOKQPTtLKdNIPr36bAVIPritNO7QbS9zsMdIPqFcf0+EP0rbe1eEN3zxnI/iD73lhFMM93PBNNPdkeD6VXzS8H0dPMiMB0aWQGmmReUg+kbfyB85TsQukijEkx3y1ai8EfpUXqUHqXRaQDWuLPdAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAuAEAAAAAAAA=eF6zYACDhldQuhtCM1RD+blQ/moo/Q9VPYOJ2paupXsYWFmh4gyB3zw6p5XO/PDwzUZh9Tz+Gb+vi8aVK21Uud6wZPdT9iVmX2ULquT32D3eLOvTUxRgLdZQU/Mn6ojdmg3piZedD9vd+Oh9407dEbsnF1dUpJYesTvhISH8seWIXVRGxMnM5iN2e6OyVK3aj9idPL78Z0fbEbuGfV2SaZ1H7NSeLO5f1HHEjiH2uMrU7iN2ibbfXp/sOmLX6/0ppLb3iJ2S5gS/0z1H7P7a5Dp+6DtiJ7rx/K9MID2pPyT+RT9Qvst3Vz+QXnwrdmP1hCN2lz2DPrgD6Ts3vZ+fBNInxNo/rAPSM15MCXkPpGdlVq95BKQVzl2f9x9IN/heMPsFpNeVP1Vgm3jETpF7OSMLkC69w9bJAaQl8lMOswPpqQn1sZxAekF9rBSIPihgvQhEX9O+2AOi75cYbAfTq/ZtANERB6z3guhqz7TdIHql/tN9IDrkbAqY/nQ1ZD+ITus0BtNOSwrAtE5TBph+blsOEY8rBdP5uyvAtFUKhF4qVAmm97FD6BKNUXqUHqVHaeJoAADdtxU=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAiAAAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUrzMkHoSKi6hsexH+uEf88vmHc338p+Effthr66lN8b7BjmvduprlG3xT42v7xmwUE72cWXfYI/77VbvW4Nc/XuQ3Yv9eJdVXsO2XVUXDkxzXOjHZuR8YXFmw7ZQa1hCOc2q5DL64PzR8EoGAWjYBQMbgAADhIuGg==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAiQAAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUrfdpYB05FQdQxKPddNdlUuZljA9+vSyY99OpFrTIK+bbZT4HC/JWL+bDvj1YvHu04dtDMv4RdhMd5vJ5/8+cCbdYfsvD+uqnZqP2RXwvhNuVJ3md3hD6vLHyw/ZAc1lcGI9dxveYUaOH8UjIJRMApGweAGAMPvLu8=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAaAEAAAAAAAA=eF6z6logFF1x23bVrJqtguIK9gHsl9f+tX9gx9Ks9d2V7YddyWrtPVLfOe0TgvZtnLNB1N6rpUX9yDJF+3W/hINa7snb/1u1OTbmgpz99Ys/S52nyNmzfbGYNv+tvP22mwdMNWUV7E/kCTQ8sFawX2P3IGu/r4L9rNNa1U+CFey3SUYV2wQq2C9oVTV38lSw5+njPstjr2Df1hWndM9EwV6aeUY6l46C/ZbfF1ZIqQL1lTgb8Coo2L8IzM6aLgNUV6gzNVxKwf5yxdmDHyUU7OecKNKQAdIX1uyyWgj0x1pmqWsxQPqTpuRPcyDNdJdtowqQvj/91WMpIP316U0uMSAdYnpnizCQFrv7aLYQkP5QxCYNogs+eD4EhceHE/+Og+h9f9X2guhNx2fuANGHCnZtA9FSrD/B4faodg2E3tkHprckNIPp30L1YDq0rRZM770DCeclayE0+wwIDQv/UXpw0gAwTzNyAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAeAEAAAAAAAA=eF6bUH3t1v2a27arZtVsFRRXsF9x/E/PrI137YQc595jSfxqJ7XycsDuKjZ7vWeNXyaoCdvHvtIv6uCQt0/87NP9IUPJ3qRAfX2Gi5p9esr5y4s5tew1mid4XLuuZu/1iy9erVzV/ltnQFGmpYq9jtdr5Rl6yvbFt9cGCIQr2Yvcu/1p7zpFe7kUI84CfUX7hRc+OThcUrDXfhg3w3+Cgv2MlxUO9skK9kxrIySbXRXs+Z6sye03UrBvnCSx+KGagv2BWW8YX8kr2LexdrY9klawZ2E53z1fUsHe6l/e/mQJBfv9338b7gH6x+LKk8WpQPrZllpBKyD9Z6PrWhUg/ULbOVcKSDeFmi0UBdIrXSryhYH0k5tqJUJAuuWVNBuIzrv7+hooXHz2iR8G0Sud3uwC0Tpn5m0H0S/9/LeB6KDmM+DwM82dBaY5Z7WD6VuTG8B0H0cdmNZUqgXTXOch4b1kLYQWXgmh7y2E0LD4GKUHBw0AuyUxrA==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAHgIAAAAAAAA=eF4LWvvxkm9Sgv0Uk56GmdMLbSdPOiDbrp9gr3kvMc5hXbx91adb2s/mxtn7bapaFiEfa2+V8jZwoW60PYPkw32aFyLshX5Mu6onHWbv1lCpOUM32P5c2eO0VToB9nq5Ubc0a73ss6x73/hyudonfrqcN8/Ywf4f+w42Hz0L+ynmyQc8agztH3IeS8pu1LBvcPda8HyJgv2i46YzXr0WtRdcka6ZxchjP1Fu+h7PGEb7ieGLH8k/+2L3xox3TkbhC7voFgPhC6du24nfe3G1Z+85u3+X9na39hyw67B2PqMevdHu4at9ttL7Ftrdl1wV+WXDRDsZ64/Fc8wb7ZTab+2tiS2wyztR9HZTepydpKFfQvR8bzs7/6DzR/ss7CYp22oUJ6raJb3nV/O2E7K7Z/I1Xav7n+2a0rTyEwovbTWjZ32zTr5i+8HDqPFu4QHbwwfNvjx7uMqWb7nLt0M2021n1T5b9PF9m22IjUlWMGeN7dY7UwImJJfa7giIOz/5c5HtlA6+9RFCRbYlwXtN3Q8V2gr9ulNnvazQVqe9pq5nTqFtsUV1msPMQtsFqUtvqcwotJWo2njzGDBe/AyOiKwA0hN6/ZXmAWmvld+/zALSs7abfQHFm84/kSUguioishtEJ0meaQbRZT+aG0G0s8dVcPz6dS8E00yJk8B00o9eMP1vMyT+xcIhtIQthFYwh9BW0PQxZYTSAGYjdIw=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAcQEAAAAAAAA=eF5zZ2W2XRLrZjd3sfzDxa8e2EotSbqkls9j5+ZnqVpvwmK3qzQvfrvMN9u7siyFExe9tdVSPmNw8cQT2ws6tzYvs3pmyzD7b4azwTNb7bC1T6euemqbt3BvsFz5E1tB0dr7bEce2+ao5i3Ikn9sm6Eodjth0iPbOV9SeG/LPbLdtE98F/ORh7Ya/lu1Pesf2nZpudWZ+T60jXrcmMOn9dCW58GrQjuhh7ZHdFyUVzE9tBWpDgs/9+OB7apDjTsPf3pg++W2wxWd9w9sC4TDf0u9fWCbsGe/4afXD2yLrpX1BwBpv6Kdz34B/eN3U2jjOSBtza5WsgNI5xTJdK8F0ne6ZwWuANKzmf/fXwqkk/V+HVwCpFfZiziB6LQcKS4QHVEp8wkULm5rOV6AaGP1R49B9L4rJo9ANEdTFzjc6hU8wPTkM3pg2lxKFUwvmagIpjeEKIDpWccg4XxsNoQuWQ2ho9sg9JzFo/RgpAFXEMzcAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAHAIAAAAAAAA=eF67s+rjJd+kAPspJj0NM6cX2prV99eFGgXY557WSM867G9/wWaq0pqrfvYqkVpyswp87Y/ITvV13OZtH3shUW2foYf9loOzTBVSXOz3eyYt3l7rYK9WG9J/8YWNvUUZu5v2BXP7DZ/OFHdyG9sbvS9QTWnTtT+yQ8pIQVvN3qT/quNxOQX7O/dWHeDfKmb/IFWslPEvn/33kuwjKnUs9nEf32cddv9lV6XN4MgW+cFuf3a57J+9j+2ULsSvVNS4bsf01lnTNO2U3ULu+SzarnvtzkhnXlkxdb3dn/O2hktmLLBLzNy2JvX/BLuz125O8P3VYBfyYkXHwp4Cuzcn3Tw+NsbbtVw4+DvO39dujWG+wO6VVnYKvMvfeO5Rt/vjGxmculHELvzOgvDM9Yx2h1cU1Hxc+9b2au3MdQ9bb9hKfJxn9VfkiO02Ue+K7J4NtotmL1wbnjHHdorvo6+VSt22c+ILg5ji6m2tCnvW7FhXZtstbJM4r7DYdt+Py7syHItsvx20W9BxudBW67G2RfEaYHhPujrpwLxCW1nnk0qsswptvXYq5WbMKLTt4dHdcRcYH0mXjLatB9LV76yWLgDSFdwnj84G0qse65jPAtIdyV47QfG2bfqfiSC6+cjNVhCtd0ilCUR/f8vaCKJXyW4Ax6/c5mlg+ptNP5jWEOkF07+LIfEf4AahgywhtKEphIalj5FKAwAXI3dOAQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAVwEAAAAAAAA=eF4t0k8oA2AYx3Et4yBpDrIT+XPZgbWDg5ayRC3JQin/TkorkQNrrbVmB0Ii87/YMMNsQ5MyE0PhNErbYTVxUbNdlnKi9n1Pn56353mf31uvI5mOG6sfGzKB76yaD+r4Pfo8aHPj+gku7WNQnCt8WHWIDxvYtoaxZSyZR7+Yu9jDKTt6j9DqQrMDVbsYdeKKqEfFnhdxf8iKshHRP4DOXvHOITwdxh8tNivRWCPyqfGzXuypwIAEU7+prDlSak+lyFuGt4V4rsKuRlyQo/uL+YME+iMYfUbTK/rCqLrG4jO0eFF+g+Egqv1YcImLcYw9omsbbbNoHkO9BpV1KOvA2k7MaEWeVgzpRB4Fzv3xf1ok1JEisScXS6XYlKZP/4Tj4t8l3KgQ7u7g4CZ2r2LePuYfY3IP311ouMOeK+yzo2kG201YbsG3aVQbMDqB/ZOo28J/aQLFrw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAABgMAAAAAAAA=eF4N01tME2YUB3DiEjKE7FKhDKZWucigcldXQcu/AoLiZFRYAVvR4XR2QztdQAlqGAIiRsaUiBFEO53gDHuQMpEs4/LAstnKA9GtJkbHpSaLy6RR/8UCO9/LLyc5Od8558vn4yPHTCQEJvW8GUIc1ajuJOwjXlu0X+/6nKjsMsV9t5eYeVL1dGiPxKFtVvdnhEffb4wQjzQ6lQW7CQ57RmtLicPe9xp7PyVerdJkunYRFWWG+WDx5bXyvuydEj9qOXSkRPKVttgbOyQ/d8zlNEm9evcVf7FyQGFcZ5T7mKgs205UJeaNthcT3n2WU44i6c/alDFfSMw6u+fixWOLHLd3Goj5nGcHmz8hqk8ExA4VEG/8onZN5xN1LzZfCRf94szb87cRp/c0BNXqibc6Ou/Z8ojmByMNUx8TQe+40oPFC9m+c1m5xNLqyNuHtxLWvoyDXR8RUdOlK51biJsxNVMLxeRS6+XUHKLv4mDxl5uJtLHHge2biOEAn3v2bGJTpqphLotwHNWmx4v6XtNsyUbi/r9VP3+bSRij2r4azCDGS/rV0+mEudU5GSZOj3o6tm2QvfiFFJ/Qydw6TaANxPFKg2MyjVhwq/ykUqz7p2VDllbmirB5K9bLXMax3s51xNstbstfqcRZu0K9UFT6Jk2mpBCt2ryOL9YSiyssRW0a4vJPTYvsHxKRT7vts2uIH5c56uPE5KJnupLVxJ3mAG/TKkL3u7p3IJn4bUGO5XkSkZtqjgkT7x9qmNAnEqabnZdqEoipiZHCnnjiwBKXYjJO3rfA1x4kVp+JrN8YS/iPZOgqVhLn5ktfX1cTKk2N7c8Y4obFesBPXNM1GJ0STQw9eTxu/oDYGupz6WIU4dSrCu+uIPY2ahWzkYR72HQ3VjzurarbESF1V7ehKZw4X9Y/82sYEf6Ds+e/5UT3I8/+5WJKcEi0fhkxkqsZ/0ZF5J80tN9aSvw9UG6YWCL9elreDRK9ibY/MhcTp8xjteXvE8Hfu9OuhxJXHypmHsi/+B9SmpVCAQAAAAAAAAAAgAAAAAAAACADAAAAAAAAzwAAAAAAAAA=eF5jYCAPGNSd/8ScxGP/OiyX70CGpn3Ey8Un4/i17VXW7TlpLq9tP3FxvNRTJW37Ip4jem2q2va+OXt/eWho2/ulmp5+rKVt31DR9YNFV9te7AqbQ4+etr3Jd3/HVn1t+ztVj08LGWjb67mc1YsG0iG79hyqAdKFu/4ydQPpdzeOTp8ApNffk1gwCUhzfSgOnAykE2aEdoFoN/dJs0G0alj7AhB95/bjhSC69pXdIhD9pjsUTFu8iQLTbpoxYDqkEEJ3VozSozSCBgDRVid6AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAPQAAAAAAAAA=eF5jYMAP9Ne9upF88KfdgXPCi9l3Stif2mTvIOUqY29ZosZ/ZreI/dzIt3eTha/bETBmFIyCIQkA78oSEg==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAATwEAAAAAAAA=eF4tzk0oA3AYx3EHWWs5KCvJaW0XJWlOK1qU1MqiuCi1klaUhiiNKLJhkyGvw+RtZTOvS63UlrJwkFg7LLSDmm2EwxyU2vd/+vR7np4XV/IjNqgKV8ejGDnCiT0ccePpBfYtYG4A9Vv4K+Zmx1Hqwf0RHHNi7AQdi2iew+cQqqziD3Fv2o4KP8oc+LSNdxb0mbHHgFVtKGnF134cbUdjBSrL8b4MbTrUlOJOMV4l0lnlUnIiB7vr0J6h/6MgD9SgWoupL/qGF/wUmoKofUT/Pe6e4eU1djhx6RKXxVzyHCfdOHWLK2ks2kH1Gj6aUWlEXwl69djbgu+N+Fcr/m7G4gaU5GHBbyrrkEzsLcT4D/UGOVkTJq/eYNMGdh3jmxMNbjTOYL0LQweo86BtHdWHGIhiZhODVqy0iTud+DCPkWF0mTBfZK8Fv8/wHxrTw/4=AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_84_t_8400.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_84_t_8400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9a93b2cb21bd71970d3f14fcef7b10516e52b9dc
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange_ts_84_t_8400.000000.vtu
@@ -0,0 +1,45 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+ <UnstructuredGrid>
+ <FieldData>
+ <DataArray type="Float64" Name="Calcite" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.0002074413507" offset="0" />
+ <DataArray type="Float64" Name="Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-2.4987393862e-184" RangeMax="0.0002074413507" offset="712" />
+ <DataArray type="Float64" Name="Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="0" RangeMax="0.0001128411771" offset="1512" />
+ <DataArray type="Float64" Name="Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-3.4393213979e-89" RangeMax="0.00011275530434" offset="1756" />
+ <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="21" format="appended" RangeMin="45" RangeMax="103" offset="2036" />
+ <DataArray type="Float64" Name="pe" NumberOfTuples="200" format="appended" RangeMin="-5.4749392229" RangeMax="8.8599764761" offset="2116" />
+ <DataArray type="Float64" Name="phi_Calcite" NumberOfTuples="200" format="appended" RangeMin="-8.8991937337e-24" RangeMax="2.4517244091e-06" offset="3772" />
+ <DataArray type="Float64" Name="phi_Calcite_prev" NumberOfTuples="200" format="appended" RangeMin="-2.953231563e-186" RangeMax="2.4517244091e-06" offset="4620" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)" NumberOfTuples="200" format="appended" RangeMin="-2.2055955777e-106" RangeMax="2.3250707982e-06" offset="5464" />
+ <DataArray type="Float64" Name="phi_Dolomite(dis)_prev" NumberOfTuples="200" format="appended" RangeMin="-7.086707186e-91" RangeMax="2.3233014065e-06" offset="5780" />
+ </FieldData>
+ <Piece NumberOfPoints="101" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="C(4)" format="appended" RangeMin="9.0683858079e-11" RangeMax="0.00014320496112" offset="6096" />
+ <DataArray type="Float64" Name="Ca" format="appended" RangeMin="1.0000000003e-10" RangeMax="0.00022019861114" offset="6824" />
+ <DataArray type="Float64" Name="Cl" format="appended" RangeMin="9.9999999786e-13" RangeMax="0.002" offset="7584" />
+ <DataArray type="Float64" Name="H" format="appended" RangeMin="1.230841961e-10" RangeMax="4.9650778322e-08" offset="8628" />
+ <DataArray type="Float64" Name="Mg" format="appended" RangeMin="9.9999999786e-13" RangeMax="0.001" offset="9384" />
+ <DataArray type="Float64" Name="darcy_velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="10420" />
+ <DataArray type="Float64" Name="pressure" format="appended" RangeMin="98704.543993" RangeMax="100000" offset="10924" />
+ </PointData>
+ <CellData>
+ <DataArray type="Float64" Name="Calcite_avg" format="appended" RangeMin="0" RangeMax="0.0002074413507" offset="12000" />
+ <DataArray type="Float64" Name="Dolomite(dis)_avg" format="appended" RangeMin="0" RangeMax="0.00011259946251" offset="12400" />
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="12560" />
+ <DataArray type="Float64" Name="velocity" format="appended" RangeMin="2.9976852e-06" RangeMax="2.9976852001e-06" offset="12624" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.5" offset="13116" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="14044" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="14344" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="14656" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAA9AEAAAAAAAA=eF5jYBhe4ObL0rWss+XtYfxzwUv6475p2XtySJ/6cV7d/kEXsxunrrb9/0/hLx2Vte3PbV9+h9NA2/4C0yW1s/ra9o6CU9Ok9bTtC19HTpQF8v9wP/SW19G2VzdRLGcB6gs18tp+R0vbvta3UqBAW9v+Rc6twFgg/63e+sb5QLpuf9ncBiAd+vT6nHogrTf7xgsJoDoZIYns50C+z//oj1+B/DeVM7qXA+lVxp3cb4DmG5hbJU4C0q9y4gVuAu3RUf19txFIxx//vHsZ0D0fpp15HQWkI0/uafACuuv0/CIPPiCdvvhExGGQux9mqk4H0hG/WxvEgP5J1NXy+wFS1/BG3gfIV2h4eMECSG/1vLs6B0j7znvGlwikuXffC68H0sVPnPZVAOlHrPFNnUC6aaLA51Ygbf/Y+EQfkC7refC5B0hvO1KYOhFIS2kH1k8A0rd1Ji+eBKSPFcibg2hmhw8mk4H0S57TzCCawfJ1CYjmtS2LB9GnlolPANEBG3jaQXTfDLvZILrgSP90EK0UtGY+iL4/y2IeiG41S1gIol+kX1oAov+eYloEoi00zoHF506wBfPfTdcC0xdlgsD0zUxvCC0eCabPPA0D04ql0WDa8EQUmH5iFwOmmW9BxBWyIHxfXwjdWTFKj9KjNC4aAAgMIy0=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAANwIAAAAAAAA=eF5jYACDBgZUGgbQ+bjEG9BoXOK41MEAunjD1Za/F9d9WdIEE7gyZ1v+4z8XZpaLvdW+79nJoyjyJNt7XZHhz3Ue3/ttdylvy0+0OySsbH/wKKfSNLHbdief1h2YKqRtn92Z+WyGgab99wjtP5F62vZrHJq9Jqlr27efPHBmib62vVRH/PSHQPrY/LQPnbra9v+3ZD+aBFR3fP5hrbXa2vZzvj2c16yjbR/Ydq1rmpa2Pfvasv3iQPG34qs/OQD5Uk4/FIqAtMy9VdvbgDT/yZlby4F0kd72Q+ZAdbprlhv9AfLnrVOzVgaacyLU5tU5oPiOcNZcOaB9i+6ZTjsEFOeVP+PBCbT3vAv/2U1A8WgLRe0LQH6Jf9OTNiD93cVXLBvozqQcVRMjIG0cPIHnNpCObLXS3wSka86GTFAz0LYvfnHUjBtIr9cOCQoD0mc2Suq6AmmXm3ZHioD0g2dhyzKA9E9trVNNQPoKt5JqLZBmeOMd0w2kzYwWx3UAaeHq77v6gbThN0OLPiDNdyP3zkQg/deiVRFEz67suD0JSH/6ea8FRG/TcPCbDKT/MM9SB9HCdz/Xg+ja1eH5IJpRdc8UEK11+3oviN7zSHEuiL6YEzcLRP+ozloAolkuhc8H0VyrViwE0aIdjmA6u0N2EYjeeOE7mN+v5A7mz9ljCaYl8kPB9IdDAWBa2CIKTH/fFAGmmW5Fg+ntaRB6W3QMmF4uDKF7KmJQxGH8UXqUHqUxaQAh5ENQAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAlQAAAAAAAAA=eF5jYKAvKCqL2Hn+wic7GJ/3vMqEdjEpe+3YpFLvVmH7LG+eiwoBsvbK8yY0fbwiY9/4yiDTuETW3jzKceWSKbL2gqlX3ycYy9r7sOzzWJ8ta78q+9q02Sdl7EP3y88255e11/2oKtTdy25/N9hBg2uhjL3k779nvjyfaPdCV+ZWtOpruL2jYBSMglEwCigDABswMWg=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAArwAAAAAAAAA=eF5jYEABDQzYAbHi6HyGCNXqX+5lUbO2z8svVCze2aj188BX4bp/drO+7e651TxxLVujyYOrtVL2V3bxXlK9LGKve6VCRDRQ1r7mVc+tgkcy9hqKHacnZ8vafxG9nfCwT9Y+NKZXikFJ1j4zwDS1LA6ors87Z8JWGftHzyPUFX/L2C/xK3DumMtgP8lhanXbOml7mDt+VbqHPVS4b4fuvlEwCkbBKBgF5AEAPVw+qQ==AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9A1tDTVTbcwtUhLMzEBkgAyrgUDAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAuAQAAAAAAAA=eF7t0ulT00cYB/BikZZDYvLbFZAqvxQPpkYQ5VAut1YxSq3CAMqloICK4lE5C0UJqDVCAIEWClNBLVqRIyKYMiAPEENECUYFLyxyyAxWEEqhHoDd5G/oy776zPfZZ/fZ2dmNumZX4415cHogssTInAtxgmecjocMnDhcbvhXMQ8eL/OpS5QgMH6d4N1fycAT9ajJjKM84ORn9FhUcEHhujVn53YMkvBlnd/I+MQtubAlvoBPhg5EJ5k78eBs3SpP63M82Dwc0Fl6BoPUL+6013qWeM3TkTeFIKitE9hs80AwEMOW9P/Okp3Oio7sBh5MkWGVexdLbjYEVnpc48HTVbV5xndY0n4tfcTwFg92u71/kHWPC6IkxxtGJ3nQvMCkSe8RSwJCW3cWXmfJapGyCvT5RBgY1WB2iSWCwvvqg7UsMcldHhZYzpLQa6byuUqWnHljMTF4imphJc24whJH6Qnd1W9YskZ5Qbe4mCU1mY/aZeks6a6cY3HgCEvmWuWsf6ZgSUGjWdDHUgRiPeXddHrucKLzhVraF7Urp083niXBOmkxenTOpDynNzuLgaq9nalbauj7BhmKh/0ZeMW1LCgIYMlur6EHj3JZ4i9wStWj97FXqox6mhlwlzyLvBPJkph7idNLs1jiuujniYRSBtLfTtczxxl47FK657c22h8imvz2KgMfNgf6/uiHwOly08x6EQNLOuut/X5lQJQzvq47GEFO+xyJQx8D0vw4xbAVgg67JQtLehjQE/gMnlmH4IZcpPz8FQNpZVFjy9UMvCgykjWeR1AjGrdvWYugSDCU7BuA4CuB0RXVGIIZU2UFC6YYMFxp+0mZDIHl/rDGm3EIVtzdahlZjMBmjVur+hIC/ds6qj//RpASsK+6uQGBsw0/d8NLBIFNHgnvJhCo01m5Zw8Cvnxia18TAtOQp9kjbxHMVsnCNj6l+zZNX3X5A4Fyw4LnfjRHp/raCq0wBKWOJvRyMTSq8ZduZhj4zckvha8RHLTNXJ7yBQb7+Rle1iYY7jh0rQyzxyDu8axzt8Eg357/fZM/hosGaTEVjhjyPnXKK9qG4aC4S+FB60mmiaZ21ENd5Z2ZezAEzw/puxiBIb6rWq2g//whiRDDbgzJmzgXHaIw6Bu4ZhyLxbBMUDLw7hCGdod/mu0kGHx22FlVf4chJLW2qCIbQ2rCe862LHrfsVozj58wyJaGKm4nYxg1Gvf2K8dgHh081ZqLAXF27eOfxTD/o1arjjIMI60+KriO4Xm2bNDxPIbj5isOT9ZhWLInLeKyFEOc57vrurfo/MH+SN0bGFzCf2mIoDlPXFNh04Lh6Adv5yEVBr/4MdPqNgwPuJn3z3Vi6B3YkZ79EMORIImeZTeGu5VC9/4nGE56yWJ9BzAYvE57ntRHXVmmXETzloJxoSt16uvFpxDV2dZg0Jwamr82l0NVLHos5FI9czfvNabOthC3ahzJC19nSLXZHxU9i1p/6WSFPrU73NvZgCoVfHZMk3XaXI5oHDs0U6JxoXXsaY1eLf4nNCo5i3/QmHeuUKzVZES7Pk/IpmkdqErR6LHDXVuvmNWrrb+MfqHtL3Vp0yqpmvrf/9B/AbFiF5Q=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAWgIAAAAAAAA=eF6zYACDhldQuhtCM1RD+blQ/moo/Q9VPUMslM8KpfUZIYzZELohACquB+V/h/IPQdWlQyiGz1BxLyj/C1T9MWFRGxCjHSrO8u7gctNQvmkMpu9asmrrOBgEq1Iu5OeuVdnNUMP99ae8aKSfodyL7XZT1J91hLao7tYLlKxnSj1id9Guz4+b+ZDdpxWqOjF9R+wE5yfM+9Z2xM42sr1z4YQjdk6a5XOLgfS8shlt2v1H7Hz8z599B6SPSW83XN17xK5U43RnBFDfrqvVhtt7jtjNvP5iqRpQfLte55UwIL/h46snHUB6vhJ3eyGI3vTlSAGQzvHS3PkESHek7PPcD6TfMp9/2Q3Udz78/5wAIL264pyAG9Dc1kcydkxAerXIm4ifQDpH1zRlO5BOdtrxagrQHXr84ZOigfTdf/tOCQDdOV/mku1jIL/oyh2jHCCfQ+LRYS8gPdfnp8RqIP3k6xbHKUD6oG3OistAeqZPLe9RIF1zfpfbGyBdn/2j9jGQbu0KT/sFpGd8Ccj4CqSDfz8UY5p4xG6Lw8TU/0C+fcYzYTYg/7XOio0sQFrxtSwfB5CWN+dJYAfSLvUH/oH4b59vXwOi1WdvyeEE6VfucAfRRjsl54Jo5fP3ekD0renLN4LoG90XVoHo2aLiu0D09Te+20F00eKCvSB6bkX4HhAdG7lrH4h+I+oPpp2CTfeD6EdNvGD66vFEMG2eFwymVwcVgulSxxwwnTKrDEzbfywB09YpFWA6+2E5mL7xB8Lv3AChu2UrwfR+dghdpDFKj9KjNC4aAK1tnEo=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAWAIAAAAAAAA=eF6zYACDhldQuhtCM1RD+blQ/moo/Q9VPUMslM8KpfUZIYzZELohACquB+V/h/IPQdWlQyiGz1BxLyj/C1T9Oz3tbyc22ze1Q8W73GO+CG3OnDmBJbzpiYIyz4n5GbZ7nXgMd9xh4bC8GawcJ79mu/WdPXaLfOMvJdaW2E0wWqQbVXXEzneR5eOjvIftvvGp3Y3qP2IXu4H1t2XXEbue/nPioROO2L0pzhRuBtLcQcH1e/uO2K0SlnzeClS3rGnp3ZDeI3bbZGxEvgNp33ussk09R+xaFFxuPAbSV8xmxTkC6dP6HyamAeknAltVyoD0nzd9p7KA9KwXOvx/gHTZ3wbG80A663eR9jqgOXFdE1LygLSm2ySdHKB9G2Lfb1ED0nF/Ii7JAu1d5Xs6/Q6Qf2G7sNMWIP+XL9PVGiCd+S6+QBPoTiOWtsx/QP6hr1MU6oB8A6ZZ5vFAeudnXfvtQPrOS73+RUA6PP9L3G0gvTb/Xdo5IN17hfnxByD97ymDzUsgrdC7nfEvkA78s2HPDyB9RSW2lnniEbuCZeZWjED6f0z7aTYgvd/z5zJWIO2gI5fBAaRfR+29xQ6kH7uKaHICaTnp9ucgcTWBW/Ug/hrpeUkgWt2vaCmIzn7/cAaIdp2itRVEKwt+WQ+ixZe17AbRUm25O0H0qSqGfSBaY6fBXhB9R+0jmB+jshhMG/d67gfRU1h1wfQ/lkwwPfFRHJi+fLUYTPP4FoDpjJJyMC3iVwamOzdUgOmQSAi9VKgSTCvfgfBLNCpRxGH8UXqUHqUxaQAEH5kTAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAzAAAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNIJkoHTpx9aNG3RSYjAypJ+DYYLM+xStPS2f9JnWHWub25ti/FBOwb5UzZ/9u60+7zrhM+D4MN2sSEOjcuZDtupOSoprcw+bFf9MtmAt/Gw3Q2//OmvHQ7bzd0ybw132mE7mzdRWnq/D9mFBL04LqwJFJfNMzfhX293OWtVY96DQ3arJ4cXPp1jbpf6feoJ0WcT7RhGwSgYBaNgFFAFAAAi7Up1AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAygAAAAAAAAA=eF5TY2QAg0cQuuEChMuwHUpvhopHQtXtgdLBEKphC5Q/BaouGcqXgNIJBZN91hUxzmLYGHvs3PGExkaWv+f/SS2xOyC2s/y9hPnaVTvyfy8sP2hXMDF+oW3Zbjs/L/noxpDDdkv7V57p5T5sx23e3PEz9bCd4otFJQ61h+0ab6XZlJketus4ue5qcOxhO3dOiZff3h6y+8G7VvSq9GG7JfMjrLVeLrXb8VUsud3hkB0DFIjETs/1b2uB80fBKBgFo2AUUAYAjnhKiQ==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAAQIAAAAAAAA=eF6LaE/T8H97w3bVrJqtguIK9rOeGdeZuW2ze2KQyPJc/YLdH9+j++dPemiXVqvXFsz10S5q5STl9B3/7NaJ1mhs7+K0v1G5aaFTj7C94p2L9ppHlew/h1S+ThJTtJ9ybvVF/RYF+5n5+87MPSdvz7WO/866Cjl7K+vpr5XiZO1Dc007ktVl7S88uvYlQVXOviymTNuuVs5+b79I6vu7cvZt1q7ZgRby9qkKeiF27fL2stZcBq9PyNs//c55+vcfeXvbA4cSPgoq2Et0LHZOlFewX1eVE7VUXcFe4WTsqgPaCvZGcydvOKSjYC90QnzxaiD/yiWRa+GaCvZl9fcMS9QU7FkWvJ1craxgn+rt4LxbQcHeTVtkQ7Scgv2yGqb1WTIK9v9c7LifSinYd1xrsdsnqWDvVbjyy1MJBfsDOqY93kB61+TtLA+A4SVRn3F1GpDW1fEWTQPSTpI9zm5A+sdUj6eGQFrndIKkKpBuOSlrIgukLycvypIA0o4yOhNEgbTS3XePhYG0i97zN0JAWk1ocQKIbtppJgSiLezZnoPihenEtPMgOkKG7zCIjgst2gOiNeYw7QTRc2Ptt4PorWzR20B06Y8X4Ph03LkBTMeETgXT7cYdYPrbv0YwfU60HkybbqkF05cEIfTai5D04L4BQrPPgNCwdDLSaAARryo4AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAGAIAAAAAAAA=eF7z/P7z5v2a27arZtVsFRRXsP9wa3HChnMb7Q5uC9u8cvlJu7iYrMWdvbftPFac/Ohl9spuDvsJq9PK3+2uZLYHvmJitv/1WH3PAxVe+/x1rjdFn0va+xW230hNlrbvi22YHs0ob5/GW/7Q7rOivcMq7QvS3ar2v9IP+bdJadkn8u5dtfiOjn1F82Lm6Ula9sq+L3Ru62jatyopbU//qG4f9mfmUcHjavZ2N21mp29Vta9x4DIJ269iP+PjPIkHT5Xt57zKrXmnoWy/MUKHo71dyd5l2vc9P/4r2r86c/0O3wSgfUvP/3U0VbTPuHTrwaM3CvYSc7Zo221RsNfZKR+6s0PBfm2nnsLLdAX7F55p+xb6K9hfSTg1icFOwd7B+sXZfEMFe7YisXlLNBTsT4lvMpyppGD/a50XX4ucgv2qKmHTc9IK9jmnvdw4pRTs38ROLXkmoWDfJpO4Ww9Id8XlNOwBhmNTdOHUMiDNumRmoQ+QZvLYkmsMpPdsWn5GBUgnsfXFyQBpsTV3pogD6Ve1c5+JAOk3S49mCwNpntvNc4WA9Dbe/Xogeqk63xdQ/DQZxV0F0Sc+tR8F0VF1ZXtBdEpP+k4QvSgqeTuIbpRP2waiA0Jfg+NVqHMDmDZbNwVM/1reDqa9FjaCaYZndWDaamEtmJ7DBaHXXoSki5rtELpgGYQWXgmh7y2E0LD0M1xpAA08O4M=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAA7gIAAAAAAAA=eF7btezjJd+kBPspJj0NM6cX2q7Kf8cnGptgf2/yse7e4AT7ayeX/NzmnGBfe4DznaR2gr3flpTka2wJ9orMfT4dy+PtXfTqJX1k4+0XvGXyu5MRZ2+SeF9HrDXWPv8gF3tVXow957IPX3cqR9vLCdxU11oWac/Vfvrct5/h9h71izNF1cLsl95nEnxnHGJ/Iegy0xG9IPsfT00c3ykF2Be9WLf+raWP/b9VBVnlgR72isIVi4IXutjfE7xi6mvlaM+uJXjyfICN/dbvH156PDWz35+ysyjfzMhelS89tLJBx37T4enb975XtV/s7HfQYraCfaXz74SnuyTsFVaHHTssI2R/Xvdckpgkl317HAOHej2T/bR9Vq6tp3/aKZWY3Vuw6aOd30qWwNNPn9sZx3V8SD103y6oTsXt+dErdoYxsiorDpyy+/aiAkgcsNslruNWqbfVLmPJ9C19vKvsvjBMOpAxcb7dAxOWBwI3JtnxdYZ/nKXdZvctlG/DrK2Vdp0by29xpOTave9YItI6IcHO0PbZifuZQXbv2K65bzJwsVvDH1wZft3ULugg16kNAep2k6abvGaeKGH3Lqr+4/s+Tjuefdqq/S6/bFPEeE0YNr6y1cg1Mrx965btgT35a5dfOWVbff6jZPvRXbZHnjuer/21yjYiTbRdj3OOLWvt06JVov22z74e6fDQb7Zdc0xa/qtWte2MFP4ZD76U2j6aLlLk2V9s2/R8/S3riiJbHqf6vzt+Fdo6TFNRenSo0Dbnr3K836pC2wlGV2xmzS+03ZLoavt3VqGt37V12awzC21Ft3xYajCj0LZg1sW114Hp6MD+59c2AelX7j67FwPpb/66M+YC6VN5j+RnA+nDvgeUZwHpOYpyB0Dp7swKiZng9Ff2qRNEJzyVbQHRJ8yfNILooE08YPrb+d3gdPpl01wwbTB1EpgWaeoD08LsvWC6pw6SnsV8IfQ9awhdZQ6hraDpfQqRNADRzXb3AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAFQIAAAAAAAA=eF7TelVf+nl5lt3cxfIPF796YBto92Sq+jdzO6vLPds5/2jZvSvTPhDDKG234FrENsUtvHY675a9lz/CaPfk955zs22+23per070VvtgW1Tw93Wl9RPbCnHPl9e5ntluORliuz3mme0kjTf9fJuf2SokstbyCj63bd5uyDXV/bnt0Y9HHz///cw2QnN5zi6xZ7avbh2pjmh7ats+MVHlw4cntjI9ojIpCU9sL36wu3XpxmNbvyib3qqYx7abXhQ6bnvzyLaxZL6xc9cjW3H7zzVepo9s1+75aZ/39qHt5etH/B9ufGh7peRkTUvjQ9uTm4s1WaMf2i6TXsr/y+ah7eLT4l9/qz60XWebdG2HyEPbudHGHxk4HtrWcUx/d/bfA9sL5sZCrj8e2NZd01299fMD24zUZ4JTPzywnVL3Mbzy3QPbe87l/cJvH9jmXb3sofUGSLNLLZvy+oHtJNdjhmZAOnKfStNvYDhaMurMuQ6ko2pS/hwE0rOOvFbfCqT7DPlXrAXS+1RtwlYC6Z6LJY+WAWm/+XLNS4H0GX67g0uA9KNFu8JBtBu3mDyIvrNClBFEM+imfwTFU7J1wEsQHdVr+BREf5CWfwyiH5wQfwSiH9VMA8enyKVwMD3vjSWYnrBTG0w3KaiC6SweJTBdtkgBTN9ig9BF5yDp4dhsCF2yGkJ/KIDQRn0Qes7ikUEDAIurymU=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAA5wIAAAAAAAA=eF77vPzjJd+kAPspJj0NM6cX2s7qXfrBKC7APuSw+7yHYQH2J66fWe/nHWDPpN6w7JpVgH0Lv91lDdUA+21mWr1MHAH2C86bJbXt8beXYmV+7HzX3/5B6lUdpzg/+9Omyzqi1vvYv4h2+753gZf9vVmZ/JVCHvZLPNNnSMq72t85JfrhxmQn+/1Vl/K+xDvaV7AG/cpMs7dPK0lvWMJiY6+xnUWjcomFfYWeudHkCFN7s7ZrIY+NjezVcuIsLXr17Ou0nNddstS09zgboKZcq2I/86GCoEycgr3Dt4xPd7ZL2lscVZbff0/Yfr3s+v4phnz2YRxr+Xz62exnPJtgu34yg/25yCjOVV+/2+lcLzb/zP7RbhWHubb5xed2GUzejP5JD+yM4n9K3NS9Zhe22Dt7ZfdZu/NdLCqJcYft2iNLHEP+77A7d+2plO3h9XZ1IZONbt9bbDd3ctlul5oZdnoTmZZW7ui1m2k0cVHI0Ua7PXsavj38UGLXrhpzXygow+7QN4e419ej7aK4PgrsyPW3S3ywd+urWke7U0yPbyWbGNsZWPxP8TurYveral21p42YnU0L5422Mna7Bu/fYowlP2zXpQm8UTN8aZui41xXs/Gm7azMPGUXllO2l+5Nk4313WXLp1xnsu/aKluDHVXMMcZzbG/t5JCWau23/frnl3RDS7OtHo/MVr4Z1bZSH6MzKkzKbKce5nPQ21RsW/Pnp1N/V5HtovXX37KyFNm6LJV4WHG80JZnIYN2y5pC2y0XnrIsXVBoe6TWennF7EJbo5yeNtuZhbaB9q/3ec4otC1ne77nETD98Lw/vmQ7kPZuzTm1FEj/O2BePg9I36qVqp0NpEsuOHjNAtJrHl47D0pv6/8XzgfRu15X94Do3yvlWkH01DMeTSB620rjRhA988QpcPrck74YTC/cPQVMewT0g+n5ob1g2m4ONB1HQOg1jhB6mwWENjSF0LD0TiwNAOxFb1U=AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAAVwEAAAAAAAA=eF4t0k8oA2AYx3Et4yBpDrIT+XPZgbWDg5ayRC3JQin/TkorkQNrrbVmB0Ii87/YMMNsQ5MyE0PhNErbYTVxUbNdlnKi9n1Pn56353mf31uvI5mOG6sfGzKB76yaD+r4Pfo8aHPj+gku7WNQnCt8WHWIDxvYtoaxZSyZR7+Yu9jDKTt6j9DqQrMDVbsYdeKKqEfFnhdxf8iKshHRP4DOXvHOITwdxh8tNivRWCPyqfGzXuypwIAEU7+prDlSak+lyFuGt4V4rsKuRlyQo/uL+YME+iMYfUbTK/rCqLrG4jO0eFF+g+Egqv1YcImLcYw9omsbbbNoHkO9BpV1KOvA2k7MaEWeVgzpRB4Fzv3xf1ok1JEisScXS6XYlKZP/4Tj4t8l3KgQ7u7g4CZ2r2LePuYfY3IP311ouMOeK+yzo2kG201YbsG3aVQbMDqB/ZOo28J/aQLFrw==AQAAAAAAAAAAgAAAAAAAACgDAAAAAAAABgMAAAAAAAA=eF4N01tME2YUB3DiEjKE7FKhDKZWucigcldXQcu/AoLiZFRYAVvR4XR2QztdQAlqGAIiRsaUiBFEO53gDHuQMpEs4/LAstnKA9GtJkbHpSaLy6RR/8UCO9/LLyc5Od8558vn4yPHTCQEJvW8GUIc1ajuJOwjXlu0X+/6nKjsMsV9t5eYeVL1dGiPxKFtVvdnhEffb4wQjzQ6lQW7CQ57RmtLicPe9xp7PyVerdJkunYRFWWG+WDx5bXyvuydEj9qOXSkRPKVttgbOyQ/d8zlNEm9evcVf7FyQGFcZ5T7mKgs205UJeaNthcT3n2WU44i6c/alDFfSMw6u+fixWOLHLd3Goj5nGcHmz8hqk8ExA4VEG/8onZN5xN1LzZfCRf94szb87cRp/c0BNXqibc6Ou/Z8ojmByMNUx8TQe+40oPFC9m+c1m5xNLqyNuHtxLWvoyDXR8RUdOlK51biJsxNVMLxeRS6+XUHKLv4mDxl5uJtLHHge2biOEAn3v2bGJTpqphLotwHNWmx4v6XtNsyUbi/r9VP3+bSRij2r4azCDGS/rV0+mEudU5GSZOj3o6tm2QvfiFFJ/Qydw6TaANxPFKg2MyjVhwq/ykUqz7p2VDllbmirB5K9bLXMax3s51xNstbstfqcRZu0K9UFT6Jk2mpBCt2ryOL9YSiyssRW0a4vJPTYvsHxKRT7vts2uIH5c56uPE5KJnupLVxJ3mAG/TKkL3u7p3IJn4bUGO5XkSkZtqjgkT7x9qmNAnEqabnZdqEoipiZHCnnjiwBKXYjJO3rfA1x4kVp+JrN8YS/iPZOgqVhLn5ktfX1cTKk2N7c8Y4obFesBPXNM1GJ0STQw9eTxu/oDYGupz6WIU4dSrCu+uIPY2ahWzkYR72HQ3VjzurarbESF1V7ehKZw4X9Y/82sYEf6Ds+e/5UT3I8/+5WJKcEi0fhkxkqsZ/0ZF5J80tN9aSvw9UG6YWCL9elreDRK9ibY/MhcTp8xjteXvE8Hfu9OuhxJXHypmHsi/+B9SmpVCAQAAAAAAAAAAgAAAAAAAACADAAAAAAAACwEAAAAAAAA=eF5jYKAPuPmydC3rbH57nnX312gna9p/2Ld2uqqGtr3/3T2Nr/W17SPrvlfP0dO2F1JofDhRR9v+5f6HyWra2vavzh7JrNPSts8KUpdtANJLw35oMgDF59hU7TsHpHO3JprtB6rn9tQtX6urba+h/+ByA9CcSVump+gAzQ1xkrRfC6RjRCLS2Q207csY0iucgDTjl8S7aUB6/3kB7hogven09l3tQPqWOde9XiD9ojzx0wQgfcEvMncSSN+yWJnJQJptt0EWiO7VntQNouv10meC6IayT/NAtJm/5EIQPavyBZh2EzFZBKJDD/uB6Uv/wsH0SrloMF3CEwOmt0VD6M6KUXow0QDW5BB6AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAVgAAAAAAAAA=eF5jYCAPFJVF7Dx/4ZGdzLTDTS/PitkfKzmn2Skkax+jZMmZ0yJrv26Lp0yHv6x9g+jvGZlvZOytlq288dxYxP6gbEvD7C237QiZPwpGwWAAAFtSHNk=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAATwEAAAAAAAA=eF4tzk0oA3AYx3EHWWs5KCvJaW0XJWlOK1qU1MqiuCi1klaUhiiNKLJhkyGvw+RtZTOvS63UlrJwkFg7LLSDmm2EwxyU2vd/+vR7np4XV/IjNqgKV8ejGDnCiT0ccePpBfYtYG4A9Vv4K+Zmx1Hqwf0RHHNi7AQdi2iew+cQqqziD3Fv2o4KP8oc+LSNdxb0mbHHgFVtKGnF134cbUdjBSrL8b4MbTrUlOJOMV4l0lnlUnIiB7vr0J6h/6MgD9SgWoupL/qGF/wUmoKofUT/Pe6e4eU1djhx6RKXxVzyHCfdOHWLK2ks2kH1Gj6aUWlEXwl69djbgu+N+Fcr/m7G4gaU5GHBbyrrkEzsLcT4D/UGOVkTJq/eYNMGdh3jmxMNbjTOYL0LQweo86BtHdWHGIhiZhODVqy0iTud+DCPkWF0mTBfZK8Fv8/wHxrTw/4=AQAAAAAAAAAAgAAAAAAAAHgJAAAAAAAAlwIAAAAAAAA=eF591EtI1FEcxXErJGcKEdOgBBdJYNEDKqQgqcBXRo8BpWBc5MJoIWYgLaYIIyWL1OxhRRjVRBaoWPZQ0TRn0jQnH6njK2sEZWIWQjELF2n3uOwL3eWH4X/vPb9zJyTkf+vHvn+l/2dtqs9RAF9Y8iJ43lRRwNZYBs/26/e34cGT97Xga8b1nYfwss2NscWBx/Dupe8/hYdVO/MTrdXwtBxt+wJeElfeEYyvgV/26Tx18NkjiVavsx6esEXnfMkcLDvNaoCXzur8r+HNrvis8o438JlHutc7eORFs62tEZ5o132b4Lm7o+pSfc3we9HKoQXu/mXxOvNb4XN9yuc9PKZ20aw2eFyycmuHNzU4zIE+wI9uUJ70mRtBBQQ/v7jD5OyCR+SZ2Lz0Z5PK3w3fe8hcK/YjfLBJ29JPx58yqxO+UKl5dcFvhfrMhelbCzTHT3DXtN0Mhp5p03y74YG2EbNxD7xwm+ZOX19lrlv8GV6/Sn3ohac4erXgE371xAM/e9zEHPUFvrJT/aFX7TLbZvXB553qFX3smnLrh69LUt/oJ/4oT3rlW/VwAD5yRjnT125SP+mZ08p/EH7ngXpLH87QXOjR4frf+ArP6NK86DcL1XP60B7NcQge+Vv9p9tqNF56RY7exTB8INbEnE+PGNV7oR+rUB9G4GXpekf0/hX6e6OHt+p9eeGHz6k/9NLtend0j1+9oq9+ol6NwtPt6hX9epS5bpDe41GvxuCWK+oV/eB+9Yp+dV69God3v1Kv6GG56hU9baN6NQEvmVKv6J13FTM91KZeTcKTrWZbB73YpV7R3RfUq2/w5QnqFT1pTr2iFz1Xr6bgrmz1ir4sRr2iHxhSr77DL5WqV/T2FPWK/hfn/pddAQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAvgAAAAAAAAA=eF5dxcciAgAAANAohOyZUhq2bLJ3yB4hiaz//waHnN67vECgrpGDHOImbuYWDnMrt3E7R7iDO7mLu7mHe7mP+3mAB3mIhznKIxzjOI9ygpM8xilOc4azPM4TPMlTPM0zPMtznON5XuBFXuJlXuFVXuN1zvMGb/IWb/MO7/Ie7/MBH/IRH/MJF/iUz/ici3zBl3zF13zDt3zH9/zAj1ziJ37mFy7zK1f4jav8zh9c40/+4m/+4V9u+P8PcvwndA==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxwAAAAAAAAA=eF4txdFGAwAAAMAkmUwymUySSZIkmWQySZJJkiQzmSTJJEmSTJKZZJJJkkmSmSRJkkkm6bN62N3LNTc1tLjVAbc56HZ3OOROh93liLvd415H3ed+D3jQQx72iEcd85jHHfeEE570lKc941knPed5L3jRS172ilNOe9UZr3ndG970lrPe9o53ved9H/jQR8752Cc+dd4Fn/ncRV/40iVf+do3vnXZd773gx9dcdVPfvaLX/3md3/40zV/+dt1//jXf/4H8McndQ==AQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF5jZqY9AAA7kgEt
+ </AppendedData>
+</VTKFile>