allowing no_obj plotting of extiction data

Ilyes Choubani
1 parent 4da8e20e
Showing 2 changed files with 217 additions and 322 deletions Show diff stats
src/idl/dustem_plugin_modify_dust_polx.pro
src/idl/dustemwrap_plot_noobj.pro
@@ -52,7 +52,10 @@ IF keyword_set(paramtag) THEN BEGIN
 ENDIF 
  
 ;Retrieving this system variable (dustem output) 
-IF isa(!dustem_current) then st=*!dustem_current else out=0 & goto, the_end
+if ~isa(!dustem_current) then begin
+    out =0 
+    goto, the_end 
+ENDIF else st=(*!dustem_current)
  
 ;below are the default values for the plugin parameters
 smallp_fact=1.     ;This is the default multiplicative factore to the dust polarization
@@ -672,8 +672,7 @@ if test_m then begin
  
  
     ENDIF
-      
-    cgcontrol, winid_m, execute=1
+    
  
 ENDIF 
  
@@ -682,171 +681,159 @@ 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
  
-    if iswinext then begin
-        
-        if !run_pol then cgwindow,wxsize=main_wxsize,wysize=main_wysize,wtitle='DUSTEMWRAP '+!dustem_version.version+' (EXTINCTION)' else cgwindow,wxsize=main_wxsize,wysize=main_wysize,wtitle='DUSTEMWRAP '+!dustem_version.version+' (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 
+    ;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
  
-        
-        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
-        
-          
-        if !dustem_end EQ 1 then begin
-            ;replacing my count - 3
-            tit=!dustem_plot_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
+    ;Plotting of the number of iterations completed by the fitter
  
-    endelse
+    ;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 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
+    widthtot = thiswidth
+    xxpos = (1 - thiswidth)/2
+    yypos = 0.96
+    
+    cgtext, 0.3, 0.95,tit,color=0,/normal,charsize=-1,width=thiswidth
+      
+    cgtext, xxpos, yypos,tit,color=0,/normal,charsize = 1.6
+    xxpos+=thiswidth*1.6
+    cgtext, xxpos, yypos,' (RUNNING)',color=0,/normal,charsize = 1.6
+    
+    
+;     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
+;     
+;       
+;     if !dustem_end EQ 1 then begin
+;         ;replacing my count - 3
+;         tit=!dustem_plot_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]
@@ -896,32 +883,18 @@ if test_x then begin
     ;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
+        dustemcgwin_dataset, st, dataset='EXT', /nodata, position=p_ext, _extra=_extra 
+        dustemcgwin_dataset, st, dataset='EXT', /nodata, /norm, position=np_ext, _extra=_extra
  
     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       
+        dustemcgwin_dataset, st, dataset='EXT', position=p_ext, _extra=_extra
+        dustemcgwin_dataset, st, dustem_ext,EXT_spec, dataset='EXT', /refresh ,position=p_ext, _extra=_extra 
+        
+        dustemcgwin_dataset, st, dataset='EXT',/norm, position=np_ext, _extra=_extra ;I think nodata can also be used too
+        
+        dustemcgwin_dataset, st, dustem_ext,EXT_spec, dataset='EXT', /norm, /refresh,position=np_ext, _extra=_extra 
+               
     endelse
  
  
@@ -929,195 +902,114 @@ if test_x 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
+            dustemcgwin_dataset, st, dataset='QEXT', /nodata, position=p_qext, _extra=_extra 
+            dustemcgwin_dataset, st, dataset='QEXT', /nodata, /norm, position=np_qext, _extra=_extra 
  
         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
+                    
+            ;Plotting of postive values if they exist (if hidden data is present it is plotted if not an empty plot is rendered)
+            dustemcgwin_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
  
-                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       
+            dustemcgwin_dataset, st, dustem_qext,QEXT_spec, dataset='QEXT',/positive_only, /refresh ,position=p_qext, _extra=_extra 
+            
+            dustemcgwin_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
+            dustemcgwin_dataset, st, dustem_qext,QEXT_spec,dataset='QEXT',/negative_only, /refresh ,position=p_qext, _extra=_extra
+        
+            dustemcgwin_dataset, st, dataset='QEXT',/norm, position=np_qext, _extra=_extra ;I think nodata can also be uqsed too
+            dustemcgwin_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
  
-            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
+            dustemcgwin_dataset, st, dataset='UEXT', /nodata, position=p_uext, _extra=_extra
+            dustemcgwin_dataset, st, dataset='UEXT', /nodata, /norm, position=np_uext, _extra=_extra 
  
         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
