diff --git a/src/idl/dustem_fit_sed_readme.pro b/src/idl/dustem_fit_sed_readme.pro
index 1b74de1..9d6b3db 100755
--- a/src/idl/dustem_fit_sed_readme.pro
+++ b/src/idl/dustem_fit_sed_readme.pro
@@ -80,12 +80,10 @@ IF keyword_set(png) THEN BEGIN
 ENDIF
 
 polarization=0.   ;default is no polarization in models
+;dustem_init,model=use_model,pol=polarization
 
 ;COMMENT: for some reson polarization=1 doesn't work for ALL of the models except for the two last A and B models
 
-dustem_init,model=use_model,pol=polarization
-!dustem_nocatch=1
-
 ;=== Set which parameters you want to fit
 CASE use_model OF
       'DBP90':BEGIN
@@ -349,6 +347,9 @@ CASE use_model OF
 
 ENDCASE
 
+dustem_init,model=use_model,pol=polarization
+!dustem_nocatch=1
+
 ;dustem_init,model=use_model,pol=polarization
 
 !dustem_verbose=1
@@ -362,10 +363,6 @@ file=dir+'Gal_composite_spectrum.xcat'
 ;file=dir+'Gal_composite_spectrum.xcat.old'
 
 spec=read_xcat(file,/silent)
-stop
-;spec.stokesQ=0.
-;spec.stokesU=0.
-;spec.stokesU=0.
 
 ;ind=where(spec.error EQ 0.,count)
 ;IF count NE 0 THEN spec(ind).error=0.2*spec(ind).spec
--
libgit2 0.21.2