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