code: def pix2ang(nside, pix, nest=False):
      
      theta, phi =  hp.pix2ang(nside, pix, nest=nest)
      lon = phi2lon(phi)
      lat = theta2lat(theta)
      return lon, lat
docs: Return (lon, lat) in degrees instead of (theta, phi) in radians