diff --git a/opm/core/fluid/RockCompressibility.cpp b/opm/core/fluid/RockCompressibility.cpp index ca4e3ec6..bb6863c2 100644 --- a/opm/core/fluid/RockCompressibility.cpp +++ b/opm/core/fluid/RockCompressibility.cpp @@ -19,12 +19,22 @@ #include #include +#include +#include #include #include namespace Opm { + RockCompressibility::RockCompressibility(const parameter::ParameterGroup& param) + : pref_(0.0), + rock_comp_(0.0) + { + pref_ = param.getDefault("rock_compressibility_pref", 100*unit::barsa); + pref_ = param.getDefault("rock_compressibility", 0.0); + } + RockCompressibility::RockCompressibility(const EclipseGridParser& deck) : pref_(0.0), rock_comp_(0.0) diff --git a/opm/core/fluid/RockCompressibility.hpp b/opm/core/fluid/RockCompressibility.hpp index 746ffeac..bcba5aa5 100644 --- a/opm/core/fluid/RockCompressibility.hpp +++ b/opm/core/fluid/RockCompressibility.hpp @@ -26,6 +26,7 @@ namespace Opm { class EclipseGridParser; + namespace parameter { class ParameterGroup; } class RockCompressibility { @@ -33,6 +34,9 @@ namespace Opm /// Construct from input deck. RockCompressibility(const EclipseGridParser& deck); + /// Construct from parameters. + RockCompressibility(const parameter::ParameterGroup& param); + /// Porosity multiplier. double poroMult(double pressure);