;;; main routine checkreko.  Check the FITS files from the automatic
  ;;; reconstructuion.  Write a list with problematic cases, move
  ;;; good ones to a folder 'OK' and bad ones to 'Bad'

filemask = 'Quiet*'
k = replicate(1b, 5, 5)
spawn, ' [ -d OK ] || mkdir OK ; [ -d Bad ] || mkdir Bad'

;window, 1, xsi=324, ysi=256, tit='Array gt 0'
;window, 2, xsi=324, ysi=256, tit='Image'
;window, 3, xsi=324, ysi=256, tit='Difference'
;window, 4, xsi=1280, ysi=1024, tit='Blink'

f = findfile(filemask, count=n)
close, 1
openw, 1, 'Logfile', /append

FOR i=0, n-1 DO BEGIN
    rdfits, p, f(i)
    m = rescale(dilate(p(*, *, 0) LE 0, k), scale=0.25)
    wset, 1 & erase & tvscl, m
    wset, 2 & erase & tvscl, rescale(p(*, *, 0), sc=.25)
    wset, 3 & erase & tvscl, rescale(p(*, *, 0)-p(*, *, 2), sc=.25)
    wset, 4 & erase
    im = 0
    REPEAT BEGIN
        tvscl, p(*, *, im)
        yorn = byte(get_kbrd(1))
        IF yorn EQ 32 THEN im = (im EQ 2 ? 0 : 2)
    ENDREP UNTIL yorn NE 32
    
    IF yorn NE 10b THEN target=' Bad/' $
    ELSE target = ' OK/'
    printf, 1, f(i)+' -->' +target
    print, f(i)+' -->' +target
    spawn, 'mv '+f(i)+target
ENDFOR

close, 1

END