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

only update diffusivity when DIFFUSIVE is in the deck

Browse Source
This commit is contained in:
Tor Harald Sandve 2021-02-25 11:14:12 +01:00
parent 60534d66b6
commit 48ad676010
1 changed files with 3 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
ebos/ecltransmissibility.hh
View File

@ -127,7 +127,7 @@ public:
// get the ntg values, the ntg values are modified for the cells merged with minpv
const std::vector<double>& ntg = eclState.fieldProps().get_double("NTG");
const bool updateDiffusivity = eclState.getSimulationConfig().isDiffusive() && enableDiffusion;
unsigned numElements = elemMapper.size();
extractPermeability_();
@ -183,7 +183,7 @@ public:
}
// if diffusion is enabled, let's do the same for the "diffusivity"
if (enableDiffusion) {
if (updateDiffusivity) {
diffusivity_->clear();
diffusivity_->reserve(numElements*3*1.05);
extractPorosity_();
@ -396,7 +396,7 @@ public:
trans_[isId_(elemIdx, outsideElemIdx)] = trans;
// update the "thermal half transmissibility" for the intersection
if (enableDiffusion) {
if (updateDiffusivity) {
Scalar halfDiffusivity1;
Scalar halfDiffusivity2;