docs: def get_equivalent_distance_inslab(mag, repi, hslab):
      """
      :param float mag:
          Magnitude
      :param repi:
          A :class:`numpy.ndarray` instance containing repi values
      :param float hslab:
          Depth of the slab
      """
      area = 10**(-3.225+0.89*mag)
      radius = (area / scipy.constants.pi)**0.5
      rjb = np.max([repi-radius, np.zeros_like(repi)], axis=0)
      rrup = (rjb**2+hslab**2)**0.5
      return rjb, rrup