diff --git a/opm/simulators/wells/WellState.hpp b/opm/simulators/wells/WellState.hpp index 3ec965b82..15e5fc068 100644 --- a/opm/simulators/wells/WellState.hpp +++ b/opm/simulators/wells/WellState.hpp @@ -44,7 +44,7 @@ namespace Opm public: typedef std::array< int, 3 > mapentry_t; typedef std::map< std::string, mapentry_t > WellMapType; - + /// Allocate and initialize if wells is non-null. /// Also tries to give useful initial values to the bhp() and /// wellRates() fields, depending on controls. The diff --git a/opm/simulators/wells/WellStateFullyImplicitBlackoil.hpp b/opm/simulators/wells/WellStateFullyImplicitBlackoil.hpp index fb6dd70fd..5b74f6e77 100644 --- a/opm/simulators/wells/WellStateFullyImplicitBlackoil.hpp +++ b/opm/simulators/wells/WellStateFullyImplicitBlackoil.hpp @@ -23,7 +23,6 @@ #include #include -#include #include #include @@ -85,7 +84,6 @@ namespace Opm const SummaryState& summary_state, const int globalNumberOfWells) { - //Todo: make it reservoir perforation rates std::vector perforationRates; perforationRates = prevState->perfphaserates_; // call init on base class @@ -187,6 +185,7 @@ namespace Opm perfRateBrine_.clear(); perfRateBrine_.resize(nperf, 0.0); + // intialize wells that have been there before // order may change so the mapping is based on the well name if (prevState && !prevState->wellMap().empty()) {