2009 May 21

Made a few changes so code works in gfortran on a mac.  This nice little program
has been sorely neglected by me.  I'll try to update more often.

2004 December 8

Modified most of the plotting routines.  The voigt point-spread-function (PSF)
profile had a bug in it.  Now it properly has two variables, the gaussian 
width and the wing component.  (see the IDL help pages for "voigt" or p. 291 
of "Radiative Processes in Astrophysics" by Rybicki & Lightman (1979).  
Also I decided it is better to plot the contours from the peak downward 
rather than a minimum value upward.  This makes it easier to compare to 
observations I think.  so you'll have a couple of extra variables to modify: 
number of contour levels, gaussian width, voigt wing component, and number 
of gaussian widths in the psf.   And there is an extra routine called 
contpsf which plots a point source and shows how to vary the voigt routine
to match your data.   This is important especially if you are doing 
polarization maps where the wings make a difference (see Whitney et al. 1997).

Finally, the plotdens routine has been updated and is easier to use, doesn't
require running a fortran code to regrid the density.

2004 June 20

Modified the IDL programs to use the astrolib convolve function
instead of IDL's convol.  This is much faster.  Also improved
the voigt convolution routines; and included a gaussian PSF 
generator.  thanks to Marshall Perrin for the convolve suggestion.

2003 March 17

1.  forgot to change one of the input files (rundisk/ttsscat.par) to
reflect changes in variables read in (CSHAPE).

2003 March 01

no new bugs.  

1.  changed the definition of EXF so density goes as  
	rho=rhoconst*r^(-exf); previously went as r^exf.
2.  changed description of EX1 and EX2 in input files.