docs: def _compute_distance(self, rup, dists, C):
          """
          Compute the distance function, equation (5).
          """
          rval = np.sqrt(dists.repi ** 2 + C['h'] ** 2)
          return C['c1'] * np.log10(rval)