diff --git a/Applications/ApplicationsLib/ProjectData.cpp b/Applications/ApplicationsLib/ProjectData.cpp
index 7189bc26a34b93e2df8434e99bfa1ddec3c444d6..a6f3693eebf2c5a6063319f285208cd00dd6c0b9 100644
--- a/Applications/ApplicationsLib/ProjectData.cpp
+++ b/Applications/ApplicationsLib/ProjectData.cpp
@@ -1043,8 +1043,8 @@ void ProjectData::parseChemicalSystem(
component_transport_process->getProcessIDToComponentNameMap();
_chemical_system = ChemistryLib::createPhreeqcIO(
- _mesh_vec[0]->getNumberOfBaseNodes(),
- process_id_to_component_name_map, *config, output_directory);
+ *_mesh_vec[0], process_id_to_component_name_map, *config,
+ output_directory);
}
else
{
diff --git a/ChemistryLib/CreatePhreeqcIO.cpp b/ChemistryLib/CreatePhreeqcIO.cpp
index 61b542d0f3e6dcf0ba41713b13a47cafca5905d5..cf95af42b8914180043f5c672f283fce901ee7e8 100644
--- a/ChemistryLib/CreatePhreeqcIO.cpp
+++ b/ChemistryLib/CreatePhreeqcIO.cpp
@@ -11,6 +11,7 @@
#include "BaseLib/Error.h"
#include "CreateOutput.h"
#include "CreatePhreeqcIO.h"
+#include "MeshLib/Mesh.h"
#include "PhreeqcIO.h"
#include "PhreeqcIOData/AqueousSolution.h"
#include "PhreeqcIOData/CreateAqueousSolution.h"
@@ -24,13 +25,13 @@
namespace ChemistryLib
{
std::unique_ptr<PhreeqcIO> createPhreeqcIO(
- std::size_t const num_nodes,
+ MeshLib::Mesh const& mesh,
std::vector<std::pair<int, std::string>> const&
process_id_to_component_name_map,
BaseLib::ConfigTree const& config,
std::string const& output_directory)
{
- auto const& num_chemical_systems = num_nodes;
+ auto const num_chemical_systems = mesh.getNumberOfBaseNodes();
// database
//! \ogs_file_param{prj__chemical_system__database}
@@ -88,18 +89,14 @@ std::unique_ptr<PhreeqcIO> createPhreeqcIO(
num_chemical_systems, aqueous_solution_per_chem_sys);
// equilibrium phases
- auto const equilibrium_phases_per_chem_sys = createEquilibriumPhases(
+ auto equilibrium_phases = createEquilibriumPhases(
//! \ogs_file_param{prj__chemical_system__equilibrium_phases}
- config.getConfigSubtreeOptional("equilibrium_phases"));
- std::vector<std::vector<EquilibriumPhase>> equilibrium_phases(
- num_chemical_systems, equilibrium_phases_per_chem_sys);
+ config.getConfigSubtreeOptional("equilibrium_phases"), mesh);
// kinetic reactants
- auto const kinetic_reactants_per_chem_sys = createKineticReactants(
+ auto kinetic_reactants = createKineticReactants(
//! \ogs_file_param{prj__chemical_system__kinetic_reactants}
- config.getConfigSubtreeOptional("kinetic_reactants"));
- std::vector<std::vector<KineticReactant>> kinetic_reactants(
- num_chemical_systems, kinetic_reactants_per_chem_sys);
+ config.getConfigSubtreeOptional("kinetic_reactants"), mesh);
// rates
auto reaction_rates = createReactionRates(
@@ -110,9 +107,8 @@ std::unique_ptr<PhreeqcIO> createPhreeqcIO(
auto const project_file_name = BaseLib::joinPaths(
output_directory,
BaseLib::extractBaseNameWithoutExtension(config.getProjectFileName()));
- auto output =
- createOutput(components_per_chem_sys, equilibrium_phases_per_chem_sys,
- kinetic_reactants_per_chem_sys, project_file_name);
+ auto output = createOutput(components_per_chem_sys, equilibrium_phases,
+ kinetic_reactants, project_file_name);
return std::make_unique<PhreeqcIO>(
project_file_name, std::move(path_to_database),
diff --git a/ChemistryLib/CreatePhreeqcIO.h b/ChemistryLib/CreatePhreeqcIO.h
index 46652a3d1bf48296d30c58bcd5413c5dcb346eb1..addc2bfabef1468684a11d39060efbe99cea6823 100644
--- a/ChemistryLib/CreatePhreeqcIO.h
+++ b/ChemistryLib/CreatePhreeqcIO.h
@@ -13,6 +13,11 @@
#include "PhreeqcIO.h"
+namespace MeshLib
+{
+class Mesh;
+}
+
namespace BaseLib
{
class ConfigTree;
@@ -26,7 +31,7 @@ class Process;
namespace ChemistryLib
{
std::unique_ptr<PhreeqcIO> createPhreeqcIO(
- std::size_t const num_nodes,
+ MeshLib::Mesh const& mesh,
std::vector<std::pair<int, std::string>> const&
process_id_to_component_name_map,
BaseLib::ConfigTree const& config,
diff --git a/ChemistryLib/PhreeqcIO.cpp b/ChemistryLib/PhreeqcIO.cpp
index 430fda2dc851686d06fe3c4ecccdca80aac1b4ce..5f3db46374ec52f11f893f5863295c0ac46303e9 100644
--- a/ChemistryLib/PhreeqcIO.cpp
+++ b/ChemistryLib/PhreeqcIO.cpp
@@ -31,16 +31,15 @@ std::ostream& operator<<(std::ostream& os,
}
} // namespace
-PhreeqcIO::PhreeqcIO(
- std::string const& project_file_name,
- std::string&& database,
- std::vector<AqueousSolution>&& aqueous_solutions,
- std::vector<std::vector<EquilibriumPhase>>&& equilibrium_phases,
- std::vector<std::vector<KineticReactant>>&& kinetic_reactants,
- std::vector<ReactionRate>&& reaction_rates,
- std::unique_ptr<Output>&& output,
- std::vector<std::pair<int, std::string>> const&
- process_id_to_component_name_map)
+PhreeqcIO::PhreeqcIO(std::string const& project_file_name,
+ std::string&& database,
+ std::vector<AqueousSolution>&& aqueous_solutions,
+ std::vector<EquilibriumPhase>&& equilibrium_phases,
+ std::vector<KineticReactant>&& kinetic_reactants,
+ std::vector<ReactionRate>&& reaction_rates,
+ std::unique_ptr<Output>&& output,
+ std::vector<std::pair<int, std::string>> const&
+ process_id_to_component_name_map)
: _phreeqc_input_file(project_file_name + "_phreeqc.inp"),
_database(std::move(database)),
_aqueous_solutions(std::move(aqueous_solutions)),
@@ -213,20 +212,24 @@ std::ostream& operator<<(std::ostream& os, PhreeqcIO const& phreeqc_io)
os << "SOLUTION " << chemical_system_id + 1 << "\n";
os << aqueous_solution << "\n";
- auto const& equilibrium_phases =
- phreeqc_io._equilibrium_phases[chemical_system_id];
+ auto const& equilibrium_phases = phreeqc_io._equilibrium_phases;
if (!equilibrium_phases.empty())
{
os << "EQUILIBRIUM_PHASES " << chemical_system_id + 1 << "\n";
- os << equilibrium_phases << "\n";
+ for (auto const& equilibrium_phase : equilibrium_phases)
+ {
+ equilibrium_phase.print(os, chemical_system_id);
+ }
}
- auto const& kinetic_reactants =
- phreeqc_io._kinetic_reactants[chemical_system_id];
+ auto const& kinetic_reactants = phreeqc_io._kinetic_reactants;
if (!kinetic_reactants.empty())
{
os << "KINETICS " << chemical_system_id + 1 << "\n";
- os << kinetic_reactants;
+ for (auto const& kinetic_reactant : kinetic_reactants)
+ {
+ kinetic_reactant.print(os, chemical_system_id);
+ }
os << "-steps " << phreeqc_io._dt << "\n" << "\n";
}
@@ -342,10 +345,8 @@ std::istream& operator>>(std::istream& in, PhreeqcIO& phreeqc_io)
auto& aqueous_solution =
phreeqc_io._aqueous_solutions[chemical_system_id];
auto& components = aqueous_solution.components;
- auto& equilibrium_phases =
- phreeqc_io._equilibrium_phases[chemical_system_id];
- auto& kinetic_reactants =
- phreeqc_io._kinetic_reactants[chemical_system_id];
+ auto& equilibrium_phases = phreeqc_io._equilibrium_phases;
+ auto& kinetic_reactants = phreeqc_io._kinetic_reactants;
for (int item_id = 0; item_id < static_cast<int>(accepted_items.size());
++item_id)
{
@@ -387,7 +388,8 @@ std::istream& operator>>(std::istream& in, PhreeqcIO& phreeqc_io)
compare_by_name,
"Could not find equilibrium phase '" + item_name +
"'.");
- equilibrium_phase.amount = accepted_items[item_id];
+ (*equilibrium_phase.amount)[chemical_system_id] =
+ accepted_items[item_id];
break;
}
case ItemType::KineticReactant:
@@ -397,7 +399,8 @@ std::istream& operator>>(std::istream& in, PhreeqcIO& phreeqc_io)
kinetic_reactants.begin(), kinetic_reactants.end(),
compare_by_name,
"Could not find kinetic reactant '" + item_name + "'.");
- kinetic_reactant.amount = accepted_items[item_id];
+ (*kinetic_reactant.amount)[chemical_system_id] =
+ accepted_items[item_id];
break;
}
}
diff --git a/ChemistryLib/PhreeqcIO.h b/ChemistryLib/PhreeqcIO.h
index 5ec9dbb5598d008fa2d6547eebecdd03a4e2f1ff..32165268bcc6150e0911835e13f013b3db574a7a 100644
--- a/ChemistryLib/PhreeqcIO.h
+++ b/ChemistryLib/PhreeqcIO.h
@@ -32,8 +32,8 @@ public:
PhreeqcIO(std::string const& project_file_name,
std::string&& database,
std::vector<AqueousSolution>&& aqueous_solutions,
- std::vector<std::vector<EquilibriumPhase>>&& equilibrium_phases,
- std::vector<std::vector<KineticReactant>>&& kinetic_reactants,
+ std::vector<EquilibriumPhase>&& equilibrium_phases,
+ std::vector<KineticReactant>&& kinetic_reactants,
std::vector<ReactionRate>&& reaction_rates,
std::unique_ptr<Output>&& output,
std::vector<std::pair<int, std::string>> const&
@@ -64,8 +64,8 @@ public:
private:
std::string const _database;
std::vector<AqueousSolution> _aqueous_solutions;
- std::vector<std::vector<EquilibriumPhase>> _equilibrium_phases;
- std::vector<std::vector<KineticReactant>> _kinetic_reactants;
+ std::vector<EquilibriumPhase> _equilibrium_phases;
+ std::vector<KineticReactant> _kinetic_reactants;
std::vector<ReactionRate> const _reaction_rates;
std::unique_ptr<Output> const _output;
std::vector<std::pair<int, std::string>> const&
diff --git a/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.cpp b/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.cpp
index 1f0f031a2eae100de2912fa82f5128f5ea7ff208..44a6b7f77d4b8702a436668e9c17e58497dad84d 100644
--- a/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.cpp
+++ b/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.cpp
@@ -7,13 +7,18 @@
*
*/
+#include <boost/optional/optional.hpp>
+
#include "BaseLib/ConfigTree.h"
+#include "CreateEquilibriumPhase.h"
#include "EquilibriumPhase.h"
+#include "MeshLib/Mesh.h"
namespace ChemistryLib
{
std::vector<EquilibriumPhase> createEquilibriumPhases(
- boost::optional<BaseLib::ConfigTree> const& config)
+ boost::optional<BaseLib::ConfigTree> const& config,
+ MeshLib::Mesh const& mesh)
{
if (!config)
{
@@ -40,8 +45,15 @@ std::vector<EquilibriumPhase> createEquilibriumPhases(
equilibrium_phase_config.getConfigParameter<double>(
"saturation_index");
+ auto amount = MeshLib::getOrCreateMeshProperty<double>(
+ const_cast<MeshLib::Mesh&>(mesh),
+ name,
+ MeshLib::MeshItemType::Node,
+ 1);
+ std::fill(amount->begin(), amount->end(), initial_amount);
+
equilibrium_phases.emplace_back(
- std::move(name), initial_amount, saturation_index);
+ std::move(name), amount, saturation_index);
}
return equilibrium_phases;
diff --git a/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.h b/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.h
index 36e6b087f95a1af1b678923e5158f24c7b92543c..2e7451f00b0d100df62ea5470d94e2f6752d2230 100644
--- a/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.h
+++ b/ChemistryLib/PhreeqcIOData/CreateEquilibriumPhase.h
@@ -9,7 +9,7 @@
#pragma once
-#include <boost/optional/optional.hpp>
+#include <boost/optional/optional_fwd.hpp>
#include <vector>
namespace BaseLib
@@ -17,6 +17,11 @@ namespace BaseLib
class ConfigTree;
}
+namespace MeshLib
+{
+class Mesh;
+}
+
namespace ChemistryLib
{
struct EquilibriumPhase;
@@ -25,5 +30,6 @@ struct EquilibriumPhase;
namespace ChemistryLib
{
std::vector<EquilibriumPhase> createEquilibriumPhases(
- boost::optional<BaseLib::ConfigTree> const& config);
+ boost::optional<BaseLib::ConfigTree> const& config,
+ MeshLib::Mesh const& mesh);
} // namespace ChemistryLib
diff --git a/ChemistryLib/PhreeqcIOData/CreateKineticReactant.cpp b/ChemistryLib/PhreeqcIOData/CreateKineticReactant.cpp
index ac9be9045e42071c1070a325fc6f8d23d23b830e..54d312c284a24807e4b983079f7efe3aaf3f0ff6 100644
--- a/ChemistryLib/PhreeqcIOData/CreateKineticReactant.cpp
+++ b/ChemistryLib/PhreeqcIOData/CreateKineticReactant.cpp
@@ -7,13 +7,18 @@
*
*/
+#include <boost/optional/optional.hpp>
+
#include "BaseLib/ConfigTree.h"
+#include "CreateKineticReactant.h"
#include "KineticReactant.h"
+#include "MeshLib/Mesh.h"
namespace ChemistryLib
{
std::vector<KineticReactant> createKineticReactants(
- boost::optional<BaseLib::ConfigTree> const& config)
+ boost::optional<BaseLib::ConfigTree> const& config,
+ MeshLib::Mesh const& mesh)
{
if (!config)
{
@@ -43,9 +48,16 @@ std::vector<KineticReactant> createKineticReactants(
reactant_config.getConfigParameter<std::vector<double>>(
"parameters", {});
+ auto amount = MeshLib::getOrCreateMeshProperty<double>(
+ const_cast<MeshLib::Mesh&>(mesh),
+ name,
+ MeshLib::MeshItemType::Node,
+ 1);
+ std::fill(std::begin(*amount), std::end(*amount), initial_amount);
+
kinetic_reactants.emplace_back(std::move(name),
std::move(chemical_formula),
- initial_amount,
+ amount,
std::move(parameters));
}
diff --git a/ChemistryLib/PhreeqcIOData/CreateKineticReactant.h b/ChemistryLib/PhreeqcIOData/CreateKineticReactant.h
index 8aaa80ac9e7fadfd46046126d15d3c36eb4ab666..65811f4b5b6b6613378e3876eda5eb4b8cfab38a 100644
--- a/ChemistryLib/PhreeqcIOData/CreateKineticReactant.h
+++ b/ChemistryLib/PhreeqcIOData/CreateKineticReactant.h
@@ -9,7 +9,7 @@
#pragma once
-#include <boost/optional/optional.hpp>
+#include <boost/optional/optional_fwd.hpp>
#include <vector>
namespace BaseLib
@@ -17,6 +17,11 @@ namespace BaseLib
class ConfigTree;
}
+namespace MeshLib
+{
+class Mesh;
+}
+
namespace ChemistryLib
{
struct KineticReactant;
@@ -25,5 +30,6 @@ struct KineticReactant;
namespace ChemistryLib
{
std::vector<KineticReactant> createKineticReactants(
- boost::optional<BaseLib::ConfigTree> const& config);
+ boost::optional<BaseLib::ConfigTree> const& config,
+ MeshLib::Mesh const& mesh);
} // namespace ChemistryLib
diff --git a/ChemistryLib/PhreeqcIOData/EquilibriumPhase.cpp b/ChemistryLib/PhreeqcIOData/EquilibriumPhase.cpp
index 378eb2a4cdbc1779015a01c08e0f243f2f00c706..52253ab3c5c221ac11dcebefa1af93c298b8b453 100644
--- a/ChemistryLib/PhreeqcIOData/EquilibriumPhase.cpp
+++ b/ChemistryLib/PhreeqcIOData/EquilibriumPhase.cpp
@@ -13,15 +13,10 @@
namespace ChemistryLib
{
-std::ostream& operator<<(std::ostream& os,
- EquilibriumPhase const& equilibrium_phase)
+void EquilibriumPhase::print(std::ostream& os,
+ std::size_t const chemical_system_id) const
{
- os << equilibrium_phase.name;
-
- os << " " << equilibrium_phase.saturation_index;
-
- os << " " << equilibrium_phase.amount << "\n";
-
- return os;
+ os << name << " " << saturation_index << " "
+ << (*amount)[chemical_system_id] << "\n";
}
} // namespace ChemistryLib
diff --git a/ChemistryLib/PhreeqcIOData/EquilibriumPhase.h b/ChemistryLib/PhreeqcIOData/EquilibriumPhase.h
index 7851377236efa1172fc25e6c6b59f56d5b6fb59d..19fe180834bc34bdda5bd579439dc4914fc7c15b 100644
--- a/ChemistryLib/PhreeqcIOData/EquilibriumPhase.h
+++ b/ChemistryLib/PhreeqcIOData/EquilibriumPhase.h
@@ -14,6 +14,7 @@
#include <vector>
#include "ChemistryLib/Output.h"
+#include "MeshLib/PropertyVector.h"
namespace BaseLib
{
@@ -25,7 +26,7 @@ namespace ChemistryLib
struct EquilibriumPhase
{
EquilibriumPhase(std::string name_,
- double amount_,
+ MeshLib::PropertyVector<double>* amount_,
double saturation_index_)
: name(std::move(name_)),
amount(amount_),
@@ -33,11 +34,10 @@ struct EquilibriumPhase
{
}
- friend std::ostream& operator<<(std::ostream& os,
- EquilibriumPhase const& equilibrium_phase);
+ void print(std::ostream& os, std::size_t const chemical_system_id) const;
std::string const name;
- double amount;
+ MeshLib::PropertyVector<double>* amount;
double const saturation_index;
static const ItemType item_type = ItemType::EquilibriumPhase;
};
diff --git a/ChemistryLib/PhreeqcIOData/KineticReactant.cpp b/ChemistryLib/PhreeqcIOData/KineticReactant.cpp
index 8f2ab34fd521a6c1ffd6a5daafc6772d40393348..f2d7402139b386aff44b6e2e25461b50d69986f0 100644
--- a/ChemistryLib/PhreeqcIOData/KineticReactant.cpp
+++ b/ChemistryLib/PhreeqcIOData/KineticReactant.cpp
@@ -13,28 +13,26 @@
namespace ChemistryLib
{
-std::ostream& operator<<(std::ostream& os,
- KineticReactant const& kinetic_reactant)
+void KineticReactant::print(std::ostream& os,
+ std::size_t const chemical_system_id) const
{
- os << kinetic_reactant.name << "\n";
+ os << name << "\n";
- if (!kinetic_reactant.chemical_formula.empty())
+ if (!chemical_formula.empty())
{
- os << "-formula " << kinetic_reactant.chemical_formula << "\n";
+ os << "-formula " << chemical_formula << "\n";
}
- os << "-m " << kinetic_reactant.amount << "\n";
+ os << "-m " << (*amount)[chemical_system_id] << "\n";
- if (!kinetic_reactant.parameters.empty())
+ if (!parameters.empty())
{
os << "-parms";
- for (auto const& parameter : kinetic_reactant.parameters)
+ for (auto const& parameter : parameters)
{
os << " " << parameter;
}
os << "\n";
}
-
- return os;
}
} // namespace ChemistryLib
diff --git a/ChemistryLib/PhreeqcIOData/KineticReactant.h b/ChemistryLib/PhreeqcIOData/KineticReactant.h
index ed8a983c5569220d11cc2ac02971e04257974019..dc993e2f31d615b6feee7e75a5b9207b597785c3 100644
--- a/ChemistryLib/PhreeqcIOData/KineticReactant.h
+++ b/ChemistryLib/PhreeqcIOData/KineticReactant.h
@@ -15,6 +15,7 @@
#include <vector>
#include "ChemistryLib/Output.h"
+#include "MeshLib/PropertyVector.h"
namespace ChemistryLib
{
@@ -22,7 +23,7 @@ struct KineticReactant
{
KineticReactant(std::string&& name_,
std::string&& chemical_formula_,
- double amount_,
+ MeshLib::PropertyVector<double>* amount_,
std::vector<double>&& parameters_)
: name(std::move(name_)),
chemical_formula(std::move(chemical_formula_)),
@@ -31,12 +32,11 @@ struct KineticReactant
{
}
- friend std::ostream& operator<<(std::ostream& os,
- KineticReactant const& kinetic_reactant);
+ void print(std::ostream& os, std::size_t const chemical_system_id) const;
std::string const name;
std::string const chemical_formula;
- double amount;
+ MeshLib::PropertyVector<double>* amount;
std::vector<double> const parameters;
static const ItemType item_type = ItemType::KineticReactant;
};
diff --git a/ProcessLib/ComponentTransport/Tests.cmake b/ProcessLib/ComponentTransport/Tests.cmake
index 3d95092f9c42841bc5344f92005d017669ec69b1..b730807279bc26810aafa70f83fad658b855b1fa 100644
--- a/ProcessLib/ComponentTransport/Tests.cmake
+++ b/ProcessLib/ComponentTransport/Tests.cmake
@@ -959,6 +959,11 @@ AddTest(
1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu 1d_isofrac_pcs_3_ts_126_t_12600.000000.vtu Syntheticb Syntheticb 1e-10 1e-16
1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu 1d_isofrac_pcs_3_ts_168_t_16800.000000.vtu Syntheticb Syntheticb 1e-10 1e-16
1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu 1d_isofrac_pcs_3_ts_210_t_21000.000000.vtu Syntheticb Syntheticb 1e-10 1e-16
+ 1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu 1d_isofrac_pcs_3_ts_42_t_4200.000000.vtu Productc Productc 1e-10 1e-16
+ 1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu 1d_isofrac_pcs_3_ts_84_t_8400.000000.vtu Productc Productc 1e-10 1e-16
+ 1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu 1d_isofrac_pcs_3_ts_126_t_12600.000000.vtu Productc Productc 1e-10 1e-16
+ 1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu 1d_isofrac_pcs_3_ts_168_t_16800.000000.vtu Productc Productc 1e-10 1e-16
+ 1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu 1d_isofrac_pcs_3_ts_210_t_21000.000000.vtu Productc Productc 1e-10 1e-16
1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu 1d_isofrac_pcs_3_ts_42_t_4200.000000.vtu H H 1e-10 1e-16
1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu 1d_isofrac_pcs_3_ts_84_t_8400.000000.vtu H H 1e-10 1e-16
1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu 1d_isofrac_pcs_3_ts_126_t_12600.000000.vtu H H 1e-10 1e-16
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu
index 5ee9b878a679ce7ed133c81da8487be7a9636541..47457c9bd9c609a132bc3a3f65a3430e8a8ee9f4 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_126_t_12600.000000_expected.vtu
@@ -1,3 +1,3 @@
version https://git-lfs.github.com/spec/v1
-oid sha256:20f88d0c2944c96400faf5632c5047cf3adbd05567f1174534dcdf01b8b400fe
-size 13290
+oid sha256:ea93899ac9913361c225bffcd3d5913aa51179c73c922bba20af82199c22fd08
+size 16506
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu
index 4e71f815d8a965cefccf8412b99361b37c7bd4bf..e18b33406a8d2ea52388f52db7540d236eec5d47 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_168_t_16800.000000_expected.vtu
@@ -1,3 +1,3 @@
version https://git-lfs.github.com/spec/v1
-oid sha256:bcc3fb0e16cd83f3d481355883f8c462ad8298524a34aae75140d4acb512fb76
-size 13358
+oid sha256:4882f6b969e527b680dce43f3b0e8befdd7ca01b42b139dcd651827d1f4489cc
+size 16570
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu
index b56011225f41ce029a7595261da367498420bf20..69e9bed01c63457552cadc0aec406d10ebf402dc 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_210_t_21000.000000_expected.vtu
@@ -1,3 +1,3 @@
version https://git-lfs.github.com/spec/v1
-oid sha256:fa28c23d46ab3d4d5ffb3474b1e737898df91934b2e9b0315131aaea01e006ce
-size 13405
+oid sha256:3734eb6db83e05196ab74170d559b4aecfa15f50b5ed4e7da4b39b6eccd14c10
+size 16594
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu
index 6b936e690eb41753293a02dddc1566f4de8d46a8..5600021c4781085de914393dac632dbe11ca4d2e 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_42_t_4200.000000_expected.vtu
@@ -1,3 +1,3 @@
version https://git-lfs.github.com/spec/v1
-oid sha256:25aea348b8aa1ab9ec6cec88b7f711e63f0951927c24f14e90fca37a7a7cae72
-size 12650
+oid sha256:a8058816d07d3088f52776b86ebb2b653f0df9bb11018d94e39839679e553f69
+size 15334
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu
index b4079f4cc85c59c30e8be76144c82aad4b8e4881..7f2bbadfb52a3eb1ff0d5af5d4b1deec0402e8ee 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_pcs_3_ts_84_t_8400.000000_expected.vtu
@@ -1,3 +1,3 @@
version https://git-lfs.github.com/spec/v1
-oid sha256:748b19cd344ebd9593d42087abdc885cbd1e16bb1a58ba0eccb43c398e4f5a99
-size 13148
+oid sha256:31ec793bc717e9910eb48b0914ff7e9427dd160777d0407c6e0b2fbbcf0ba8d1
+size 16274