RL85_MWGC¶

class dust_extinction.averages.RL85_MWGC(*args, meta=None, name=None, **kwargs)[source]¶

Bases: BaseExtModel

Reike & Lebofsky (1985) MW Average Extinction Curve

Parameters:

None

Raises:

None

Notes

From Rieke & Lebofsky (1985, ApJ,288, 618)

Example showing the average curve

import numpy as np
import matplotlib.pyplot as plt
import astropy.units as u

from dust_extinction.averages import RL85_MWGC

fig, ax = plt.subplots()

# define the extinction model
ext_model = RL85_MWGC()

# generate the curves and plot them
x = np.arange(1.0/ext_model.x_range[1], 1.0/ext_model.x_range[0], 0.1) * u.micron

ax.plot(x,ext_model(x),label='RL85_MWGC')
ax.plot(1.0/ext_model.obsdata_x, ext_model.obsdata_axav, 'ko',
        label='obsdata')

ax.set_xlabel(r'$\lambda$ [$\mu m$]')
ax.set_ylabel(r'$A(x)/A(V)$')

ax.legend(loc='best')
plt.show()

(Source code, png, hires.png, pdf)

../_images/dust_extinction-averages-RL85_MWGC-1.png

Attributes Summary

`Rv`
`obsdata_axav`
`obsdata_tolerance`
`obsdata_x`
`x_range`

Methods Summary

evaluate(x)

RL85 MWGC function

Attributes Documentation

Rv = 3.09¶

obsdata_axav = array([0.027, 0.03 , 0.037, 0.047, 0.06 , 0.074, 0.083, 0.087, 0.074, 0.043, 0.02 , 0.023, 0.058, 0.112, 0.175, 0.282])¶

obsdata_tolerance = 1e-06¶

obsdata_x = array([0.07692308, 0.08 , 0.08333333, 0.08695652, 0.09090909, 0.0952381 , 0.1 , 0.10526316, 0.11111111, 0.11764706, 0.125 , 0.20833333, 0.28571429, 0.45045045, 0.60606061, 0.8 ])¶

x_range = [0.07692307692307693, 0.8]¶

Methods Documentation

classmethod evaluate(x)[source]¶

RL85 MWGC function

Parameters:

x: float

expects either x in units of wavelengths or frequency or assumes wavelengths in wavenumbers [1/micron]

internally wavenumbers are used

Returns:

axav: np array (float): A(x)/A(V) extinction curve [mag]

Raises:

ValueError: Input x values outside of defined range

Navigation

RL85_MWGC¶