OV-GSO-DC / dustem-wrapper_idl

Commit 283845dadec7d4b6c754c8cc72ecaabb536c3526

Authored by Ilyes Choubani
1 parent b9bf16a2
Exists in master

Small corrections to test on scopes in write_isrf_lv + starting to update help s… 

…ections for "compute" procs
Showing 8 changed files with 21 additions and 17 deletions   Show diff stats
Inline Side-by-side
src/idl/dustem_compute_ext.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_ext,p_dim,st,$
  1 +FUNCTION dustem_compute_ext,p_dim,$
  2 +st=st,$
2 3 EXT_spec=EXT_spec,$
3 4 out_st=out_st,$
4 5 wave0=wave0,$
5 6 _extra=extra,$
6 7 help=help
7 8  
8   -
9 9 ;+
10 10 ; NAME:
11   -; dustem_compute_sed
  11 +; dustem_compute_ext
12 12 ; PURPOSE:
13   -; Computes an EXT curve from a given Dustem spectrum
  13 +; Computes an extinction curve from a given Dustem spectrum
14 14 ; CATEGORY:
15 15 ; Dustem
16 16 ; CALLING SEQUENCE:
... ...
src/idl/dustem_compute_polext.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_polext ,p_dim,sti,$
  1 +FUNCTION dustem_compute_polext ,p_dim,$
  2 +sti=sti,$
2 3 POLEXT_spec=POLEXT_spec,$
3 4 SPEXT_spec=SPEXT_spec,$
4 5 dustem_fpolext=dustem_fpolext,$
... ...
src/idl/dustem_compute_polsed.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_polsed,p_dim,st,$
  1 +FUNCTION dustem_compute_polsed,p_dim,$
  2 +st=st,$
2 3 P_spec=P_spec,SP_spec=SP_spec,$
3 4 dustem_polfrac=dustem_polfrac,$
4 5 out_st=out_st,$
... ...
src/idl/dustem_compute_sed.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_sed,p_dim,st,$
  1 +FUNCTION dustem_compute_sed,p_dim,$
  2 +st=st,$
2 3 SED_spec=SED_spec,out_st=out_st,$
3 4 wave0=wave0,$
4 5 _extra=extra,$
... ... @@ -50,7 +51,6 @@ ENDIF
50 51 IF not keyword_set(st) THEN BEGIN
51 52  
52 53 dustem_activate_plugins,p_dim/(*(*!dustem_fit).param_init_values),st
53   - ;st=dustem_run(p_dim)
54 54  
55 55 ENDIF
56 56  
... ...
src/idl/dustem_compute_stokes.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_stokes,p_dim,st,$
  1 +FUNCTION dustem_compute_stokes,p_dim,$
  2 +st=st,$
2 3 Q_spec=Q_spec,$
3 4 U_spec=U_spec,$
4 5 PSI_spec=PSI_spec,$
... ...
src/idl/dustem_compute_stokext.pro
  Diff comments   View file @283845d
1   -FUNCTION dustem_compute_stokext,p_dim,st,$
  1 +FUNCTION dustem_compute_stokext,p_dim,$
  2 +st=st,$
2 3 QEXT_spec=QEXT_spec,$
3 4 UEXT_spec=UEXT_spec,$
4 5 PSIEXT_spec=PSIEXT_spec,$
... ...
src/idl/dustem_mpfit_run.pro
  Diff comments   View file @283845d
... ... @@ -99,7 +99,7 @@ FOR i_tag=0,count_data_tag-1 DO BEGIN
99 99 CASE tagnames(ind_data_tag(i_tag)) OF
100 100  
101 101 'SED' : BEGIN
102   - dustem_sed = dustem_compute_sed(p_dim,st,SED_spec=SED_spec)
  102 + dustem_sed = dustem_compute_sed(p_dim,st=st,SED_spec=SED_spec)
103 103 chi2_sed = total(((*(*!dustem_data).sed).values-dustem_sed)^2/(*(*!dustem_data).sed).sigma^2)
104 104 n_sed = n_elements((*(*!dustem_data).sed).values)
105 105 rchi2_sed = chi2_sed / n_sed
... ... @@ -126,7 +126,7 @@ FOR i_tag=0,count_data_tag-1 DO BEGIN
126 126 END
127 127  
128 128 'EXT' : BEGIN
129   - dustem_ext = dustem_compute_ext(p_dim,st,EXT_spec=EXT_spec)
  129 + dustem_ext = dustem_compute_ext(p_dim,st=st,EXT_spec=EXT_spec)
