Commit ae0273b187c60b28b12a33338d7ebe7491a17920
1 parent
eafc1f8f
Exists in
master
removed redundant file
Showing
1 changed file
with
0 additions
and
98 deletions
Show diff stats
src/idl/dustem_plot_ext_notgood.pro deleted
... | ... | @@ -1,98 +0,0 @@ |
1 | -PRO dustem_plot_ext,st,st_model,_Extra=extra,help=help,UV=UV | |
2 | - | |
3 | -;+ | |
4 | -; NAME: | |
5 | -; dustem_plot_ext | |
6 | -; PURPOSE: | |
7 | -; Plots a Dustem extinction | |
8 | -; CATEGORY: | |
9 | -; Dustem | |
10 | -; CALLING SEQUENCE: | |
11 | -; dustem_plot_ext,st,st_model[,/help][,_extra=] | |
12 | -; INPUTS: | |
13 | -; st = Dustem model output structure | |
14 | -; st_model = Dustem model parameters structure | |
15 | -; OPTIONAL INPUT PARAMETERS: | |
16 | -; _extra = extra parameters for the plot routine | |
17 | -; OUTPUTS: | |
18 | -; None | |
19 | -; OPTIONAL OUTPUT PARAMETERS: | |
20 | -; None | |
21 | -; ACCEPTED KEY-WORDS: | |
22 | -; help = If set, print this help | |
23 | -; COMMON BLOCKS: | |
24 | -; None | |
25 | -; SIDE EFFECTS: | |
26 | -; model is plotted | |
27 | -; RESTRICTIONS: | |
28 | -; The dustem idl wrapper must be installed | |
29 | -; PROCEDURE: | |
30 | -; None | |
31 | -; EXAMPLES | |
32 | -; dustem_init,/wrap | |
33 | -; dir_in=getenv('DUSTEM_SOFT_DIR')+'/src/dustem3.8_web/' | |
34 | -; st_model=dustem_read_all(dir_in) | |
35 | -; dir=getenv('DUSTEM_RES') | |
36 | -; st=dustem_read_all_res(dir_in) | |
37 | -; loadct,13 | |
38 | -; dustem_plot_ext,st,st_model,yr=[1.e-6,1.e0],/ysty,xr=[1.,400],/xsty | |
39 | -; MODIFICATION HISTORY: | |
40 | -; | |
41 | -; Modified by VG with flags /UV,/IR | |
42 | -; | |
43 | -; Written by J.-Ph. Bernard | |
44 | -; see evolution details on the dustem cvs maintained at CESR | |
45 | -; Contact J.-Ph. Bernard (Jean-Philippe.Bernard@cesr.fr) in case of problems. | |
46 | -;- | |
47 | - | |
48 | -IF keyword_set(help) THEN BEGIN | |
49 | - doc_library,'dustem_plot_ext' | |
50 | - goto,the_end | |
51 | -ENDIF | |
52 | - | |
53 | -Ngrains=n_tags(st.sed)-2 | |
54 | -colors=dustem_grains_colors(Ngrains) | |
55 | - | |
56 | -;ytit=textoidl('extinction (cm^2/g)') & fact=1. | |
57 | -ytit=textoidl('\sigma_{ext} (1e-21 cm^2/H)') & mH=1.66e-24 & fact=mH*1.e21 | |
58 | -tit='Dustem extinction' | |
59 | - | |
60 | -;VG | |
61 | -if keyword_set(UV) then begin | |
62 | - xtit=textoidl('wave number (\mum^-1)') | |
63 | - xrange=1./st.ext.wav | |
64 | - plot,xrange,/nodata,_Extra=extra,xtit=xtit,ytit=ytit,tit=tit | |
65 | -endif else begin | |
66 | - xtit=textoidl('\lambda (\mum)') | |
67 | - xrange=st.ext.wav | |
68 | - plot_oo,xrange,/nodata,_Extra=extra,xtit=xtit,ytit=ytit,tit=tit | |
69 | -endelse | |
70 | -;VG | |
71 | - | |
72 | -abs_tot=st.ext.abs_grain(0)*0. | |
73 | -sca_tot=st.ext.sca_grain(0)*0. | |
74 | -FOR i=0L,Ngrains-1 DO BEGIN | |
75 | - oplot,xrange,(st.ext.abs_grain(i)+st.ext.sca_grain(i))*fact,color=colors(i),psym=psym | |
76 | -; oplot,xrange,st.ext.abs_grain(i)*fact,color=colors(i-1),psym=psym | |
77 | -; oplot,xrange,st.ext.sca_grain(i)*fact,color=colors(i-1),psym=psym,linestyle=2 | |
78 | - abs_tot=abs_tot+st.ext.abs_grain(i)*fact | |
79 | - sca_tot=sca_tot+st.ext.sca_grain(i)*fact | |
80 | -ENDFOR | |
81 | -oplot,xrange,(abs_tot+sca_tot),color=255,psym=psym | |
82 | -;oplot,xrange,abs_tot*fact,color=255,psym=psym | |
83 | -;oplot,xrange,sca_tot*fact,color=255,psym=psym,linestyle=2 | |
84 | - | |
85 | -defsysv,'!dustem_data',exists=exists | |
86 | -if exists then begin | |
87 | - if tag_exist(!dustem_data,'ext') then begin | |
88 | - if keyword_set(UV) then begin | |
89 | - oploterror,1./(*!dustem_data.ext).wav,(*!dustem_data.ext).values,(*!dustem_data.ext).wav*0.,(*!dustem_data.ext).sigma,psym=4 | |
90 | - endif else begin | |
91 | - oploterror,(*!dustem_data.ext).wav,(*!dustem_data.ext).values,(*!dustem_data.ext).wav*0.,(*!dustem_data.ext).sigma,psym=4 | |
92 | - endelse | |
93 | - endif | |
94 | -endif | |
95 | - | |
96 | -the_end: | |
97 | - | |
98 | -END |