Merge pull request #4841 from GitPaean/set_bhp_update_well_state_2

setting the bhp to be bhp_limit when bhp value is not initialized with rate target control when updateWellStateWithTarget()

opm/simulators/wells/WellInterface_impl.hpp

@@ -874,6 +874,10 @@ namespace Opm
             // within the wellbore from the previous result, and hopefully it is a good
             // initial guess for the zero rate target.
             ws.surface_rates[phasePos] = std::max(1.e-7, ws.surface_rates[phasePos]);
+
+            if (ws.bhp == 0.) {
+                ws.bhp = controls.bhp_limit;
+            }
         }
         //Producer
         else
@@ -1096,6 +1100,10 @@ namespace Opm
 
                 break;
             } // end of switch
+
+            if (ws.bhp == 0.) {
+                ws.bhp = controls.bhp_limit;
+            }
         }
     }