;;;;$ls -1 *.air > filelist ;;;doesn't seem to work in .r mode.

phyle=''
phyleout=''
openr,1,'filelist'
while not eof(1) do begin
  readf,1,phyle
  print,phyle
  openr,2,phyle
  phyleout= strmid(phyle,0,9)+'vac'
  openw,3,phyleout

  while not eof(2) do begin
    readf,2,airlam,w

    lam= airlam*10.  ; convert to Angstroms

    sigma2 = 1.d8/lam/lam
    fact = 1.D0 + 6.4328D-5 + 2.94981D-2/(146.D0 - sigma2) + $
                            2.5540D-4/( 41.D0 - sigma2)

    vaclam=fact*airlam
    printf,3,vaclam, w, format='(2f10.6)'
  endwhile
  close,2
  close,3

endwhile
close,1
end