PRO dustem_old2new_sed_format

dustem_init

dir=!dustem_wrap_soft_dir+'/Data/SEDs/'
file=dir+'Gal_composite_spectrum.xcat'

str='cp '+file+' '+file+'.old'
spawn,str

old_sed=read_xcat(file,/silent,comments=old_comments)
Nsed=n_elements(old_sed)

new_sed=dustem_initialize_sed(Nsed,comments=new_comments,help=help)
comments=['Written by dustem_old2new_sed_format.pro on '+systime(0),old_comments,new_comments]


new_sed.instru=old_sed.instru
new_sed.filter=old_sed.filter
new_sed.wave=old_sed.wave
new_sed.StokesI=old_sed.spec
new_sed.sigmaII=(old_sed.error)^2

write_xcat,new_sed,file,comments=comments

END