WD01

class dust_extinction.grain_models.WD01(modelname='MWRV31', **kwargs)[source]

Bases: GMBase

Weingartner & Draine (2001) Grain Models

Parameters:
None
Raises:
None

Notes

From Weingartner & Draine (2001, ApJ, 548, 296)

Example showing the possible curves

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

from dust_extinction.grain_models import WD01

fig, ax = plt.subplots()

tmod = WD01()
possmodels = tmod.possnames.keys()

# generate the curves and plot them
lam = np.logspace(-2.0, 4.0, num=1000)
x = (1.0 / lam) / u.micron

for cmodel in possmodels:
    # define the extinction model
    ext_model = WD01(cmodel)
    ax.plot(lam,ext_model(x),label=cmodel)

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

ax.set_xscale('log')
ax.set_yscale('log')

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

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

../_images/dust_extinction-grain_models-WD01-1.png

Attributes Summary

possnames

x_range

Attributes Documentation

possnames = {'LMC2': ('kext_albedo_WD_LMC2_10.dat', 3.1), 'LMCAvg': ('kext_albedo_WD_LMCavg_20.dat', 3.1), 'MWRV31': ('kext_albedo_WD_MW_3.1_60.dat', 3.1), 'MWRV40': ('kext_albedo_WD_MW_4.0B_40.dat', 4.0), 'MWRV55': ('kext_albedo_WD_MW_5.5B_30.dat', 5.5), 'SMCBar': ('kext_albedo_WD_SMCbar_0.dat', 3.1)}
x_range = [0.0001, 100.0]