Commit 474c4b97746a6ec7e6efc2761d07ef87ee164b39
1 parent
608a8312
Exists in
master
modified to cope for new parameters/keyword convention
Showing
3 changed files
with
12 additions
and
14 deletions
Show diff stats
src/idl/dustem_compute_polsed.pro
| 1 | 1 | FUNCTION dustem_compute_polsed,p_dim,$ |
| 2 | -st=st,$ | |
| 3 | -P_spec=P_spec,$ | |
| 4 | -SP_spec=SP_spec,$ | |
| 5 | -dustem_polfrac=dustem_polfrac,$ | |
| 6 | -out_st=out_st,$ | |
| 7 | -dustem_sed=dustem_sed,$ | |
| 8 | -_extra=extra | |
| 9 | - | |
| 2 | + st=st,$ | |
| 3 | + P_spec=P_spec,$ | |
| 4 | + SP_spec=SP_spec,$ | |
| 5 | + dustem_polfrac=dustem_polfrac,$ | |
| 6 | + dustem_sed=dustem_sed | |
| 10 | 7 | |
| 11 | 8 | ;+ |
| 12 | 9 | ; NAME: |
| ... | ... | @@ -16,7 +13,7 @@ _extra=extra |
| 16 | 13 | ; CATEGORY: |
| 17 | 14 | ; Dustem |
| 18 | 15 | ; CALLING SEQUENCE: |
| 19 | -; sed=dustem_compute_polsed(p_dim[,st=][,P_spec=],[SP_spec=][dustem_polfrac=][out_st=][dustem_sed=][,_extra=][,/help]) | |
| 16 | +; sed=dustem_compute_polsed(p_dim[,st=][,P_spec=],[SP_spec=][dustem_polfrac=][dustem_sed=][,/help]) | |
| 20 | 17 | ; INPUTS: |
| 21 | 18 | ; p_dim = parameter values |
| 22 | 19 | ; OPTIONAL INPUT PARAMETERS: |
| ... | ... | @@ -29,7 +26,8 @@ _extra=extra |
| 29 | 26 | ; SP_spec = Dustem polarized emission fraction |
| 30 | 27 | ; dustem_polfrac = Computed polarization emission fraction for spectrum points in !dustem_data |
| 31 | 28 | ; OPTIONAL OUTPUT PARAMETERS: |
| 32 | -; out_st = Dustem output structure | |
| 29 | +; st = Dustem structure | |
| 30 | +; dustem_sed = Computed emission for spectrum points in !dustem_data | |
| 33 | 31 | ; ACCEPTED KEY-WORDS: |
| 34 | 32 | ; help = If set, print this help |
| 35 | 33 | ; COMMON BLOCKS: |
| ... | ... | @@ -162,7 +160,7 @@ out_st=st |
| 162 | 160 | ;GENERATING THE INTERPOLATES FOR POLFRAC (dustem_polfrac) |
| 163 | 161 | if !run_lin then begin |
| 164 | 162 | |
| 165 | - if not keyword_set(dustem_sed) then dustem_sed = dustem_compute_sed(p_dim,st,SED_spec) | |
| 163 | + if not keyword_set(dustem_sed) then dustem_sed = dustem_compute_sed(p_dim,st=st,SED_spec=SED_spec) | |
| 166 | 164 | |
| 167 | 165 | If n_elements(dustem_sed) ne n_elements(dustem_polsed) then begin |
| 168 | 166 | if n_elements(dustem_sed) gt n_elements(dustem_polsed) then begin | ... | ... |
src/idl/dustem_fit_sed_ext_pol_example.pro
| ... | ... | @@ -418,7 +418,7 @@ ext.EXT_Q = toto[0] |
| 418 | 418 | ext.EXT_U = toto[1] |
| 419 | 419 | |
| 420 | 420 | sed.StokesI = dustem_compute_sed(rv,st=st) |
| 421 | -toto = dustem_compute_stokes(rv,st=st,Qspec=dustem_qsed,Uspec=dustem_used) | |
| 421 | +toto = dustem_compute_stokes(rv,st=st,Q_spec=dustem_qsed,U_spec=dustem_used) | |
| 422 | 422 | sed.StokesQ = toto[0] |
| 423 | 423 | sed.StokesU = toto[1] |
| 424 | 424 | ... | ... |
src/idl/dustem_write_fits_table.pro
| ... | ... | @@ -60,9 +60,9 @@ IF ptr_valid((*!dustem_data).sed) THEN $ |
| 60 | 60 | |
| 61 | 61 | IF !run_pol THEN BEGIN |
| 62 | 62 | ;This is needed to compute Q,U of both SED and spectra |
| 63 | - dustem_predicted_polsed=dustem_compute_stokes(*(*!dustem_fit).current_param_values,bidon, $ | |
| 63 | + dustem_predicted_polsed=dustem_compute_stokes(*(*!dustem_fit).current_param_values, $ | |
| 64 | 64 | Q_spec=dustem_predicted_Q_spec,U_spec=dustem_predicted_U_spec,PSI_spec=dustem_predicted_PSI_spec, $ |
| 65 | - dustem_psi_em,out_st=out_st) | |
| 65 | + dustem_psi_em=dustem_psi_em) | |
| 66 | 66 | dustem_predicted_Qsed=dustem_predicted_polsed[0] |
| 67 | 67 | dustem_predicted_Used=dustem_predicted_polsed[1] |
| 68 | 68 | ;Same in extinction, but does not work in polarization extinction was not fitted. | ... | ... |