D03

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

Bases: GMBase

Draine (2003) Grain Models

Parameters:
None
Raises:
None

Notes

From Draine (2003, ARA&A, 41, 241; 2003, ApJ, 598, 1017). Using Weingartner & Draine (2001, ApJ, 548, 296) size distributions

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 D03

fig, ax = plt.subplots()

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

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

for cmodel in possmodels:
    # define the extinction model
    ext_model = D03(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-D03-1.png

Attributes Summary

possnames

x_range

Attributes Documentation

possnames = {'MWRV31': ('kext_albedo_WD_MW_3.1_60_D03.dat', 3.1), 'MWRV40': ('kext_albedo_WD_MW_4.0A_40_D03.dat', 4.0), 'MWRV55': ('kext_albedo_WD_MW_5.5A_30_D03.dat', 5.5)}
x_range = [0.0001, 10000.0]