From f933f113a818d99e69595aca5e2e481909daff6e Mon Sep 17 00:00:00 2001
From: Markus Wolff <markus.wolff@twt-gmbh.de>
Date: Tue, 26 Feb 2013 10:09:22 +0000
Subject: [PATCH] Bugfix: corrected an index in the pn-formulation
reviewed by Melanie
git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@10274 2fb0f335-1f38-0410-981e-8018bf24f1b0
---
.../fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh b/dumux/decoupled/2p/diffusion/fvmpfa/lmethod/fvmpfal2pfaboundpressure2padaptive.hh
index d23ccf3834..d974ebe174 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]);
--
GitLab