if (n_elements(infil) eq 0) then begin
    infil = file_list('$DIR_SXT_SFT', 'sft94*')
    rd_roadmap, infil, rmap
    cen = sxt_cen(rmap, /cmd)
    dx = cen(0,*) - 482
    dy = cen(1,*) - 585
    ;
    ss = where( (gt_res(rmap) eq 1) and (gt_dpe(rmap) eq 25) and $
	(abs(dx) lt 2) and (abs(dy) lt 2) and (gt_filtb(rmap) eq 2) )
    ;
    rr = rmap(ss)
    tim2orbit, rr, tim2night=tim2night, /simple
end
;
if (n_elements(bkg3) eq 0) then begin
    ss2 = where((tim2night*60 lt 20) and (tim2night*60 gt 16))	;15 images
    rd_xda, infil, ss(ss2), index, data, rmap
    temp = data(10:500, 200:360, *)
    temp = sxt_decomp(temp)
    bkg3 = total(temp, 3) / n_elements(ss2)
end
if (n_elements(avg1) eq 0) then begin
    n = n_elements(ss)
    avg0 = fltarr(n)	;raw no BKG subtraction
    avg1 = fltarr(n)	;full frame, BKG subtracted
    avg2 = fltarr(n)	;box around leak spot
    avg3 = fltarr(n)	;box around central lattitudes
    ;
    for i=0,n-1 do begin
	print, i+1, ' of ', n, '  ', fmt_tim(rmap(ss(i)))
	;
	rd_xda, infil, ss(i), index, data, rmap
	avg0(i) = total(sxt_decomp(data)) / n_elements(data)
	;
	data2 = dark_sub(index, data)
	avg1(i) = total(data2) / n_elements(data)
	;
	temp = data2(70:270, 30:200)
	avg2(i) = total(temp) / n_elements(temp)
	;
	temp = sxt_decomp( data(10:500, 200:360) )
	avg3(i) = total(temp - bkg3) / n_elements(temp)
    end
end
;
wdef, win, 1000, 700, /already
!p.multi=[0,2,2]
;
x = tim2night*60
xrange0 = [max(x), min(x)]
xrange1 = [20, 0]
xtit='FEM derived TIM2NIGHT (seconds)'
ytit0 = 'Full Image Raw average'
ytit1 = 'Full Image Bkg Subtracted avg'
ytit2 = '200x170 box around leak spot'
ytit3 = '500x160 box around equator (xrays) - Average BKG'
mtit = '1994 SFT Data (DN)'
plot, x, avg0, psym=2, xrange=xrange0, xtit=xtit, tit=ytit0, ytit=mtit
plot, x, avg1, psym=2, xrange=xrange0, xtit=xtit, tit=ytit1, ytit=mtit
plot, x, avg2, psym=2, xrange=xrange0, xtit=xtit, tit=ytit2, ytit=mtit
plot, x, avg3, psym=2, xrange=xrange0, xtit=xtit, tit=ytit3, ytit=mtit

;plot, x, avg0, psym=2, xrange=xrange1, xtit=xtit, ytit=ytit0, tit=mtit
;plot, x, avg,  psym=2, xrange=xrange1, xtit=xtit, ytit=ytit1, tit=mtit
plottime, 0, 0, 'SFT_PLOT  Ver 1.0' & plottime
;
end