import matplotlib.pyplot as plt
from skycalc_ipy import SkyCalc

tbl = SkyCalc().get_sky_spectrum()

plt.figure(figsize=(8,8))
plt.subplot(211)
plt.plot(tbl["lam"], tbl["trans"])
plt.xlabel("Wavelength [um]")
plt.ylabel("Transmission")
plt.subplot(212)
plt.plot(tbl["lam"], tbl["flux"])
plt.xlabel("Wavelength [um]")
plt.ylabel("Emission [ph s-1 m-2 um-1 arcsec-2]")
plt.semilogy()