diff --git a/src/idl/dustem_run_readme.pro b/src/idl/dustem_run_readme.pro index 5ff6beb..1ed075a 100644 --- a/src/idl/dustem_run_readme.pro +++ b/src/idl/dustem_run_readme.pro @@ -113,7 +113,8 @@ st=dustem_run() ;stt=dustem_set_data(sed=sed) ;== The following plots the resulting emission SED -window,0 +win=0L +window,win & win=win+2 xtit=textoidl('\lambda (\mum)') ;=== same plot as in Compiegne et al. 2010 ytit=textoidl('\nuI_\nu^{em} (W/m^2/sr for N_H=1.e20 H/cm^2)') @@ -121,24 +122,38 @@ tit='Dustem Emitted Intensity '+use_model yr=[1e-11,6.e-7] xr=[1,5e3] dustem_plot_nuinu_em,st,yr=yr,/ysty,xr=xr,/xsty,/xlog,/ylog,title=tit,xtit=xtit,ytit=ytit -IF keyword_set(png) THEN BEGIN - file_png=dir_png+'Last_dustem_run_readme_nuinuem.png' - write_png,file_png,tvrd(/true) - message,'Wrote '+file_png,/info -ENDIF ;== The following plots the resulting extinction -window,1 +window,win & win=win+1 yr=[0,2.5] ;range of 1/lambda xr=[1,10] ;range of sigma xtit=textoidl('1/\lambda (\mum^{-1})') ytit=textoidl('\sigma_{ext} (1e-21 cm^2/H)') tit='Dustem extinction '+use_model dustem_plot_extinction,st,st_model,xr=xr,yr=yr,/xsty,/ysty,xtit=xtit,ytit=ytit,title=tit + IF keyword_set(png) THEN BEGIN - file_png=dir_png+'Last_dustem_run_readme_extinction.png' + win=10L + window,win ;& win=win+2 + file_png=dir_png+'Last_dustem_run_readme_nuinuem.png' + xtit=textoidl('\lambda (\mum)') + ytit=textoidl('\nuI_\nu^{em} (W/m^2/sr for N_H=1.e20 H/cm^2)') + tit='Dustem Emitted Intensity '+use_model + yr=[1e-11,6.e-7] + xr=[1,5e3] + dustem_plot_nuinu_em,st,yr=yr,/ysty,xr=xr,/xsty,/xlog,/ylog,title=tit,xtit=xtit,ytit=ytit write_png,file_png,tvrd(/true) message,'Wrote '+file_png,/info + ;===== + file_png=dir_png+'Last_dustem_run_readme_extinction.png' + yr=[0,2.5] ;range of 1/lambda + xr=[1,10] ;range of sigma + xtit=textoidl('1/\lambda (\mum^{-1})') + ytit=textoidl('\sigma_{ext} (1e-21 cm^2/H)') + tit='Dustem extinction '+use_model + dustem_plot_extinction,st,st_model,xr=xr,yr=yr,/xsty,/ysty,xtit=xtit,ytit=ytit,title=tit + write_png,file_png ,tvrd(/true),/verbose + message,'Wrote '+file_png,/info ENDIF the_end: -- libgit2 0.21.2