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

Use Density from EclipseState

Browse Source
This commit is contained in:
Jørgen Kvalsvik
2016-12-15 16:03:53 +01:00
parent f8b1bef974
commit 5278a1f193
2 changed files with 7 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
opm/core/props/IncompPropertiesSinglePhase.cpp
View File

@@ -37,9 +37,9 @@ namespace Opm
{
rock_.init(eclState, grid.number_of_cells, grid.global_cell, grid.cartdims);
if (deck.hasKeyword("DENSITY")) {
const auto& densityRecord = deck.getKeyword("DENSITY").getRecord(0);
surface_density_ = densityRecord.getItem("OIL").getSIDouble(0);
const auto& densities = eclState.getTableManager().getDensityTable();
if( !densities.empty() ) {
surface_density_ = densities[0].oil;
} else {
surface_density_ = 1000.0;
OPM_MESSAGE("Input is missing DENSITY -- using a standard density of "

8
opm/core/props/pvt/PvtPropertiesIncompFromDeck.cpp
View File

@@ -47,10 +47,10 @@ namespace Opm
}
// Surface densities. Accounting for different orders in eclipse and our code.
if (deck.hasKeyword("DENSITY")) {
const auto& densityRecord = deck.getKeyword("DENSITY").getRecord(region_number);
surface_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] = densityRecord.getItem("OIL").getSIDouble(0);
surface_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] = densityRecord.getItem("WATER").getSIDouble(0);
const auto& densities = es.getTableManager().getDensityTable();
if (!densities.empty()) {
surface_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] = densities[region_number].water;
surface_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] = densities[region_number].oil;
} else {
OPM_THROW(std::runtime_error, "Input is missing DENSITY\n");
}