Now accepts both BHP and SURFACE_RATE wells.
Change example to use a surface rate constraint for one well.