barak.absorb.calc_iontau¶

barak.absorb.calc_iontau(wa, ion, zp1, logN, b, debug=False, ticks=False, maxdv=1000.0, label_tau_threshold=0.01, vpad=500.0, verbose=True)[source]¶

Returns tau values at each wavelength for transitions in ion.

Parameters :

wa : array of floats wavelength array ion : atom.data entry ion entry from readatom output dictionary zp1 : float redshift + 1 logN : float log10(column density in cm**-2) b : float b parameter (km/s). Assumes thermal broadening. maxdv : float (default 1000) For performance reasons, only calculate the Voigt profile for a single line to +/- maxdv. Increase this if you expect DLA-type extended wings. None for no maximum. vpad : float (default 500) Include transitions that are within vpad km/s of either edge of the wavelength array.

Returns :

tau : array of floats Array of optical depth values.