diff --git a/opm/core/pressure/tpfa/cfs_tpfa_residual.c b/opm/core/pressure/tpfa/cfs_tpfa_residual.c
index 02c6d254..cb3d8173 100644
--- a/opm/core/pressure/tpfa/cfs_tpfa_residual.c
+++ b/opm/core/pressure/tpfa/cfs_tpfa_residual.c
@@ -1090,7 +1090,7 @@ compute_wflux(int                        np    ,
 {
     int           w, c, i, p;
     double        pw, dp, t;
-    const double *pmob;
+    const double *pmob, *wdp;
 
     struct Wells          *W;
     struct CompletionData *cdata;
@@ -1101,13 +1101,14 @@ compute_wflux(int                        np    ,
 
     W     = wells->W;
     cdata = wells->data;
+    wdp   = cdata->wdp;
 
     for (w = i = 0; w < W->number_of_wells; w++) {
         pw = wpress[w];
 
         for (; i < W->well_connpos[w + 1]; i++) {
             c  = W->well_cells[ i ];
-            dp = pw - cpress[c];
+            dp = pw + wdp[ i ] - cpress[c];
 
             if (dp > 0) { pmob = cdata->phasemob + (i * np); } /* w->c */
             else        { pmob = pmobc           + (c * np); } /* c->w */