;+
; NAME: eis_status_sts1_conversions.pro
;
; PURPOSE: To convert EIS status 1 parameters from telemetered values
;          to engineering values
;      
; CATEGORY:
;
; CALLING SEQUENCE: 
;
; INPUTS: 
;
; KEYWORD PARAMETERS: None
;
; OUTPUTS: 
;
; CALLS:
;
; COMMON BLOCKS: None
;
; PROCEDURE:
;
; RESTRICTIONS: Uses SolarSoft
;
; MODIFICATION HISTORY:
;	03/June/2006	mcrw wrote
;
;-------------------------------------------------------------------------------------------------
; files - string array of filenames in time order
; output - string pathname to where to write ps files
;
pro icu_temperatures_check, files, log=logfile

    p1 = reform(eis_get_daily_data(files,'CCDA_TEMP',1,3,maxp2,minp2,stim,etim))
    p2 = reform(eis_get_daily_data(files,'CCDB_TEMP',2,3,maxp3,minp3))

    ccda_min = eis_convert_ccda_temperature([minp2])
    ccda_max = eis_convert_ccda_temperature([maxp2])
    ccdb_min = eis_convert_ccdb_temperature([minp3])
    ccdb_max = eis_convert_ccdb_temperature([maxp3])

    print,format='(A29,A29)', $
      stim,etim

;    print, 'CCD-A temperature : min ' + $
;      strtrim(string(eis_convert_ccda_temperature([minp2])),2) + $
;      ' degC	max ' + $
;      strtrim(string(eis_convert_ccda_temperature([maxp2])),2) + $
;      ' degC'
;    
;    print, 'CCD-B temperature : min ' + $
;      strtrim(string(eis_convert_ccdb_temperature([minp3])),2) + $
;      ' degC	max ' + $
;      strtrim(string(eis_convert_ccdb_temperature([maxp3])),2) + $
;      ' degC'

;    print
    

end

pro eis_ccd_temp_trend, dir, log=log_file

    files = findfile(dir + 'eis_sts1_*fits*', count=count)
    if count eq 0 then begin
        print, 'No files'
        return
    endif

    p1 = reform(eis_get_daily_data(files,'CCDA_TEMP',1,3,maxp2,minp2,stim,etim))
    p2 = reform(eis_get_daily_data(files,'CCDB_TEMP',2,3,maxp3,minp3))

    ccda_min = eis_convert_ccda_temperature([minp2])
    ccda_max = eis_convert_ccda_temperature([maxp2])
    ccdb_min = eis_convert_ccdb_temperature([minp3])
    ccdb_max = eis_convert_ccdb_temperature([maxp3])

    mom = moment(p1)
    ccda_avg = eis_convert_ccda_temperature([mom[0]])
    ccda_std = sqrt(mom[1])
    mom = moment(p2)
    ccdb_avg = eis_convert_ccda_temperature([mom[0]])
    ccdb_std = sqrt(mom[1])

;    print,format='(2(A23," "),6(F8.2), 2(F8.2))', $
;      stim,etim,ccda_min,ccda_max,ccda_avg,ccdb_min,ccdb_max,ccdb_avg, ccda_std, ccdb_std

    print,format='(2(A23," "),6(F8.2))', $
      stim,etim,ccda_min,ccda_max,ccda_avg,ccdb_min,ccdb_max,ccdb_avg

    if n_elements(log_file) ne 0 then begin
        openw, lu, log_file, /append, /get_lun
        printf,lu,format='(2(A23,","),5(F6.2,","),F6.2)', $
          stim,etim,ccda_min,ccda_max,ccda_avg,ccdb_min,ccdb_max,ccdb_avg
        close, lu
        free_lun, lu
    endif
    

end