From eca28a8ade00252547c7e09bd78735c0d2be88ca Mon Sep 17 00:00:00 2001
From: Kai Bao <kai.bao@sintef.no>
Date: Thu, 29 Jun 2017 14:31:17 +0200
Subject: [PATCH] adding a function outputWellState for StandardWellsDense

to output the well state for debugging purpose.
---
 opm/autodiff/StandardWellsDense.hpp      |  2 ++
 opm/autodiff/StandardWellsDense_impl.hpp | 27 ++++++++++++++++++++++++
 2 files changed, 29 insertions(+)

diff --git a/opm/autodiff/StandardWellsDense.hpp b/opm/autodiff/StandardWellsDense.hpp
index 94555f1e9..b5f0f010d 100644
--- a/opm/autodiff/StandardWellsDense.hpp
+++ b/opm/autodiff/StandardWellsDense.hpp
@@ -425,6 +425,8 @@ enum WellVariablePositions {
 
             void computeAverageFormationFactor(Simulator& ebosSimulator,
                                                std::vector<double>& B_avg) const;
+
+            void outputWellState(const WellState& well_state) const;
         };
 
 
diff --git a/opm/autodiff/StandardWellsDense_impl.hpp b/opm/autodiff/StandardWellsDense_impl.hpp
index 1eb7d4d7e..7279a5617 100644
--- a/opm/autodiff/StandardWellsDense_impl.hpp
+++ b/opm/autodiff/StandardWellsDense_impl.hpp
@@ -2511,4 +2511,31 @@ namespace Opm {
 
 
 
+    template<typename TypeTag>
+    void
+    StandardWellsDense<TypeTag>::
+    outputWellState(const WellState& well_state) const
+    {
+        std::cout << " output the bhp " << std::endl;
+        for (const double bhp : well_state.bhp()) {
+            std::cout << bhp << " ";
+        }
+        std::cout << std::endl;
+
+        std::cout << " output the well rates " << std::endl;
+        for (const double rate : well_state.wellRates()) {
+            std::cout << rate << " ";
+        }
+        std::cout << std::endl;
+
+        std::cout << " output the wellSolutions " << std::endl;
+        for (const double solution : well_state.wellSolutions()) {
+            std::cout << solution << " ";
+        }
+        std::cout << std::endl;
+    }
+
+
+
+
 } // namespace Opm