ctools

import gammalib as gl

import argparse

#__________________________________________________________

# MAIN FUNCTION

if __name__=="__main__":

    # Get the command line options

    parser = argparse.ArgumentParser()

    parser.add_argument("--incoordsys", help="Input coordinate system <GAL|CEL>",

                        type=str)

    parser.add_argument("--xref", help="X Coordinate (RA or Galactic longitude)",

                        type=float)

    parser.add_argument("--yref", help="Y Coordinate (DEC ro GLAT)",

                        type=float)

    opts = parser.parse_args()

    # Create coordinate object

    coords = gl.GSkyDir()

    # Convert GAL -> CEL

    if opts.incoordsys == "GAL":

        coords.lb_deg( opts.xref, opts.yref )

        print("\nConverting GAL -> CEL:")

        print("   RA : {0}\n   DEC: {1}\n".format(coords.ra_deg(), coords.dec_deg()))

    # Convert CEL -> GAL

    elif opts.incoordsys == "CEL":

        coords.radec_deg( opts.xref, opts.yref )

        print("\nConverting CEL -> GAL:")

        print("   glon: {0}\n   glat: {1}\n".format(coords.l_deg(), coords.b_deg()))

    # Unknown coordinate system (or none) supplied

    else:

        print("\nERROR: incoordsys param must be either GAL or CEL.\n")