pro make_ascii, file ;read-in spectrum f = readfits(file,header) crval1 = strtrim(sxpar(header,'CRVAL1'),2) ;readout coordinate at reference pixel cdelt1 = strtrim(sxpar(header,'CDELT1'),2) ;readout coord. incr. per pixel crpix1 = strtrim(sxpar(header,'CRPIX1'),2) ;readout coord. at reference pixel exptime = strtrim(sxpar(header,'EXPTIME'),2) ;readout exposure time l = dblarr(n_elements(f)) l[0] = double(crval1) for j=0L,n_elements(f)-2 do begin l[j+1] = l[j] + double(cdelt1) endfor l = l-((double(crpix1)-1.d0)*double(cdelt1)) ;correct wavelength for reference pixel data = dblarr(2,n_elements(l)) data[0,*] = l data[1,*] = f name = strmid(file,0,strpos(file,'.')) openw, iunit, name+'.ascii', /get_lun printf, iunit, data close, iunit free_lun, iunit end