The program double.c runs using an input file ad.dat with the format:

---------------------------
table name with mu/muo data
number of points in table
number of layers
for each layer:
  layer type
 *g 
  w_haze
  tau_rayleigh
  tau_haze
  tau_gas
---------------------------

*For layer type 2 (double H-G) g is replaced by g1, g2, and f

I've included all tables, named xxxx.yydat, xxxx.yy dat, where xxxx is the
wavelength and yy the south latitude.

The programs makexxxx.ind.pro are used to create the format files to run the
code and integrate over the filter bandpasses xxxx.table.new while
integrating the result into a file called outfile.  They use precomputed
optical depth files tau.dat.xxxx and haze files HAZEXXXX that were computed
with f_hze=20.  You'll need to change the word "double" in each of these to
whatever you call the executable.

The programs slicexxxx.pro are used to create figures 4-7 using the outfiles.

nrutil.c, nrutil.h, and spline.c all come from Numerical Recipes in C
(http://www.nrbook.com/a/bookcpdf.php).