diff --git a/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp b/MaterialLib/PorousMedium/UnsaturatedProperty/RelativePermeability/CreateRelativePermeabilityModel.cpp index 242033299c5fd0bce2bfc8b3d7ae22894eca5244..95275cfe9a274e6b888cdc07cd933838400fda34 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"