code: def ang2pix(nside, lon, lat, nest=False):
      
      theta = np.radians(90. - lat)
      phi = np.radians(lon)
      return healpy.ang2pix(nside, theta, phi, nest=nest)
docs: Input (lon, lat) in degrees instead of (theta, phi) in radians