Adds has_disgas and has_vapoil flags to hydrocarbonstate initialization

Gas saturation must be used as primalvariable if DISGAS and VAPOIL is
not used.
This commit is contained in:
Tor Harald Sandve 2016-06-06 15:45:50 +02:00
parent 87c7b2c720
commit d035881a0a
2 changed files with 4 additions and 2 deletions

View File

@ -552,7 +552,9 @@ namespace Opm
props.capPress(numCells, state_->saturation().data(), cells.data(), pc.data(), nullptr);
fluidprops_->setSwatInitScaling(state_->saturation(), pc);
}
initHydroCarbonState(*state_, pu, Opm::UgGridHelpers::numCells(grid));
initHydroCarbonState(*state_, pu, Opm::UgGridHelpers::numCells(grid), deck_->hasKeyword("DISGAS"), deck_->hasKeyword("VAPOIL"));
}

View File

@ -87,7 +87,7 @@ namespace Opm
if (output_writer_.isRestart()) {
// This is a restart, populate WellState and ReservoirState state objects from restart file
output_writer_.initFromRestartFile(props_.phaseUsage(), props_.permeability(), grid_, state, prev_well_state);
initHydroCarbonState(state, props_.phaseUsage(), Opm::UgGridHelpers::numCells(grid_));
initHydroCarbonState(state, props_.phaseUsage(), Opm::UgGridHelpers::numCells(grid_), has_disgas_, has_vapoil_);
}
// Create timers and file for writing timing info.