code: def nside2pixarea(nside, degrees=False):
      
      area = nside_to_pixel_area(nside)
      if degrees:
          return area.to(u.deg ** 2).value
      else:
          return area.to(u.sr).value
docs: Drop-in replacement for healpy `~healpy.pixelfunc.nside2pixarea`.