docs: def _get_pga_bc(self, f0, f1, f2, SC, mag, rrup, vs30, scale_fac):
          """
          Compute and return PGA on BC boundary
          """
          pga_bc = np.zeros_like(vs30)
          self._compute_mean(self.COEFFS_BC[PGA()], f0, f1, f2, SC, mag,
                             rrup, vs30 < 2000.0, pga_bc, scale_fac)
  
          return (10 ** pga_bc) * 1e-2 / g