docs: def _compute_anelestic_attenuation_term(self, C, dists):
          """
          Compute and return anelastic attenuation term in equation 5,
          page 970.
          """
          f_aat = np.zeros_like(dists.rjb)
          idx = dists.rjb > 80.0
          f_aat[idx] = C["b10"] * (dists.rjb[idx] - 80.0)
          return f_aat