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)
# skycalc["msolflux"] = 130       # [sfu] For dates after 2019-01-31
tbl = skycalc.get_sky_spectrum()

plt.semilogy(tbl["lam"], tbl["flux"])
plt.xlabel(f"Wavelength [{tbl['lam'].unit}]")
plt.ylabel(f"Flux [{tbl['flux'].unit}]")
plt.grid(True, ls=":")