adding variableStateWellIndices to StandardWells

to handle different types of variables later for different types of
wells.

opm/autodiff/BlackoilModelBase_impl.hpp

@@ -604,8 +604,7 @@ namespace detail {
         if (active_[Gas]) {
             indices[Xvar] = next++;
         }
-        indices[Qs] = next++;
+        asImpl().stdWells().variableStateWellIndices(indices, next);
-        indices[Bhp] = next++;
         assert(next == fluid_.numPhases() + 2);
         return indices;
     }

opm/autodiff/StandardWells.hpp

@@ -180,6 +180,11 @@ namespace Opm {
                                   WellState& well_state);
 
 
+            void
+            variableStateWellIndices(std::vector<int>& indices,
+                                     int& next) const;
+
+
 
         protected:
             bool wells_active_;

opm/autodiff/StandardWells_impl.hpp

@@ -1118,4 +1118,16 @@ namespace Opm
 
     }
 
+
+
+
+
+    void
+    StandardWells::variableStateWellIndices(std::vector<int>& indices,
+                                            int& next) const
+    {
+        indices[Qs] = next++;
+        indices[Bhp] = next++;
+    }
+
 }