+                dustemcgwin_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
-                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
+                dustemcgwin_dataset, st, dustem_uext,UEXT_spec, dataset='UEXT',/positive_only, /refresh ,position=p_uext, _extra=_extra 
  
-                cgwindow,'dustemcgwin_dataset', st, dataset='UEXT', position=p_uext, /addcmd, winid=winid_x,/negative_only, _extra=_extra & cmdind_x+=1
+              
+                dustemcgwin_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
-                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
-                
+                dustemcgwin_dataset, st, dustem_uext,UEXT_spec, dataset='UEXT',/negative_only, /refresh ,position=p_uext, _extra=_extra
  
-                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       
+                dustemcgwin_dataset, st, dataset='UEXT',/norm, position=np_uext, _extra=_extra ;I think nodata can also be uused too
+                dustemcgwin_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
  
-          cgwindow,'dustemcgwin_dataset', st, dataset='PSI_EXT', /nodata, position=p_psiext, /addcmd, winid=winid_x, _extra=_extra  & cmdind_x+=1
+          dustemcgwin_dataset, st, dataset='PSI_EXT', /nodata, position=p_psiext, _extra=_extra  
  
       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
+              dustemcgwin_dataset, st, dataset='psi_ext', position=p_psiext, _extra=_extra 
+              dustemcgwin_dataset, st, dustem_psi_ext,PSIEXT_spec, dataset='psi_ext', /refresh ,position=p_psiext, _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
+          dustemcgwin_dataset, st, dataset='polext', /nodata, position=p_pext, _extra=_extra 
+          dustemcgwin_dataset, st, dataset='polext', /nodata, /norm, position=np_pext, _extra=_extra 
  
       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
+        dustemcgwin_dataset, st, dataset='polext', position=p_pext, _extra=_extra 
+        dustemcgwin_dataset, st, dustem_polext,POLEXT_spec, dataset='polext', /refresh ,position=p_pext, _extra=_extra 
+         
+        dustemcgwin_dataset, st, dataset='polext',/norm, position=np_pext, _extra=_extra  ;I think nodata can also be upolsed too
+        dustemcgwin_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
  
-          cgwindow,'dustemcgwin_dataset', st, dataset='fpolext', /nodata, position=p_spext, /addcmd, winid=winid_x, _extra=_extra  & cmdind_x+=1
+          dustemcgwin_dataset, st, dataset='fpolext', /nodata, position=p_spext, _extra=_extra  
+          
  
       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
+        
+        dustemcgwin_dataset, st, dataset='fpolext', position=p_spext, _extra=_extra 
+        dustemcgwin_dataset, st, dustem_fpolext,SPEXT_spec,EXT_spec, dataset='fpolext', /refresh ,position=p_spext, _extra=_extra
+
  
-          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)
-    window,1,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? 
-;    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
-    dustemcgwin_dataset, st, res, errors, p_dim,dataset='PARAMETERS',_extra=_extra ;& cmdind_prms+=1
-    dustemcgwin_dataset, st, res, errors, p_dim,dataset='PARAMETERS',/refresh, _extra=_extra 
-    ;!dustemcgwin_ncmds.prms.pl = cmdind_prms
-    cgcontrol, winid_prms ,execute=1        
-; ENDELSE 
-
+;PLOTTING OF THE PARAMETERS
+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? 
+dustemcgwin_dataset, st, res, errors, p_dim,dataset='PARAMETERS',_extra=_extra ;& cmdind_prms+=1
+dustemcgwin_dataset, st, res, errors, p_dim,dataset='PARAMETERS',/refresh, _extra=_extra 
  
+;PLOTTING OF THE PLUGINS
 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)
-        window,2,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? 
-        ;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
-        dustemcgwin_dataset, st, res, errors, p_dim,dataset='PLUGINS',_extra=_extra ;& cmdind_prms+=1
-        dustemcgwin_dataset, st, res, errors, p_dim,dataset='PLUGINS', /refresh , winid=winid_plgns, _extra=_extra 
-        ;!dustemcgwin_ncmds.plgns.pl = cmdind_plgns
-        cgcontrol, winid_plgns ,execute=1        
-;    ENDELSE
+    
+    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? 
+    dustemcgwin_dataset, st, res, errors, p_dim,dataset='PLUGINS',_extra=_extra ;& cmdind_prms+=1
+    dustemcgwin_dataset, st, res, errors, p_dim,dataset='PLUGINS', /refresh , winid=winid_plgns, _extra=_extra 
  
 ENDIF