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

Bugfixes in RockCompressibility::isActive() and constructor.

Browse Source
This commit is contained in:
Atgeirr Flø Rasmussen 2012-03-20 21:36:31 +01:00
parent faa334a4d0
commit 3fe00d0ea5
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
opm/core/fluid/RockCompressibility.cpp
View File

@ -41,12 +41,12 @@ namespace Opm
{
if (deck.hasField("ROCKTAB")) {
const table_t& rt = deck.getROCKTAB().rocktab_;
int n = rt[0].size();
int n = rt[0][0].size();
p_.resize(n);
poromult_.resize(n);
for (int i = 0; i < n; ++i) {
p_[i] = rt[0][i][0];
poromult_[i] = rt[0][i][1];
p_[i] = rt[0][0][i];
poromult_[i] = rt[0][1][i];
}
} else if (deck.hasField("ROCK")) {
const ROCK& r = deck.getROCK();
@ -59,7 +59,7 @@ namespace Opm
bool RockCompressibility::isActive() const
{
return p_.empty() && (rock_comp_ == 0.0);
return !p_.empty() || (rock_comp_ != 0.0);
}
double RockCompressibility::poroMult(double pressure) const