From 7e1e1e795138e58efea48f75e20ae7ee92dc2175 Mon Sep 17 00:00:00 2001 From: Xavier Raynaud Date: Tue, 21 Aug 2012 10:20:59 +0200 Subject: [PATCH] Fixed sign error and replaced porevol->porevol0. --- opm/core/pressure/tpfa/cfs_tpfa_residual.c | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/opm/core/pressure/tpfa/cfs_tpfa_residual.c b/opm/core/pressure/tpfa/cfs_tpfa_residual.c index ab79dbce7..b3ca1f9d8 100644 --- a/opm/core/pressure/tpfa/cfs_tpfa_residual.c +++ b/opm/core/pressure/tpfa/cfs_tpfa_residual.c @@ -1247,7 +1247,7 @@ cfs_tpfa_res_comprock_assemble( /* Assemble usual system (without rock compressibility). */ cfs_tpfa_res_assemble(G, dt, forces, zc, cq, trans, gravcap_f, - cpress, wpress, porevol, h); + cpress, wpress, porevol0, h); /* Check if we have only Neumann wells. */ well_is_neumann = 1; @@ -1276,7 +1276,7 @@ cfs_tpfa_res_comprock_assemble( } h->J->sa[j] += porevol[c] * rock_comp[c] / dt; - h->F[c] -= dpvdt; + h->F[c] += dpvdt; } /* Re-do the singularity-removing adjustment if necessary */