from msaexp import spectrum
import msaexp.utils

result = spectrum.multiplot_spectra(
    file="smacs0723-ero-v4_g395m-f290lp_2736_6355.spec.fits",
    sx=8, sy=2,
    ny=4,
)

# Overplot line list
li = msaexp.utils.lines.LineList()
for ax in result['axes']:
    li.add_to_axis(ax, alpha=0.5)

result['axes'][-1].set_xlabel(
    f"rest wavelength, z={result['z']:.3f}"
)

result['fig'].tight_layout(pad=1)