diff --git a/dumux/implicit/2p2c/2p2cmodel.hh b/dumux/implicit/2p2c/2p2cmodel.hh index 1d97824668ef57193726b7ebb6351a24d8cd64ef..47c8a6dd7ec3e9d338dc64c5291ac9e07abdc8bc 100644 --- a/dumux/implicit/2p2c/2p2cmodel.hh +++ b/dumux/implicit/2p2c/2p2cmodel.hh @@ -219,12 +219,16 @@ public: ElementIterator elemIt = this->gridView_().template begin<0>(); const ElementIterator elemEndIt = this->gridView_().template end<0>(); for (; elemIt != elemEndIt; ++elemIt) { - this->localResidual().evalPhaseStorage(*elemIt, phaseIdx); + if(elemIt->partitionType() == Dune::InteriorEntity) + { + + + this->localResidual().evalPhaseStorage(*elemIt, phaseIdx); for (unsigned int i = 0; i < this->localResidual().storageTerm().size(); ++i) storage += this->localResidual().storageTerm()[i]; - } - + } + } if (this->gridView_().comm().size() > 1) storage = this->gridView_().comm().sum(storage); } diff --git a/dumux/implicit/3p3c/3p3cmodel.hh b/dumux/implicit/3p3c/3p3cmodel.hh index 800680b03a395dbcc64f05a52e2558f2969445f9..75486f1dc9db9a4da334649c9b169ee5815e9f4a 100644 --- a/dumux/implicit/3p3c/3p3cmodel.hh +++ b/dumux/implicit/3p3c/3p3cmodel.hh @@ -207,12 +207,14 @@ public: ElementIterator elemIt = this->gridView_().template begin<0>(); const ElementIterator elemEndIt = this->gridView_().template end<0>(); for (; elemIt != elemEndIt; ++elemIt) { + if(elemIt->partitionType() == Dune::InteriorEntity) + { this->localResidual().evalPhaseStorage(*elemIt, phaseIdx); for (unsigned int i = 0; i < this->localResidual().storageTerm().size(); ++i) storage += this->localResidual().storageTerm()[i]; - } - + } + } if (this->gridView_().comm().size() > 1) storage = this->gridView_().comm().sum(storage); }