dustem_read_mix.pro 570 Bytes
FUNCTION dustem_read_mix,file,silent=silent
  
  st=''
  Nmax=10000L
  sts=strarr(Nmax)
  openr,unit,file,/get_lun
  
;==read comments
  str='' & first_char='#'
  WHILE first_char EQ '#' DO BEGIN
     readf,unit,str
     first_char=strmid(str,0,1)
  ENDWHILE
  
;==read FMIX
;last line in previous WHILE loop is first FMIX
  tmp = [double(str)]
  WHILE NOT EOF(unit) DO BEGIN
     READF, unit, newline
     tmp = [tmp, newline]
  ENDWHILE
  fmix = tmp

;==save and return structure
  full_st={file:file,fmix:fmix}
  
  close,unit
  free_lun,unit
  return,full_st
  
END