dustem_read_mix.pro 701 Bytes
FUNCTION dustem_read_mix,file,silent=silent
  
  st=''
  Nmax=10000L
  sts=strarr(Nmax)
  openr,unit,file,/get_lun
  
;==read comments
  str=''
  first_char='#'
  comments=['']
  Ncomments=0L
  WHILE first_char EQ '#' DO BEGIN
     readf,unit,str
     first_char=strmid(str,0,1)
     comments=[comments,str]
     IF first_char EQ '#' THEN Ncomments=Ncomments+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,comments:comments}
  
  close,unit
  free_lun,unit
  return,full_st
  
END