OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests 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 e53b0a58fa
commit 59cc0a1635
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); rock_.init(eclState, grid.number_of_cells, grid.global_cell, grid.cartdims);
if (deck.hasKeyword("DENSITY")) { const auto& densities = eclState.getTableManager().getDensityTable();
const auto& densityRecord = deck.getKeyword("DENSITY").getRecord(0); if( !densities.empty() ) {
surface_density_ = densityRecord.getItem("OIL").getSIDouble(0); surface_density_ = densities[0].oil;
} else { } else {
surface_density_ = 1000.0; surface_density_ = 1000.0;
OPM_MESSAGE("Input is missing DENSITY -- using a standard density of " 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. // Surface densities. Accounting for different orders in eclipse and our code.
if (deck.hasKeyword("DENSITY")) { const auto& densities = es.getTableManager().getDensityTable();
const auto& densityRecord = deck.getKeyword("DENSITY").getRecord(region_number); if (!densities.empty()) {
surface_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] = densityRecord.getItem("OIL").getSIDouble(0); surface_density_[phase_usage.phase_pos[PhaseUsage::Aqua]] = densities[region_number].water;
surface_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] = densityRecord.getItem("WATER").getSIDouble(0); surface_density_[phase_usage.phase_pos[PhaseUsage::Liquid]] = densities[region_number].oil;
} else { } else {
OPM_THROW(std::runtime_error, "Input is missing DENSITY\n"); OPM_THROW(std::runtime_error, "Input is missing DENSITY\n");
} }