diff --git a/ebos/ecltransmissibility.hh b/ebos/ecltransmissibility.hh index 9d2bcaf3b..8e2f5b69c 100644 --- a/ebos/ecltransmissibility.hh +++ b/ebos/ecltransmissibility.hh @@ -127,7 +127,7 @@ public: // get the ntg values, the ntg values are modified for the cells merged with minpv const std::vector& 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;