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

Merge pull request #269 from rolk/269_sink

Browse Source 
Use a more realistic measure for source and show use of units
This commit is contained in:
Bård Skaflestad 2013-06-28 04:56:36 -07:00
commit d3065fb3f6
1 changed files with 9 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
tutorials/tutorial2.cpp
View File

@ -87,8 +87,8 @@ int main()
using namespace Opm::unit;
using namespace Opm::prefix;
int num_phases = 1;
std::vector<double> mu(num_phases, 1.0*centi*Poise);
std::vector<double> rho(num_phases, 1000.0*kilogram/cubic(meter));
std::vector<double> viscosities(num_phases, 1.0*centi*Poise);
std::vector<double> densities(num_phases, 1000.0*kilogram/cubic(meter));
/// \internal [fluid]
/// \endinternal
/// \page tutorial2
@ -96,7 +96,7 @@ int main()
/// We define a permeability equal to 100 mD.
/// \snippet tutorial2.cpp perm
/// \internal [perm]
double k = 100.0*milli*darcy;
double permeability = 100.0*milli*darcy;
/// \internal [perm]
/// \endinternal
@ -105,8 +105,10 @@ int main()
/// We set up a simple property object for a single-phase situation.
/// \snippet tutorial2.cpp single-phase property
/// \internal [single-phase property]
Opm::IncompPropertiesBasic props(1, Opm::SaturationPropsBasic::Constant, rho,
mu, 1.0, k, dim, num_cells);
const double porosity = 1.;
Opm::IncompPropertiesBasic props(1, Opm::SaturationPropsBasic::Constant,
densities, viscosities, porosity,
permeability, dim, num_cells);
/// \internal [single-phase property]
/// /endinternal
@ -125,8 +127,8 @@ int main()
/// \snippet tutorial2.cpp source
/// \internal [source]
std::vector<double> src(num_cells, 0.0);
src[0] = 100.;
src[num_cells-1] = -100.;
src[0] = 150.*cubic(meter)/day;
src[num_cells-1] = -src[0];
/// \internal [source]
/// \endinternal