; jtm_rd_ccd.pro - read and process some RAL data for CMB
;   (6-Jul-04) JTM

;    infile = "md_data.028" ; Ne source, translator 110, short wave
    infile = "/Users/oivindw/solarb/eis/data/science/ccsds/md_data.029" ; Ne source, translator 110, short wave
    dark = "/Users/oivindw/solarb/eis/data/science/ccsds/md_data.027"

; Create the data object
    data = obj_new('eis_data', infile, datasource='ccsds', hdr=hdr)
    data_u = obj_new('eis_data', infile, datasource='ccsds', hdr=hdr)

; Correct with dark
    eis_cal_dc, dark, data

; New data for CMB
    ccd_size = data->getccd_sz()
    ccd = fltarr(ccd_size[0], ccd_size[1])
    nwin = data->getnwin()
    for i = 0, nwin-1 do begin
        data->getwin, i, wd, pos
        ccd[pos[0]:pos[0]+pos[1]-1, pos[2]:pos[2]+pos[3]-1] = reform(wd[*,*,0])
    endfor

    ccd_size_u = data_u->getccd_sz()
    ccd_u = fltarr(ccd_size_u[0], ccd_size_u[1])
    nwin_u = data_u->getnwin()
    for i = 0, nwin_u-1 do begin
        data_u->getwin, i, wd_u, pos_u
        ccd_u[pos_u[0]:pos_u[0]+pos_u[1]-1, pos_u[2]:pos_u[2]+pos_u[3]-1] = reform(wd_u[*,*,0])
    endfor

; Put together data for CMB
    d028 = total(ccd[*, 700:749], 2)/50.0
    d028_u = total(ccd_u[*, 700:749], 2)/50.0

; Save as text file for CMB
    openw, unit, 'd029-027.dat', /get_lun
    for i = 0, n_elements(d028)-1 do begin
        printf, unit, i, d028[i], d028_u[i]
    endfor
    free_lun, unit

; Do a plot
;    set_plot, 'ps'
;    ps_long
;    device, file='image029.ps'
    !p.multi= [0, 1, 2]
    plot, d028, xrange=[0, 2147], /xstyle, psym=10, $
        xtitle='CCD Pixel', ytitle='DN', $
        title='Fsrun 154, image 029[700:749,*], dark subtracted'
    plot, d028, xrange=[2148, 4295], /xstyle, psym=10, $
        xtitle='CCD Pixel', ytitle='DN'
;    device, /close
;    set_plot, 'win'

end