OilfieldToolsNet/opm-simulators
4
0
Fork 0
mirror of https://github.com/OPM/opm-simulators.git synced 2025-02-25 18:55:30 -06:00
Code Issues Packages Projects Releases Wiki Activity

fix the threshold pressure implementation

Browse Source 
before the variable for threshold pressure is used, it should be
initialized!
This commit is contained in:
Andreas Lauser 2016-02-25 15:56:30 +01:00
parent 0012aac5bd
commit 34d0184f7d
1 changed files with 4 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
applications/ebos/eclfluxmodule.hh
View File

@ -213,9 +213,11 @@ protected:
exteriorDofIdx_ = scvf.exteriorIndex();
assert(interiorDofIdx_ != exteriorDofIdx_);
trans_ = problem.transmissibility(stencil.globalSpaceIndex(interiorDofIdx_),
stencil.globalSpaceIndex(exteriorDofIdx_));
unsigned I = stencil.globalSpaceIndex(interiorDofIdx_);
unsigned J = stencil.globalSpaceIndex(exteriorDofIdx_);
trans_ = problem.transmissibility(I, J);
faceArea_ = scvf.area();
thpres_ = problem.thresholdPressure(I, J);
// estimate the gravity correction: for performance reasons we use a simplified
// approach for this flux module that assumes that gravity is constant and always