diff --git a/examples/problems/groundwaterproblem.hh b/examples/problems/groundwaterproblem.hh index bb019f5fa..6f1fad393 100644 --- a/examples/problems/groundwaterproblem.hh +++ b/examples/problems/groundwaterproblem.hh @@ -56,23 +56,6 @@ namespace TTag { struct GroundWaterBaseProblem {}; } -template -struct LensLowerLeftX { using type = UndefinedProperty; }; -template -struct LensLowerLeftY { using type = UndefinedProperty; }; -template -struct LensLowerLeftZ { using type = UndefinedProperty; }; -template -struct LensUpperRightX { using type = UndefinedProperty; }; -template -struct LensUpperRightY { using type = UndefinedProperty; }; -template -struct LensUpperRightZ { using type = UndefinedProperty; }; -template -struct Permeability { using type = UndefinedProperty; }; -template -struct PermeabilityLens { using type = UndefinedProperty; }; - template struct Fluid { @@ -92,55 +75,6 @@ template struct Problem { using type = Opm::GroundWaterProblem; }; -template -struct LensLowerLeftX -{ - using type = GetPropType; - static constexpr type value = 0.25; -}; -template -struct LensLowerLeftY -{ - using type = GetPropType; - static constexpr type value = 0.25; -}; -template -struct LensLowerLeftZ -{ - using type = GetPropType; - static constexpr type value = 0.25; -}; -template -struct LensUpperRightX -{ - using type = GetPropType; - static constexpr type value = 0.75; -}; -template -struct LensUpperRightY -{ - using type = GetPropType; - static constexpr type value = 0.75; -}; -template -struct LensUpperRightZ -{ - using type = GetPropType; - static constexpr type value = 0.75; -}; -template -struct Permeability -{ - using type = GetPropType; - static constexpr type value = 1e-10; -}; -template -struct PermeabilityLens -{ - using type = GetPropType; - static constexpr type value = 1e-12; -}; - // Use the conjugated gradient linear solver with the default preconditioner (i.e., // ILU-0) from dune-istl template @@ -154,6 +88,30 @@ struct LinearSolverWrapper namespace Opm::Parameters { +template +struct LensLowerLeftX { using type = Properties::UndefinedProperty; }; + +template +struct LensLowerLeftY { using type = Properties::UndefinedProperty; }; + +template +struct LensLowerLeftZ { using type = Properties::UndefinedProperty; }; + +template +struct LensUpperRightX { using type = Properties::UndefinedProperty; }; + +template +struct LensUpperRightY { using type = Properties::UndefinedProperty; }; + +template +struct LensUpperRightZ { using type = Properties::UndefinedProperty; }; + +template +struct Permeability { using type = Properties::UndefinedProperty; }; + +template +struct PermeabilityLens { using type = Properties::UndefinedProperty; }; + // Enable gravity template struct EnableGravity @@ -180,6 +138,62 @@ struct InitialTimeStepSize static constexpr type value = 1; }; +template +struct LensLowerLeftX +{ + using type = GetPropType; + static constexpr type value = 0.25; +}; + +template +struct LensLowerLeftY +{ + using type = GetPropType; + static constexpr type value = 0.25; +}; + +template +struct LensLowerLeftZ +{ + using type = GetPropType; + static constexpr type value = 0.25; +}; + +template +struct LensUpperRightX +{ + using type = GetPropType; + static constexpr type value = 0.75; +}; + +template +struct LensUpperRightY +{ + using type = GetPropType; + static constexpr type value = 0.75; +}; + +template +struct LensUpperRightZ +{ + using type = GetPropType; + static constexpr type value = 0.75; +}; + +template +struct Permeability +{ + using type = GetPropType; + static constexpr type value = 1e-10; +}; + +template +struct PermeabilityLens +{ + using type = GetPropType; + static constexpr type value = 1e-12; +}; + } // namespace Opm::Parameters namespace Opm { @@ -246,20 +260,20 @@ public: eps_ = 1.0e-3; - lensLowerLeft_[0] = Parameters::get(); + lensLowerLeft_[0] = Parameters::get(); if (dim > 1) - lensLowerLeft_[1] = Parameters::get(); + lensLowerLeft_[1] = Parameters::get(); if (dim > 2) - lensLowerLeft_[2] = Parameters::get(); + lensLowerLeft_[2] = Parameters::get(); - lensUpperRight_[0] = Parameters::get(); + lensUpperRight_[0] = Parameters::get(); if (dim > 1) - lensUpperRight_[1] = Parameters::get(); + lensUpperRight_[1] = Parameters::get(); if (dim > 2) - lensUpperRight_[2] = Parameters::get(); + lensUpperRight_[2] = Parameters::get(); - intrinsicPerm_ = this->toDimMatrix_(Parameters::get()); - intrinsicPermLens_ = this->toDimMatrix_(Parameters::get()); + intrinsicPerm_ = this->toDimMatrix_(Parameters::get()); + intrinsicPermLens_ = this->toDimMatrix_(Parameters::get()); } /*! @@ -269,28 +283,28 @@ public: { ParentType::registerParameters(); - Parameters::registerParam + Parameters::registerParam ("The x-coordinate of the lens' lower-left corner [m]."); - Parameters::registerParam + Parameters::registerParam ("The x-coordinate of the lens' upper-right corner [m]."); if (dimWorld > 1) { - Parameters::registerParam + Parameters::registerParam ("The y-coordinate of the lens' lower-left corner [m]."); - Parameters::registerParam + Parameters::registerParam ("The y-coordinate of the lens' upper-right corner [m]."); } if (dimWorld > 2) { - Parameters::registerParam + Parameters::registerParam ("The z-coordinate of the lens' lower-left corner [m]."); - Parameters::registerParam + Parameters::registerParam ("The z-coordinate of the lens' upper-right corner [m]."); } - Parameters::registerParam + Parameters::registerParam ("The intrinsic permeability [m^2] of the ambient material."); - Parameters::registerParam + Parameters::registerParam ("The intrinsic permeability [m^2] of the lens."); }