130 130 chi2_ext = total(((*(*!dustem_data).ext).values-dustem_ext)^2/(*(*!dustem_data).ext).sigma^2)
131 131 n_ext = n_elements((*(*!dustem_data).ext).values)
132 132 rchi2_ext = chi2_ext / n_ext
... ... @@ -153,14 +153,14 @@ FOR i_tag=0,count_data_tag-1 DO BEGIN
153 153  
154 154 IF !run_lin then begin
155 155  
156   - dustem_polext = dustem_compute_polext(p_dim,st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext,dustem_ext=dustem_ext)
  156 + dustem_polext = dustem_compute_polext(p_dim,st=st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext,dustem_ext=dustem_ext)
157 157  
158 158 ENDIF
159 159 END
160 160  
161 161 'QSED': BEGIN
162 162  
163   - stokes = dustem_compute_stokes(p_dim,st,Q_spec=Q_spec,U_spec=U_spec,PSI_spec=PSI_spec,dustem_psi_em=dustem_psi_em)
  163 + stokes = dustem_compute_stokes(p_dim,st=st,Q_spec=Q_spec,U_spec=U_spec,PSI_spec=PSI_spec,dustem_psi_em=dustem_psi_em)
164 164 dustem_qsed = stokes[0]
165 165 ;taking care of the color correction for spass ; SHOULD NOT LEAVE THIS ON
166 166 ;dustem_qsed(-1)=dustem_qsed(-2)
... ... @@ -194,7 +194,7 @@ FOR i_tag=0,count_data_tag-1 DO BEGIN
194 194 IF !run_pol and !run_lin THEN BEGIN
195 195 ;the keyword for dustem_sed here is used so that the function uses it to compute the polarization fraction instead of
196 196 ;running dustem_compute_sed again in the procedure. Which it used to do.
197   - dustem_polsed = dustem_compute_polsed(p_dim,st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac,dustem_sed=dustem_sed)
  197 + dustem_polsed = dustem_compute_polsed(p_dim,st=st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac,dustem_sed=dustem_sed)
198 198  
199 199 ENDIF
200 200 END
... ... @@ -202,7 +202,7 @@ FOR i_tag=0,count_data_tag-1 DO BEGIN
202 202  
203 203 'QEXT': BEGIN ;UPDATED
204 204  
205   - Stokext = dustem_compute_stokext(p_dim,st,QEXT_spec=QEXT_spec,UEXT_spec=UEXT_spec,PSIEXT_spec=PSIEXT_spec,dustem_psi_ext=dustem_psi_ext)
  205 + Stokext = dustem_compute_stokext(p_dim,st=st,QEXT_spec=QEXT_spec,UEXT_spec=UEXT_spec,PSIEXT_spec=PSIEXT_spec,dustem_psi_ext=dustem_psi_ext)
206 206 dustem_qext = stokext[0]
207 207 chi2_qext =total(((*(*!dustem_data).qext).values-dustem_qext)^2/(*(*!dustem_data).qext).sigma^2)
208 208 n_qext = n_elements((*(*!dustem_data).qext).values)
... ...
src/idl/dustem_write_isrf_lv.pro
  Diff comments   View file @283845d
... ... @@ -42,7 +42,7 @@ scopes = strarr(n_tags(*!dustem_plugin))
42 42 IF (tag_names(*!dustem_plugin))[0] EQ 'NONE' then goto, the_after
43 43 ;Getting the actual scopes
44 44 FOR i=0L, n_tags(*!dustem_plugin)-1 DO BEGIN
45   - if isa((*(*!dustem_plugin).(i).scope)) THEN BEGIN
  45 + if isa((*!dustem_plugin).(i).scope) THEN BEGIN
46 46 scopes[i] = (*(*!dustem_plugin).(i).scope)
47 47 ENDIF ELSE goto, the_after
48 48 ENDFOR
... ...