Parametric bootstrap goodness-of-fit testing for Wehrly–Johnson bivariate circular distributions