diff --git a/opm/polymer/fullyimplicit/SimulatorFullyImplicitBlackoilPolymer_impl.hpp b/opm/polymer/fullyimplicit/SimulatorFullyImplicitBlackoilPolymer_impl.hpp index cd275cd3f..7a745ef41 100644 --- a/opm/polymer/fullyimplicit/SimulatorFullyImplicitBlackoilPolymer_impl.hpp +++ b/opm/polymer/fullyimplicit/SimulatorFullyImplicitBlackoilPolymer_impl.hpp @@ -328,12 +328,7 @@ namespace Opm props_.permeability()); const Wells* wells = wells_manager.c_wells(); WellStateFullyImplicitBlackoil well_state; - well_state.init(wells, state.blackoilState()); - if (timer.currentStepNum() != 0) { - // Transfer previous well state to current. - well_state.partialCopy(prev_well_state, *wells, prev_well_state.numWells()); - } - + well_state.init(wells, state.blackoilState(), prev_well_state); // compute polymer inflow std::unique_ptr polymer_inflow_ptr; if (deck_->hasKeyword("WPOLYMER")) { diff --git a/opm/polymer/fullyimplicit/SimulatorFullyImplicitCompressiblePolymer.cpp b/opm/polymer/fullyimplicit/SimulatorFullyImplicitCompressiblePolymer.cpp index edb555ae2..f4d35ed51 100644 --- a/opm/polymer/fullyimplicit/SimulatorFullyImplicitCompressiblePolymer.cpp +++ b/opm/polymer/fullyimplicit/SimulatorFullyImplicitCompressiblePolymer.cpp @@ -278,11 +278,7 @@ namespace Opm props_.permeability()); const Wells* wells = wells_manager.c_wells(); WellStateFullyImplicitBlackoil well_state; - well_state.init(wells, state.blackoilState()); - if (timer.currentStepNum() != 0) { - // Transfer previous well state to current. - well_state.partialCopy(prev_well_state, *wells, prev_well_state.numWells()); - } + well_state.init(wells, state.blackoilState(), prev_well_state); //Compute polymer inflow. std::unique_ptr polymer_inflow_ptr; if (deck_->hasKeyword("WPOLYMER")) {