;== Initialize DustemWrap
dustem_init,mode='COMPIEGNE_ETAL2010'

;== set DustemWrap behavior
!dustem_verbose=1
!dustem_show_plot=1

;== The following structure contains the inputs to the model
st_model=dustem_read_all(!dustem_soft_dir)

;== The following line writes inputs to the Fortran
dustem_write_all,st_model,!dustem_res
;== This runs the Fortran. The ouput structure contains the results
st=dustem_run()

;== The following just plots the resulting emission SED
loadct,13
window,0
dustem_plot_nuinu_em,st,yr=[1e-11,6.e-7],/ysty,xr=[1,5e3],/xsty
;== The following just plots the resulting extinction
window,1
dustem_plot_extinction,st,st_model,yr=[1.e-6,1.e0],/ysty,xr=[1.,400],/xsty