diff --git a/.hgsubstate b/.hgsubstate
index ecfdd2f..7fc450b 100644
--- a/.hgsubstate
+++ b/.hgsubstate
@@ -1 +1 @@
-c35871a1f64c0704b3f0fa7d706f1b0cc0810800 dune/porsol/opmpressure
+c05905415adff332a3df8533ebe055acd115443d dune/porsol/opmpressure
diff --git a/dune/porsol/mimetic/TpfaCompressible.hpp b/dune/porsol/mimetic/TpfaCompressible.hpp
index bd8073a..e8e544a 100644
--- a/dune/porsol/mimetic/TpfaCompressible.hpp
+++ b/dune/porsol/mimetic/TpfaCompressible.hpp
@@ -283,10 +283,13 @@ namespace Dune
                     }
                 }
 
+                std::vector<double> wellperfA, phasemobwellperf, wellperf_gpot;
+
                 // Assemble system matrix and rhs.
                 psolver_.assemble(src, bctypes_, bcvalues_, dt,
-                                  fp_.totcompr, initial_voldiscr, fp_.cellA, fp_.faceA, fp_.phasemobf,
-                                  cell_pressure_scalar, &(pfluid_->surfaceDensities()[0]));
+                                  fp_.totcompr, initial_voldiscr, fp_.cellA, fp_.faceA,
+                                  wellperfA, fp_.phasemobf, phasemobwellperf,
+                                  cell_pressure_scalar, wellperf_gpot, &(pfluid_->surfaceDensities()[0]));
                 // Solve system.
                 PressureSolver::LinearSystem s;
                 psolver_.linearSystem(s);