PRO dustemwrap_plot_noobj,p_dim,st, $ dustem_sed,SED_spec, $ dustem_qsed,Q_spec, $ dustem_used,U_spec, $ dustem_polsed,P_spec, $ dustem_polfrac,SP_spec, $ dustem_psi_em,PSI_spec, $ dustem_ext,EXT_spec, $ dustem_qext,QEXT_spec, $ dustem_uext,UEXT_spec, $ dustem_polext,POLEXT_spec, $ dustem_fpolext,SPEXT_spec, $ dustem_psi_ext,PSIEXT_spec, $ _extra=_extra,help=help ;+ ; NAME: ; dustemwrap_plot_noobj ; PURPOSE: ; plots dustemwrap observed sed, predicted sed, predicted spectra for total Intensity, polarization in emission and extinction without OOP routines ; CATEGORY: ; DustEMWrap ; CALLING SEQUENCE: ; dustemwrap_plot_noobj,p_dim,st,dustem_sed,SED_spec,dustem_qsed,Q_spec,dustem_used,U_spec,dustem_polsed,P_spec,dustem_polfrac,SP_spec,dustem_psi_em,PSI_spec,dustem_ext,EXT_spec,dustem_qext,QEXT_spec,dustem_uext,UEXT_spec,dustem_polext,POLEXT_spec,dustem_fpolext,SPEXT_spec,dustem_psi_ext,PSIEXT_spec,_extra=_extra ; INPUTS: ; p_dim ; OPTIONAL INPUT PARAMETERS: ; st ; dustem_sed ; dustem_qsed ; dustem_used ; dustem_polsed ; dustem_polfrac ; dustem_psi_em ; dustem_ext ; dustem_qext ; dustem_uext ; dustem_polext ; dustem_fpolext ; dustem_psi_ext ; OUTPUTS: ; dustem_sed ; dustem_qsed ; dustem_used ; dustem_polsed ; dustem_polfrac ; dustem_psi_em ; dustem_ext ; dustem_qext ; dustem_uext ; dustem_polext ; dustem_fpolext ; dustem_psi_ext ; SED_spec ; Q_spec ; U_spec ; P_spec ; SP_spec ; PSI_spec ; EXT_spec ; QEXT_spec ; UEXT_spec ; POLEXT_spec ; SPEXT_spec ; PSIEXT_spec ; OPTIONAL OUTPUT PARAMETERS: ; ; ACCEPTED KEY-WORDS: ; help ; COMMON BLOCKS: ; None ; SIDE EFFECTS: ; Does a plot ; RESTRICTIONS: ; The DustEM fortran code must be installed ; The DustEMWrap IDL code must be installed ; PROCEDURES AND SUBROUTINES USED: ; ; EXAMPLES ; ; MODIFICATION HISTORY: ; Written by Ilyes Choubani ; Evolution details on the DustEMWrap gitlab. ;- IF keyword_set(help) THEN BEGIN doc_library,'dustemwrap_plot' goto,the_end END ;This is also necessary for the plotting of the results of the fit (Last iteration) IF not keyword_set(st) THEN BEGIN ;Activation of the plugins is needed because of the plotting of the total emission model that includes them. dustem_activate_plugins,p_dim/(*(*!dustem_fit).param_init_values),st=st ; 0/0 division case? ENDIF ;# Emission if not keyword_set(dustem_sed) and isa((*!dustem_data).sed) then begin ;better test? dustem_sed = dustem_compute_sed(p_dim,st=st,SED_spec=SED_spec) 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,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] endif endif ;# Extinction if not keyword_set(dustem_ext) and isa((*!dustem_data).ext) then begin dustem_ext = dustem_compute_ext(p_dim,st=st,EXT_spec=EXT_spec) 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,QEXT_spec=QEXT_spec,UEXT_spec=UEXT_spec,PSIEXT_spe=PSIEXT_spec,dustem_psi_ext=dustem_psi_ext) dustem_qext = toto[0] dustem_uext = toto[1] endif endif mouchard=*!iteration_mouchard new_mouchard=mouchard[0] new_mouchard.iteration=!dustem_iter.act new_mouchard.chi2=(*!dustem_fit).chi2 new_mouchard.rchi2=(*!dustem_fit).rchi2 ;JPB: need here to store values of current parameter values (including plugins) into !iteration_mouchard ;help,(*!dustem_fit).CURRENT_PARAM_VALUES ;This is not currently filled up ... ;Below will work only for non plugin parameters ;goto,skip1 ;param_descs=*(*!dustem_fit).param_descs ;params_current=[0.D0] ;FOR i=0L,n_elements(param_descs)-1 DO BEGIN ; str='params_current=[params_current,'+param_descs[i]+']' ; print,str ; toto=execute(str) ;ENDFOR ;params_current=params_current[1:*] ;print,params_current ; ;stop ;skip1: mouchard=[mouchard,new_mouchard] !iteration_mouchard=ptr_new(mouchard) ;LIST OF TESTS THAT NEED TO BE RAN SO THAT THE PLOTTING OF THE MODEL SPECTRA OCCURS test_sed = isa((*!dustem_show).sed) test_ext = isa((*!dustem_show).ext) if !run_pol then begin test_polext = isa((*!dustem_show).polext) test_polsed = isa((*!dustem_show).polsed) test_polfrac = isa((*!dustem_show).polfrac) test_fpolext = isa((*!dustem_show).fpolext) test_psi_em = isa((*!dustem_show).psi_em) test_psi_ext = isa((*!dustem_show).psi_ext) test_qsed = isa((*!dustem_show).qsed) test_used = isa((*!dustem_show).used) test_qext = isa((*!dustem_show).qext) test_uext = isa((*!dustem_show).uext) test_m = test_sed or test_polsed or test_polfrac or test_qsed or test_used or test_psi_em test_x = test_ext or test_polext or test_fpolext or test_qext or test_uext or test_psi_ext endif else begin test_m = test_sed test_x = test_ext endelse ;testing on the fitting of emission vs extinction data n_plgns = n_tags(*!dustem_plugin) frmt0='(A36)' frmt1='(1E10.2)' frmt2='(10F10.2)' ;Instead of setting the extra tags for plotting in dustemcgwin, we'll place them in this (parent) procedure. tgnms_extra = tag_names(_extra) dustem_set_plot_range,test_sed, test_ext,_extra= _extra if isa((*!dustem_fit).current_param_errors) then begin errors = (*(*!dustem_fit).current_param_errors) endif else begin errors = (*(*!dustem_fit).param_init_values)*0. ENDELSE if isa((*!dustem_fit).chi2) then begin chi2 = (*!dustem_fit).chi2 endif else begin chi2 = la_undef() endelse if isa((*!dustem_fit).rchi2) then begin rchi2 = (*!dustem_fit).rchi2 endif else begin rchi2 = la_undef() endelse ;==== These are arameters for window positioning if !run_pol then BEGIN wdelta_x=60 main_wxsize=910 main_wysize=658 param_wxpos=main_wxsize+wdelta_x param_wypos=15 param_wxsize=350 param_wysize=200 plugin_wxpos=main_wxsize+wdelta_x plugin_wypos=300 plugin_wxsize=350 plugin_wysize=200 ENDIF ELSE BEGIN wdelta_x=60 main_wxsize=640 ;default of cgwindow main_wysize=512 ;default of cgwindow param_wxpos=main_wxsize+wdelta_x;700 param_wypos=15 param_wxsize=350 param_wysize=200 plugin_wxpos=main_wxsize+wdelta_x;700 plugin_wypos=300 plugin_wxsize=350 plugin_wysize=200 ENDELSE ;JPB: I moved these 3 lines out of the if test, so that they apply also ;when only extinction is present TVLCT, 255, 255, 255, 254 ; White color !P.Color = 253 !P.Background = 254 if test_m then begin ;Generating emission window and saving window ID IF !d.name NE 'PS' THEN BEGIN IF !run_pol THEN BEGIN window,0, xsize=main_wxsize,ysize=main_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (EMISSION)',xpos=0,ypos=0 ENDIF else BEGIN window,0,xsize=main_wxsize,ysize=main_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (EMISSION)',xpos=0, ypos=0;, wback='grey';,wobject=winobj_m ENDELSE ENDIF ;Plotting model string (legend) xxpos = 0.015 yypos = 0.97 ;Plotting Hydrogen column density string (legend) cgtext, xxpos, yypos,'Model : '+!dustem_model,color=0,/normal,charsize = 1.2;,/addcmd ;& cmdind_m+=1 This is the first command in the window yypos = 0.95 cgtext, xxpos, yypos,textoidl('N_{H}')+textoidl(' (H/cm^2) : ')+strmid(string(*!dustem_HCD,format='(1E10.2)'),2),color=0,/normal,charsize = 1.2;,/addcmd & cmdind_m+=1 ;Plotting of the number of iterations completed by the fitter ;Plotting of the number of runs ;Getting the width of the unrefreshed string widthtot=0. ;cgtext, xxpos, yypos,'Runs : ',color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_m+=1 cgtext, xxpos, yypos,'Iteration : ',color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width !dustemcgwin_ncmds.runs.txtwdth = thiswidth ;Getting the width of the refreshed string cgtext, xxpos, yypos,strtrim(!dustem_iter.act,2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;These two lines will also have to be modified when plotting the refreshed values xxpos = 0.015 ;1.0 - widthtot - 0.03 yypos = 0.015 ;Displaying fixed part ;cgtext, xxpos, yypos,'Runs : ',color=0,/normal,charsize = 1.35,/addcmd & cmdind_m+=1 cgtext, xxpos, yypos,'Iteration : ',color=0,/normal,charsize = 1.35 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext*1.35 + 0.015 , yypos,strtrim(!dustem_iter.act,2),color=0,/normal,charsize = 1.0 ;plotting of (mostly) fixed data for rchi2 ;Since the format of the strings is specified, the length/width of the string in thex direction is fixed. ;Getting the width of the unrefreshed string widthtot=0. cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width ;Getting the width of the refresh.ed string cgtext, xxpos, yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;saving the command string width to replace it xxpos = 1.0 - widthtot - 0.03 yypos = 0.97 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize = 1.0 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext , yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize = 1.0 ;plotting of (mostly) fixed data for rchi2 ;Getting the width of the unrefreshed string widthtot=0. cgtext, xxpos, yypos,textoidl('r_{\chi^{2}} : '),color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width ;!dustemcgwin_ncmds.rchi2.txtwdth1 = thiswidth ;Getting the width of the refreshed string cgtext, xxpos, yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;These two lines will also have to be modified when plotting the refreshed values xxpos = 1.0 - widthtot - 0.03 yypos = 0.95 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('r_{\chi^{2}} : '),color=0,/normal,charsize = 1.35 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext , yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize = 1.0 ;saving the command id to replace it ;Plotting of the title of the dashboard (plot(s)) ;NB issue here because the user can't modify the title after the dustem_mpfit_data call at the main procedure level (user proc) ;Getting the width of the title tit=!dustem_plot_range.title_m xxpos = 0.5 yypos =0.96 ;the runs inbetween if !dustem_end EQ 0 then cgtext, xxpos , yypos,tit+' (RUNNING)',color=0,/normal,charsize = 1.6,alignment=0.5 ;subsequent iterations to the end of the fit if !dustem_end EQ -1 then cgtext, xxpos , yypos,tit,color=0,/normal,charsize = 1.6,alignment=0.5 ;last iteration of the fit if !dustem_end EQ 1 then cgtext, xxpos , yypos,tit+' (FINAL RUN)',color=0,/normal,charsize = 1.6,alignment=0.5 if !run_pol then begin ;,wback='grey' ;position arrays for plot and normalized graph ;THIS IS HARD-CODED IN A VERY SILLY WAY p_sed = [0.05,0.67,0.475,0.90] p_psed = [0.525,0.67,0.95,0.90] p_spsed = [0.05,0.375,0.475,0.575] p_psised = [0.525,0.375,0.95,0.575] p_qsed = [0.05,0.12,0.475,0.35] p_used = [0.525,0.12,0.95,0.35] np_sed = [0.05,0.60,0.475,0.67] np_psed = [0.525,0.60,0.95,0.67] np_qsed = [0.05,0.05,0.475,0.12] np_used = [0.525,0.05,0.95,0.12] p_sed[1]+= 0.0125 & p_sed[3]+= 0.0125 p_psed[1]+= 0.0125 & p_psed[3]+= 0.0125 p_spsed[1]+= 0.0125 & p_spsed[3]+= 0.0125 p_psised[1]+= 0.0125 & p_psised[3]+= 0.0125 p_qsed[1]+= 0.0125 & p_qsed[3]+= 0.0125 p_used[1]+= 0.0125 & p_used[3]+= 0.0125 np_sed[1]+= 0.0125 & np_sed[3]+= 0.0125 np_psed[1]+= 0.0125 & np_psed[3]+= 0.0125 np_qsed[1]+= 0.0125 & np_qsed[3]+= 0.0125 np_used[1]+= 0.0125 & np_used[3]+= 0.0125 endif else begin p_sed = [0.07,0.30,0.97,0.85] np_sed = [0.07,0.10,0.97,0.30] endelse ;Plotting of sed data axes (whether or not data is present) plotsym,0,/fill ;you might need to execute this again. hmm... ;stop ;This is not the right test the user might want to show the prediction for sed and fit the stokes parameters if ~test_sed then begin ; SED data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='SED', /nodata, position=p_sed, _extra=_extra dustem_plot_dataset, st, dataset='SED', /nodata, /norm, position=np_sed, _extra=_extra endif else begin ;SED exists ;this plots the observed SED dustem_plot_dataset, st, dataset='SED', position=p_sed, _extra=_extra ;this plots the predicted spectrum and predicted SED dustem_plot_dataset, st, dustem_sed,SED_spec, dataset='SED', /refresh ,position=p_sed, _extra=_extra ;This plots an empty normalized plot dustem_plot_dataset, st, dataset='SED',/norm, position=np_sed, _extra=_extra ;I think nodata can also be used too ;This plots the normalized observed SED - best fit dustem_plot_dataset, st, dustem_sed,SED_spec, dataset='SED', /norm, /refresh,position=np_sed, _extra=_extra endelse if !run_pol then begin ;then goto, end_m if ~test_qsed then begin ;Qsed is not present dustem_plot_dataset, st, dataset='QSED', /nodata, position=p_qsed, _extra=_extra dustem_plot_dataset, st, dataset='QSED', /nodata, /norm, position=np_qsed, _extra=_extra endif else begin ;Qsed data is present ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendered) dustem_plot_dataset, st, dataset='QSED', position=p_qsed,/positive_only, _extra=_extra ;Refreshing the positive_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_qsed,Q_spec, dataset='QSED',/positive_only, /refresh ,position=p_qsed, _extra=_extra dustem_plot_dataset, st, dataset='QSED', position=p_qsed, /addcmd,/negative_only, _extra=_extra ;Refreshing the negative_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_qsed,Q_spec,dataset='QSED',/negative_only, /refresh ,position=p_qsed, _extra=_extra dustem_plot_dataset, st, dataset='QSED',/norm, position=np_qsed, _extra=_extra ;I think nodata can also be uqsed too dustem_plot_dataset, st, dustem_qsed,Q_spec,dataset='QSED', /norm, /refresh,position=np_qsed, _extra=_extra endelse if ~test_used then begin ;used is not present dustem_plot_dataset, st, dataset='USED', /nodata, position=p_used, _extra=_extra dustem_plot_dataset, st, dataset='USED', /nodata, /norm, position=np_used, _extra=_extra endif else begin ;used data is present ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendred) dustem_plot_dataset, st, dataset='USED', position=p_used,/positive_only, _extra=_extra ;Refreshing the positive_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_used,U_spec, dataset='USED',/positive_only, /refresh ,position=p_used, _extra=_extra dustem_plot_dataset, st, dataset='USED', position=p_used, winid=winid_m,/negative_only, _extra=_extra ;Refreshing the negative_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_used,U_spec, dataset='USED',/negative_only, /refresh ,position=p_used, _extra=_extra dustem_plot_dataset, st, dataset='USED',/norm, position=np_used, _extra=_extra ;I think nodata can also be uused too dustem_plot_dataset, st, dustem_used,U_spec, dataset='USED', /norm, /refresh,position=np_used, _extra=_extra endelse if ~test_psi_em then begin ; psi_em data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='psi_em', /nodata, position=p_psised, /addcmd, winid=winid_m, _extra=_extra endif else begin ;psi_em exists dustem_plot_dataset, st, dataset='psi_em', position=p_psised, _extra=_extra dustem_plot_dataset, st, dustem_psi_em,PSI_spec, dataset='psi_em', /refresh ,position=p_psised, _extra=_extra endelse if ~test_polsed then begin ; polsed data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='polsed', /nodata, position=p_psed, _extra=_extra dustem_plot_dataset, st, dataset='polsed', /nodata, /norm, position=np_psed, _extra=_extra endif else begin dustem_plot_dataset, st, dataset='polsed', position=p_psed, _extra=_extra dustem_plot_dataset, st, dustem_polsed,P_spec, dataset='polsed', /refresh ,position=p_psed, _extra=_extra dustem_plot_dataset, st, dataset='polsed',/norm, position=np_psed, _extra=_extra ;I think nodata can also be upolsed too dustem_plot_dataset, st, dustem_polsed,P_spec, dataset='polsed', /norm, /refresh,position=np_psed, _extra=_extra endelse if ~test_polfrac then begin ; polfrac data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='polfrac', /nodata, position=p_spsed, _extra=_extra endif else begin ;polfrac exists dustem_plot_dataset, st, dataset='polfrac', position=p_spsed, _extra=_extra dustem_plot_dataset, st, dustem_polfrac,SP_spec,SED_spec, dataset='polfrac', /refresh ,position=p_spsed, _extra=_extra endelse ENDIF ENDIF ;stop if test_x then begin ;Generating extinction window and saving window ID if !run_pol then window,1,xsize=main_wxsize,ysize=main_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (EXTINCTION)' else window,1,xsize=main_wxsize,ysize=main_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (EXTINCTION)';, wback='grey';,wobject=winobj_m ;Plotting model string (legend) xxpos = 0.015 yypos = 0.97 ;Plotting Hydrogen column density string (legend) cgtext, xxpos, yypos,'Model : '+!dustem_model,color=0,/normal,charsize = 1.2;,/addcmd ;& cmdind_m+=1 This is the first command in the window yypos = 0.95 cgtext, xxpos, yypos,textoidl('N_{H}')+textoidl(' (H/cm^2) : ')+strmid(string(*!dustem_HCD,format='(1E10.2)'),2),color=0,/normal,charsize = 1.2 ;Plotting of the number of iterations completed by the fitter ;Plotting of the number of runs ;Getting the width of the unrefreshed string widthtot=0. ;cgtext, xxpos, yypos,'Runs : ',color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_x+=1 cgtext, xxpos, yypos,'Iteration : ',color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width ;Getting the width of the refreshed string cgtext, xxpos, yypos,strtrim(!dustem_iter.act,2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;These two lines will also have to be modified when plotting the refreshed values xxpos = 0.015 ;1.0 - widthtot - 0.03 yypos = 0.015 ;Displaying fixed part ;cgtext, xxpos, yypos,'Runs : ',color=0,/normal,charsize = 1.35,/addcmd & cmdind_x+=1 cgtext, xxpos, yypos,'Iteration : ',color=0,/normal,charsize = 1.35 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext + 0.015 , yypos,strtrim(!dustem_iter.act,2),color=0,/normal,charsize = 1.0 ;saving the command id to replace it ;plotting of (mostly) fixed data for rchi2 ;Since the format of the strings is specified, the length/width of the string in thex direction is fixed. ;Getting the width of the unrefreshed string widthtot=0. cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width ;Getting the width of the refresh.ed string cgtext, xxpos, yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;saving the command string width to replace it xxpos = 1.0 - widthtot - 0.03 yypos = 0.97 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize = 1.0 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext , yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize = 1.0 ;saving the command id to replace it ;plotting of (mostly) fixed data for rchi2 ;Getting the width of the unrefreshed string widthtot=0. cgtext, xxpos, yypos,textoidl('r_{\chi^{2}} : '),color=0,/normal,charsize=-1,width=thiswidth widthtext = thiswidth widthtot+= thiswidth ;fixed width ;!dustemcgwin_ncmds.rchi2.txtwdth1 = thiswidth ;Getting the width of the refreshed string cgtext, xxpos, yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize=-1,width=thiswidth widthtot+= thiswidth;varying width ;These two lines will also have to be modified when plotting the refreshed values xxpos = 1.0 - widthtot - 0.03 yypos = 0.95 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('r_{\chi^{2}} : '),color=0,/normal,charsize = 1.35 ;Displaying refreshed part (first run . This is the command that will be replaced) cgtext, xxpos + widthtext , yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize = 1.0 ;saving the command id to replace it ;plotting of title tit=!dustem_plot_range.title_x xxpos = 0.5 yypos =0.96 ;the runs inbetween if !dustem_end EQ 0 then cgtext, xxpos , yypos,tit+' (RUNNING)',color=0,/normal,charsize = 1.6,alignment=0.5 ;subsequent iterations to the end of the fit if !dustem_end EQ -1 then cgtext, xxpos , yypos,tit,color=0,/normal,charsize = 1.6,alignment=0.5 ;last iteration of the fit if !dustem_end EQ 1 then cgtext, xxpos , yypos,tit+' (FINAL RUN)',color=0,/normal,charsize = 1.6,alignment=0.5 if !run_pol then begin ;position arrays for plot and normalized graph ;THIS IS HARD-CODED IN A VERY SILLY WAY p_ext = [0.05,0.67,0.475,0.90] p_pext = [0.525,0.67,0.95,0.90] p_spext = [0.05,0.375,0.475,0.575] ; Should we plot this? tau_polext/tau_ext_tot. YES FOR NOW... p_psiext = [0.525,0.375,0.95,0.575] p_qext = [0.05,0.12,0.475,0.35] p_uext = [0.525,0.12,0.95,0.35] np_ext = [0.05,0.60,0.475,0.67] np_pext = [0.525,0.60,0.95,0.67] np_qext = [0.05,0.05,0.475,0.12] np_uext = [0.525,0.05,0.95,0.12] p_ext[1]+= 0.0125 & p_ext[3]+= 0.0125 p_pext[1]+= 0.0125 & p_pext[3]+= 0.0125 p_spext[1]+= 0.0125 & p_spext[3]+= 0.0125 p_psiext[1]+= 0.0125 & p_psiext[3]+= 0.0125 p_qext[1]+= 0.0125 & p_qext[3]+= 0.0125 p_uext[1]+= 0.0125 & p_uext[3]+= 0.0125 np_ext[1]+= 0.0125 & np_ext[3]+= 0.0125 np_pext[1]+= 0.0125 & np_pext[3]+= 0.0125 np_qext[1]+= 0.0125 & np_qext[3]+= 0.0125 np_uext[1]+= 0.0125 & np_uext[3]+= 0.0125 endif else begin p_ext = [0.07,0.30,0.97,0.85] np_ext = [0.07,0.10,0.97,0.30] endelse ;Plotting of sed data axes (whether or not data is present) plotsym,0,/fill ;you might need to execute this again. hmm... ;This is not the right test the user might want to show the prediction for sed and fit the stokes parameters if ~test_ext then begin ; EXT data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='EXT', /nodata, position=p_ext, _extra=_extra dustem_plot_dataset, st, dataset='EXT', /nodata, /norm, position=np_ext, _extra=_extra endif else begin ;EXT exists dustem_plot_dataset, st, dataset='EXT', position=p_ext, _extra=_extra dustem_plot_dataset, st, dustem_ext,EXT_spec, dataset='EXT', /refresh ,position=p_ext, _extra=_extra dustem_plot_dataset, st, dataset='EXT',/norm, position=np_ext, _extra=_extra ;I think nodata can also be used too dustem_plot_dataset, st, dustem_ext,EXT_spec, dataset='EXT', /norm, /refresh,position=np_ext, _extra=_extra endelse if !run_pol then begin if ~test_qext then begin ;Qext is not present dustem_plot_dataset, st, dataset='QEXT', /nodata, position=p_qext, _extra=_extra dustem_plot_dataset, st, dataset='QEXT', /nodata, /norm, position=np_qext, _extra=_extra endif else begin ;Qext data is present ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendered) dustem_plot_dataset, st, dataset='QEXT', position=p_qext,/positive_only, _extra=_extra ;Refreshing the positive_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_qext,QEXT_spec, dataset='QEXT',/positive_only, /refresh ,position=p_qext, _extra=_extra dustem_plot_dataset, st, dataset='QEXT', position=p_qext,/negative_only, _extra=_extra ;Refreshing the negative_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_qext,QEXT_spec,dataset='QEXT',/negative_only, /refresh ,position=p_qext, _extra=_extra dustem_plot_dataset, st, dataset='QEXT',/norm, position=np_qext, _extra=_extra ;I think nodata can also be uqsed too dustem_plot_dataset, st, dustem_qext,QEXT_spec,dataset='QEXT', /norm, /refresh,position=np_qext, _extra=_extra endelse ;stop if ~test_uext then begin ;uext is not present dustem_plot_dataset, st, dataset='UEXT', /nodata, position=p_uext, _extra=_extra dustem_plot_dataset, st, dataset='UEXT', /nodata, /norm, position=np_uext, _extra=_extra endif else begin ;uext data is present ;stop ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendred) dustem_plot_dataset, st, dataset='UEXT', position=p_uext,/positive_only, _extra=_extra ;Refreshing the positive_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_uext,UEXT_spec, dataset='UEXT',/positive_only, /refresh ,position=p_uext, _extra=_extra dustem_plot_dataset, st, dataset='UEXT', position=p_uext,/negative_only, _extra=_extra ;Refreshing the negative_only plot ONLY if there is data. If not do nothing dustem_plot_dataset, st, dustem_uext,UEXT_spec, dataset='UEXT',/negative_only, /refresh ,position=p_uext, _extra=_extra dustem_plot_dataset, st, dataset='UEXT',/norm, position=np_uext, _extra=_extra ;I think nodata can also be uused too dustem_plot_dataset, st, dustem_uext,UEXT_spec, dataset='UEXT', /norm, /refresh,position=np_uext, _extra=_extra endelse if ~test_psi_ext then begin ; psi_ext data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='PSI_EXT', /nodata, position=p_psiext, _extra=_extra endif else begin ;psi_ext exists dustem_plot_dataset, st, dataset='psi_ext', position=p_psiext, _extra=_extra dustem_plot_dataset, st, dustem_psi_ext,PSIEXT_spec, dataset='psi_ext', /refresh ,position=p_psiext, _extra=_extra endelse if ~test_polext then begin ; polext data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='polext', /nodata, position=p_pext, _extra=_extra dustem_plot_dataset, st, dataset='polext', /nodata, /norm, position=np_pext, _extra=_extra endif else begin ;polext exists dustem_plot_dataset, st, dataset='polext', position=p_pext, _extra=_extra dustem_plot_dataset, st, dustem_polext,POLEXT_spec, dataset='polext', /refresh ,position=p_pext, _extra=_extra dustem_plot_dataset, st, dataset='polext',/norm, position=np_pext, _extra=_extra ;I think nodata can also be upolsed too dustem_plot_dataset, st, dustem_polext,POLEXT_spec, dataset='polext', /norm, /refresh,position=np_pext, _extra=_extra endelse if ~test_fpolext then begin ; fpolext data is not present (For completeness) # This is only valid when !run_pol=1 dustem_plot_dataset, st, dataset='fpolext', /nodata, position=p_spext, _extra=_extra endif else begin ;fpolext exists dustem_plot_dataset, st, dataset='fpolext', position=p_spext, _extra=_extra dustem_plot_dataset, st, dustem_fpolext,SPEXT_spec,EXT_spec, dataset='fpolext', /refresh ,position=p_spext, _extra=_extra endelse ENDIF ENDIF if !dustem_end then !dustem_end = -1 ;so that the last iteration is only plotted once. ;PLOTTING OF THE PARAMETERS IF !d.name NE 'PS' THEN BEGIN window,2,xpos=param_wxpos,ypos=param_wypos,xsize=param_wxsize,ysize=param_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (PARAMETERS)' ; keeping the same y dimension of the emission/extinction plot? ENDIF dustem_plot_dataset, st, p_dim, errors,dataset='PARAMETERS',_extra=_extra ;& cmdind_prms+=1 dustem_plot_dataset, st, p_dim,errors,dataset='PARAMETERS',/refresh, _extra=_extra ;PLOTTING OF THE PLUGINS IF (tag_names(*!dustem_plugin))(0) NE 'NONE' THEN BEGIN IF !d.name NE 'PS' THEN BEGIN window,3,xpos=plugin_wxpos,ypos=plugin_wypos,xsize=plugin_wxsize,ysize=plugin_wysize,title='DUSTEMWRAP '+!dustem_version.version+' (PLUGINS)' ; keeping the same y dimension of the emission/extinction plot? ENDIF dustem_plot_dataset, st, p_dim, errors,dataset='PLUGINS',_extra=_extra ;& cmdind_prms+=1 dustem_plot_dataset, st, p_dim, errors,dataset='PLUGINS', /refresh , winid=winid_plgns, _extra=_extra ENDIF the_end: END