OV-GSO-DC / dustem-wrapper_idl

Commit 54e685fcd8fcfb844fdd9437beb66e9a8d48e82a

Authored by Ilyes Choubani
1 parent 91fe3103
Exists in master

small corrections to dustemwrap_plot

Showing 1 changed file with 12 additions and 5 deletions   Show diff stats
Inline Side-by-side
src/idl/dustemwrap_plot.pro
  Diff comments   View file @54e685f
... ... @@ -101,7 +101,9 @@ IF not keyword_set(st) THEN BEGIN
101 101 if !run_pol && !run_lin then begin
102 102  
103 103 dustem_polsed = dustem_compute_polsed(p_dim,st=st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac)
104   - toto = dustem_compute_stokes(p_dim,st=st,dustem_qsed=dustem_qsed,dustem_used=dustem_used,Q_spec=Q_spec,U_spec=U_spec,PSI_spec=PSI_spec,dustem_psi_em=dustem_psi_em)
  104 + toto = 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)
  105 + dustem_qsed = toto[0]
  106 + dustem_used = toto[1]
105 107 ;stop
106 108  
107 109 endif
... ... @@ -116,8 +118,9 @@ IF not keyword_set(st) THEN BEGIN
116 118  
117 119 ;stop
118 120 dustem_polext = dustem_compute_polext(p_dim,st=st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext)
119   - toto = dustem_compute_stokext(p_dim,st=st,dustem_qext=dustem_qext,dustem_uext=dustem_uext,QEXT_spec=QEXT_spec,UEXT_spec=UEXT_spec,PSIEXT_spec=PSIEXT_spec,dustem_psi_ext=dustem_psi_ext)
120   -
  121 + toto = 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)
  122 + dustem_qext = toto[0]
  123 + dustem_uext = toto[1]
121 124 endif
122 125 endif
123 126  
... ... @@ -132,7 +135,9 @@ ENDIF ELSE BEGIN ;st is provided
132 135 if !run_pol && !run_lin then begin
133 136  
134 137 dustem_polsed = dustem_compute_polsed(p_dim,st=st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac)
135   - toto = dustem_compute_stokes(p_dim,st=st,dustem_qsed=dustem_qsed,dustem_used=dustem_used,Q_spec=Q_spec,U_spec=U_spec,PSI_spec=PSI_spec,dustem_psi_em=dustem_psi_em)
  138 + toto = 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)
  139 + dustem_qsed = toto[0]
  140 + dustem_used = toto[1]
136 141 ;stop
137 142 endif
138 143 endif
... ... @@ -147,7 +152,9 @@ ENDIF ELSE BEGIN ;st is provided
147 152 if !run_pol && !run_lin then begin
148 153  
149 154 dustem_polext = dustem_compute_polext(p_dim,st=st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext)
150   - toto = dustem_compute_stokext(p_dim,st=st,dustem_qext=dustem_qext,dustem_uext=dustem_uext,QEXT_spec=QEXT_spec,UEXT_spec=UEXT_spec,PSIEXT_spec=PSIEXT_spec,dustem_psi_ext=dustem_psi_ext)
  155 + toto = 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)
  156 + dustem_qsed = toto[0]
  157 + dustem_used = toto[1]
151 158  
152 159 endif
153 160 endif
... ...