Commit 8709270cbcb2fd926e904d2e5f92605099b06843

Authored by Jean-Philippe Bernard
1 parent 9c7cd75f
Exists in master

modified for new form of dustem_data

Showing 1 changed file with 10 additions and 7 deletions   Show diff stats
src/idl/dustem_plot_nuinu_em.pro
... ... @@ -67,10 +67,13 @@ cgoplot,st.sed.wav,st.sed.(i)*fact,color='black',psym=psym
67 67 ;Remarks: not sure that the plot of the sed below actually works. See dustem_sed_plot.pro
68 68  
69 69 defsysv,'!dustem_data',exists=exists
  70 +
  71 +;stop
  72 +
70 73 IF exists THEN BEGIN
71   - IF ptr_valid(!dustem_data.sed) THEN BEGIN
72   - ind_filt=where((*!dustem_data.sed).filt_names NE 'SPECTRUM',count_filt)
73   - ind_spec=where((*!dustem_data.sed).filt_names EQ 'SPECTRUM',count_spec)
  74 + IF ptr_valid((*!dustem_data).sed) THEN BEGIN
  75 + ind_filt=where((*!dustem_data).sed.filt_names NE 'SPECTRUM',count_filt)
  76 + ind_spec=where((*!dustem_data).sed.filt_names EQ 'SPECTRUM',count_spec)
74 77  
75 78 ;=== Plot the data
76 79 ;use_col_data_spec=255
... ... @@ -80,14 +83,14 @@ IF exists THEN BEGIN
80 83 IF count_spec NE 0 THEN BEGIN
81 84 ;plotsym,0,/fill
82 85 ;oplot,((*!dustem_data.sed).wav)(ind_spec),((*!dustem_data.sed).values)(ind_spec)*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_spec))/1.d20,psym=8,_extra=extra,color=use_col_data_spec,syms=0.5
83   - cgoplot,((*!dustem_data.sed).wav)[ind_spec],((*!dustem_data.sed).values)[ind_spec]*(3.e8/1.e-6/(*!dustem_data.sed).wav[ind_spec])/1.d20,psym='Filled Circle',_extra=extra,color=use_col_data_spec
84   - err_bar,((*!dustem_data.sed).wav)(ind_spec),((*!dustem_data.sed).values)(ind_spec)*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_spec))/1.d20,yrms=3.*((*!dustem_data.sed).sigma)(ind_spec)/2.*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_spec))/1.d20,color=use_col_data_spec
  86 + cgoplot,((*!dustem_data).sed.wav)[ind_spec],((*!dustem_data).sed.values)[ind_spec]*(3.e8/1.e-6/(*!dustem_data).sed.wav[ind_spec])/1.d20,psym='Filled Circle',_extra=extra,color=use_col_data_spec
  87 + err_bar,((*!dustem_data).sed.wav)(ind_spec),((*!dustem_data).sed.values)(ind_spec)*(3.e8/1.e-6/(*!dustem_data).sed.wav(ind_spec))/1.d20,yrms=3.*((*!dustem_data).sed.sigma)(ind_spec)/2.*(3.e8/1.e-6/(*!dustem_data).sed.wav(ind_spec))/1.d20,color=use_col_data_spec
85 88 ENDIF
86 89 IF count_filt NE 0 THEN BEGIN
87 90 plotsym,0,/fill
88 91 ;oplot,((*!dustem_data.sed).wav)(ind_filt),((*!dustem_data.sed).values)(ind_filt)*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_filt))/1.d20,psym=8,_extra=extra,color=use_col_data_filt
89   - cgoplot,((*!dustem_data.sed).wav)(ind_filt),((*!dustem_data.sed).values)(ind_filt)*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_filt))/1.d20,psym=8,_extra=extra,color=use_col_data_filt
90   - err_bar,((*!dustem_data.sed).wav)(ind_filt),((*!dustem_data.sed).values)(ind_filt)*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_filt))/1.d20,yrms=3.*((*!dustem_data.sed).sigma)(ind_filt)/2*(3.e8/1.e-6/(*!dustem_data.sed).wav(ind_filt))/1.d20,color=use_col_data_filt
  92 + cgoplot,((*!dustem_data).sed.wav)(ind_filt),((*!dustem_data).sed.values)(ind_filt)*(3.e8/1.e-6/(*!dustem_data).sed.wav(ind_filt))/1.d20,psym=8,_extra=extra,color=use_col_data_filt
  93 + err_bar,((*!dustem_data).sed.wav)(ind_filt),((*!dustem_data).sed.values)(ind_filt)*(3.e8/1.e-6/(*!dustem_data).sed.wav(ind_filt))/1.d20,yrms=3.*((*!dustem_data).sed.sigma)(ind_filt)/2*(3.e8/1.e-6/(*!dustem_data).sed.wav(ind_filt))/1.d20,color=use_col_data_filt
91 94 ENDIF
92 95 ENDIF
93 96 ENDIF
... ...