PRO get_data_from_hk_database
start_time = 8.6408641e+14
stop_time = 9.2197442e+14
  hk_data = get_aim_hk_telemetry( start_time, stop_time, tlmId=22330, /gps )		; PX
  hkdata = (*hk_data.cips$px_ccd_temp)[*]
  px_ccd_temp = hkdata.housekeeping.eu[*]
  px_ccd_time = hkdata.housekeeping.timetag[*]

  hk_data = get_aim_hk_telemetry( start_time, stop_time, tlmId=22176, /gps )		; PY
  hkdata = (*hk_data.cips$py_ccd_temp)[*]
  py_ccd_temp = hkdata.housekeeping.eu[*]
  py_ccd_time = hkdata.housekeeping.timetag[*]

  hk_data = get_aim_hk_telemetry( start_time, stop_time, tlmId=22237, /gps )		; MX
  hkdata = (*hk_data.cips$mx_ccd_temp)[*]
  mx_ccd_temp = hkdata.housekeeping.eu[*]
  mx_ccd_time = hkdata.housekeeping.timetag[*]

  hk_data = get_aim_hk_telemetry( start_time, stop_time, tlmId=22213, /gps )		; MY
  hkdata = (*hk_data.cips$my_ccd_temp)[*]
  my_ccd_temp = hkdata.housekeeping.eu[*]
  my_ccd_time = hkdata.housekeeping.timetag[*]

; Now save the data by hour
  min_usec_time = MIN( px_ccd_time )
  max_usec_time = MAX( px_ccd_time )
  start_date = LONG( usec2yd( min_usec_time ) )
  stop_date = LONG( usec2yd( max_usec_time ) )
  openw, unit, "/home/frodo/btemple/AIM/webstuff/ccd_temps.txt", /get_lun

  current_date = start_date
  WHILE current_date LE stop_date DO BEGIN

  num_temp = N_ELEMENTS( px_ccd_temp )
   mid_time = yd2usec( current_date )
   FOR i = 0UL, 23 DO BEGIN
      px_index = WHERE( px_ccd_time GE (mid_time - 1800000000) AND  px_ccd_time LT (mid_time + 1800000000) )
      py_index = WHERE( py_ccd_time GE (mid_time - 1800000000) AND  py_ccd_time LT (mid_time + 1800000000) )
      mx_index = WHERE( mx_ccd_time GE (mid_time - 1800000000) AND  mx_ccd_time LT (mid_time + 1800000000) )
      my_index = WHERE( my_ccd_time GE (mid_time - 1800000000) AND  my_ccd_time LT (mid_time + 1800000000) )
      IF px_index[0] NE -1 AND py_index[0] NE -1 AND mx_index[0] NE -1 AND my_index[0] NE -1 THEN BEGIN
         date_str = STRING( DOUBLE(current_date) + (DOUBLE(i)/12.0), format='(%"%13.4f")' )
         PRINTF, unit, format='(%"%s, %f, %f, %f, %f")', date_str, MEDIAN(px_ccd_temp[px_index]), MEDIAN(py_ccd_temp[py_index]), MEDIAN(mx_ccd_temp[mx_index]), MEDIAN(my_ccd_temp[my_index])
;         PRINT, format='(%"%s, %f, %f, %f, %f")', date_str, MEDIAN(px_ccd_temp[px_index]), MEDIAN(py_ccd_temp[py_index]), MEDIAN(mx_ccd_temp[mx_index]), MEDIAN(my_ccd_temp[my_index])
      ENDIF
      mid_time += (i * 60 * 60 * 1000000)
   ENDFOR
   current_date += 1
   if current_date EQ 2007366 THEN current_date = 2008001
   if current_date EQ 2008366 THEN current_date = 2009001
   ENDWHILE
   CLOSE, unit
   FREE_LUN, unit
END