import matplotlib.pyplot as plt
from skycalc_ipy import SkyCalc

tbl = SkyCalc().get_sky_spectrum()

fig, ax = plt.subplots(2, 1, sharex=True, figsize=(8,8))
fig.subplots_adjust(hspace=0)
ax[0].plot(tbl["lam"], tbl["trans"], lw=1)
ax[1].semilogy(tbl["lam"], tbl["flux"], lw=1)
ax[0].set_ylabel("Transmission")
ax[1].set_ylabel(f"Emission [{tbl['flux'].unit}]")
ax[1].set_xlabel(f"Wavelength [{tbl['lam'].unit}]")
ax[0].grid(True, ls=":")
ax[1].grid(True, ls=":")