pro visgencoll
@c_system.inc
@c_header.inc

tab=string(9B)

exptim=fltarr(500)
wvl=fltarr(500)
shutter=intarr(500)
sisig=fltarr(500)
vis=fltarr(500)
spsi=fltarr(500)
spflsi=fltarr(500)
spingas=fltarr(500)
meastype=strarr(500)

spawn,'lsf',files,/noshell ; get the file names
s=size(files)
fcount=s(1)-1

for i=0,fcount do begin
 d_read,files(i),h,p
 exptim(i)=d_value(h,h_exptime)
 x=d_value(h,h_meas_type)
 if x eq 6 then meastype(i)='DLV'
 if x eq 7 then meastype(i)='ULV'
 if x eq 8 then meastype(i)='DLIS'
 if x eq 9 then meastype(i)='ULIS'
 if x eq 12 then meastype(i)='SLI strip'
 if x eq 14 then meastype(i)='SA'
 if x eq 16 then meastype(i)='DLVS'
 if x eq 17 then meastype(i)='ULVS'
 if x eq 18 then meastype(i)='Dark'
 if x eq 19 then meastype(i)='CCD'
 if x eq 21 then meastype(i)='DLI2'
 if x eq 22 then meastype(i)='SLI'
 if x eq 23 then meastype(i)='DLI1'
 wvl(i)=d_value(h,h_mono_wavl)
 shutter(i)=d_value(h,h_cal_shut)
 sisig(i)=d_value(h,h_mono_si)
 spsi(i)=d_value(h,h_spsi)
 spflsi(i)=d_value(h,h_spflsi)
 spingas(i)=d_value(h,h_spingas)
 sdev=stdev(p(*,0:199),visavg)
 vis(i)=visavg
endfor

cd,'$RED/Linearity', current=olddir

openw,16,'junk.dat'
for i=0,fcount do begin
 printf,16,i,tab,meastype(i),tab,files(i),tab,wvl(i),tab,$
  shutter(i),tab,sisig(i),tab,spsi(i),tab,spflsi(i),tab,$
  spingas(i),tab,exptim(i),tab,vis(i)
 endfor
close,16

cd, olddir

return
end