spawn, 'ls gof/OC*/FH0e_???????-???????', acsfiles
restore, 'qqmaster.sav'

tt1 = ttmaster
qq1 = qinertial
t0 = tt1(0)
for i = 0, n_elements(acsfiles)-1 do begin
    acsfile = acsfiles(i)
    print, acsfile, '...', format='(A,A,$)'
    newfile = 'new'+acsfile
    spawn, 'cp '+acsfile+' '+newfile
    errmsg = ''
    fxbopen, unit, newfile, 1, hh, errmsg=errmsg, access='RW'
    if errmsg NE '' then goto, NEXT_FILE
    t0 = fxpar(hh, 'TSTART')
    t1 = fxpar(hh, 'TSTOP')
    tmid = 0.5*(t0+t1)

    if tmid LT 210870664L OR tmid GT 211413960L then begin
        fxbfinish, unit
        spawn, 'rm '+newfile
        print, '(unaffected) ', format='(A,$)'
        goto, NEXT_FILE
    endif

    fxbreadm, unit, ['time', 'acsestq'+['1','2','3','4']], $
      ftt, q1, q2, q3, q4, errmsg=errmsg
    if errmsg NE '' then begin
        fxbfinish, unit
        goto, NEXT_FILE
    endif
        
    fqq = transpose([[q1],[q2],[q3],[q4]])

    nq = n_elements(ftt)
    qint = qterp(tt1-t0, qq1, ftt-t0, /slerp)
    fqq1 = qtmult(qint, fqq)
    fnrm = sqrt(total(fqq1^2,1))
    fqq1 = float(fqq1/rebin(reform([fnrm],1,nq),4,nq))

    fxbwritm, unit, 'acsestq'+['1','2','3','4'], $
      fqq1(0,*), fqq1(1,*), fqq1(2,*), fqq1(3,*)
    fxbfinish, unit
    print, 'H ', format='(A,$)'
    spawn, 'fmodhead '+newfile+' attkeys.dat'
    print, strtrim(nq,2), format='(A," ",$)'

    NEXT_FILE:
    if errmsg NE '' then print, 'FAILED' $
    else print, 'done'
endfor

end