fixed for polarization

Jean-Philippe Bernard
1 parent 378717b3
Showing 2 changed files with 3 additions and 1 deletions Show diff stats
src/idl/dustem_fit_intensity_example.pro
src/idl/dustem_read_fits_table.pro
@@ -382,6 +382,7 @@ IF keyword_set(fits_save_and_restore) THEN BEGIN
    dustem_write_fits_table,filename=fits_save_and_restore,help=help
 ;=== At this point, you could erase all dustem system variables, or exit idl... all the
 ;=== information needed to recover the results and remake the plots has been saved in the FITS table
+  ;stop
   dustem_read_fits_table,filename=fits_save_and_restore,dustem_spectra_st=dustem_spectra_st
   ;==== plot result taken from the saved fits table
   res=*(*!dustem_fit).CURRENT_PARAM_VALUES
@@ -61,7 +61,7 @@ str_predicted_SED=mrdfits(file,unit_predicted_sed,header_predicted_SED)
 ;str_predicted_spectrum_tot=mrdfits(file,3,header_predicted_spectrum_tot)
  
 used_model=strtrim(sxpar(header_predicted_SED,'MODEL'),2)
-used_pol=strtrim(sxpar(header_predicted_SED,'POL'),2)
+used_pol=long(strtrim(sxpar(header_predicted_SED,'POL'),2))
 IF used_model NE 'DEFAULT' THEN BEGIN
 	dustem_init,model=used_model,pol=used_pol
 ENDIF ELSE BEGIN
@@ -134,6 +134,7 @@ ENDFOR
  
 ;=== Below is only to get the same output form as dustem_compute_sed
 ;Note: If the same model is used, the outputs should be the same as what is in the fits file ....
+;stop
 dustem_predicted_sed=dustem_compute_sed(*(*!dustem_fit).current_param_values,out_st=dustem_spectra_st)
  
 ;==== get the predicted emission spectra from extension 4