From a53bcb11bb859e00f99bed69979e6dd7726de2e7 Mon Sep 17 00:00:00 2001
From: Wenqing Wang <wenqing.wang@ufz.de>
Date: Fri, 11 Nov 2016 12:22:40 +0100
Subject: [PATCH] [Unsat] Check m<0 for van Genuchten relative permeability
model
---
.../RelativePermeability/CreateRelativePermeabilityModel.cpp | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp b/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp
index 242033299c5..95275cfe9a2 100644
--- a/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp
+++ b/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp
@@ -52,7 +52,7 @@ std::unique_ptr<RelativePermeability> createWettingPhaseVanGenuchten(
//! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__WettingPhaseVanGenuchten__m}
const double m = config.getConfigParameter<double>("m");
- if (m > 1.0) // m <= 1
+ if (m < 0. || m > 1.0) // m <= 1
{
OGS_FATAL(
"The exponent parameter of WettingPhaseVanGenuchten relative\n"
@@ -84,7 +84,7 @@ std::unique_ptr<RelativePermeability> createNonWettingPhaseVanGenuchten(
//! \ogs_file_param{material_property__porous_medium__porous_medium__relative_permeability__NonWettingPhaseVanGenuchten__m}
const double m = config.getConfigParameter<double>("m");
- if (m > 1.0) // m <= 1
+ if (m < 0.|| m > 1.0) // m <= 1
{
OGS_FATAL(
"The exponent parameter of NonWettingPhaseVanGenuchten relative\n"
--
GitLab