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

communicate global size setting up blackoil extension boundary conditions

Browse Source
This commit is contained in:
Arne Morten Kvarving 2020-01-29 16:46:27 +01:00
parent 55b5ebd3df
commit ef89d656fb
1 changed files with 8 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
ebos/eclproblem.hh
View File

@ -2753,11 +2753,16 @@ private:
const auto& simulator = this->simulator(); const auto& simulator = this->simulator();
const auto& vanguard = simulator.vanguard(); const auto& vanguard = simulator.vanguard();
const auto& eclState = vanguard.eclState(); const auto& eclState = vanguard.eclState();
const auto& comm = vanguard.gridView().comm();
size_t numDof = this->model().numGridDof(); size_t numDof = this->model().numGridDof();
size_t globalSize;
if (comm.rank() == 0)
globalSize = eclState.getInputGrid().getCartesianSize();
comm.broadcast(&globalSize, 1, 0);
if (enableSolvent) { if (enableSolvent) {
std::vector<double> solventSaturationData(eclState.getInputGrid().getCartesianSize(), 0.0); std::vector<double> solventSaturationData(globalSize, 0.0);
if (eclState.fieldProps().has_double("SSOL")) if (eclState.fieldProps().has_double("SSOL"))
solventSaturationData = eclState.fieldProps().get_global_double("SSOL"); solventSaturationData = eclState.fieldProps().get_global_double("SSOL");
@ -2771,7 +2776,7 @@ private:
} }
if (enablePolymer) { if (enablePolymer) {
std::vector<double> polyConcentrationData(eclState.getInputGrid().getCartesianSize(), 0.0); std::vector<double> polyConcentrationData(globalSize, 0.0);
if (eclState.fieldProps().has_double("SPOLY")) if (eclState.fieldProps().has_double("SPOLY"))
polyConcentrationData = eclState.fieldProps().get_global_double("SPOLY"); polyConcentrationData = eclState.fieldProps().get_global_double("SPOLY");
@ -2785,7 +2790,7 @@ private:
} }
if (enablePolymerMolarWeight) { if (enablePolymerMolarWeight) {
std::vector<double> polyMoleWeightData(eclState.getInputGrid().getCartesianSize(), 0.0); std::vector<double> polyMoleWeightData(globalSize, 0.0);
if (eclState.fieldProps().has_double("SPOLYMW")) if (eclState.fieldProps().has_double("SPOLYMW"))
polyMoleWeightData = eclState.fieldProps().get_global_double("SPOLYMW"); polyMoleWeightData = eclState.fieldProps().get_global_double("SPOLYMW");
polymerMoleWeight_.resize(numDof, 0.0); polymerMoleWeight_.resize(numDof, 0.0);