docs: def ang2disc(nside, lon, lat, radius, inclusive=False, fact=4, nest=False):
      """
      Wrap `query_disc` to use lon, lat, and radius in degrees.
      """
      vec = ang2vec(lon,lat)
      return query_disc(nside,vec,radius,inclusive,fact,nest)