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