FUNCTION dustem_read_ionfrac,file,silent=silent readcol,file,frac_cat, frac_neu, frac_ani,silent=silent nlines=n_elements(frac_cat) one_st={frac_cat:0., frac_neu:0., frac_ani:0.} st=replicate(one_st,Nlines) st.frac_cat = frac_cat st.frac_neu = frac_neu st.frac_ani = frac_ani return,st END