From 5b9a35a2df4669ad4a51fd1e5beaa353a172048c Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Thu, 7 Jul 2022 12:03:34 +0200
Subject: [PATCH] [LiquidFlow] Added an input for aperture size
---
.../LiquidFlow/CreateLiquidFlowProcess.cpp | 18 +++++++++++++++++-
1 file changed, 17 insertions(+), 1 deletion(-)
diff --git a/ProcessLib/LiquidFlow/CreateLiquidFlowProcess.cpp b/ProcessLib/LiquidFlow/CreateLiquidFlowProcess.cpp
index 61c646866cc..f083b608e69 100644
--- a/ProcessLib/LiquidFlow/CreateLiquidFlowProcess.cpp
+++ b/ProcessLib/LiquidFlow/CreateLiquidFlowProcess.cpp
@@ -124,11 +124,27 @@ std::unique_ptr<Process> createLiquidFlowProcess(
checkMPLProperties(mesh, *media_map);
DBUG("Media properties verified.");
+ auto const* aperture_size_parameter = &ParameterLib::findParameter<double>(
+ ProcessLib::Process::constant_one_parameter_name, parameters, 1);
+
+ auto const aperture_config =
+ //! \ogs_file_param{prj__processes__process__LIQUID_FLOW__aperture_size}
+ config.getConfigSubtreeOptional("aperture_size");
+ if (aperture_config)
+ {
+ aperture_size_parameter = &ParameterLib::findParameter<double>(
+ //! \ogs_file_param_special{prj__processes__process__LIQUID_FLOW__aperture_size__parameter}
+ *aperture_config, "parameter", parameters, 1);
+ }
+
LiquidFlowData process_data{
std::move(media_map),
MeshLib::getElementRotationMatrices(
mesh_space_dimension, mesh.getDimension(), mesh.getElements()),
- mesh_space_dimension, std::move(specific_body_force), has_gravity};
+ mesh_space_dimension,
+ std::move(specific_body_force),
+ has_gravity,
+ *aperture_size_parameter};
return std::make_unique<LiquidFlowProcess>(
std::move(name), mesh, std::move(jacobian_assembler), parameters,
--
GitLab