FUNCTION dustem_read_specem,file,silent=silent ; READ(UNIT = fspecem%unit,FMT = *) COEFFISRF, NTEMP,NSPINC,NIT readcol,file,coeffisrf,ntemp,nspinc,nit,silent=silent st={coeffisrf:0.,ntemp:0,nspinc:0,nit:0} st.coeffisrf=coeffisrf st.ntemp=ntemp st.nspinc=nspinc st.nit=nit return,st END