diff --git a/opm/autodiff/WellStateMultiSegment.hpp b/opm/autodiff/WellStateMultiSegment.hpp index 5ce218bb5..e80d2402f 100644 --- a/opm/autodiff/WellStateMultiSegment.hpp +++ b/opm/autodiff/WellStateMultiSegment.hpp @@ -208,7 +208,8 @@ namespace Opm for (int p = 0; p < np; ++p) { perfPhaseRates()[np * (i + start_perforation) + p] = wellRates()[np * w + p] / double(number_of_perforations); } - perfPress()[i + start_perforation] = state.pressure()[wells[w]->wellCells()[i]]; + const double safety_factor = (wells[w]->wellType() == INJECTOR) ? 1.01 : 0.99; + perfPress()[i + start_perforation] = safety_factor * state.pressure()[wells[w]->wellCells()[i]]; } // 5. Segment rates and pressures