Dust Emissivity toolkit

Only API docs for now

Reference/API

dust_emissivity Package

This is an Astropy affiliated package.

Functions

test([package, test_path, args, plugins, ...]) Run the tests using py.test.

dust_emissivity.dust Module

Functions

colofsnu(nu, snu_per_beam[, temperature, muh2])
kappa(nu[, nu0, kappa0, beta]) Compute the opacity $kappa$ given a reference frequency (or wavelength)
kappa_table(nu[, url]) Use an opacity table to compute kappa
massofsnu(nu, snu[, distance, temperature, ...])
planck_average_kappa(temperature[, kappa]) Compute the Blackbody-weighted opacity assuming a power law dust
planck_average_kappa_table(temperature[, ...])
quad(func, a, b[, args, full_output, ...]) Compute a definite integral.
snu(nu, column, kappa, temperature) Compute the flux density for a given column of gas assuming some opacity kappa
snudnu(nu, column, kappa, temperature, bandwidth)
snuofmass(nu, mass[, beamomega, distance, ...]) Flux density for a given mass and beam area
tauofsnu(nu, snu_per_beam[, temperature]) nu in GHz

dust_emissivity.blackbody Module

Includes both wavelength and frequency blackbody functions. Has flexible units. Also allows for a few varieties of modified blackbody.

Functions

blackbody(nu, temperature[, outunit]) Planck’s Law Blackbody (Frequency units)
blackbody_wavelength(lam, temperature[, outunit])
integrate_sed(vmin, vmax[, function]) Integrate one of the SED functions over frequency
modified_blackbody(nu, temperature[, beta, ...]) Snu = 2hnu^3 c^-2 (e^(hnu/kT) - 1)^-1 (1 - e^(-tau_nu) ) Kappa0 and Nu0 are set as per http://arxiv.org/abs/1101.4654 which uses OH94 values.

dust_emissivity.fit_sed Module

Functions

blackbody(nu, temperature[, outunit]) Planck’s Law Blackbody (Frequency units)
fit_modified_bb(xdata, flux, error, guesses) Wrapper SED fitter that uses astropy’s units and returns the fitted
fit_modifiedbb_montecarlo(frequency, flux[, ...]) An MCMC version of the fitter.
fit_sed_lmfit_hz(xdata, flux[, guesses, ...])
Parameters:
fit_sed_mpfit_hz(xdata, flux[, guesses, ...])
Parameters:
modified_blackbody(nu, temperature[, beta, ...]) Snu = 2hnu^3 c^-2 (e^(hnu/kT) - 1)^-1 (1 - e^(-tau_nu) ) Kappa0 and Nu0 are set as per http://arxiv.org/abs/1101.4654 which uses OH94 values.