updating the interface of function processEclipseFormat

This commit is contained in:
Kai Bao 2021-03-22 15:20:08 +01:00
parent 41fd301d16
commit 6061b5f5b0

View File

@ -343,12 +343,14 @@ protected:
{
const EclipseGrid * input_grid = nullptr;
std::vector<double> porv;
std::unordered_map<size_t, double> aquifer_cell_volumes;
// At this stage the ParallelEclipseState instance is still in global
// view; on rank 0 we have undistributed data for the entire grid, on
// the other ranks the EclipseState is empty.
if (mpiRank == 0) {
input_grid = &this->eclState().getInputGrid();
porv = this->eclState().fieldProps().porv(true);
aquifer_cell_volumes = this->eclState().aquifer().numericalAquifers().aquiferCellVolumes();
}
grid_.reset(new Dune::CpGrid());
@ -357,7 +359,8 @@ protected:
/*flipNormals=*/false,
/*clipZ=*/false,
porv,
this->eclState().getInputNNC());
this->eclState().getInputNNC(),
aquifer_cell_volumes);
// we use separate grid objects: one for the calculation of the initial condition
// via EQUIL and one for the actual simulation. The reason is that the EQUIL code