pro make_tables
; generate table of densities
ntab=501
nbar=10.^(findgen(ntab)/(ntab-1)*10)
dir = 'TOOLS/krumholz/'

; set gas temperature
t=10.0

; c18o 2->1
jc18o=fltarr(ntab)
dummy=emissivity(100., 10., 2, 1, dir + 'c18o.dat')
for n=0,ntab-1 do jc18o[n]=emissivity(nbar[n], t, 2, 1, /same_mol)
absorptivityc18o=fltarr(ntab)
for n=0,ntab-1 do absorptivityc18o[n]=absorptivity(nbar[n], t, 2, 1, /same_mol)

; n2h+ 1->0
jn2hplus=fltarr(ntab)
dummy=emissivity(100., 10., 1, 0, dir + 'n2hplus.dat')
for n=0,ntab-1 do jn2hplus[n]=emissivity(nbar[n], t, 1, 0, /same_mol)
absorptivityn2hplus=fltarr(ntab)
for n=0,ntab-1 do absorptivityn2hplus[n]=absorptivity(nbar[n], t, 2, 1, /same_mol)

; nh3 1-1
jnh3_1=fltarr(ntab)
dummy=emissivity(100., 10., 1, 0, dir + 'p-nh3.dat')
for n=0,ntab-1 do jnh3_1[n]=emissivity(nbar[n], t, 1, 0, /same_mol)
absorptivitynh3_1=fltarr(ntab)
for n=0,ntab-1 do absorptivitynh3_1[n]=absorptivity(nbar[n], t, 1, 0, /same_mol)

openw, fp, 'moltab.txt', /get_lun
for n=0,ntab-1 do printf, fp, nbar[n], jc18o[n], jn2hplus[n], $
	jnh3_1[n]
free_lun, fp

openw, fp, 'kappatab.txt', /get_lun
for n=0,ntab-1 do printf, fp, nbar[n], absorptivityc18o[n], $
	absorptivityn2hplus[n], absorptivitynh3_1[n]
free_lun, fp

end