extinction.odonnell94¶
-
extinction.
odonnell94
(wave, a_v, r_v, unit='aa', out=None)¶ O’Donnell (1994) extinction function.
Like Cardelli, Clayton, & Mathis (1989) [1] but using the O’Donnell (1994) [2] optical coefficients between 3030 A and 9091 A.
Parameters: - wave : numpy.ndarray (1-d)
Wavelengths or wavenumbers.
- a_v : float
Scaling parameter, A_V: extinction in magnitudes at characteristic V band wavelength.
- r_v : float
Ratio of total to selective extinction, A_V / E(B-V).
- unit : {‘aa’, ‘invum’}, optional
Unit of wave: ‘aa’ (Angstroms) or ‘invum’ (inverse microns).
- out : np.ndarray, optional
If specified, store output values in this array.
Returns: - Extinction in magnitudes at each input wavelength.
Notes
This function matches the Goddard IDL astrolib routine CCM_UNRED. From the documentation for that routine:
- The CCM curve shows good agreement with the Savage & Mathis (1979) [3] ultraviolet curve shortward of 1400 A, but is probably preferable between 1200 and 1400 A.
- Curve is extrapolated between 912 and 1000 A as suggested by Longo et al. (1989) [4]
- Valencic et al. (2004) [5] revise the ultraviolet CCM curve (3.3 – 8.0 um^-1). But since their revised curve does not connect smoothly with longer and shorter wavelengths, it is not included here.
References
[1] Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [2] O’Donnell, J. E. 1994, ApJ, 422, 158O [3] Savage & Mathis 1979, ARA&A, 17, 73 [4] Longo et al. 1989, ApJ, 339,474 [5] Valencic et al. 2004, ApJ, 616, 912