On the computation of planetary boundary-layer height using the bulk Richardson number method