close, /all

openr, 1, 'jultimes'
jultimes = strarr(307)
readf, 1, jultimes
close, 1

bb = [33, 40, 81, 99, 102]

FOR ib=0, n_elements(bb)-1 DO BEGIN
;    IF ib EQ 0 THEN GOTO, Weiter
    burst = bb(ib)
     ;;; convert input data to integer
    infile = 'RAW/file01.'+nnumber(burst, 3)
    openr, 1, infile
    bild = assoc(1, bytarr(768, 572))
    openw, 2, 'raw/rawdata1'
    FOR i=0, 99 DO writeu, 2, fix(bild(i))
    close, 2
    close, 1
    spawn, 'rm '+infile
Weiter:    
      ;;; Preprocessing
    
    pre_speckle, 'allvars1.dat' 
    
    sp_mask, /all, /no
    
    caldat, jultimes(burst-1), month, day, year, hrs, mins, secs
    tim = nnumber(hrs, 2)+nnumber(mins, 2)+nnumber(fix(secs+.5), 2)
    
    name = 'AR8739_23oct99_'+tim+'.fits'
    
      ;;; put together the image
      ;;; with bad seeing, this fails sometimes.  Catch that
      ;;; In any case, save the reconstructions, maybe connect does
      ;;; not work correctly.
    
    cmd = 'tar cfI '+tim+'.tbz results1/schrittweite results1/teilbild*'+$
      ' results1/langzeitbild results1/alpha results1/referenzbild '+$
      ' support1/rawsum support1/flaeche* support1/ruhdimensionen '+$
      ' support1/basisdimensionen'
    spawn, cmd
    
    catch, Error_Status
    IF Error_Status NE 0 THEN BEGIN
        print, 'An error occured during connect!'
    ENDIF ELSE BEGIN
        p = connect(/no_ini, /int, /safe)
        IF min(p) EQ 0 THEN BEGIN
            print, 'Min data = 0!'
        ENDIF ELSE $
          save_result, p, 'FITS/'+name, /no
    ENDELSE
    
    catch, /cancel
    
      ;;; clean up disk
    
    spawn, 'rm -rf results1 support1 raw/rawdata1'
    
ENDFOR

END