;+ ; Project : SOHO - CDS ; ; Name : CHIANTI_NE ; ; Purpose : Calculate and plot CHIANTI density sensitive line ratios. ; ; Explanation : CHIANTI_NE (density ratios) ; calculates and plots density sensitive line ratios based on ; the CHIANTI atomic database of Dere et. al. ; ; Use : IDL> chianti_ne ; ; Inputs : None ; ; Opt. Inputs : None ; ; Outputs : None ; ; Opt. Outputs: None. ; ; Keywords : None ; ; Calls : CALC_DMM_DR, plot_dmm_dr_fig , make_ion_list ; ; Common : dmm_dr_com dmm_lines(with plot_dmm_dr_fig), ; ; Restrictions: None ; ; Side effects: None ; ; Category : Spectral ; ; Prev. Hist. : Started life as 'density_ratios' by Ken Dere ; ; Written : C D Pike, RAL, 13-Jan-96 ; ; Modified : Added selection of els/ions from master file. CDP, 21-Jan-96 ; Added hardcopy of line list, refs etc. CDP, 22-Jan-96 ; Include multiple lines in ratio. CDP, 27-Jan-96 ; General upgrade and added temp/unit selections. CDP, 6-Jun-97 ; Added intensity ratio selection. CDP, 17-Jul-97 ; Added wavelength ranges. CDP, 18-Jul-97 ; Fixed typos introduced in version 7. CDP, 22-Jul-97 ; Float the user-supplied Log Temperature. CDP, 1-Aug-97 ; Update for IDL v5.2. CDP, 20-Apr-99 ; v. 11 Update list of elements. CDP, 18-Jun-99 ; ; V.12. Added ratio plots and hardcopies in linear scale, added various ; checks and minor things. Added a few tags to the ; ratio output structure (temperature, units, comment). ; Removed optional output structure. ; Updated to be CHIANTI v.3-compatible. ; Giulio Del Zanna, DAMTP, 7-Oct-2000 ; Version 13, 21-Dec-2000, William Thompson, GSFC ; Modified for better cross-platform graphics capability ; ; V. 14, 1-May-02, GDZ, commented out a few refs. ; V. 15, 21-May-2002, GDZ, changed the way to deal with fonts. ; V.16, 17-Sep-2002 (GDZ) ; added !p.multi = 0 upon exit and added X-label. ; ; ; Version : Version 16, 17-Sep-2002 ;- pro chianti_ne_event,ev ; ; event processor for density ratio program (chianti_ne) ; ; ; pick up common definitions ; common dmm_dr_com, dmm_dr_base, image_area,$ mess_win, plot_title, $ xuvtop, wavmin_id, wavmax_id, densmin_id, densmax_id,$ wslider_min, wslider_max, dslider_min, dslider_max,$ el_info, ion_info, temp_info, inten_info, unit_info,$ ion_base, els, ions, phot_mod, wind_keep common dmm_lines, list_wvl, list_int, list_descr1, list_descr2, species,$ density, nden, ratio, description, nlist, savetext,$ temper, intens common wavmm, wminw1, wminw2, wminw3, wminw4, $ wmaxw1, wmaxw2, wmaxw3, wmaxw4, $ minw1, minw2, minw3, minw4, $ maxw1, maxw2, maxw3, maxw4 common dmm_refs, ioneq_ref, wgfaref, eref, upsref ; ; list of elements ; ;element=['H','He','Li','Be','B ','C ','N ','O ','F ','Ne','Na',$ ; 'Mg','Al','Si','P ','S ','Cl','Ar','K ','Ca','Sc','Ti',$ ; 'V ','Cr','Mn','Fe','Co','Ni','Cu','Zn'] ; ; get the values, user values and type of widget involved ; widget_control,ev.id,get_uvalue=uvalue name = strmid(tag_names(ev,/str),7,1000) if n_elements(uvalue) eq 0 then uvalue = '' if n_elements(value) eq 0 then value = '' ;print,'name: ',name,' value: ',value,' uvalue: ',uvalue ;help,ev,/str ;print,'**',name,'**' ; ; act according to type of widget ; case name of ; ; Button will be the type of most of the user input from various wwidgets ; 'BUTTON': begin case strmid(uvalue,0,4) of ; ; Exit from the system (read text widgets first) ; 'EXIT': begin widget_control,/destroy,ev.top DEVICE, WINDOW_STATE=STATE IF N_ELEMENTS(STATE) GT 32 THEN FOR I_WIN=32,N_ELEMENTS(STATE)-1 DO $ IF STATE(I_WIN) THEN WDELETE, I_WIN !p.multi = 0 end ; ; HELP will put the help file to the screen (not implemented yet..) ; 'HELP': begin widg_help,'dmm_ne.hlp' end ; ; Change UNITS ; 'UNIT': begin if phot_mod then begin widget_control, unit_info, set_val='Units for the ratio plot: ERGS' phot_mod = 0 endif else begin widget_control, unit_info, set_val='Units for the ratio plot: PHOTONS' phot_mod = 1 endelse end ; ; Set the calculation going ; 'CALC': begin widget_control,wminw1, get_v=minw1 minw1 = float(minw1(0)) widget_control,wminw2, get_v=minw2 minw2 = float(minw2(0)) widget_control,wminw3, get_v=minw3 minw3 = float(minw3(0)) widget_control,wminw4, get_v=minw4 minw4 = float(minw4(0)) widget_control,wmaxw1, get_v=maxw1 maxw1 = float(maxw1(0)) widget_control,wmaxw2, get_v=maxw2 maxw2 = float(maxw2(0)) widget_control,wmaxw3, get_v=maxw3 maxw3 = float(maxw3(0)) widget_control,wmaxw4, get_v=maxw4 maxw4 = float(maxw4(0)) widget_control, el_info,get_val=el el = strtrim(el(0),2) widget_control, ion_info,get_val=ion ion = strtrim(ion(0),2) widget_control, temp_info,get_val=temper temper = float(strtrim(temper(0),2)) widget_control, inten_info,get_val=intens intens = float(strtrim(intens(0),2)) widget_control, image_area, get_v=wind wset, wind !p.multi=[0,2,1] widget_control, dmm_dr_base,/hour calc_dmm_dr,el,ion,dslider_min,$ dslider_max,mess_win,$ temperature=temper !p.multi=[0,2,1] if nlist gt 0 then plot_dmm_dr_fig DEVICE, WINDOW_STATE=STATE wind_keep = where(state eq 1) wind_keep = wind_keep(0) end ; ; Delete plot windows to tidy up ; 'DELP': begin DEVICE, WINDOW_STATE=STATE if wind_keep ge 0 then begin state(wind_keep) = 0 IF N_ELEMENTS(STATE) GT 32 THEN FOR I_WIN=32,N_ELEMENTS(STATE)-1 DO $ IF STATE(I_WIN) THEN WDELETE, I_WIN endif end ; ; plot ratio after user supplies lines involved ; 'PRAT': BEGIN CASE strmid(uvalue,5,3) OF 'LIN': plot_log = 0 'LOG': plot_log = 1 END xinput,text,['Enter lines involved in ratio',$ 'eg. enter 2,3 4 / 1+5',$ 'to calculate ratio from sum of lines',$ '2, 3 and 4 all divided by lines 1+5'],$ group=dmm_dr_base,xoff=400,yoff=50 stext = text if strpos(text,'/') gt 0 then begin text = repchar(text,',',' ') text = repchar(text,'+',' ') text = repchar(text,'.',' ') text = strtrim(text,2) text = str2arr(text,del='/') text(0) = strtrim(text(0),2) text(1) = strtrim(text(1),2) numer = fix(str2arr(text(0),del=' ')) denom = fix(str2arr(text(1),del=' ')) numval = 0.0 denval = 0.0 nlines = indgen(n_elements(list_wvl)) for kk=0,n_elements(numer)-1 do begin hit = where(nlines eq numer(kk),count) if count eq 0 then begin bell widget_control, mess_win, /append,$ set_val='Invalid line number in numerator --> ' + stext return endif if phot_mod then begin numval = numval + list_int(numer(kk),*) endif else begin numval = numval + list_int(numer(kk),*)/$ list_wvl(numer(kk)) endelse endfor for kk=0,n_elements(denom)-1 do begin hit = where(nlines eq denom(kk),count) if count eq 0 then begin bell widget_control, mess_win, /append,$ set_val='Invalid line number in denominator --> ' + stext return endif if phot_mod then begin denval = denval + list_int(denom(kk),*) endif else begin denval = denval + list_int(denom(kk),*)/$ list_wvl(denom(kk)) endelse endfor ratio = numval/denval ratio=reform(ratio) if n_elements(numer) gt 1 then begin numdesc = fmt_vect(list_wvl(numer),$ form='(f10.4)',delim='+') endif else numdesc = string(list_wvl(numer(0)),'(f10.4)') numdesc = strcompress(numdesc,/rem) if n_elements(denom) gt 1 then begin dendesc = fmt_vect(list_wvl(denom),$ form='(f10.4)',delim='+') endif else dendesc = string(list_wvl(denom(0)),'(f10.4)') dendesc = strcompress(dendesc,/rem) description=species+' '+numdesc+'/'+dendesc window,/free wset,!d.window !p.multi=0 if phot_mod then ytit = 'Intensity (photons) Ratio' else $ ytit = 'Intensity (ergs) Ratio' IF plot_log THEN $ plot_oo,density,ratio,xr=[density(0),density(nden-1)],$ title=description,xtitle='Log!d10!n (Electron density [cm!u-3!n]) ', $ ytitle=ytit,chars=1.3 ELSE $ plot,density,ratio,xr=[density(0),density(nden-1)],$ title=description,xtitle='Electron density [cm!u-3!n] ',$ ytitle=ytit,chars=1.3 endif else begin bell widget_control, mess_win, /append,$ set_val='Invalid ratio specification -->' + stext endelse end ; ; SAVE the calculated ratio ; 'DUMP': BEGIN ;check that we have defined the ratio: IF n_elements(description) EQ 0 THEN BEGIN bell widget_control, mess_win, /append,$ set_val='Still no ratio to SAVE!' ENDIF ELSE BEGIN xinput,text,'Enter file name to store ratio data (a .CH_NE will be appended)',$ group=dmm_dr_base,/modal,xoff=400,yoff=50 if text ne '' then BEGIN if phot_mod then ytit = 'Intensity (photons) Ratio' else $ ytit = 'Intensity (ergs) Ratio' ne_ratio = {temperature:temper, $ units:ytit, $ comment:'Calculated with CHIANTI at constant temperature', $ desc:description,density:density,ratio:ratio} save,ne_ratio,file=text+'.CH_NE',/xdr bell widget_control, mess_win, /append,$ set_val='Ratio saved in file '+text+'.CH_NE' ENDIF ENDELSE end ; ; HARDCOPY requested ; 'HARD': begin case strmid(uvalue,5,3) of 'RAT': BEGIN ;check that we have a ratio to plot: IF n_elements(description) EQ 0 THEN BEGIN bell widget_control, mess_win, /append,$ set_val='Still no ratio defined to plot!' ENDIF ELSE BEGIN IF strmid(uvalue,8,1) EQ '1' THEN plot_log = 0 ELSE $ IF strmid(uvalue,8,1) EQ '2' THEN plot_log = 1 ps !p.multi=0 if phot_mod then ytit = 'Intensity (photons) Ratio' else $ ytit = 'Intensity (ergs) Ratio' IF plot_log THEN $ plot_oo,density,ratio,$ xr=[density(0),density(nden-1)],$ title='CHIANTI: '+description,$ xtit='Log!d10!n (Electron density [cm!u-3!n])', $ ytitle=ytit,chars=1.3 ELSE $ plot, density,ratio,$ xr=[density(0),density(nden-1)],$ title='CHIANTI: '+description,$ xtit='Electron density [cm!u-3!n]', $ ytitle=ytit,chars=1.3 get_utc,utc utc = anytim2cal(utc) xyouts, 0.05,0.0,'Printed: '+utc,/norm,$ chars=0.7 prin = getenv('PRINTER') IF prin EQ '' THEN $ xsel_printer,prin,group=dmm_dr_base psplot,queu=prin bell widget_control, mess_win, /append,$ set_val='Plot sent to printer '+prin END END 'INT': begin if nlist gt 0 then begin ps plot_dmm_dr_fig prin = getenv('PRINTER') IF prin EQ '' THEN $ xsel_printer,prin,group=dmm_dr_base,$ instruct='Select printer for hardcopy' psplot,queu=prin bell widget_control, mess_win, /append,$ set_val='Plot sent to printer '+prin endif else BEGIN bell widget_control, mess_win, /append,$ set_val='Still nothing to plot!' endelse end 'LIN': begin get_utc,utc utc = anytim2cal(utc) text = 'CHIANTI '+species+' Printed: '+$ utc savetext = [text,$ 'Calculation at constant T='+string(temper,format='(e10.3)')+' (K)', $ ' ',$ '# Wavelength Intensity Transition ', $ ' (A) ', $ savetext] ; savetext = [savetext,' ','Ion. Eq. refs:',$ ; ioneq_ref] ; savetext = [savetext,' ','Wgfa refs:',$ ; wgfaref] ; savetext = [savetext,' ','Elvl refs:',$ ; eref] ; savetext = [savetext,' ','Upsilon refs:',$ ; upsref] print_str,savetext,/hard,/quiet end else: endcase end else: print,'Button not recognised' endcase end ; ; wavelength sliders ; 'SLIDER' : begin case strmid(uvalue,0,4) of 'WAVE': begin widget_control,ev.id,get_uvalue=uvalue,get_value=value if strmid(uvalue,4,3) eq 'MIN' then begin if value(0) ge wslider_max then value(0) = value(0)-2 wslider_min = value(0) if value(0) gt wslider_max then begin widget_control, wavmax_id,set_val=value(0)+2 wslider_max = value(0)+2 endif endif else begin if value(0) le wslider_min then value(0) = value(0)+2 wslider_max = value(0) if value(0) lt wslider_min then begin widget_control, wavmin_id,set_val=value(0)-2 wslider_min = value(0)-2 endif endelse widget_control,wavmin_id,set_val=string(wslider_min) widget_control,wavmax_id,set_val=string(wslider_max) end 'DENS': begin widget_control,ev.id,get_uvalue=uvalue,get_value=value if strmid(uvalue,4,3) eq 'MIN' then begin if value(0) ge dslider_max then value(0) = value(0)-1 dslider_min = value(0) if value(0) ge dslider_max then begin widget_control, densmax_id,set_val=value(0)+1 dslider_max = value(0)+1 endif endif else begin if value(0) le dslider_min then value(0) = value(0)+1 dslider_max = value(0) if value(0) le dslider_min then begin widget_control, densmin_id,set_val=value(0)-1 dslider_min = value(0)-1 endif endelse widget_control,densmin_id,set_val=string(dslider_min) widget_control,densmax_id,set_val=string(dslider_max) end else: endcase end ; ; must have been from cw_pdmenu ; '': begin case uvalue of 'ELEMENT': begin widget_control, el_info ,set_val=ev.value nit = where(strupcase(els) eq strupcase(ev.value)) for i=0,n_elements(els) do begin widget_control, ion_base(i), map=0 endfor widget_control, ion_base(nit(0)+1), map=1 widget_control, ion_info ,set_val='' end 'ION' : widget_control, ion_info ,set_val=ev.value else: endcase end 'TEXT_CH': begin end else: print,'name not recognised' endcase end ; ; -------------------------------- ; The main routine ; -------------------------------- ; pro chianti_ne, modal=modal, group_leader=group_leader,$ font=font ; ; pick up common definitions ; common dmm_dr_com, dmm_dr_base, image_area,$ mess_win, plot_title, $ xuvtop, wavmin_id, wavmax_id, densmin_id, densmax_id,$ wslider_min, wslider_max, dslider_min, dslider_max,$ el_info, ion_info, temp_info, inten_info, unit_info,$ ion_base, els, ions, phot_mod, wind_keep common dmm_lines, list_wvl, list_int, list_descr1, list_descr2, species,$ density, nden, ratio, description, nlist, savetext, $ temper, intens common wavmm, wminw1, wminw2, wminw3, wminw4, $ wmaxw1, wmaxw2, wmaxw3, wmaxw4, $ minw1, minw2, minw3, minw4, $ maxw1, maxw2, maxw3, maxw4 defsysv,'!xuvtop', EXISTS = EXISTS IF NOT EXISTS THEN $ message, 'system variable !xuvtop must be set ' xuvtop = !xuvtop nlist = 0 wind_keep = -1 ; ; initial checks ; if (!d.flags and 65536) eq 0 then message,'Widgets are NOT available.' ; ; check not already in use ; if xregistered('chianti_ne') then begin print,'CDS chianti_ne program is already registered.' return endif ; ; set only possible graphics device and load LUT ; set_x loadct,0 ;clear all ratio arrays delvarx, list_wvl, list_int, list_descr1, list_descr2, species,$ density, nden, ratio, description, nlist ; ; the base widget is all-screen and column-orientated ; if not keyword_set(font) then begin ; font = '-adobe-helvetica-bold-r-normal--0-120-75-75-p-0-iso8859-1' font = get_dfont('-adobe-helvetica-bold-r-normal*120*75*') IF font(0) NE '' THEN font = font(0) ELSE font = 'fixed' endif device, get_screen_size = sz sz = 0.9*sz if (sz(0) ge 1280) and (sz(1) ge 1024) then sz(*) = 0 sz = sz < [1280,1024] dmm_dr_base = widget_base(title='CHIANTI Density Dependent Ratios',$ /column,/frame,x_scroll=sz(0),y_scroll=sz(1)) ; ; the TOP widget is where all the user/numerical/text input takes place ; top = widget_base(dmm_dr_base,/row,space=10) ; ; the left hand side of the top contains... ; top_left = widget_base(top,/column,space=5) wltit = widget_label(top_left,$ value='Select up to 4 wavelength ranges.',font=font) ss1 = widget_base(top_left,/row) wlss = widget_label(ss1,value='Minimum Wavelength.',font=font) sss1 = widget_base(ss1,/row) wminw1 = widget_text(sss1,/edit,xsize=5,ysize=1,uvalue='MINW1',$ value='1', font=font) wminw2 = widget_text(sss1,/edit,xsize=5,ysize=1,uvalue='MINW2',$ value='0', font=font) wminw3 = widget_text(sss1,/edit,xsize=5,ysize=1,uvalue='MINW3',$ value='0', font=font) wminw4 = widget_text(sss1,/edit,xsize=5,ysize=1,uvalue='MINW4',$ value='0', font=font) ss2 = widget_base(top_left,/row) wlss = widget_label(ss2,value='Maximum Wavelength.',font=font) sss2 = widget_base(ss2,/row) wmaxw1 = widget_text(sss2,/edit,xsize=5,ysize=1,uvalue='MAXW1',$ value='1700', font=font) wmaxw2 = widget_text(sss2,/edit,xsize=5,ysize=1,uvalue='MAXW2',$ value='0', font=font) wmaxw3 = widget_text(sss2,/edit,xsize=5,ysize=1,uvalue='MAXW3',$ value='0', font=font) wmaxw4 = widget_text(sss2,/edit,xsize=5,ysize=1,uvalue='MAXW4',$ value='0', font=font) ; ; density range sliders ; densmin_id = widget_slider(top_left,title='Log minimum density',$ uvalue='DENSMIN',xsize=500,min=6,max=14) dslider_min = 6 densmax_id = widget_slider(top_left,title='Log maximum density',$ uvalue='DENSMAX',xsize=500,min=6,max=14) widget_control, densmax_id,set_val='14' dslider_max = 14 ; ; the top middle widget contains... ; top_middle= widget_base(top,/column,space=10) ; ; selection of element ; cs3 = widget_base(top_middle,/row) junk = {CW_PDMENU_S,FLAGS:0,NAME:''} ;element=['H','He','Li','Be','B ','C ','N ','O ','F ','Ne','Na',$ ; 'Mg','Al','Si','P ','S ','Cl','Ar','K ','Ca','Sc','Ti',$ ; 'V ','Cr','Mn','Fe','Co','Ni','Cu','Zn'] ; ; large font ; bfont = get_dfont('-adobe-helvetica-bold-r-normal*180*') ; bfont = '-adobe-helvetica-bold-r-normal--0-120-75-75-p-0-iso8859-1' IF bfont(0) NE '' THEN bfont = bfont(0) ELSE bfont = 'fixed' make_ion_list,els,ions desc = [{CW_PDMENU_S,FLAGS:1,NAME:'Select element'}] for i=0,n_elements(els)-1 do begin desc = [desc,{CW_PDMENU_S,FLAGS:0,NAME:els(i)}] endfor desc(n_elements(desc)-1).flags=2 menu = cw_pdmenu(cs3,desc,/return_name,uvalue='ELEMENT',font=font) cs4 = widget_base(cs3,/column) el_info = cw_field(cs4,title='Current element:',value=' ',$ /row,xsize=4,fieldfont=bfont) ; ; selection of ion ; junk = {CW_PDMENU_S,FLAGS:0,NAME:''} ionstage=['I','II','III','IV','V','VI','VII','VIII','IX','X','XI','XII',$ 'XIII','XIV','XV','XVI','XVII','XVIII','XIX','XX','XXI','XXII',$ 'XXIII','XXIV','XXV','XXVI','XXVII','XXVIII'] cs51 = widget_base(top_middle,/row) cs5 = widget_base(cs51) nel = n_elements(els) ion_base = lonarr(nel+1) ion_base(0) = widget_base(cs5,map=0) for k=1,nel do begin desc = [{CW_PDMENU_S,FLAGS:1,NAME:'Select ionization'}] ionstage = ions(k-1,*) ionstage = ionstage(where(ionstage ne '')) for i=0,n_elements(ionstage)-1 do begin desc = [desc,{CW_PDMENU_S,FLAGS:0,NAME:ionstage(i)}] endfor desc(n_elements(desc)-1).flags=2 ion_base(k) = cw_pdmenu(cs5,desc,/return_name,uvalue='ION',font=font) widget_control, ion_base(k), map=0 endfor cs6 = widget_base(cs51,/column) ion_info = cw_field(cs6,title='Current ion:',value=' ',$ /row,xsize=5,fieldfont=bfont) unit_info = widget_button(top_middle,value='Units for the ratio plot: ERGS',uvalue='UNIT',$ frame=2) ;define the default: phot_mod = 0 widget_control,unit_info,set_val='Units for the ratio plot: ERGS' ; ; temperature selection ; eb2 = widget_base(top_middle,/column) temp_info = cw_field(eb2,title='Temperature (K):',value='0',$ /row,xsize=8,/return,uvalue='USERT') ;define default: ;temper = 1.e6 ; ; minimun line intensity selection ; eb3 = widget_base(top_middle,/column) inten_info = cw_field(eb3,title='Minimum intensity ratio (for the calculation):',$ value='0.0001',$ /row,xsize=10,/return,uvalue='USERR') widget_control, ion_base(1), map=1 ; ; MIDDLE "command" line has the options... ; middle = widget_base(dmm_dr_base,/row,space=10,/frame) ; ; ...Quit without any further ado ; m1 = widget_base(middle,/column) dummy = widget_button(m1,value='Quit',uvalue='EXIT',/no_rel,font=font) ; ; ...calculate new data with current instrument settings ; m2 = widget_base(middle,/column) dummy = widget_button(m2,value='CALCULATE line intensities',uvalue='CALC',$ /no_rel,font=font) ; ; replot graph with revised wavelength range ; m9 = widget_base(middle,/column) XPDmenu, ['"PLOT ratio" {','"in linear scale" PRAT_LIN',$ '" in log scale" PRAT_LOG','}'], m9, font=font ; ; request HARDCOPY ; m4 = widget_base(middle,/column) XPDmenu, ['"HARDCOPY" {','"Of ratio plot (linear scale)" HARD_RAT1',$ '"Of ratio plot (log scale)" HARD_RAT2',$ '"Line details (+ refs)" HARD_LIN',$ '"Of intensity plots" HARD_INT','}'], m4, font=font ; ; dump details ; m6 = widget_base(middle,/column) XPDmenu, ['"SAVE" {','"Ratio values" DUMP_DR','}'],m6, font=font ;dummy = widget_button(m9,value='PLOT ratio',uvalue='PRAT',$ ; /no_rel,font=font) ; ; delete plot windows to tidy up ; m10 = widget_base(middle,/column) dummy = widget_button(m10,value='Delete plot windows',uvalue='DELP',$ /no_rel,font=font) ; ; BOTTOM display area ; bottom = widget_base(dmm_dr_base,/frame,/column) ; ; ...plotting area ; image_area = widget_draw(bottom,xsize=900,ysize=250,/button_events) have_vert = 0 widget_control,image_area,sensitive=0 ; ; messages ; mess_win = widget_text(dmm_dr_base,ysize=8,xsize=80,/scroll) ; ; make the whole thing happen ; widget_control,dmm_dr_base,/realize ; ; find out window id for plotting ; widget_control,get_value=window,image_area wset,window ; ; make the boss do some work ; xmanager,'chianti_ne',dmm_dr_base, modal=modal, group_leader=group_leader end