From fb95c6da191c7a9ec09c28b557e91da8ba0a4a53 Mon Sep 17 00:00:00 2001 From: Kjetil Olsen Lye Date: Fri, 13 Apr 2012 14:22:44 +0200 Subject: [PATCH] Added computations for total flow for each well --- examples/wells_example.cpp | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/examples/wells_example.cpp b/examples/wells_example.cpp index 09329be4e..44201ff54 100644 --- a/examples/wells_example.cpp +++ b/examples/wells_example.cpp @@ -73,9 +73,12 @@ int main(int argc, char** argv) { std::vector face_flux; std::vector well_bhp; - std::vector well_rate; - pressure_solver.solve(totmob, omega, src, wdp, bcs.c_bcs(), pressure, face_flux, well_bhp, well_rate); + std::vector well_rate_per_cell; + pressure_solver.solve(totmob, omega, src, wdp, bcs.c_bcs(), pressure, face_flux, well_bhp, well_rate_per_cell); std::cout << "Solved" << std::endl; + + std::vector well_rate; + computeFlowRatePerWell(*wells.c_wells(), well_rate_per_cell, well_rate); if(wells.wellCollection().conditionsMet(well_bhp, well_rate, *grid.c_grid(), state.saturation() )) { std::cout << "Conditions met for wells!" << std::endl; }