pro ted7b, seq, item, ipos
;
;
infil = map_filnam(seq, 0, 0, item=item)
;
if (n_elements(item) eq 1) then begin
    img = rfits(infil, h=h)
end else begin
    amp = rfits(infil(0), h=h)
    offset = rfits(infil(1))
    img = 2*amp / (offset+amp)
end
;
printf, 1, 'Input File(s): ' + arr2str(infil)
printf, 1, '  '
;
nx = 10
ny = 20
;;lasx = sxpar(h, 'laspos_x')
;;lasy = sxpar(h, 'laspos_y')
if (ipos eq 0) then begin
    lasx = 32
    lasy = 32
end else begin
    lasx = 32 + cos(-2*!pi*(ipos-1)/4.)*20
    lasy = 32 + sin(-2*!pi*(ipos-1)/4.)*20
end
xst = (lasx-nx/2)>0<(n_elements(img(*,0))-nx)
yst = (lasy-ny/2)>0<(n_elements(img(0,*))-ny)
;
ss = indgen(nx)+xst
printf, 1, string(ss, format='(5x, 100i12)')
for i=yst+ny-1,yst,-1 do printf, 1, string(i, img(ss, i), format='(i5, 100f12.5)')
;
printf, 1, string(12b)
end
;------------------------------------------------------------------------------
;
;
  dir = '/md66/log/sci160k/'
  ff = file_list(dir, '*.tfr', /cd, file=file)
  ss = where((file ge '940316_041320_01.tfr') and (file le '940316_043822_01.tfr'))
seq = file(ss)
seq = strmid(seq, 0, strlen(seq(0))-7)
nseq = n_elements(seq)
;
close, 1
openw, 1, 'ted7b.txt'
printf, 1, 'TED7B  Program Run ' + !stime
;
for i=0,nseq-1 do begin
    ted7b, seq(i), '_m1', i
    ted7b, seq(i), '_m1s', i
    ted7b, seq(i), '_a1', i
    ted7b, seq(i), '_o1', i
    ted7b, seq(i), '_p1', i
    ted7b, seq(i), ['_a1', '_o1'], i

    ted7b, seq(i), '_m2', i
    ted7b, seq(i), '_m2s', i
    ted7b, seq(i), '_a2', i
    ted7b, seq(i), '_o2', i
    ted7b, seq(i), '_p2', i
    ted7b, seq(i), ['_a2', '_o2'], i
end

close, 1
end