From 5664fcd7d7ba22faf125f6641929e056be4c1cd2 Mon Sep 17 00:00:00 2001
From: Kjetil Olsen Lye <kjetil.olsen.lye@sintef.no>
Date: Thu, 12 Apr 2012 16:39:12 +0200
Subject: [PATCH] Made us of wells output arguments in IncompTpfa

---
 opm/core/pressure/IncompTpfa.cpp | 14 ++++++++++++++
 1 file changed, 14 insertions(+)

diff --git a/opm/core/pressure/IncompTpfa.cpp b/opm/core/pressure/IncompTpfa.cpp
index 540c1ebd..cd371bc2 100644
--- a/opm/core/pressure/IncompTpfa.cpp
+++ b/opm/core/pressure/IncompTpfa.cpp
@@ -25,6 +25,7 @@
 #include <opm/core/linalg/LinearSolverInterface.hpp>
 #include <opm/core/linalg/sparse_sys.h>
 #include <opm/core/utility/ErrorMacros.hpp>
+#include <opm/core/newwells.h>
 
 namespace Opm
 {
@@ -135,6 +136,13 @@ namespace Opm
         soln.cell_press = &pressure[0];
         soln.face_flux  = &faceflux[0];
 
+        if(wells_ != NULL) {
+            well_bhp.resize(wells_->number_of_wells);
+            well_rate.resize(wells_->number_of_wells);
+            soln.well_flux = &well_rate[0];
+            soln.well_press = &well_bhp[0];
+        }
+        
         ifs_tpfa_press_flux(gg, &F, &trans_[0], h_, &soln);
     }
 
@@ -224,6 +232,12 @@ namespace Opm
         soln.cell_press = &pressure[0];
         soln.face_flux  = &faceflux[0];
 
+        if(wells_ != NULL) {
+            well_bhp.resize(wells_->number_of_wells);
+            well_rate.resize(wells_->number_of_wells);
+            soln.well_flux = &well_rate[0];
+            soln.well_press = &well_bhp[0];
+        }
         ifs_tpfa_press_flux(gg, &F, &trans_[0], h_, &soln);
     }