import astropy
import astropy.units as u
# this is to convert the flux denisty from erg / s / cm^2 / Hz (that we have) to Jy, by using astropy.units to define the 10^23 "annoying constant factor" that appears somewhere in the code
SpectralFluxDensity2Jy = u.Jy.to(u.erg / u.cm**2 / u.s / u.Hz, 1., equivalencies=u.spectral_density) ** -1