diff --git a/src/idl/dustemwrap_plot.pro b/src/idl/dustemwrap_plot.pro
index 5b72fb2..51da7ef 100755
--- a/src/idl/dustemwrap_plot.pro
+++ b/src/idl/dustemwrap_plot.pro
@@ -101,7 +101,9 @@ IF not keyword_set(st) THEN BEGIN
       if !run_pol && !run_lin then begin
           
           dustem_polsed = dustem_compute_polsed(p_dim,st=st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac)
-          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) 
+          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) 
+          dustem_qsed = toto[0]
+          dustem_used = toto[1]
           ;stop
           
       endif         
@@ -116,8 +118,9 @@ IF not keyword_set(st) THEN BEGIN
 
           ;stop                              
           dustem_polext = dustem_compute_polext(p_dim,st=st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext)
-          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) 
-          
+          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) 
+          dustem_qext = toto[0]
+          dustem_uext = toto[1]
       endif
   endif
   
@@ -132,7 +135,9 @@ ENDIF ELSE BEGIN ;st is provided
         if !run_pol && !run_lin then begin
             
             dustem_polsed = dustem_compute_polsed(p_dim,st=st,P_spec=P_spec,SP_spec=SP_spec,dustem_polfrac=dustem_polfrac)
-            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) 
+            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) 
+            dustem_qsed = toto[0]
+            dustem_used = toto[1]
             ;stop    
         endif
     endif
@@ -147,7 +152,9 @@ ENDIF ELSE BEGIN ;st is provided
         if !run_pol && !run_lin then begin
             
             dustem_polext = dustem_compute_polext(p_dim,st=st,POLEXT_spec=POLEXT_spec,SPEXT_spec=SPEXT_spec,dustem_fpolext=dustem_fpolext)
-            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) 
+            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) 
+            dustem_qsed = toto[0]
+            dustem_used = toto[1]
             
         endif
     endif
--
libgit2 0.21.2