;Programme utilise pour creer les flat contenus dans flat19980304.save
;la liste 1 a servi a creer flat1, la liste 2 flat2
;150 iterations dans les deux cas.


list = ['efr19980304.133517', 'efr19980304.134020', 'efr19980304.134512', 'efr19980304.135008', 'efr19980304.135502',$
        'efr19980304.135957', 'efr19980304.140439', 'efr19980304.140929', 'efr19980304.141419']
off = [[252.30740, 294.59860], [231.54613, 330.23659], [270.58000, 261.73000], [241.44000, 271.54000], $
       [239.21000, 263.58000], [239.37000, 323.42000], [247.28000, 285.29000], [241.75342, 320.43182], $
       [252.01762, 295.36643]]

;list = ['efr19980304.193551', 'efr19980304.194040', 'efr19980304.194527', 'efr19980304.195022', 'efr19980304.205436']
;off = [[270.64000, 282.00000], [261.98000, 298.80000], $
;       [251.13687, 297.90150], [252.20539, 294.83624], $
;       [251.76989, 295.90868]]


n = n_elements(list)
cube = fltarr(512, 512, n)

for i = 0, n-1 do begin
  img = (readfits(eit_file2path(list(i)), hdr)-852)>0
  sat = where(img ge 32768)
  if sat ne -1 then img(sat) = 0
  img = eit_flat(eit_degrid(img, hdr), hdr) / (eit_fxpar(hdr, 'EXPTIME') + eit_fxpar(hdr, 'SHUTTER CLOSE TIME'))
  cube(*, *, i) = rebin(img, 512, 512)
  
endfor

mk_kuhn_flat, cube, 1, flat, offsets=off, niter=150

end