Commit 0c7104bf61fded7c830fdadb2c103edb294b5e9e

Authored by Jean-Philippe Bernard
1 parent 5a7a4415
Exists in master

adapted to new fileds of sed sturcture

src/idl/dustem_initialize_sed.pro
... ... @@ -62,7 +62,8 @@ ENDIF
62 62  
63 63 ;one_st={instru:'',filter:'',wave:0.,spec:0.,error:0.}
64 64 iv=la_undef(4)
65   -one_st={instru:'',filter:'',wave:iv,spec:iv,error:iv,StokesQ:iv,StokesU:iv,largeP:iv,smallp:iv,psi:iv,sigmaII:iv,sigmaQQ:iv,sigmaUU:iv,sigmaIQ:iv,sigmaIU:iv,sigmaQU:iv,sigma_largep:iv,sigma_smallp:iv,sigma_psi:iv}
  65 +one_st={instru:'',filter:'',wave:iv,StokesI:iv,StokesQ:iv,StokesU:iv,largeP:iv,smallp:iv,psi:iv,sigmaII:iv,sigmaQQ:iv,sigmaUU:iv,sigmaIQ:iv,sigmaIU:iv,sigmaQU:iv,sigma_largep:iv,sigma_smallp:iv,sigma_psi:iv}
  66 +;one_st={instru:'',filter:'',wave:iv,spec:iv,error:iv,StokesQ:iv,StokesU:iv,largeP:iv,smallp:iv,psi:iv,sigmaII:iv,sigmaQQ:iv,sigmaUU:iv,sigmaIQ:iv,sigmaIU:iv,sigmaQU:iv,sigma_largep:iv,sigma_smallp:iv,sigma_psi:iv}
66 67 comments=['Dustem-Wrap SED']
67 68 comments=[comments,'instru: Instrument name']
68 69 comments=[comments,'filter: Instrument filter']
... ...
src/idl/dustem_set_data.pro
... ... @@ -76,8 +76,9 @@ IF keyword_set(sed) THEN BEGIN
76 76 wavs(ind)=dustem_filter2wav(sed(ind).filter)
77 77 ENDIF
78 78 ;=== define observations structure
79   - obs_st={instru_names:sed.instru,filt_names:sed.filter,wav:wavs, values:sed.spec,sigma:sed.error}
80   - ;=== fill !dustem_data
  79 + ;obs_st={instru_names:sed.instru,filt_names:sed.filter,wav:wavs, values:sed.spec,sigma:sed.error}
  80 + obs_st={instru_names:sed.instru,filt_names:sed.filter,wav:wavs, values:sed.StokesI,sigma:sqrt(sed.sigmaII)}
  81 + ;=== fill !dustem_data
81 82 !dustem_data.sed = ptr_new(obs_st)
82 83 ;VGb
83 84 ; If f_HI is specified, multiply the data by f_HI
... ... @@ -88,6 +89,7 @@ IF keyword_set(sed) THEN BEGIN
88 89 ;VGe
89 90 ENDIF
90 91  
  92 +;===== not sure if this below is vere happening anymore
91 93 IF keyword_set(polsed) THEN BEGIN
92 94 wavs=polsed.wave
93 95 ;=== Impose central wavelengths for photometric channels
... ... @@ -96,7 +98,8 @@ IF keyword_set(polsed) THEN BEGIN
96 98 wavs(ind)=dustem_filter2wav(polsed(ind).filter)
97 99 ENDIF
98 100 ;=== define observations structure
99   - obs_st={instru_names:polsed.instru,filt_names:polsed.filter,wav:wavs, values:polsed.spec,sigma:polsed.error}
  101 + ;obs_st={instru_names:polsed.instru,filt_names:polsed.filter,wav:wavs, values:polsed.spec,sigma:polsed.error}
  102 + obs_st={instru_names:polsed.instru,filt_names:polsed.filter,wav:wavs, values:polsed.largeP,sigma:sqrt(sed.sigma_largeP)}
100 103 ;=== fill !dustem_data
101 104 !dustem_data.polsed = ptr_new(obs_st)
102 105 ;VGe
... ... @@ -183,6 +186,7 @@ IF keyword_set(uext) THEN BEGIN
183 186 ;VGe
184 187 ENDIF
185 188  
  189 +;===== not sure if this below is vere happening anymore
186 190 IF keyword_set(polfrac) THEN BEGIN
187 191 wavs=polfrac.wave
188 192 ;=== Impose central wavelengths for photometric channels
... ...