;
;
;

dpath='/usr2/mpenn/02sep17/'

dark=fltarr(256,256)
suncent=fltarr(256,256)
for i=0,63 do begin
 fname=dpath+'suncent02.d.i.'+string(i+1,format='(i3.3)')
 temp=swap_endian(rfits(fname))
 dark=dark+float(temp)/64.
 fname=dpath+'suncent02.o.i.'+string(i+1,format='(i3.3)')
 temp=swap_endian(rfits(fname))
 suncent=suncent+float(temp)/64.
endfor

suncent=suncent-dark
suncent=rotate(suncent,1)
flat=csun_irflat(suncent)
flat=rotate(flat,3)

k=0
fname=dpath+'sunlimb01.o.i.'+string(k+1,format='(i3.3)')
limb=swap_endian(rfits(fname))
callimb=fd_correct2(limb,flat,dark)

plot,alog(findgen(256))/alog(10.),alog(callimb(100,*))/alog(10.),/ynoz


;limb=swap_endian(rfits(dpath+'sunlimb01_64coad.o.i.001'))
;callimb=fd_correct2(limb,flat,dark)
;window,0,xs=256,ys=256
;tvhist,callimb





the_end:
end