OV-GSO-DC / dustem-wrapper_idl

Commit 9cb769264545bcfaeb9eb74f401b5d80b2898c4e

Authored by Jean-Philippe Bernard
2 parents d4dc9d03 506ec75f
Exists in master

?? 

Merge branch 'master' of https://gitlab.irap.omp.eu/OV-GSO-DC/dustem-wrapper_idl
Showing 3 changed files with 74 additions and 73 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_plugin_continuum.pro
  Diff comments   View file @9cb7692
1 -FUNCTION dustem_plugin_continuum,key=key,val=val,scope=scope,paramtag=paramtag,help=help 1 +FUNCTION dustem_plugin_continuum,key=key,val=val,scope=scope,paramtag=paramtag,paramdefault=paramdefault,help=help
2 2
3 ;+ 3 ;+
4 ; NAME: 4 ; NAME:
5 ; dustem_plugin_continuum 5 ; dustem_plugin_continuum
6 ; PURPOSE: 6 ; PURPOSE:
7 -; Make products for a given flight 7 +; DUSTEM plugin to compute a blackbody continuum
8 ; CATEGORY: 8 ; CATEGORY:
9 ; DUSTEM Wrapper 9 ; DUSTEM Wrapper
10 ; CALLING SEQUENCE: 10 ; CALLING SEQUENCE:
@@ -17,7 +17,9 @@ FUNCTION dustem_plugin_continuum,key=key,val=val,scope=scope,paramtag=paramtag,h @@ -17,7 +17,9 @@ FUNCTION dustem_plugin_continuum,key=key,val=val,scope=scope,paramtag=paramtag,h
17 ; OUTPUTS: 17 ; OUTPUTS:
18 ; cont = continuum spectrum (on dustem wavelengths) 18 ; cont = continuum spectrum (on dustem wavelengths)
19 ; OPTIONAL OUTPUT PARAMETERS: 19 ; OPTIONAL OUTPUT PARAMETERS:
20 -; None 20 +; scope = scope of the plugin
  21 +; paramdefault = default values of parameters
  22 +; paramtag = plugin parameter names as strings
21 ; ACCEPTED KEY-WORDS: 23 ; ACCEPTED KEY-WORDS:
22 ; help = if set, print this help 24 ; help = if set, print this help
23 ; COMMON BLOCKS: 25 ; COMMON BLOCKS:
@@ -48,6 +50,7 @@ temp=1000. @@ -48,6 +50,7 @@ temp=1000.
48 ampl=1.d-2 50 ampl=1.d-2
49 scope='ADD_SED' 51 scope='ADD_SED'
50 paramtag=['T_BB (K)','Amp'] 52 paramtag=['T_BB (K)','Amp']
  53 +paramdefault=[temp,ampl]
51 54
52 IF keyword_set(key) THEN BEGIN 55 IF keyword_set(key) THEN BEGIN
53 a=where(key EQ 1,count1) 56 a=where(key EQ 1,count1)
src/idl/dustem_plugin_freefree.pro
  Diff comments   View file @9cb7692
1 -FUNCTION dustem_plugin_freefree,key=key,val=val,scope=scope,paramtag=paramtag,help=help 1 +FUNCTION dustem_plugin_freefree,key=key,val=val,scope=scope,paramtag=paramtag,paramdefault=paramdefault,help=help
2 2
3 ;+ 3 ;+
4 ; NAME: 4 ; NAME:
@@ -18,6 +18,8 @@ FUNCTION dustem_plugin_freefree,key=key,val=val,scope=scope,paramtag=paramtag,he @@ -18,6 +18,8 @@ FUNCTION dustem_plugin_freefree,key=key,val=val,scope=scope,paramtag=paramtag,he
18 ; freefree = free-free spectrum (on dustem wavelengths) 18 ; freefree = free-free spectrum (on dustem wavelengths)
19 ; OPTIONAL OUTPUT PARAMETERS: 19 ; OPTIONAL OUTPUT PARAMETERS:
20 ; scope = scope of the plugin 20 ; scope = scope of the plugin
  21 +; paramdefault = default values of parameters
  22 +; paramtag = plugin parameter names as strings
21 ; ACCEPTED KEY-WORDS: 23 ; ACCEPTED KEY-WORDS:
22 ; help = if set, print this help 24 ; help = if set, print this help
23 ; COMMON BLOCKS: 25 ; COMMON BLOCKS:
@@ -46,7 +48,9 @@ Tgas=10000. ;default gas temperature @@ -46,7 +48,9 @@ Tgas=10000. ;default gas temperature
46 Amplitude=1. ;Amplitude 48 Amplitude=1. ;Amplitude
47 smallp=0.0 ;default polarization fraction 49 smallp=0.0 ;default polarization fraction
48 psi=0. ;default polarization angle 50 psi=0. ;default polarization angle
  51 +scope='ADD_SED'
