docs: def galToCel(ll, bb):
      """
      Converts Galactic (deg) to Celestial J2000 (deg) coordinates
      """
      bb = np.radians(bb)
      sin_bb = np.sin(bb)
      cos_bb = np.cos(bb)
  
      ll = np.radians(ll)
      ra_gp = np.radians(192.85948)
      de_gp = np.radians(27.12825)
      lcp = np.radians(122.932)
  
      sin_lcp_ll = np.sin(lcp - ll)
      cos_lcp_ll = np.cos(lcp - ll)
  
      sin_d = (np.sin(de_gp) * sin_bb) \
              + (np.cos(de_gp) * cos_bb * cos_lcp_ll)
      ramragp = np.arctan2(cos_bb * sin_lcp_ll,
                              (np.cos(de_gp) * sin_bb) \
                              - (np.sin(de_gp) * cos_bb * cos_lcp_ll))
      dec = np.arcsin(sin_d)
      ra = (ramragp + ra_gp + (2. * np.pi)) % (2. * np.pi)
      return np.degrees(ra), np.degrees(dec)