The files for the the 3D Cartesian grid code illuminated by point sources are
available in this
gzipped tar file.
The current setup reproduces the example
circumbinary disk scattered light and polarization images described in the
booklet and shown below.
The FORTRAN files are:
density.f
gridset.f
iarray.f
mcpolar.f
peeloff.f
ran2.f
scatt1.f
scattp.f
stokes.f
sourceph.f
sources.f
taufind1.f
tauint.f
tauint2.f
Include files are:
grid.txt
images.txt
photon.txt
sources.txt
Input parameters are in:
params.par
The file that compiles the code and creates the executable file 'mcgrid' is:
Makefile
The code produces 2D images in unformatted f77 files which can be viewed
with a simple IDL routine included with the codes. If you do not use IDL,
Tom Robitaille has provided a Python script that will convert the images into
FITS files. The script is
available here.
Below is example ouput from this code. The sets of three
images are the scattered
flux (stretched to bring up the faint scattered light disk), and the Stokes
Q and U flux images.
Kenny Wood
Code and gzipped tar file last updated: 8 July 2013