diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh index d23ccf38344a8bbefe568244df360df72d3181d6..d974ebe174e048d9982e2e778b523f6f0b7bfdb5 100644 --- a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh +++ b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh @@ -2252,7 +2252,7 @@ void FVMPFAL2PFABoundPressure2PAdaptive<TypeTag>::assemble() if (i == wPhaseIdx) { //add capillary pressure term to right hand side - this->f_[globalIdx1] += (pcFluxReal[0] - pcFluxReal[1]); + this->f_[globalIdx1] += (pcFluxReal[0] - pcFluxReal[3]); this->f_[globalIdx2] += (pcFluxReal[1] - pcFluxReal[0]); this->f_[globalIdx3] += (pcFluxReal[2] - pcFluxReal[1]); this->f_[globalIdx4] += (pcFluxReal[3] - pcFluxReal[2]);