we definitely should use quantities here: wave = Quantity(wave, u.angstrom, u.spectral()), ebv = Quantity(ebv, u.mag), ...