PRO dustemwrap_plot,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 ;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 ; 0/0 division case? ;# Emission ;if isa(!dustem_data.sed) then begin ;using !dustem_data instead of !dustem_show because in this release the user is compelled to fit IQU whether in M or X if not keyword_set(dustem_sed) and isa(!dustem_data.sed) then begin ;better test? dustem_sed = dustem_compute_sed(p_dim,st,SED_spec) if !run_pol && !run_lin then begin dustem_polsed = dustem_compute_polsed(p_dim,st,P_spec,SP_spec,dustem_polfrac) toto = dustem_compute_stokes(p_dim,st,dustem_qsed,dustem_used,Q_spec,U_spec,PSI_spec,dustem_psi_em) ;stop endif endif ;# Extinction ; if isa(!dustem_data.ext) then begin if not keyword_set(dustem_ext) and isa(!dustem_data.ext) then begin dustem_ext = dustem_compute_ext(p_dim,st,EXT_spec) if !run_pol && !run_lin then begin dustem_polext = dustem_compute_polext(p_dim,st,POLEXT_spec,SPEXT_spec,dustem_fpolext) toto = dustem_compute_stokext(p_dim,st,dustem_qext,dustem_uext,QEXT_spec,UEXT_spec,PSIEXT_spec,dustem_psi_ext) endif endif ENDIF ELSE BEGIN ;st is provided ;stop ;# Emission ;if isa(!dustem_data.sed) then begin if not keyword_set(dustem_sed) and isa(!dustem_data.sed) then begin dustem_sed = dustem_compute_sed(p_dim,st,SED_spec) if !run_pol && !run_lin then begin dustem_polsed = dustem_compute_polsed(p_dim,st,P_spec,SP_spec,dustem_polfrac) toto = dustem_compute_stokes(p_dim,st,dustem_qsed,dustem_used,Q_spec,U_spec,PSI_spec,dustem_psi_em) ;stop endif endif ;# Extinction ;if isa(!dustem_data.ext) then begin if not keyword_set(dustem_ext) and isa(!dustem_data.ext) then begin dustem_ext = dustem_compute_ext(p_dim,st,EXT_spec) if !run_pol && !run_lin then begin dustem_polext = dustem_compute_polext(p_dim,st,POLEXT_spec,SPEXT_spec,dustem_fpolext) toto = dustem_compute_stokext(p_dim,st,dustem_qext,dustem_uext,QEXT_spec,UEXT_spec,PSIEXT_spec,dustem_psi_ext) endif endif ENDELSE ;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)' ;I have decided to keep track of the command index manually as I did not find any command online that does this. cmdind_m = 0 ;for emission cmdind_x = 0 ;for extinction ;trying to do the same thing for parameters and plugins cmdind_prms = 0 ;for parameters cmdind_plgns = 0 ;for plugins ;ADDING PLUGIN(S) TO SPECTRUM---------------- ;these following indices will be used to change the keep the same information in the _extra structure ;this is an initial solution until I finish coding the _extra-filtering procedure. iswinsed = !dustemcgwin_id.sed EQ la_undef() iswinext = !dustemcgwin_id.ext EQ la_undef() ;trying to do the same thing for parameters and plugins ;But we need to know if the extinction or the emission mode is on. ;so the use of ths pointer will be different than the others iswinprms = !dustemcgwin_id.prms EQ la_undef() iswinplgns = !dustemcgwin_id.plgns EQ la_undef() ;plotstrct = !p degtorad = !pi/180 ;factor that the arctan will be devided by in order to have an axis in degrees. Talk to JP maybe he wants angle in radians. ;I think I will test over the presence of the errors because they are not computed at that point. ;errors = ((*!dustem_fit).current_param_errors) if isa((*!dustem_fit).current_param_errors) then begin errors = (*(*!dustem_fit).current_param_errors) * (*(*!dustem_fit).param_init_values) endif else begin errors = (*(*!dustem_fit).param_init_values)*0.;+la_undef();Initialized to 0 instead of la_undef() ENDELSE if isa((*!dustem_fit).chi2) then begin chi2 = (*!dustem_fit).chi2 endif else begin chi2 = la_undef();maybe not the best initialization endelse if isa((*!dustem_fit).rchi2) then begin rchi2 = (*!dustem_fit).rchi2 endif else begin rchi2 = la_undef();maybe not the best initialization endelse if isa(p_dim) then begin if isa((*!dustem_fit).current_param_values) then begin res = (*(*!dustem_fit).current_param_values) endif else begin res = p_dim endelse endif else begin res = (*(*!dustem_fit).param_init_values)*0.;+la_undef() endelse ;stop if test_m then begin ;Generating emission window and saving window ID if iswinsed then begin ;This is when the cgwindow does not exist already if !run_pol then cgwindow,wxsize=910,wysize=658,wtitle='DUSTEMWRAP v2.0 (EMISSION)' else cgwindow, wtitle='DUSTEMWRAP v2.0 (EMISSION)';, wback='grey';,wobject=winobj_m winid_m = cgquery(dimensions=dim_m,/current) ;this has to be changed because the dimensions of this array change with the data sets present (extinction/emission). !dustemcgwin_id.sed = winid_m cgcontrol, winid_m, execute=0 ;cgset, winid_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,/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,/addcmd & cmdind_m+=1 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,/addcmd & cmdind_m+=1 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,/addcmd & cmdind_m+=1 ;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,/addcmd & cmdind_m+=1 ;saving the command id to replace it !dustemcgwin_ncmds.runs.pl = cmdind_m ;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,/addcmd & cmdind_m+=1 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,/addcmd & cmdind_m+=1 widthtot+= thiswidth;varying width ;saving the command string width to replace it !dustemcgwin_ncmds.chi2.txtwdth = thiswidth xxpos = 1.0 - widthtot - 0.03 yypos = 0.97 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize = 1.0,/addcmd & cmdind_m+=1 ;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,/addcmd & cmdind_m+=1 ;saving the command id to replace it !dustemcgwin_ncmds.chi2.pl = cmdind_m ;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,/addcmd & cmdind_m+=1 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,/addcmd & cmdind_m+=1 widthtot+= thiswidth;varying width !dustemcgwin_ncmds.rchi2.txtwdth = thiswidth ;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,/addcmd & cmdind_m+=1 ;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,/addcmd & cmdind_m+=1 ;saving the command id to replace it !dustemcgwin_ncmds.rchi2.pl = cmdind_m ;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=string(_extra.title); cgtext, 0.3, 0.95,tit+' (RUNNING)',color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_m+=1 widthtot = thiswidth xxpos = (1 - thiswidth)/2 yypos = 0.96 cgtext, 0.3, 0.95,tit,color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_m+=1 !dustemcgwin_ncmds.pltit.txtwdth = thiswidth*1.6 + (1 - widthtot)/2 cgtext, xxpos, yypos,tit,color=0,/normal,charsize = 1.6,/addcmd & cmdind_m+=1 xxpos+=thiswidth*1.6 cgtext, xxpos, yypos,' (RUNNING)',color=0,/normal,charsize = 1.6,/addcmd & cmdind_m+=1 !dustemcgwin_ncmds.pltit.pl = cmdind_m endif else begin ;This is when the cgwindow exists already winid_m = !dustemcgwin_id.sed cgcontrol, winid_m, execute=0 cgset, winid_m ;apparently specifying the window index isn't enough. This might be because of the use of cgcontrol ;Refreshing chi2 data: cmdin_m = !dustemcgwin_ncmds.chi2.pl xxpos = 1 - !dustemcgwin_ncmds.chi2.txtwdth - 0.03 yypos = 0.97 cgwindow,'cgtext', xxpos , yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdin_m, winid=winid_m ;& cmdind_m+=1 ;Refreshing rchi2 data: cmdin_m = !dustemcgwin_ncmds.rchi2.pl xxpos = 1 - !dustemcgwin_ncmds.rchi2.txtwdth - 0.03 yypos = 0.95 cgwindow,'cgtext', xxpos , yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdin_m, winid=winid_m ;& cmdind_m+=1 ;Refreshing number of runs data: cmdin_m = !dustemcgwin_ncmds.runs.pl xxpos = 0.03 + !dustemcgwin_ncmds.runs.txtwdth*1.35 ;(because of charsize) yypos = 0.015 if !dustem_iter.act eq 0 then txt = 'Fit failed. ' else txt = strtrim(!dustem_iter.act,2) cgwindow,'cgtext', xxpos , yypos,txt,color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdin_m, winid=winid_m ;& cmdind_m+=1 ;update of the title here ? if isa((*!dustem_fit).current_param_values) then begin ;replacing my count - 3 tit=string(_extra.title) cmdin_m = !dustemcgwin_ncmds.pltit.pl xxpos = !dustemcgwin_ncmds.pltit.txtwdth yypos = 0.96 cgwindow,'cgtext', xxpos , yypos,' (Final Run)',color=0,/normal,charsize = 1.6,/replacecmd,cmdindex=cmdin_m, winid=winid_m endif endelse ;cgWindow_GetDefs,xsize=xsize,ysize=ysize ;cgcontrol, resize=[xsize,ysize] 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... ;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 cgwindow,'dustemcgwin_dataset', st, dataset='SED', /nodata, position=p_sed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='SED', /nodata, /norm, position=np_sed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;SED exists if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.sed.pl cgwindow,'dustemcgwin_dataset', st, dustem_sed,SED_spec,dataset='SED', /refresh ,position=p_sed, winid=winid_m,cmdindex=cmdind_m,/replacecmd,_extra=_extra cmdind_m = !dustemcgwin_ncmds.sed.nrm cgwindow,'dustemcgwin_dataset', st, dustem_sed,SED_spec,dataset='SED', /norm, /refresh,position=np_sed, winid=winid_m,cmdindex=cmdind_m,/replacecmd,_extra=_extra endif else begin cgwindow,'dustemcgwin_dataset', st, dataset='SED', position=p_sed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dustem_sed,SED_spec, dataset='SED', /refresh ,position=p_sed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.sed.pl = cmdind_m cgwindow,'dustemcgwin_dataset', st, dataset='SED',/norm, position=np_sed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 ;I think nodata can also be used too cgwindow,'dustemcgwin_dataset', st, dustem_sed,SED_spec, dataset='SED', /norm, /refresh,position=np_sed, /addcmd ,winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.sed.nrm = cmdind_m endelse endelse ;if ~!run_pol then begin;I don't think the order of the ;goto, end_m if !run_pol then begin ;then goto, end_m if ~test_qsed then begin ;Qsed is not present cgwindow,'dustemcgwin_dataset', st, dataset='QSED', /nodata, position=p_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='QSED', /nodata, /norm, position=np_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;Qsed data is present if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.qsed.pl ;refreshing the positive_only and negative_only plots if data is present in them cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec,dataset='QSED',/positive_only, nodata=0, /refresh ,position=p_qsed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec, dataset='QSED',/negative_only, nodata=0, /refresh ,position=p_qsed, winid=winid_m,/replacecmd,cmdindex=cmdind_m+2, _extra=_extra cmdind_m = !dustemcgwin_ncmds.qsed.nrm cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec, dataset='QSED', /refresh, nodata=0, /norm,position=np_qsed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra endif else begin ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendered) cgwindow,'dustemcgwin_dataset', st, dataset='QSED', position=p_qsed, /addcmd, winid=winid_m,/positive_only, _extra=_extra & cmdind_m+=1 ;Refreshing the positive_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec, dataset='QSED',/positive_only, /refresh ,position=p_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.qsed.pl = cmdind_m cgwindow,'dustemcgwin_dataset', st, dataset='QSED', position=p_qsed, /addcmd, winid=winid_m,/negative_only, _extra=_extra & cmdind_m+=1 ;Refreshing the negative_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec,dataset='QSED',/negative_only, /refresh ,position=p_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='QSED',/norm, position=np_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 ;I think nodata can also be uqsed too cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec,dataset='QSED', /norm, /refresh,position=np_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.qsed.nrm = cmdind_m endelse endelse ;stop if ~test_used then begin ;used is not present cgwindow,'dustemcgwin_dataset', st, dataset='USED', /nodata, position=p_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='USED', /nodata, /norm, position=np_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;used data is present if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.used.pl ;refreshing the positive_only and negative_only plots if data is present in them cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED',/positive_only, /refresh, nodata=0 ,position=p_used, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED',/negative_only, /refresh, nodata=0 ,position=p_used, winid=winid_m,/replacecmd,cmdindex=cmdind_m+2, _extra=_extra cmdind_m = !dustemcgwin_ncmds.used.nrm cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED', /refresh, nodata=0, /norm,position=np_used, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra endif else begin ;stop ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendred) cgwindow,'dustemcgwin_dataset', st, dataset='USED', position=p_used, /addcmd, winid=winid_m,/positive_only, _extra=_extra & cmdind_m+=1 ;Refreshing the positive_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED',/positive_only, /refresh ,position=p_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.used.pl = cmdind_m cgwindow,'dustemcgwin_dataset', st, dataset='USED', position=p_used, /addcmd, winid=winid_m,/negative_only, _extra=_extra & cmdind_m+=1 ;Refreshing the negative_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED',/negative_only, /refresh ,position=p_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='USED',/norm, position=np_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 ;I think nodata can also be uused too cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec, dataset='USED', /norm, /refresh,position=np_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.used.nrm = cmdind_m endelse endelse if ~test_psi_em then begin ; psi_em data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='psi_em', /nodata, position=p_psised, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;psi_em exists ;stop if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.psi_em.pl cgwindow,'dustemcgwin_dataset', st, dustem_psi_em,PSI_spec, dataset='psi_em', /refresh ,position=p_psised, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra endif else begin ;stop cgwindow,'dustemcgwin_dataset', st, dataset='psi_em', position=p_psised, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dustem_psi_em,PSI_spec, dataset='psi_em', /refresh ,position=p_psised, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.psi_em.pl = cmdind_m endelse endelse if ~test_polsed then begin ; polsed data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='polsed', /nodata, position=p_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dataset='polsed', /nodata, /norm, position=np_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;polsed exists if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.polsed.pl cgwindow,'dustemcgwin_dataset', st, dustem_polsed,P_spec, dataset='polsed', /refresh ,position=p_psed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra cmdind_m = !dustemcgwin_ncmds.polsed.nrm cgwindow,'dustemcgwin_dataset', st, dustem_polsed,P_spec, dataset='polsed', /norm, /refresh,position=np_psed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra endif else begin ;stop cgwindow,'dustemcgwin_dataset', st, dataset='polsed', position=p_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dustem_polsed,P_spec, dataset='polsed', /refresh ,position=p_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.polsed.pl = cmdind_m cgwindow,'dustemcgwin_dataset', st, dataset='polsed',/norm, position=np_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 ;I think nodata can also be upolsed too cgwindow,'dustemcgwin_dataset', st, dustem_polsed,P_spec, dataset='polsed', /norm, /refresh,position=np_psed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.polsed.nrm = cmdind_m endelse endelse if ~test_polfrac then begin ; polfrac data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='polfrac', /nodata, position=p_spsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 endif else begin ;polfrac exists if ~iswinsed then begin cmdind_m = !dustemcgwin_ncmds.polfrac.pl cgwindow,'dustemcgwin_dataset', st, dustem_polfrac,SP_spec,SED_spec, dataset='polfrac', /refresh ,position=p_spsed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra endif else begin ;stop cgwindow,'dustemcgwin_dataset', st, dataset='polfrac', position=p_spsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 cgwindow,'dustemcgwin_dataset', st, dustem_polfrac,SP_spec,SED_spec, dataset='polfrac', /refresh ,position=p_spsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1 !dustemcgwin_ncmds.polfrac.pl = cmdind_m endelse endelse ENDIF cgcontrol, winid_m, execute=1 ENDIF ;stop if test_x then begin ;Generating extinction window and saving window ID if iswinext then begin if !run_pol then cgwindow,wxsize=910,wysize=658,wtitle='DUSTEMWRAP v2.0 (EXTINCTION)' else cgwindow, wtitle='DUSTEMWRAP v2.0 (EXTINCTION)';, wback='grey';,wobject=winobj_m winid_x = cgquery(dimensions=dim_x,/current) ;this has to be changed because the dimensions of this array change with the data sets present (extinction/emission). !dustemcgwin_id.ext = winid_x cgcontrol, winid_x, execute=0 ;cgset, winid_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,/addcmd & cmdind_x+=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_x+=1 cgtext, xxpos, yypos,'Iteration : ',color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_x+=1 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,/addcmd & cmdind_x+=1 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,/addcmd & cmdind_x+=1 ;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,/addcmd & cmdind_x+=1 ;saving the command id to replace it !dustemcgwin_ncmds.runs.pl = cmdind_x ;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,/addcmd & cmdind_x+=1 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,/addcmd & cmdind_x+=1 widthtot+= thiswidth;varying width ;saving the command string width to replace it !dustemcgwin_ncmds.chi2.txtwdth = thiswidth xxpos = 1.0 - widthtot - 0.03 yypos = 0.97 ;Displaying fixed part cgtext, xxpos, yypos,textoidl('\chi^{2} : '),color=0,/normal,charsize = 1.0,/addcmd & cmdind_x+=1 ;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,/addcmd & cmdind_x+=1 ;saving the command id to replace it !dustemcgwin_ncmds.chi2.pl = cmdind_x ;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,/addcmd & cmdind_x+=1 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,/addcmd & cmdind_x+=1 widthtot+= thiswidth;varying width !dustemcgwin_ncmds.rchi2.txtwdth = thiswidth ;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,/addcmd & cmdind_x+=1 ;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,/addcmd & cmdind_x+=1 ;saving the command id to replace it !dustemcgwin_ncmds.rchi2.pl = cmdind_x ;Plotting of the title of the dashboard (plot(s)) ;NB: JP is right, the user should be able to set the entirety of the title. ;Make this block aware of the saving of the data (from the fits table) ;Getting the width of the title tit=string(_extra.title); cgtext, 0.3, 0.95,tit+' (RUNNING)',color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_x+=1 widthtot = thiswidth xxpos = (1 - thiswidth)/2 yypos = 0.96 cgtext, 0.3, 0.95,tit,color=0,/normal,charsize=-1,width=thiswidth,/addcmd & cmdind_x+=1 !dustemcgwin_ncmds.pltit.txtwdth = thiswidth*1.6 + (1 - widthtot)/2 cgtext, xxpos, yypos,tit,color=0,/normal,charsize = 1.6,/addcmd & cmdind_x+=1 xxpos+=thiswidth*1.6 cgtext, xxpos, yypos,' (RUNNING)',color=0,/normal,charsize = 1.6,/addcmd & cmdind_x+=1 !dustemcgwin_ncmds.pltit.pl = cmdind_x endif else begin winid_x = !dustemcgwin_id.ext cgcontrol, winid_x, execute=0 cgset, winid_x ;apparently specifying the window index isn't enough. This might be because of the use of cgcontrol ;Refreshing chi2 data: cmdind_x = !dustemcgwin_ncmds.chi2.pl xxpos = 1 - !dustemcgwin_ncmds.chi2.txtwdth - 0.03 yypos = 0.97 cgwindow,'cgtext', xxpos , yypos,strtrim(string(chi2,format=frmt2),2),color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdind_x, winid=winid_x ;& cmdind_m+=1 ;Refreshing rchi2 data: cmdind_x = !dustemcgwin_ncmds.rchi2.pl xxpos = 1 - !dustemcgwin_ncmds.rchi2.txtwdth - 0.03 yypos = 0.95 cgwindow,'cgtext', xxpos , yypos,strtrim(string(rchi2,format=frmt2),2),color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdind_x, winid=winid_x ;& cmdind_m+=1 ;Refreshing number of runs data: cmdind_x = !dustemcgwin_ncmds.runs.pl xxpos = 0.03 + !dustemcgwin_ncmds.runs.txtwdth*1.35 ;(because of charsize) yypos = 0.015 cgwindow,'cgtext', xxpos , yypos,strtrim(!dustem_iter.act,2),color=0,/normal,charsize = 1.0,/replacecmd,cmdindex=cmdind_x, winid=winid_x ;& cmdind_m+=1 ;update of the title here ? ;NB: JP is right, the user should be able to set the entirety of the title. ;Make this block aware of the saving of the data (from the fits table) if isa((*!dustem_fit).current_param_values) then begin ;replacing my count - 3 tit=string(_extra.title) cmdind_x = !dustemcgwin_ncmds.pltit.pl xxpos = !dustemcgwin_ncmds.pltit.txtwdth yypos = 0.96 cgwindow,'cgtext', xxpos , yypos,' (Final Run)',color=0,/normal,charsize = 1.6,/replacecmd,cmdindex=cmdind_x, winid=winid_x endif endelse ;cgWindow_GetDefs,xsize=xsize,ysize=ysize ;cgcontrol, resize=[xsize,ysize] if !run_pol then begin ;,wback='grey' ;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 cgwindow,'dustemcgwin_dataset', st, dataset='EXT', /nodata, position=p_ext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='EXT', /nodata, /norm, position=np_ext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;EXT exists if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.ext.pl cgwindow,'dustemcgwin_dataset', st, dustem_ext,EXT_spec,dataset='EXT', /refresh ,position=p_ext, winid=winid_x,cmdindex=cmdind_x,/replacecmd,_extra=_extra cmdind_x = !dustemcgwin_ncmds.ext.nrm cgwindow,'dustemcgwin_dataset', st, dustem_ext,EXT_spec,dataset='EXT', /norm, /refresh,position=np_ext, winid=winid_x,cmdindex=cmdind_x,/replacecmd,_extra=_extra endif else begin cgwindow,'dustemcgwin_dataset', st, dataset='EXT', position=p_ext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dustem_ext,EXT_spec, dataset='EXT', /refresh ,position=p_ext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.ext.pl = cmdind_x cgwindow,'dustemcgwin_dataset', st, dataset='EXT',/norm, position=np_ext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 ;I think nodata can also be used too cgwindow,'dustemcgwin_dataset', st, dustem_ext,EXT_spec, dataset='EXT', /norm, /refresh,position=np_ext, /addcmd ,winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.ext.nrm = cmdind_x endelse endelse if !run_pol then begin if ~test_qext then begin ;Qext is not present cgwindow,'dustemcgwin_dataset', st, dataset='QEXT', /nodata, position=p_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='QEXT', /nodata, /norm, position=np_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;Qext data is present if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.qext.pl ;refreshing the positive_only and negative_only plots if data is present in them cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec,dataset='QEXT',/positive_only, nodata=0, /refresh ,position=p_qext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec, dataset='QEXT',/negative_only, nodata=0, /refresh ,position=p_qext, winid=winid_x,/replacecmd,cmdindex=cmdind_x+2, _extra=_extra cmdind_x = !dustemcgwin_ncmds.qext.nrm cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec, dataset='QEXT', /refresh, nodata=0, /norm,position=np_qext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra endif else begin ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendered) cgwindow,'dustemcgwin_dataset', st, dataset='QEXT', position=p_qext, /addcmd, winid=winid_x,/positive_only, _extra=_extra & cmdind_x+=1 ;Refreshing the positive_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec, dataset='QEXT',/positive_only, /refresh ,position=p_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.qext.pl = cmdind_x cgwindow,'dustemcgwin_dataset', st, dataset='QEXT', position=p_qext, /addcmd, winid=winid_x,/negative_only, _extra=_extra & cmdind_x+=1 ;Refreshing the negative_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec,dataset='QEXT',/negative_only, /refresh ,position=p_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='QEXT',/norm, position=np_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 ;I think nodata can also be uqsed too cgwindow,'dustemcgwin_dataset', st, dustem_qext,QEXT_spec,dataset='QEXT', /norm, /refresh,position=np_qext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.qext.nrm = cmdind_x endelse endelse ;stop if ~test_uext then begin ;uext is not present cgwindow,'dustemcgwin_dataset', st, dataset='UEXT', /nodata, position=p_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='UEXT', /nodata, /norm, position=np_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;uext data is present if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.uext.pl ;refreshing the positive_only and negative_only plots if data is present in them cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT',/positive_only, nodata=0, /refresh ,position=p_uext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT',/negative_only, nodata=0, /refresh ,position=p_uext, winid=winid_x,/replacecmd,cmdindex=cmdind_x+2, _extra=_extra cmdind_x = !dustemcgwin_ncmds.uext.nrm cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT', /refresh, nodata=0, /norm,position=np_uext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra endif else begin ;stop ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendred) cgwindow,'dustemcgwin_dataset', st, dataset='UEXT', position=p_uext, /addcmd, winid=winid_x,/positive_only, _extra=_extra & cmdind_x+=1 ;Refreshing the positive_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT',/positive_only, /refresh ,position=p_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.uext.pl = cmdind_x cgwindow,'dustemcgwin_dataset', st, dataset='UEXT', position=p_uext, /addcmd, winid=winid_x,/negative_only, _extra=_extra & cmdind_x+=1 ;Refreshing the negative_only plot ONLY if there is data. If not do nothing cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT',/negative_only, /refresh ,position=p_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='UEXT',/norm, position=np_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 ;I think nodata can also be uused too cgwindow,'dustemcgwin_dataset', st, dustem_uext,UEXT_spec, dataset='UEXT', /norm, /refresh,position=np_uext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.uext.nrm = cmdind_x endelse endelse if ~test_psi_ext then begin ; psi_ext data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='PSI_EXT', /nodata, position=p_psiext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;psi_ext exists if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.psi_ext.pl cgwindow,'dustemcgwin_dataset', st, dustem_psi_ext,PSIEXT_spec, dataset='PSI_EXT', /refresh ,position=p_psiext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra endif else begin ;stop cgwindow,'dustemcgwin_dataset', st, dataset='psi_ext', position=p_psiext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dustem_psi_ext,PSIEXT_spec, dataset='psi_ext', /refresh ,position=p_psiext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.psi_ext.pl = cmdind_x endelse endelse if ~test_polext then begin ; polext data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='polext', /nodata, position=p_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dataset='polext', /nodata, /norm, position=np_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;polext exists if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.polext.pl cgwindow,'dustemcgwin_dataset', st, dustem_polext,POLEXT_spec, dataset='polext', /refresh ,position=p_pext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra cmdind_x = !dustemcgwin_ncmds.polext.nrm cgwindow,'dustemcgwin_dataset', st, dustem_polext,POLEXT_spec, dataset='polext', /norm, /refresh,position=np_pext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra endif else begin ;stop cgwindow,'dustemcgwin_dataset', st, dataset='polext', position=p_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dustem_polext,POLEXT_spec, dataset='polext', /refresh ,position=p_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.polext.pl = cmdind_x cgwindow,'dustemcgwin_dataset', st, dataset='polext',/norm, position=np_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 ;I think nodata can also be upolsed too cgwindow,'dustemcgwin_dataset', st, dustem_polext,POLEXT_spec, dataset='polext', /norm, /refresh,position=np_pext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.polext.nrm = cmdind_x endelse endelse if ~test_fpolext then begin ; fpolext data is not present (For completeness) # This is only valid when !run_pol=1 cgwindow,'dustemcgwin_dataset', st, dataset='fpolext', /nodata, position=p_spext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 endif else begin ;fpolext exists if ~iswinext then begin cmdind_x = !dustemcgwin_ncmds.fpolext.pl cgwindow,'dustemcgwin_dataset', st, dustem_fpolext,SPEXT_spec,EXT_spec, dataset='fpolext', /refresh ,position=p_spext, winid=winid_x,/replacecmd,cmdindex=cmdind_x, _extra=_extra endif else begin cgwindow,'dustemcgwin_dataset', st, dataset='fpolext', position=p_spext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 cgwindow,'dustemcgwin_dataset', st, dustem_fpolext,SPEXT_spec,EXT_spec, dataset='fpolext', /refresh ,position=p_spext, /addcmd, winid=winid_x, _extra=_extra & cmdind_x+=1 !dustemcgwin_ncmds.fpolext.pl = cmdind_x endelse endelse ENDIF cgcontrol, winid_x, execute=1 ENDIF IF ~iswinprms THEN BEGIN;(NEXT RUNS) winid_prms = !dustemcgwin_id.prms cgcontrol, winid_prms, execute=0 cgset, winid_prms cmdind_prms = !dustemcgwin_ncmds.prms.pl cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PARAMETERS', /refresh, winid=winid_prms,/replacecmd,cmdindex=cmdind_prms, _extra=_extra cgcontrol, winid_prms ,execute=1 ENDIF ELSE BEGIN ;(FIRST RUN) cgwindow,wxsize=350,wysize=200,wtitle='DUSTEMWRAP v2.0 (PARAMETERS)' ; keeping the same y dimension of the emission/extinction plot? winid_prms = cgquery(dimensions=dim_prms,/current) ;we don't need the dimension of the window so far/ I don't know why I needed it before cgcontrol, winid_prms, execute=0 !dustemcgwin_id.prms = winid_prms cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PARAMETERS', /addcmd , winid=winid_prms,_extra=_extra ;& cmdind_prms+=1 cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PARAMETERS', /addcmd, /refresh , winid=winid_prms, _extra=_extra & cmdind_prms+=1 !dustemcgwin_ncmds.prms.pl = cmdind_prms cgcontrol, winid_prms ,execute=1 ENDELSE ;DISPLAY OF PLUGINS HAS A SMALL ISSUE IF (tag_names(*!dustem_plugin))(0) NE 'NONE' THEN BEGIN IF ~iswinplgns THEN BEGIN;(NEXT RUNS) winid_plgns = !dustemcgwin_id.plgns cgcontrol, winid_plgns, execute=0 cgset, winid_plgns cmdind_plgns = !dustemcgwin_ncmds.plgns.pl cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PLUGINS', /refresh, winid=winid_plgns,/replacecmd,cmdindex=cmdind_plgns, _extra=_extra cgcontrol, winid_plgns ,execute=1 ENDIF ELSE BEGIN ;(FIRST RUN) cgwindow,wxsize=350,wysize=200,wtitle='DUSTEMWRAP v2.0 (PLUGINS)' ; keeping the same y dimension of the emission/extinction plot? winid_plgns = cgquery(dimensions=dim_plgns,/current) ;we don't need the dimension of the window so far/ I don't know why I needed it before cgcontrol, winid_plgns, execute=0 !dustemcgwin_id.plgns = winid_plgns cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PLUGINS', /addcmd , winid=winid_plgns,_extra=_extra ;& cmdind_prms+=1 cgwindow,'dustemcgwin_dataset', st, res, errors, p_dim,dataset='PLUGINS', /addcmd, /refresh , winid=winid_plgns, _extra=_extra & cmdind_plgns+=1 !dustemcgwin_ncmds.plgns.pl = cmdind_plgns cgcontrol, winid_plgns ,execute=1 ENDELSE ENDIF END