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,_extra=_extra 


;This is 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)  ; 0/0 division case?
  st=dustem_run(p_dim)
  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) 
      
  endif  
   
ENDIF



;if ~windowavailable(cgquery()) then begin

;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_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_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

;if iswinsed then begin ;I believe the iswin test should be below the case condition because I am not copying all the first 'HUGE LOOP' but rather replacing the positional arrays.            
if test_m then begin
    ;Generating emission window and saving window ID
    
    if iswinsed then begin
        
        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.0,/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.0,/addcmd & cmdind_m+=1
        
        ;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.015
        
        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.015
        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))
        
        ;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 
    
        
        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.015   
        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.015   
        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
         
         
        ;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,SED_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,SED_spec, dataset='QSED',/negative_only, nodata=0, /refresh ,position=p_qsed, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra
                
                cmdind_m = !dustemcgwin_ncmds.qsed.nrm
                cgwindow,'dustemcgwin_dataset', st, dustem_qsed,Q_spec,SED_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,SED_spec, dataset='QSED',/positive_only, /refresh ,position=p_qsed, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1
                
                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,SED_spec,dataset='QSED',/negative_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',/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,SED_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,SED_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,SED_spec, dataset='USED',/negative_only, /refresh, nodata=0 ,position=p_used, winid=winid_m,/replacecmd,cmdindex=cmdind_m, _extra=_extra
                
                cmdind_m = !dustemcgwin_ncmds.used.nrm
                cgwindow,'dustemcgwin_dataset', st, dustem_used,U_spec,SED_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,SED_spec, dataset='USED',/positive_only, /refresh ,position=p_used, /addcmd, winid=winid_m, _extra=_extra & cmdind_m+=1
                
                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,SED_spec, dataset='USED',/negative_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',/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,SED_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

         if  ~iswinsed then begin
             cmdind_m = !dustemcgwin_ncmds.psi_em.pl
             cgwindow,'dustemcgwin_dataset', st, dustem_psi_em,PSI_spec,SED_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,SED_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
        
    end_m:
    
    ;cgcontrol, update=1
;     varrr=cgquery(dimensions=dim_m) ;be careful because the dim_m array can be multi-dimensional because of the presence of two windows.
;     cgcontrol, update=1
;     cgcontrol, resize=[dim_m(0),dim_m(1)]
;     
    ;cgcontrol, resize=1
;     cgcontrol,update=1
;     cgcontrol,update=1
        
;     cgcontrol, resize=[!D.X_SIZE,!D.Y_SIZE]
    
    ;cgcontrol, update=1
    ;cgcontrol, update=0
    
    
 ;   cgcontrol, update=1 
    
    
    
    
    ;varr=cgquery(dimensions=dim_m,objectref=winobj_m)
   ;stop
    cgcontrol, winid_m, execute=1
    ;cgcontrol, winid_m, execute=0      
         
ENDIF 


;Plotting of the parameter and plugin values / This will have to be placed after the extinction block which hasn't been 
;written yet. 

;This block will be different because it does not depend on wether dustemwrap is used in polarization or not


;stop
;Plotting the parameter window 
;cgcontrol, execute=0
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 

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 


END