Internal function to optBH
getBH(pars, SSB, rec, SSBpR, mode = 1, plot = F)
an initial guess at model parameters steepness and R0
'observations' of spawning biomass
'observations' (model predictions) of recruitment
spawning stock biomass per recruit at unfished conditions
should fit or recruitment deviations be returned
should a plot of the model fit be produced?