docs: def get_ruptures_within(dstore, bbox):
      """
      Extract the ruptures within the given bounding box, a string
      minlon,minlat,maxlon,maxlat.
      Example:
      http://127.0.0.1:8800/v1/calc/30/extract/ruptures_with/8,44,10,46
      """
      minlon, minlat, maxlon, maxlat = map(float, bbox.split(','))
      hypo = dstore['ruptures']['hypo'].T  # shape (3, N)
      mask = ((minlon <= hypo[0]) * (minlat <= hypo[1]) *
              (maxlon >= hypo[0]) * (maxlat >= hypo[1]))
      return dstore['ruptures'][mask]