Remove injector crossflow from tracer production.

This commit is contained in:
Ove Sævareid 2021-06-09 16:39:01 +02:00
parent 581408c760
commit 9478e8d03d

View File

@ -447,7 +447,7 @@ protected:
const size_t cartIdx = simulator_.vanguard().cartesianIndex(cartesianCoordinate); const size_t cartIdx = simulator_.vanguard().cartesianIndex(cartesianCoordinate);
const int I = this->cartToGlobal_[cartIdx]; const int I = this->cartToGlobal_[cartIdx];
Scalar rate = simulator_.problem().wellModel().well(well.name())->volumetricSurfaceRateForConnection(I, tr.phaseIdx_); Scalar rate = simulator_.problem().wellModel().well(well.name())->volumetricSurfaceRateForConnection(I, tr.phaseIdx_);
if (rate < 0) { //Injection rates already reported during assembly if (rate < 0 && well.isProducer()) { //Injection rates already reported during assembly
for (int tIdx =0; tIdx < tr.numTracer(); ++tIdx) { for (int tIdx =0; tIdx < tr.numTracer(); ++tIdx) {
this->wellTracerRate_.at(std::make_pair(well.name(),this->tracerNames_[tr.idx_[tIdx]])) += rate*tr.concentration_[tIdx][I]; this->wellTracerRate_.at(std::make_pair(well.name(),this->tracerNames_[tr.idx_[tIdx]])) += rate*tr.concentration_[tIdx][I];
} }