diff --git a/dune/upscaling/SteadyStateUpscalerImplicit.hpp b/dune/upscaling/SteadyStateUpscalerImplicit.hpp index 18ededa..cce128b 100644 --- a/dune/upscaling/SteadyStateUpscalerImplicit.hpp +++ b/dune/upscaling/SteadyStateUpscalerImplicit.hpp @@ -114,6 +114,7 @@ namespace Dune // ------- Data members ------- std::vector last_saturation_state_; + bool use_gravity_; bool output_vtk_; bool output_ecl_; bool print_inoutflows_; diff --git a/dune/upscaling/SteadyStateUpscalerImplicit_impl.hpp b/dune/upscaling/SteadyStateUpscalerImplicit_impl.hpp index 4c70f12..10e01d8 100644 --- a/dune/upscaling/SteadyStateUpscalerImplicit_impl.hpp +++ b/dune/upscaling/SteadyStateUpscalerImplicit_impl.hpp @@ -77,7 +77,7 @@ namespace Dune inline void SteadyStateUpscalerImplicit::initImpl(const Opm::parameter::ParameterGroup& param) { Super::initImpl(param); - use_gravity_ = param.getDefault("output_vtk", use_gravity_); + use_gravity_ = param.getDefault("use_gravity", use_gravity_); output_vtk_ = param.getDefault("output_vtk", output_vtk_); output_ecl_ = param.getDefault("output_ecl", output_ecl_); if (output_ecl_) { @@ -157,8 +157,8 @@ namespace Dune Opm::SparseVector injection(num_cells); // Gravity. FieldVector gravity(0.0); - if(use_gravity_){ - gravity[2] = -Opm::unit::gravity; + if (use_gravity_) { + gravity[2] = Opm::unit::gravity; } if (gravity.two_norm() > 0.0) { diff --git a/dune/upscaling/SteadyStateUpscaler_impl.hpp b/dune/upscaling/SteadyStateUpscaler_impl.hpp index ee346d4..1dc9b4e 100644 --- a/dune/upscaling/SteadyStateUpscaler_impl.hpp +++ b/dune/upscaling/SteadyStateUpscaler_impl.hpp @@ -143,8 +143,8 @@ namespace Dune Opm::SparseVector injection(num_cells); // Gravity. FieldVector gravity(0.0); - if(use_gravity_){ - gravity[2] = -Opm::unit::gravity; + if (use_gravity_) { + gravity[2] = Opm::unit::gravity; } if (gravity.two_norm() > 0.0) { MESSAGE("Warning: Gravity is experimental for flow solver.");