OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2024-12-28 18:21:00 -06:00
Code Issues Packages Projects Releases Wiki Activity

Add gravity() method.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2013-05-23 14:48:34 +02:00
parent d171f043e5
commit e5bd6547ae
1 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
opm/autodiff/GeoProps.hpp
View File

@ -54,6 +54,8 @@ namespace Opm
tpfa_htrans_compute(ug, props.permeability(), htrans.data()); tpfa_htrans_compute(ug, props.permeability(), htrans.data());
tpfa_trans_compute (ug, htrans.data() , trans_.data()); tpfa_trans_compute (ug, htrans.data() , trans_.data());
// Gravity potential
std::fill(gravity_, gravity_ + 3, 0.0);
if (grav != 0) { if (grav != 0) {
const typename Vector::Index nd = grid.dimensions; const typename Vector::Index nd = grid.dimensions;
@ -71,17 +73,20 @@ namespace Opm
} }
} }
} }
std::copy(grav, grav + nd, gravity_);
} }
} }
const Vector& poreVolume() const { return pvol_ ; } const Vector& poreVolume() const { return pvol_ ; }
const Vector& transmissibility() const { return trans_; } const Vector& transmissibility() const { return trans_; }
const Vector& gravityPotential() const { return gpot_ ; } const Vector& gravityPotential() const { return gpot_ ; }
const double* gravity() const { return gravity_; }
private: private:
Vector pvol_ ; Vector pvol_ ;
Vector trans_; Vector trans_;
Vector gpot_ ; Vector gpot_ ;
double gravity_[3]; // Size 3 even if grid is 2-dim.
}; };
} }