diff --git a/MaterialLib/MPL/Properties/RelPermBrooksCorey.h b/MaterialLib/MPL/Properties/RelPermBrooksCorey.h index 010b5985e0016bafa71e396b2fdd7b7e3777d46b..247a9b600ec965005d7ecfb54c188556d2925bac 100644 --- a/MaterialLib/MPL/Properties/RelPermBrooksCorey.h +++ b/MaterialLib/MPL/Properties/RelPermBrooksCorey.h @@ -95,6 +95,10 @@ inline std::unique_ptr<RelPermBrooksCorey> createRelPermBrooksCorey( auto const exponent = //! \ogs_file_param{prj__media__medium__properties__property__RelPermBrooksCorey__lambda} config.getConfigParameter<double>("lambda"); + if (exponent == 0.) + { + OGS_FATAL("Exponent 'lambda' must be positive."); + } return std::make_unique<RelPermBrooksCorey>( residual_liquid_saturation,