PRO count_1min5minmag, numb_1min, numb_5min, numb_rollover, dt=dt
;
;      program for choosing good data
;
openr,5,'tmp'
srt_num=0 & readf,5,srt_num
close,5
num=srt_num
mdi_id= strcompress(string(num),/remove_all)

mdiname = strarr(60)
numb_1min=0.0 & numb_5min=0.0 & numb_rollover=0.0
; define good data
maskstring = '0x402c05f2' & qualmask = 0L
reads, maskstring, qualmask, format='(%"%x")'
; definition end

  mdiquality, mdi_qual, mdi_qualvld
  if (mdi_qual(0) eq 'NoFile') then goto, input_dat
  data_dir = '/mag/fd_M_96m_01d.00' + mdi_id + '/fd_M_96m_01d.'
  print, data_dir
  spawn,'\rm dlist_mdi'
  spawn,'ls ' + data_dir + '*.fits' +  ' > dlist_mdi'
  close,5
  mdifil=0
  log_name='dlist_mdi'
  openr,5,log_name
  while not eof(5) do begin
  data_file=''
  readf,5,data_file
  mdiname(mdifil)=data_file
  mdifil = mdifil + 1
  mdifil_end = mdifil - 1
  endwhile
  close,5
  if ( mdifil eq 0 ) then goto, input_dat

name_mdi = STRING('/mag/fd_M_96m_01d.00'+mdi_id+'/fd_M_96m_01d.'+mdi_id+'.00*.fits')
      dat=readfits(name_mdi, hd, /silent)
      pos=strpos(name_mdi, '.fits')
      file_id=fix(strmid(name_mdi, pos-3, 3))
      qualitystr=mdi_qual(file_id)
      qualvldstr=mdi_qualvld(file_id)
      quality=0L & qualvld=0L
      reads, qualitystr, quality, format='(%"%x")'
      reads, qualvldstr, qualvld, format='(%"%x")'
      ; evaluate the quality of the data

input_dat:
spawn,'\rm dlist_mdi'
END