From a3f658f3b16c001f9b16bb0cdf091cbe5cebdc20 Mon Sep 17 00:00:00 2001 From: Tor Harald Sandve Date: Wed, 7 Jan 2015 10:32:17 +0100 Subject: [PATCH] Add warning to c solvers when stopped wells occur Stopped wells are treated as shut wells in the c solvers in opm-core. I.e they are completly isolated from the domain. --- opm/core/pressure/tpfa/cfs_tpfa_residual.c | 2 +- opm/core/pressure/tpfa/ifs_tpfa.c | 3 +-- 2 files changed, 2 insertions(+), 3 deletions(-) diff --git a/opm/core/pressure/tpfa/cfs_tpfa_residual.c b/opm/core/pressure/tpfa/cfs_tpfa_residual.c index 0d209611..c3c6040b 100644 --- a/opm/core/pressure/tpfa/cfs_tpfa_residual.c +++ b/opm/core/pressure/tpfa/cfs_tpfa_residual.c @@ -870,7 +870,7 @@ assemble_completion_to_well(int i, int w, int c, int nc, int np, ctrl = W->ctrls[ w ]; if (well_controls_well_is_stopped( ctrl )) { - /* Interpreting a negative current control index to mean a shut well */ + fprintf(stderr, "Stopped well detected: will be treated as completely shut\n"); welleq_coeff_shut(np, h, &res, &w2c, &w2w); } else { diff --git a/opm/core/pressure/tpfa/ifs_tpfa.c b/opm/core/pressure/tpfa/ifs_tpfa.c index 724c5b95..84f4ac0a 100644 --- a/opm/core/pressure/tpfa/ifs_tpfa.c +++ b/opm/core/pressure/tpfa/ifs_tpfa.c @@ -364,9 +364,8 @@ assemble_well_contrib(int nc , ctrls = W->ctrls[ w ]; if (well_controls_well_is_stopped(ctrls) ) { - + fprintf(stderr, "Stopped well detected: will be treated as completely shut\n"); /* Treat this well as a shut well, isolated from the domain. */ - assemble_shut_well(nc, w, W, mt, h); } else {