49 paramtag=['T_gas (K)','Amp','p','Psi (deg)'] 52 paramtag=['T_gas (K)','Amp','p','Psi (deg)']
  53 +paramdefault=[Tgas,Amplitude,smallp,psi]
50 IF keyword_set(key) THEN BEGIN 54 IF keyword_set(key) THEN BEGIN
51 a=where(key EQ 1,count1) 55 a=where(key EQ 1,count1)
52 b=where(key EQ 2,count2) 56 b=where(key EQ 2,count2)
@@ -112,8 +116,6 @@ CASE use_method OF @@ -112,8 +116,6 @@ CASE use_method OF
112 END 116 END
113 ENDCASE 117 ENDCASE
114 118
115 -scope='ADD_SED'  
116 -  
117 ;polarization this actually need to 119 ;polarization this actually need to
118 polar_ippsi2iqu,output[*,0],Q,U,replicate(smallp,Nwavs),replicate(psi,Nwavs) 120 polar_ippsi2iqu,output[*,0],Q,U,replicate(smallp,Nwavs),replicate(psi,Nwavs)
119 121
src/idl/dustemcgwin_dataset.pro
  Diff comments   View file @9cb7692
@@ -2124,40 +2124,42 @@ if keyword_set(dataset) then begin @@ -2124,40 +2124,42 @@ if keyword_set(dataset) then begin
2124 endif 2124 endif
2125 endif else begin 2125 endif else begin
2126 ;stop 2126 ;stop
2127 - if ct_hdnpts then begin  
2128 - if ct_spec_hdn_pstv ne 0 then begin  
2129 -  
2130 - xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_spec_hdn_pstv)  
2131 - yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_spec_hdn_pstv)  
2132 - rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_spec_hdn_pstv)/2.  
2133 -  
2134 - dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr  
2135 - dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,noerase=1,color='Black',/overplot,rms=rms,yr=yr,xtickformat='(A1)',xr=xr,xs=1,/xlog,/ylog;,psym=8,syms=0.8  
2136 -  
2137 - endif  
2138 -  
2139 - if ct_filt_hdn_pstv ne 0 then begin  
2140 -  
2141 - xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_filt_hdn_pstv)  
2142 - yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_filt_hdn_pstv)  
2143 - rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_filt_hdn_pstv)/2.  
2144 -  
2145 - dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr  
2146 - dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,color='Black',/overplot,rms=rms,xr=xr,yr=yr,xtickformat='(A1)',/xlog,/ylog;,psym=8,syms=0.8  
2147 -  
2148 - endif  
2149 -  
2150 - if ct_spec_hdn_pstv eq 0 and ct_filt_hdn_pstv eq 0 then begin  
2151 - ; stop  
2152 - if ct_spec ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_spec),-((*!dustem_data.used).values)(idx_spec),ppositions=position,charsize=1.15,/nodata,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/positive_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'  
2153 -  
2154 -  
2155 - if ct_filt ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_filt),-((*!dustem_data.used).values)(idx_filt),/nodata,ppositions=position,charsize=1.15,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/positive_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'  
2156 -  
2157 -  
2158 - endif  
2159 - endif  
2160 - 2127 + dustem_plot_mlog,((*!dustem_data.used).wav)(idx_spec),-((*!dustem_data.used).values)(idx_spec),ppositions=position,charsize=1.15,/nodata,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/positive_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr
  2128 +
  2129 +; if ct_hdnpts then begin
  2130 +; if ct_spec_hdn_pstv ne 0 then begin
  2131 +;
  2132 +; xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_spec_hdn_pstv)
  2133 +; yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_spec_hdn_pstv)
  2134 +; rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_spec_hdn_pstv)/2.
  2135 +;
  2136 +; dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr
  2137 +; dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,noerase=1,color='Black',/overplot,rms=rms,yr=yr,xtickformat='(A1)',xr=xr,xs=1,/xlog,/ylog;,psym=8,syms=0.8
  2138 +;
  2139 +; endif
  2140 +;
  2141 +; if ct_filt_hdn_pstv ne 0 then begin
  2142 +;
  2143 +; xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_filt_hdn_pstv)
  2144 +; yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_filt_hdn_pstv)
  2145 +; rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_filt_hdn_pstv)/2.
  2146 +;
  2147 +; dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr
  2148 +; dustem_plot_mlog,xx,yy,ppositions=position,/positive_only,xs=1,noerase=1,color='Black',/overplot,rms=rms,xr=xr,yr=yr,xtickformat='(A1)',/xlog,/ylog;,psym=8,syms=0.8
  2149 +;
  2150 +; endif
  2151 +;
  2152 +; if ct_spec_hdn_pstv eq 0 and ct_filt_hdn_pstv eq 0 then begin
  2153 +; ; stop
  2154 +; if ct_spec ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_spec),-((*!dustem_data.used).values)(idx_spec),ppositions=position,charsize=1.15,/nodata,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/positive_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'
  2155 +;
  2156 +;
  2157 +; if ct_filt ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_filt),-((*!dustem_data.used).values)(idx_filt),/nodata,ppositions=position,charsize=1.15,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/positive_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'
  2158 +;
  2159 +;
  2160 +; endif
  2161 +; endif
  2162 +;
