pro aocs, files, output, end_time
    start_file = 1
    last_file = n_elements(files) - 1
    maxr = -99999L
    minr = 99999L
    for i = 0, last_file do begin

        p1 = reform(eis_read_fits_data(files[i], 'SAT_ROT', 1, header=p1_header))
        p2 = reform(eis_read_fits_data(files[i], 'UFSAXANG', 1))
        p3 = reform(eis_read_fits_data(files[i], 'UFSAYANG', 1))
        p4 = reform(eis_read_fits_data(files[i], 'UFSBXANG', 1))
        p5 = reform(eis_read_fits_data(files[i], 'UFSBYANG', 1))
        p6 = reform(eis_read_fits_data(files[i], 'ERRANGZ', 1))

        if i eq 0 then begin 
            base_time = fxpar(p1_header, 'UT_START',count=count)
            if count eq 0 then $
              base_time = fxpar(p1_header, 'DATE_OBS')
            p1_array = [p1]
            p2_array = [p2]
            p3_array = [p3]
            p4_array = [p4]
            p5_array = [p5]
            p6_array = [p6]
        endif else begin
            p1_array = [p1_array, p1]
            p2_array = [p2_array, p2]
            p3_array = [p3_array, p3]
            p4_array = [p4_array, p4]
            p5_array = [p5_array, p5]
            p6_array = [p6_array, p6]
        endelse
        if i eq last_file then begin
;;;            end_time = fxpar(p1_header, 'UT_END')
;;;            end_time = '01-Apr-06 23:57:03'
;;;            end_time = '24-Sep-06 00:00:00'
            print, 'LAST FILE'
        endif
        mx = max([p1, p2, p3, p4, p5, p6])
        mn = min([p1, p2, p3, p4, p5, p6])
        if mx gt maxr then maxr = mx
        if mn lt minr then minr = mn
    endfor
    x = lindgen(n_elements(p1_array))
    print, 'BASE_TIME = ', base_time
    print, 'END_TIME = ',  end_time

    cs=1.6

;    window, /free, xsize=1024
;    !p.multi=0
;    utplot, x, htr_a_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4,title='Htr A', ytitl='Temp. (Deg C)',charsize=cs
;    outplot,x, htr_b_array,color='0000FF'x
;    outplot,x, htr_c_array,color='00FFFF'x

    set_plot,'ps'
    ; Could use encapsulated here
;;;    device,xoff=2.5,yoff=6,xs=16,ys=18,file='/tmp/eis_icu_temps.ps'
    device,xoff=2.5,yoff=6,xs=16,ys=18,file=output+'/aocs1.ps'
    !p.multi=[0,1,3]
    utplot, x, p1_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='Sat Rot', ytitl='',charsize=cs
    utplot, x, p2_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='UFS A X Angle', ytitl='',charsize=cs
    utplot, x, p3_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='UFS A Y Angle', ytitl='',charsize=cs
    utplot, x, p4_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='UFS B X Angle', ytitl='',charsize=cs
    utplot, x, p5_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='UFS B Y Angle', ytitl='',charsize=cs
    utplot, x, p6_array, base_time, timerange=[base_time,end_time],yrange=[minr-1,maxr+1],xstyle=2,pstyle=4, $
      title='Error Angle Z', ytitl='',charsize=cs
    device,/close
    set_plot,'X'
end

pro plot_aocs, dir, output, end_time

    ;dir = '/tmp/'
    files = findfile(dir + 'eis_aocs1_*fits', count=count)
    if count eq 0 then begin
        print, 'No files'
        return
    endif
    
;    print, count, ' files'
;    print, 'Files: ', files
;    print, 'n ', n_elements(files)
 
    aocs, files, output, end_time
    
    !p.multi=0
end