OV-GSO-DC / dustem-wrapper_idl

Commit b36e1e5e1e83ea654185d7c0f3b602f32a931dde

Authored by Jean-Philippe Bernard
1 parent 2e1f81db
Exists in master

added png plot options

Showing 2 changed files with 20 additions and 4 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_fit_sed_readme.pro
  Diff comments   View file @b36e1e5
@@ -71,7 +71,7 @@ ENDIF ELSE BEGIN @@ -71,7 +71,7 @@ ENDIF ELSE BEGIN
71 ENDELSE 71 ENDELSE
72 72
73 IF keyword_set(png) THEN BEGIN 73 IF keyword_set(png) THEN BEGIN
74 - dir_png='/tmp/Dustem/Figures/' 74 + dir_png=!dustem_data+'/Figures/'
75 force_mkdir,dir_png 75 force_mkdir,dir_png
76 ENDIF 76 ENDIF
77 77
src/idl/dustem_run_readme.pro
  Diff comments   View file @b36e1e5
1 -PRO dustem_run_readme,postcript=postcript,model=model,help=help,plot_it=plot_it 1 +PRO dustem_run_readme,postcript=postcript,model=model,help=help,plot_it=plot_it,png=png
2 2
3 ;+ 3 ;+
4 ; NAME: 4 ; NAME:
@@ -51,6 +51,11 @@ IF keyword_set(help) THEN BEGIN @@ -51,6 +51,11 @@ IF keyword_set(help) THEN BEGIN
51 goto,the_end 51 goto,the_end
52 ENDIF 52 ENDIF
53 53
  54 +IF keyword_set(png) THEN BEGIN
  55 + dir_png=!dustem_data+'/Figures/'
  56 + force_mkdir,dir_png
  57 +ENDIF
  58 +
54 loadct,13 59 loadct,13
55 60
56 IF keyword_set(model) THEN BEGIN 61 IF keyword_set(model) THEN BEGIN
@@ -107,7 +112,7 @@ st=dustem_run() @@ -107,7 +112,7 @@ st=dustem_run()
107 ;stop 112 ;stop
108 ;stt=dustem_set_data(sed=sed) 113 ;stt=dustem_set_data(sed=sed)
109 114
110 -;== The following just plots the resulting emission SED 115 +;== The following plots the resulting emission SED
111 window,0 116 window,0
112 xtit=textoidl('\lambda (\mum)') 117 xtit=textoidl('\lambda (\mum)')
113 ;=== same plot as in Compiegne et al. 2010 118 ;=== same plot as in Compiegne et al. 2010
@@ -116,7 +121,13 @@ tit='Dustem Emitted Intensity '+use_model @@ -116,7 +121,13 @@ tit='Dustem Emitted Intensity '+use_model
116 yr=[1e-11,6.e-7] 121 yr=[1e-11,6.e-7]
117 xr=[1,5e3] 122 xr=[1,5e3]
118 dustem_plot_nuinu_em,st,yr=yr,/ysty,xr=xr,/xsty,/xlog,/ylog,title=tit,xtit=xtit,ytit=ytit 123 dustem_plot_nuinu_em,st,yr=yr,/ysty,xr=xr,/xsty,/xlog,/ylog,title=tit,xtit=xtit,ytit=ytit
119 -;== The following just plots the resulting extinction 124 +IF keyword_set(png) THEN BEGIN
  125 + file_png=dir_png+'Last_dustem_run_readme_nuinuem.png'
  126 + write_png,file_png,tvrd(/true)
  127 + message,'Wrote '+file_png,/info
  128 +ENDIF
  129 +
  130 +;== The following plots the resulting extinction
120 window,1 131 window,1
121 yr=[0,2.5] ;range of 1/lambda 132 yr=[0,2.5] ;range of 1/lambda
122 xr=[1,10] ;range of sigma 133 xr=[1,10] ;range of sigma
@@ -124,6 +135,11 @@ xtit=textoidl('1/\lambda (\mum^{-1})') @@ -124,6 +135,11 @@ xtit=textoidl('1/\lambda (\mum^{-1})')
124 ytit=textoidl('\sigma_{ext} (1e-21 cm^2/H)') 135 ytit=textoidl('\sigma_{ext} (1e-21 cm^2/H)')
125 tit='Dustem extinction '+use_model 136 tit='Dustem extinction '+use_model
126 dustem_plot_extinction,st,st_model,xr=xr,yr=yr,/xsty,/ysty,xtit=xtit,ytit=ytit,title=tit 137 dustem_plot_extinction,st,st_model,xr=xr,yr=yr,/xsty,/ysty,xtit=xtit,ytit=ytit,title=tit
  138 +IF keyword_set(png) THEN BEGIN
  139 + file_png=dir_png+'Last_dustem_run_readme_extinction.png'
  140 + write_png,file_png,tvrd(/true)
  141 + message,'Wrote '+file_png,/info
  142 +ENDIF
127 143
128 the_end: 144 the_end:
129 145