add functions to find the column density per unit wavelength/velocity interval from the optical depth