diff --git a/opm/core/pressure/CompressibleTpfa.cpp b/opm/core/pressure/CompressibleTpfa.cpp index 96e162268..cc15f7f56 100644 --- a/opm/core/pressure/CompressibleTpfa.cpp +++ b/opm/core/pressure/CompressibleTpfa.cpp @@ -181,11 +181,12 @@ namespace Opm const int nperf = wells_->well_connpos[nw]; const int dim = grid_.dimensions; const double grav = gravity_ ? gravity_[dim - 1] : 0.0; + wellperf_gpot_.clear(); + wellperf_gpot_.resize(np*nperf, 0.0); if (grav == 0.0) { - wellperf_gpot_.clear(); - wellperf_gpot_.resize(np*nperf, 0.0); return; } + // Temporary storage for perforation A matrices and densities. std::vector A(np*np, 0.0); std::vector rho(np, 0.0);