diff --git a/src/idl/dustem_plugin_modify_dust_pol.pro b/src/idl/dustem_plugin_modify_dust_pol.pro index eb1c630..31d6c15 100644 --- a/src/idl/dustem_plugin_modify_dust_pol.pro +++ b/src/idl/dustem_plugin_modify_dust_pol.pro @@ -3,33 +3,47 @@ FUNCTION dustem_plugin_modify_dust_pol, key=key, val=val, scope=scope, paramtag= ;+ ; NAME: ; dustem_plugin_modify_dust_pol +; ; PURPOSE: -; modifies stokes Q,U values for dust polarization in st according to the given keywords and values +; modifies Emission Stokes Q,U values for dust polarization in st according to the given keywords and values +; ; CATEGORY: -; DUSTEM Wrapper +; DustEMWrap, Plugin, Mid-level, Distributed +; ; CALLING SEQUENCE: ; a=dustem_plugin_modify_dust_pol(st,[key=][val=][,scope=][,paramtag=][,paramdefault=][,/help]) +; ; INPUTS: ; st = dustem structure +; ; OPTIONAL INPUT PARAMETERS: ; key = input parameter numbers (first = polarization fraction in %, default=1.%, second=polarization angle, default=0.) ; val = input parameter values +; ; OUTPUTS: ; out = array containing the stokes emission parameters associated to the dust/synchrotron component or a concatenated version for both +; ; OPTIONAL OUTPUT PARAMETERS: ; scope = if set, returns only the scope of the pluggin ; paramtag = if set, returns only the parameter tags +; ; ACCEPTED KEY-WORDS: ; help = if set, print this help +; ; COMMON BLOCKS: ; None +; ; SIDE EFFECTS: ; None +; ; RESTRICTIONS: ; The dustem fortran code must be installed ; The dustem idl wrapper must be installed -; PROCEDURE: -; This is a dustem plugin +; +; MODIFICATION HISTORY: +; Written by IC 2022 +; Evolution details on the DustEMWrap gitlab. +; See http://dustemwrap.irap.omp.eu/ for FAQ and help. ;- IF keyword_set(help) THEN BEGIN @@ -74,9 +88,9 @@ Nwaves=(size(I))[1] ;0/0 division in frac ;since we divide by I we only need to avoid its null indices +frac_model=P*0. indI = where(I ne 0, countI) -if countI ne 0 then frac_model = P*0. & frac_model[indI] = P[indI]/I[indI] ;This is the polarization fraction in the model - +if countI ne 0 then frac_model[indI] = P[indI]/I[indI] ;This is the polarization fraction in the model frac_used=frac_model*smallp_fact diff --git a/src/idl/dustem_plugin_modify_dust_polx.pro b/src/idl/dustem_plugin_modify_dust_polx.pro index 43e4c7f..ea7c14b 100644 --- a/src/idl/dustem_plugin_modify_dust_polx.pro +++ b/src/idl/dustem_plugin_modify_dust_polx.pro @@ -3,35 +3,50 @@ FUNCTION dustem_plugin_modify_dust_polx, key=key, val=val, scope=scope, paramtag ;+ ; NAME: ; dustem_plugin_modify_dust_polx +; ; PURPOSE: -; modifies stokes Q,U values for dust polarization in st according to the given keywords and values +; modifies Extinction Stokes Q,U values for dust polarization in st according to the given keywords and values +; ; CATEGORY: -; DUSTEM Wrapper +; DustEMWrap, Plugin, Mid-level, Distributed +; ; CALLING SEQUENCE: ; a=dustem_plugin_modify_dust_polx(st,[key=][val=][,scope=][,paramtag=][,/help]) +; ; INPUTS: ; st = dustem structure +; ; OPTIONAL INPUT PARAMETERS: ; key = input parameter numbers (first = polarization fraction in %, default=1.%, second=polarization angle, default=0.) ; val = input parameter values +; ; OUTPUTS: ; out = array containing the stokes emission parameters associated to the dust/synchrotron component or a concatenated version for both +; ; OPTIONAL OUTPUT PARAMETERS: ; scope = if set, returns only the scope of the pluggin ; paramtag = if set, returns only the parameter tags +; ; ACCEPTED KEY-WORDS: ; help = if set, print this help +; ; COMMON BLOCKS: ; None +; ; SIDE EFFECTS: ; None +; ; RESTRICTIONS: ; The dustem fortran code must be installed ; The dustem idl wrapper must be installed -; PROCEDURE: -; This is a dustem plugin +; +; MODIFICATION HISTORY: +; Written by IC 2022 +; Evolution details on the DustEMWrap gitlab. +; See http://dustemwrap.irap.omp.eu/ for FAQ and help. ;- + IF keyword_set(help) THEN BEGIN doc_library,'dustem_plugin_modify_dust_polx' out=0. @@ -88,7 +103,8 @@ ENDIF ;since we divide by I we only need to avoid its null indices indx = where(EXT_spec ne 0, countx) -if countx ne 0 then frac_model = POLEXT_spec*0. & frac_model[indx] = POLEXT_spec[indx]/EXT_spec[indx] ;This is the polarization fraction in the model +frac_model = POLEXT_spec*0. +if countx ne 0 then frac_model[indx] = POLEXT_spec[indx]/EXT_spec[indx] ;This is the polarization fraction in the model frac_used=frac_model*smallp_fact diff --git a/src/idl/dustem_plugin_modify_isrf.pro b/src/idl/dustem_plugin_modify_isrf.pro index 16de28f..9458b8a 100644 --- a/src/idl/dustem_plugin_modify_isrf.pro +++ b/src/idl/dustem_plugin_modify_isrf.pro @@ -83,7 +83,7 @@ ENDIF if isa(!dustem_isrf_file) then user_ISRF = dustem_read_isrf(*!dustem_isrf_file) out=amp*User_ISRF.isrf - + the_end: scope='USER_ISRF' -- libgit2 0.21.2