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

Added surfaceDensities().

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2011-01-11 16:11:50 +01:00
parent 30a7644e9f
commit 4752cecffa
1 changed files with 11 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

11
dune/porsol/blackoil/BlackoilFluid.hpp
View File

@ -26,6 +26,7 @@
#include <dune/porsol/blackoil/fluid/FluidStateBlackoil.hpp>
#include <dune/common/EclipseGridParser.hpp>
#include <dune/common/fvector.hh>
#include <dune/common/Units.hpp>
#include <vector>
@ -49,6 +50,11 @@ namespace Opm
{
fmi_params_.init(parser);
FluidSystemBlackoil<>::init(parser);
const double density_unit = parser.units().density;
const std::vector<double>& dens = parser.getDENSITY().densities_[0];
surface_densities_[Oil] = Dune::unit::convert::from(dens[0], density_unit);
surface_densities_[Water] = Dune::unit::convert::from(dens[1], density_unit);
surface_densities_[Gas] = Dune::unit::convert::from(dens[2], density_unit);
}
FluidState computeState(PhaseVec phase_pressure, CompVec z) const
{
@ -63,8 +69,13 @@ namespace Opm
}
return state;
}
const CompVec& surfaceDensities() const
{
return surface_densities_;
}
private:
FluidMatrixInteractionBlackoilParams<double> fmi_params_;
CompVec surface_densities_;
};