#!/usr/bin/env python
"""
Example Kitt Peak

./radec2azel.py 257.96295344 15.437854 31.9583 -111.5967 2014-12-25T22:00:00MST
"""
from pymap3d import radec2azel
from argparse import ArgumentParser


def main():
    p = ArgumentParser(description="RightAscension,Declination =>" "Azimuth,Elevation")
    p.add_argument("ra", help="right ascension [degrees]", type=float)
    p.add_argument("dec", help="declination [degrees]", type=float)
    p.add_argument("lat", help="WGS84 latitude of observer [degrees]", type=float)
    p.add_argument("lon", help="WGS84 latitude of observer [degrees]", type=float)
    p.add_argument("time", help="UTC time of observation YYYY-mm-ddTHH:MM:SSZ")
    P = p.parse_args()

    az_deg, el_deg = radec2azel(P.ra, P.dec, P.lat, P.lon, P.time)
    print("azimuth: [deg]", az_deg)
    print("elevation [deg]:", el_deg)


if __name__ == "__main__":
    main()