Code

Collisional Grooming Algorithm

This algorithm uses the output of a collisionless disk simulation to iteratively solve the mass flux equation for the density distribution of a collisional disk in three dimensions. The algorithm can be run on a single processor in just a few hours.  This code has not yet been polished for public release.  For details on how this algorithm works, see this paper.

dustmap v3.0

dustmap synthesizes images and histograms of debris disks using discrete data points as input.  When creating your debris disk model, simply save the positions of the dust grains to file in one of the required binary formats.  Then use dustmap to read in the data and produce scattered light and thermal emission images.

dustmap v3.0 is a hybrid of IDL and C code.  The computationally intense work is done in C, so it’s much faster, and the routine is called in IDL, so it’s still user-friendly.  A lot of new features were added, including:

• Ability to quickly create datacubes of images at many wavelengths
• Easier all-sky maps with Aitoff projections
• Faster and more accurate grain temperature calculations
• Ability to create optical depth histograms
• Expanded library of optical constants calculated on-the-fly via Mie scattering theory

Note: some outdated features were removed and the installation procedure has changed.  Please see the README file for complete details.

dustmap v2.2

An older, all-IDL version of dustmap.  This version is no longer updated and is significantly slower than v3.0.  However, if you are C-phobic and need to edit the code, this may be useful.

dmmovie v1.0

dmmovie is an IDL suite used along with dustmap to produce frames for movies/animations of debris disks.  Simply define the motion of the camera with user-friendly keywords, and sit back while the animation frames are generated.  Requires dustmap v3.0 or later.

kowalsky.pro

An IDL routine that deprojects an observed on-sky ellipse (e.g. debris disk ring, binary orbit, etc.) using the Kowalsky method.