import matplotlib.pyplot as plt
from skycalc_ipy import SkyCalc
skycalc = SkyCalc()
skycalc.get_almanac_data(ra=83.8221, dec=-5.3911,
                         date="2017-12-24T04:00:00",
                         update_values=True)
tbl = skycalc.get_sky_spectrum()
plt.plot(tbl["lam"], tbl["flux"])
plt.xlabel("Wavelength " + str(tbl["lam"].unit))
plt.ylabel("Flux " + str(tbl["flux"].unit))
plt.semilogy()