CreateChemicalSystem.cpp

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright (c) OpenGeoSys Community (opengeosys.org)
// SPDX-License-Identifier: BSD-3-Clause
 
#include "CreateChemicalSystem.h"
 
#include "BaseLib/ConfigTree.h"
#include "ChemicalSystem.h"
#include "CreateAqueousSolution.h"
#include "CreateEquilibriumReactants.h"
#include "CreateExchange.h"
#include "CreateKineticReactant.h"
#include "CreateSurface.h"
#include "EquilibriumReactant.h"
#include "Exchange.h"
#include "KineticReactant.h"
#include "Surface.h"
 
namespace ChemistryLib
{
namespace PhreeqcIOData
{
std::unique_ptr<ChemicalSystem> createChemicalSystem(
    BaseLib::ConfigTree const& config, MeshLib::Mesh& mesh)
{
    // solution
    auto aqueous_solution = createAqueousSolution(
        config.getConfigSubtree("solution"), mesh);
 
    // kinetic reactants
    auto kinetic_reactants = createKineticReactants(
        config.getConfigSubtreeOptional("kinetic_reactants"), mesh);
 
    // equilibrium reactants
    auto equilibrium_reactants = createEquilibriumReactants(
        config.getConfigSubtreeOptional("equilibrium_reactants"), mesh);
 
    // exchangers
    auto exchangers = createExchange(
        config.getConfigSubtreeOptional("exchangers"), mesh);
 
    // surface
    auto surface = createSurface(
        config.getConfigSubtreeOptional("surface"), mesh);
 
    return std::make_unique<ChemicalSystem>(std::move(aqueous_solution),
                                            std::move(kinetic_reactants),
                                            std::move(equilibrium_reactants),
                                            std::move(exchangers),
                                            std::move(surface));
}
}  // namespace PhreeqcIOData
}  // namespace ChemistryLib