pro allopacities

	openr, lun, 'temp_pelog.dat', /get_lun
	openw, 5, 'opacities.dat'
	fmt = '(1x, f7.2, 2x, f7.4, 3(2x, f8.5))'
	while (eof(lun) ne 1) do begin
		readf, lun, temp, pelog
;
;NOTE: I used 3 different procedures called from this procedure.  Each of the called
;procedures listed below then called functions that returned the alpha values
;(cross sections). Conversion to kappa per gram was done before the values were
;passed to this routine.  You should substitute your procedures(s) here.  Be sure
;to pass the variables back and forth correctly!
;
		abcoh, temp, pelog, kappa_anuh			; neutral H
		acohneg, temp, pelog, kappa_anuhng		; neg H ion
		scatter, temp, pelog, kappa_sigtot		; all scattering
;
		printf, 5, temp, pelog, alog10(kappa_anuh), alog10(kappa_anuhng), alog10(kappa_sigtot), $
		  format=fmt
	endwhile
	free_lun, lun
	free_lun, 5
	end