OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

Remove redundent check

Browse Source
This commit is contained in:
Svenn Tveit 2023-08-11 16:38:49 +02:00
parent b84837fc61
commit 161c895c63
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ebos/equil/initstateequil_impl.hh
View File

@ -1669,6 +1669,8 @@ applyNumericalAquifers_(const GridView& gridView,
ElementMapper elemMapper(gridView, Dune::mcmgElementLayout()); ElementMapper elemMapper(gridView, Dune::mcmgElementLayout());
auto elemIt = gridView.template begin</*codim=*/0>(); auto elemIt = gridView.template begin</*codim=*/0>();
const auto& elemEndIt = gridView.template end</*codim=*/0>(); const auto& elemEndIt = gridView.template end</*codim=*/0>();
const auto oilPos = FluidSystem::oilPhaseIdx;
const auto gasPos = FluidSystem::gasPhaseIdx;
for (; elemIt != elemEndIt; ++elemIt) { for (; elemIt != elemEndIt; ++elemIt) {
const Element& element = *elemIt; const Element& element = *elemIt;
const unsigned int elemIdx = elemMapper.index(element); const unsigned int elemIdx = elemMapper.index(element);
@ -1676,16 +1678,12 @@ applyNumericalAquifers_(const GridView& gridView,
const auto search = num_aqu_cells.find(cartIx); const auto search = num_aqu_cells.find(cartIx);
if (search != num_aqu_cells.end()) { if (search != num_aqu_cells.end()) {
// numerical aquifer cells are filled with water initially // numerical aquifer cells are filled with water initially
if (FluidSystem::phaseIsActive(watPos)) { this->sat_[watPos][elemIdx] = 1.;
this->sat_[watPos][elemIdx] = 1.;
}
const auto oilPos = FluidSystem::oilPhaseIdx;
if (!co2store_or_h2store && FluidSystem::phaseIsActive(oilPos)) { if (!co2store_or_h2store && FluidSystem::phaseIsActive(oilPos)) {
this->sat_[oilPos][elemIdx] = 0.; this->sat_[oilPos][elemIdx] = 0.;
} }
const auto gasPos = FluidSystem::gasPhaseIdx;
if (FluidSystem::phaseIsActive(gasPos)) { if (FluidSystem::phaseIsActive(gasPos)) {
this->sat_[gasPos][elemIdx] = 0.; this->sat_[gasPos][elemIdx] = 0.;
} }