2161 endelse 2163 endelse
2162 2164
2163 endif 2165 endif
@@ -2218,40 +2220,34 @@ if keyword_set(dataset) then begin @@ -2218,40 +2220,34 @@ if keyword_set(dataset) then begin
2218 endif 2220 endif
2219 2221
2220 endif else begin 2222 endif else begin
  2223 + dustem_plot_mlog,((*!dustem_data.used).wav)(idx_spec),-((*!dustem_data.used).values)(idx_spec),ppositions=position,charsize=1.15,/nodata,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/negative_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr
  2224 +
2221 2225
2222 - if ct_hdnpts then begin  
2223 -  
2224 - if ct_spec_hdn_ngtv ne 0 then begin  
2225 - xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_spec_hdn_ngtv)  
2226 - yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_spec_hdn_ngtv)  
2227 - rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_spec_hdn_ngtv)/2.  
2228 -  
2229 - dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr  
2230 - dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,noerase=1,color='Black',/overplot,rms=rms,yr=yr,xtickformat='(A1)',xr=xr,xs=1,/xlog,/ylog;,psym=8,syms=0.8  
2231 -  
2232 - endif  
2233 -  
2234 - if ct_filt_hdn_ngtv ne 0 then begin  
2235 -  
2236 - xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_filt_hdn_ngtv)  
2237 - yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_filt_hdn_ngtv)  
2238 - rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_filt_hdn_ngtv)/2.  
2239 -  
2240 - dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr  
2241 - dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,color='Black',/overplot,rms=rms,xr=xr,yr=yr,xtickformat='(A1)',/xlog,/ylog;,psym=8,syms=0.8  
2242 -  
2243 - endif  
2244 -  
2245 - if ct_spec_hdn_ngtv eq 0 and ct_filt_hdn_ngtv eq 0 then begin  
2246 -  
2247 - if ct_spec ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_spec),-((*!dustem_data.used).values)(idx_spec),ppositions=position,charsize=1.15,/nodata,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/negative_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'  
2248 -  
2249 -  
2250 - if ct_filt ne 0 then dustem_plot_mlog,((*!dustem_data.used).wav)(idx_filt),-((*!dustem_data.used).values)(idx_filt),/nodata,ppositions=position,charsize=1.15,color='Dodger Blue',psym=16,syms=0.8,xs=1,xr=xr,noerase=1,/negative_only,xtickformat='(A1)',xtit='',/xlog,/ylog,yr=yr;,ytickformat='(A1)'  
2251 -  
2252 -  
2253 - endif  
2254 - endif 2226 +; if ct_hdnpts then begin
  2227 +;
  2228 +; if ct_spec_hdn_ngtv ne 0 then begin
  2229 +; xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_spec_hdn_ngtv)
  2230 +; yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_spec_hdn_ngtv)
  2231 +; rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_spec_hdn_ngtv)/2.
  2232 +;
  2233 +; dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr
  2234 +; dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,noerase=1,color='Black',/overplot,rms=rms,yr=yr,xtickformat='(A1)',xr=xr,xs=1,/xlog,/ylog;,psym=8,syms=0.8
  2235 +;
  2236 +; endif
  2237 +;
  2238 +; if ct_filt_hdn_ngtv ne 0 then begin
  2239 +;
  2240 +; xx = (((*!dustem_show.used).wav)[idx_rmv_sed])(idx_filt_hdn_ngtv)
  2241 +; yy = (((*!dustem_show.used).values)[idx_rmv_sed])(idx_filt_hdn_ngtv)
  2242 +; rms=3.*(((*!dustem_show.used).sigma)[idx_rmv_sed])(idx_filt_hdn_ngtv)/2.
  2243 +;
  2244 +; dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,charsize=1.15,xtickformat='(A1)',xr=xr,color='Black',psym=16,syms=0.8,/xlog,/ylog,yr=yr
  2245 +; dustem_plot_mlog,xx,yy,ppositions=position,/negative_only,xs=1,noerase=1,color='Black',/overplot,rms=rms,xr=xr,yr=yr,xtickformat='(A1)',/xlog,/ylog;,psym=8,syms=0.8
  2246 +;
  2247 +; endif
  2248 +;
  2249 +;
  2250 +; endif
2255 ;stop 2251 ;stop
2256 endelse 2252 endelse
2257 ;stop 2253 ;stop