modified for chain run

Jean-Philippe Bernard
1 parent d0f68bb6
Showing 2 changed files with 12 additions and 11 deletions Show diff stats
src/idl/dustem_fit_intensity_mbb_example.pro
src/idl/dustem_fit_sed_ext_pol_example.pro
 PRO dustem_fit_intensity_mbb_example,model=model $
-                                ,sed_file=sed_file $
-                                ,Nitermax=Nitermax $
-                                ,postscript=postscript $
-                                ,fits_save_and_restore=fits_save_and_restore $
-                                 ,help=help $
-                                 ,wait=wait $
-                                 ,verbose=verbose
+                                    ,sed_file=sed_file $
+                                    ,Nitermax=Nitermax $
+                                    ,postscript=postscript $
+                                    ,fits_save_and_restore=fits_save_and_restore $
+                                    ,help=help $
+                                    ,wait=wait $
+                                    ,verbose=verbose
  
 ;+
 ; NAME:
@@ -242,15 +242,16 @@ end
  
 ;=== MAKE THE FINAL PLOT
 IF keyword_set(postscript) THEN BEGIN
-    dir_ps='./'
+    ;dir_ps='./'
     set_plot,'PS'
-    ps_file=dir_ps+postscript
+    ;ps_file=dir_ps+postscript
+    ps_file=postscript
     device,filename=ps_file,/color
 ENDIF
 dustemwrap_plot,*(*!dustem_fit).CURRENT_PARAM_VALUES,dummy,xr=xr,/xstyle,yr=yr,/ysty,/ylog,/xlog,title=tit+' (Final fit)'
 IF keyword_set(postscript) THEN BEGIN
-  set_plot,'X'
   device,/close
+  set_plot,'X'
   message,'Wrote '+ps_file,/info
 ENDIF
  
@@ -469,7 +469,7 @@ IF keyword_set(fits_save_and_restore) THEN BEGIN
 ;=== 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
+  dustem_read_fits_table,filename=fits_save_and_restore,dustem_st=dustem_spectra_st
   ;==== plot result taken from the saved fits table
   res=*(*!dustem_fit).CURRENT_PARAM_VALUES
   IF !dustem_noobj THEN BEGIN