diff --git a/src/idl/dustem_fit_intensity_example.pro b/src/idl/dustem_fit_intensity_example.pro
index 514f7cc..756d1d5 100644
--- a/src/idl/dustem_fit_intensity_example.pro
+++ b/src/idl/dustem_fit_intensity_example.pro
@@ -387,7 +387,7 @@ IF keyword_set(fits_save_and_restore) THEN BEGIN
   ;==== plot result taken from the saved fits table
   res=*(*!dustem_fit).CURRENT_PARAM_VALUES
   IF !dustem_noobj THEN BEGIN
-     dustemwrap_plot_noobj,res,dustem_spectra_st,xr=xr,/xstyle,yr=yr,/ysty,/ylog,/xlog,title=tit+' (From Saved FITS file)'
+    dustemwrap_plot_noobj,res,dustem_spectra_st,xr=xr,/xstyle,yr=yr,/ysty,/ylog,/xlog,title=tit+' (From Saved FITS file)'
   ENDIF ELSE BEGIN
     dustemwrap_plot,res,dustem_spectra_st,xr=xr,/xstyle,yr=yr,/ysty,/ylog,/xlog,title=tit+' (From Saved FITS file)'
   ENDELSE
diff --git a/src/idl/dustem_write_fits_table.pro b/src/idl/dustem_write_fits_table.pro
index b30fb20..1072322 100755
--- a/src/idl/dustem_write_fits_table.pro
+++ b/src/idl/dustem_write_fits_table.pro
@@ -63,9 +63,9 @@ ENDELSE
 
 IF !run_pol THEN BEGIN
 ;This is needed to compute Q,U of both SED and spectra
-   toto=dustem_compute_stokes(*(*!dustem_fit).current_param_values,bidon, $
+   toto=dustem_compute_stokes(*(*!dustem_fit).current_param_values, $
                                                  Q_spec=dustem_predicted_Q_spec,U_spec=dustem_predicted_U_spec,PSI_spec=dustem_predicted_PSI_spec, $
-                                                 dustem_psi_em,out_st=out_st)
+                                                 dustem_psi_em=dustem_psi_em)
    	dustem_predicted_Qsed=toto[0]
     dustem_predicted_Used=toto[1]
 ;Same in extinction, but does not work in polarization extinction was not fitted.
--
libgit2 0.21.2