OilfieldToolsNet/opm-core
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

use bar instead of Pascal as the unit for the pressure field

Browse Source 
Eclipse seems to write bars...
This commit is contained in:
Andreas Lauser 2013-10-28 11:23:31 +01:00 committed by Roland Kaufmann
parent b935697838
commit 07ed5b4339
1 changed files with 9 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
opm/core/io/eclipse/BlackoilEclipseOutputWriter.cpp
View File

@ -90,7 +90,15 @@ void BlackoilEclipseOutputWriter::writeReservoirState(const BlackoilState& reser
ecl_rst_file_start_solution(rst_file); ecl_rst_file_start_solution(rst_file);
{ {
ecl_kw_type* pressure_kw = newEclDoubleKeyword_("PRESSURE", reservoirState.pressure()); // convert the pressures from Pascals to bar because eclipse
// seems to write bars
std::vector<double> pressureBar(reservoirState.pressure());
auto it = pressureBar.begin();
const auto &endIt = pressureBar.end();
for (; it != endIt; ++it)
(*it) /= 1e5;
ecl_kw_type* pressure_kw = newEclDoubleKeyword_("PRESSURE", pressureBar);
ecl_rst_file_add_kw(rst_file, pressure_kw); ecl_rst_file_add_kw(rst_file, pressure_kw);
ecl_kw_free(pressure_kw); ecl_kw_free(pressure_kw);
} }