pro read_ice_index,file,wave,index
   
;-----------------------------------------------------------------
; Routine reads an ascii file of refractive index data 
;
; Input:
;    file......refractive index data file name
;
; Output:
;    wave.....wavelength (microns)
;    index....complex refractive index
;    
; Source: Mark Hervig, GATS
;-----------------------------------------------------------------

   a     = 0. 
   b     = 0. 
   c     = 0.
   headr = ' '
   nwave = 0L
   
;- read the refractive index file   
   
   openr,lu,file,/get_lun

   for i = 1,7 do readf,lu,headr 

   while (not eof(lu)) do begin
   
     readf,lu,a,b,c
     
     if nwave eq 0 then begin 
       wave = a
       dxr  = b
       dxi  = c
     endif
     if nwave gt 0 then begin 
       wave = [wave,a]
       dxr  = [dxr,b]
       dxi  = [dxi,c]
     endif
     
     nwave = nwave + 1L
   endwhile

   close,lu
   free_lun,lu
   
   index = complex(dxr,dxi)
   
;- done

   return
   end