************************************************
README - what are these files in this directory?
************************************************

oldprograms/
  Old fortran routines since the days of nk2q
  need to take a look someday...

dorschner95_olmg50.optc
  Dorschner, J.; Begemann, B.; Henning, T.; Jaeger, C.; Mutschke, H.
  1995A&A...300..503D
  Olivine Fe 50% Mg 50% 

  Note: original data: 0.2um to 500.0um
        long wavelength end extrapolated to 4000.0um

dorschner95_olmg40.optc
  Dorschner, J.; Begemann, B.; Henning, T.; Jaeger, C.; Mutschke, H.
  1995A&A...300..503D
  Olivine Fe 60% Mg 40%

  Note: original data: 0.2um to 500.0um
        long wavelength end extrapolated to 4000.0um

begemann97_porousalmina.optc
  1.2um > l         : Harman et al. 1994, J. Appl. Phys., 76, 8032
  7.8um > l > 1.2um : Harris, 1955, J. Opt. Soc. Am., 45, 27
          l > 7.8um : Begemann et al. 1997, ApJ, 476, 199

H2O_franktest.optc
  Molster, F.
  2001
  Water Ice

  Note: original data: 0.0443um to 333.33um
        long wavelength end extrapolated to 500.0um

frankfors1.optc
  Molster, F. 
  2001
  Forsterite

  Note: original data: 0.0427um to 125um
        long wavelength end extrapolated to 500.0um

frankfors2.optc
  Molster, F. 
  2001
  Forsterite

  Note: original data: 0.0427um to 125um
        long wavelength end extrapolated to 500.0um

frankfors3.optc
  Molster, F. 
  2001
  Forsterite

  Note: original data: 0.0427um to 125um
        long wavelength end extrapolated to 500.0um

ens_p_am.optc
  Molster, F.
  2001
  Enstatite

  Note: original data: 0.0427um to 98.913um
        long wavelength end extrapolated to 500.0um

ens_s1_am.optc
  Molster, F.
  2001
  Enstatite

  Note: original data: 0.0427um to 98.913um
        long wavelength end extrapolated to 500.0um

ens_s2_am.optc
  Molster, F.
  2001
  Enstatite

  Note: original data: 0.0427um to 98.913um
        long wavelength end extrapolated to 500.0um

rouleau91_ac.optc
  Rouleau, Francois; Martin, P. G.
  1991ApJ...377..526R
  Amorphous Carbon

  Note: 0.000355179um to 301.218um
        These values are for the "BE" type.

rouleau+martin.dat
  The original table from the above ref. in the electronic form.

zubko96_ac_acar.optc
zubko96_ac_ach2.optc
zubko96_ac_be.optc
  Zubko, V. G.; Mennella, V.; Colangeli, L.; Bussoletti, E.
  1996MNRAS.282.1321Z
  Amorphous Carbon

_acar is amorphous carbon in Argon
_ach2 is amorphous carbon in H2
_be is burning benzene

  Note: micron, n, and K
        ~0.04um to 1000 to 2000um

cosil.lnk    
pyrmg40.lnk  
pyrmg50.lnk
  Jena Group N&K values
  (http://www.astro.uni-jena.de/Group/Subgroups/Labor/Labor/silicates.html)

  cosil.lnk:  values usually a little higher than DL for lambda < 12um,
              and lower than DL for lambda > 12um

	      not very different than the two below

  pyrmg40.lnk:  slightly higher than DL  at short lambda,  lower at long

  pyrmg50.lnk:  slightly higher than DL at short lambda,  lower at long.,
		not very different from pyrmg40.lnk

draine
  This directory contains n&k values of 
    graphite,
    silicate 
    SiC grain constants
  by Draine and collaborators dated 1995/10/25.

  See the README file in the directory for contents.

  Bascially for our purposes, we want the 

  eps_suvSil  file  which contains the "smoothed UV" silicate grain constants

  the five columns are:
  1) wavelength in microns
  2) real dielectric constant  minus one
  3) imaginary dielectric constant
  4) real refractive index minus one
  5) imaginary refractive index

  We want to use the columns 1, 4 and 5

  Need to edit the file,  use awk to add the one
  and then use reverse to reverse the wavelength order

  awk '{print($1,$4+1,$5)}' eps_suvSil > dwsuvSil.optical

  These numbers are in better agreement with the Jena values for silicates

dwsuvSil.optical_rev.optc
  Weingartner, J.C. & Draine, B.T. 1999, ApJ in prep.

  Note: this seems to be the best "astronomical silicates" values.

preibisch93_dice.optc
  Preibisch, Th.; Ossenkopf, V.; Yorke, H. W.; Henning, Th.
  1993A&A...279..577P

  Note: dirty ice; 0.1 - 800um

..................................................................

from www.astro.uni-jena.de/Laboratory/Database/carbon.html

Schnaiter et al, ``Matrix-isolated nano-sized carbon grains as an
analog for the 217.5 nm feature carrier'', Astrophys. J. 498, 486

"from matrix-isolation spectroscopy of non-agglomerated carbon
particles (Schnaiter et al. 1998). This unique experiment combined a
molecular-beam method for extraction of carbon grains out of their
condensation zone with the isolation of the particles in rare-gas
matrices. This allowed UV spectroscopy at non-agglomerated grains and
the application of simple Mie theory for evaluation of optical data. 
Data are available for carbon particles condensed in:

pure argon atmosphere -- schnaiter_ac_ar.optc
2/3 argon, 1/3 hydrogen -- schnaiter_ac_arh2.optc

data only from 0.2 to 1.0 microns wavelength...


====================================================================
Updated by Megan in June 2007

The original files only went out to 500um, I've interpolated and 
extrapolated them out to 4000, you can cut them short if you need fewer wavelengths

Included in this directory or tar pack is the python code I used
to extrapolate them (interp.py). It's options are as follows:

> interp.py --help
usage: interp.py [options]

options:
  -h, --help            show this help message and exit
  
  -i INFILE, --infile=INFILE
                        Input filename
  -o OUTFILE, --outfile=OUTFILE
                        Output filename
  -w ENDWAVE, --wave=ENDWAVE
                        Ending wavelength
 
  -l, --nolimit         Dont limit the n&k values as  in the Chiang paper
 
  -f FITPOINTS, --fitpoints=FITPOINTS
                        Number of points at the end for the fit
  -s STARTVAL, --startval=STARTVAL
                        Wavelength to start n&k limits at

The n&k values are limited using the specification in Chiang, E. I. et al. 2001, ApJ, 547, 1077
in order to keep k above zero. The paper uses a wavelength of 500um as the cutoff,
but I've coded interp.py such that you can specify a different cutoff if you
have actual data that extends past that value.

The code also produces a PDF file with plots of the original data and the
extrapolated data for reference. If you want to produce a different format
plot (say jpg) just change the output name of the file in the code.
Also, make sure you have python and the packages specified in the top
of the code loaded on your computer in order to run the full code.



