;OBSOLETE FILE - WAS replaced by thm_gen_overplot ;+ ; Purpose: To make mission overview plots of all instruments ; ; Inputs: PROBES: spacecraft ('a','b','c','d','e') ; DATE: the date string or seconds since 1970 ('2007-03-23') ; DUR: duration (default units are days) ; DAYS: redundant keyword to set the units of duration (but its comforting to have) ; HOURS: keyword to make the duration be in units of hours ; DEVICE: sets the device (x or z) (default is x) ; MAKEPNG: keyword to generate 5 png files ; DIRECTORY: sets the directory where the above pngs are placed (default is './') ; DONT_DELETE_DATA: keyword to not delete all existing tplot variables before loading data in for ; the overview plot (sometimes old variables can interfere with overview plot) ; ERROR: Tells the calling routine whether an error that prevented overplot completion occured. (error=1 indicates an error, error=0 indicates no error) ; ;WARNING: ; Any code that requires windowStorage or loadedData should be put in the "~keyword_set(no_draw)" block. ; spd_ui_overplot is called without valid values for these parameters during document load ; ; ; Example: thm_gen_overplot,probe='a',date='2007-03-23',dur=1 ; The above example will produce a full day plot in the X window. ; ; ;Version: ; $LastChangedBy: nikos $ ; $LastChangedDate: 2014-03-20 18:58:37 -0700 (Thu, 20 Mar 2014) $ ; $LastChangedRevision: 14624 $ ; $URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/trunk/spedas/gui/utilities/spd_ui_overplot.pro $ ;- ;functionalizes repeated code, this is called after each panel is added to gui ;it should help manage memory pro spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data compile_opt hidden if ~keyword_set(dont_delete_data) then begin to_delete = ssl_set_complement([tn_before],[tnames('*'),'']) if(size(to_delete,/n_dim) ne 0 || to_delete[0] ne -1L) then $ del_data,to_delete endif end Pro quick_set_panel_labels, panel_in, string_in, colors_in = colors_in, zaxis = zaxis, zhorizontal = zhorizontal compile_opt hidden If(is_string(string_in,/blank) Eq 0) Then Return n = n_elements(string_in) If(keyword_set(colors_in)) Then Begin nc = n_elements(colors_in[0,*]) ;Colors_in should be 3xN, if there aren't the same number of colors, ;as labels, then use the first color If(nc Ne n) Then colors = rebin(colors_in[*,0], 3, n) $ Else colors = colors_in Endif Else colors = rebin([0b, 0b, 0b], 3, n) If(keyword_set(zaxis)) Then Begin panel_in -> getproperty, zaxis = zobj If(obj_valid(zobj)) Then Begin to = obj_new('spd_ui_text') to -> setproperty, value = string_in[0] ;only one label to -> setproperty, size = 8.0 zobj -> setproperty, labeltextobject = to if keyword_set(zhorizontal) then begin zobj -> setproperty, labelorientation = 0 zobj -> setproperty, labelmargin = 30.0 endif else zobj -> setproperty, labelmargin = 45.0 Endif Endif Else Begin panel_in -> getproperty, yaxis = yobj If(obj_valid(yobj)) Then Begin yobj -> setproperty, margin = 32.0 yobj -> getproperty, labels = obj2 ;added to support single simplified labels yobj -> setproperty, blacklabels=1 If obj_valid(obj2) then begin lobj = obj2 -> get(/all) ;Keep only the first object, clean out the rest lobj0 = lobj[0] If(n_elements(lobj) Gt 1) Then For j = 1, n_elements(lobj)-1 Do obj_destroy, lobj[j] obj2 -> remove, /all lobj0 -> setproperty, color = colors[*, 0] lobj0 -> setproperty, value = string_in[0] lobj0 -> setproperty, size = 10.0 lobj0 -> setproperty, show = 1 obj2 -> add, lobj0 If(n Gt 1) Then Begin For j = 1, n-1 Do Begin lobjj = lobj0 -> copy() lobjj -> setproperty, color = colors[*, j] lobjj -> setproperty, value = string_in[j] lobjj -> setproperty, show = 1 obj2 -> add, lobjj Endfor Endif Endif Endif Endelse Return End ;boolean helper function to check that time dimensions ;match y dimensions function check_keogram_dims,names compile_opt idl2,hidden if names[0] eq '' then return, '' tr = timerange(/current) for i = 0,n_elements(names)-1 do begin name = names[i] get_data,name,data=d ;checks that the quantity is well formed if is_struct(d) && $ in_set('x',strlowcase(tag_names(d))) && $ in_set('y',strlowcase(tag_names(d))) then begin dim_x = dimen(d.x) dim_y = dimen(d.y) ;this checks that the dimensions aren't messed ;It also checks that this particular quantity is ;from the date/duration requested if dim_x[0] eq dim_y[0] && $ d.x[0] le tr[1] && $ d.x[dim_x[0]-1] ge tr[0] then return,name endif endfor return,'' end ;updates the metadata of quantities being loaded into the gui ;This way the information will be properly autodetected pro spd_ui_update_dlimits,varname,observatory,instrument,coord_sys,units compile_opt idl2,hidden if tnames(varname) eq '' then return get_data,varname,data=d,dlimit=dl ;handle pseudo variables with a recursive call if is_string(d[0]) then begin if n_elements(d) eq 1 then begin d = strsplit(d,' ') endif for i = 0,n_elements(d)-1 do begin spd_ui_update_dlimits,d[i],observatory,instrument,coord_sys,units endfor return endif if ~is_struct(dl) then begin dl = {data_att:''} endif str_element,dl,'data_att',success=s if ~s then begin str_element,dl,'data_att',/add, {PROJECT:'SPEDAS',$ OBSERVATORY:observatory,$ INSTRUMENT:instrument,$ COORD_SYS:coord_sys,$ UNITS:units} endif else begin data_att = dl.data_att str_element,data_att,'PROJECT',/add,'SPEDAS' str_element,data_att,'OBSERVATORY',/add,observatory str_element,data_att,'INSTRUMENT',/add,instrument str_element,data_att,'COORD_SYS',/add,coord_sys str_element,data_att,'UNITS',/add,units str_element,dl,'data_att',/add,data_att endelse store_data,varname,dl=dl end ;WARNING: ; Any code that requires windowStorage or loadedData should be put in the "~keyword_set(no_draw)" block. ; spd_ui_overplot is called without valid values for these parameters during document load ; ;This code adds overview plot to the current window. Calling routine should create a new window for the overview plot. ;This prevents weak window menu management code from exploding. pro spd_ui_overplot, windowStorage,loadedData,drawObject,$ probes = probes, date = date, dur = dur, $ days = days, hours = hours, device = device, $ directory = directory, makepng = makepng, $ dont_delete_data = dont_delete_data, $ oplot_calls=oplot_calls,no_draw=no_draw,$ error=error ;catch errors and fail gracefully ;------------------------------------------------------- common overplot_position, load_position, error_count error=1 ;indicate error to calling routine(set to zero if the routine reaches completion !quiet=0 error_count=0 load_position = 'init' tn_before = tnames('*') ; get list of pre-existing tvars ;catch statement to allow program to recover from errors ;------------------------------------------------------- catch,error_status if error_status ne 0 then begin error_count++ if error_count ge 1000. then begin dprint, ' ' dprint, 'The program is quitting because it fears its in an infinite loop.' dprint, 'To eliminate this fear add the keyword /fearless to the call.' return endif if !error_state.name eq 'THM_SPINMODEL_POST_PROCESS_NO_TVAR' then begin ok = dialog_message('Cannot find STATE data. Check the GUI configuration'+ $ 'to verify data is being downloaded.', /center, $ title='Error creating overview plot') return endif print, '***********Catch error**************' help, /last_message, output = err_msg For j = 0, n_elements(err_msg)-1 Do print, err_msg[j] print, 'load_position: ' , load_position case load_position of 'init' : goto, SKIP_DAY 'fgm' : goto, SKIP_FGM_LOAD 'fbk' : goto, SKIP_FBK_LOAD 'sst' : goto, SKIP_SST_LOAD 'esa' : goto, SKIP_ESA_LOAD 'gmag' : goto, SKIP_GMAG_LOAD 'roi' : goto, SKIP_ROI_LOAD 'asi' : goto, SKIP_ASI_LOAD 'pos' : goto, SKIP_POS_LOAD 'mode' : goto, SKIP_SURVEY_MODE 'bound' : goto, SKIP_BOUNDS else : goto, SKIP_DAY endcase endif ;check some inputs ;------------------------------------------------------- if keyword_set(probes) then sc = strlowcase(probes) ; quick change of variable name vsc = ['a','b','c','d','e'] if not keyword_set(sc) then begin dprint, 'You did not enter a spacecraft into the program call.' dprint, "Valid inputs are: 'a','b','c','d','e' (ie, sc='b')" return endif if total(strmatch(vsc,strtrim(strlowcase(sc)))) gt 1 then begin dprint, 'This program is only designed to accept a single spacecraft as input.' dprint, "Valid inputs are: 'a','b','c','d','e' (ie, sc='b')" return endif if total(strmatch(vsc,strtrim(strlowcase(sc)))) eq 0 then begin dprint, "The input sc= '",strtrim(sc),"' is not a valid input." dprint, "Valid inputs are: 'a','b','c','d','e' (ie, sc='b')" return endif if ~keyword_set(dur) then begin dprint, 'duration not input, setting dur = 1' dur = 1 endif if (dur Lt 0) then begin dprint, 'Invalid duration, setting dur = 1' dur = 1 endif if keyword_set(hours) then dur=dur/24. if ~keyword_set(date) then begin dprint, 'You did not enter a date into the program call.' dprint, "Example: thm_gen_overplot,sc='b',date='2007-03-23'" return endif else begin t0 = time_double(date) t1 = t0+dur*60D*60D*24D if(t1 Lt time_double('2007-02-17/00:00:00')) then begin dprint, 'Invalid time entered: ', time_string(date) ok = dialog_message(time_string(date) + ' is not a valid date. '+ $ 'Overview plots cannot be generated before 2007-02-17.', $ /center, title='Generate Overview Plot') return endif if(t0 Gt systime(/seconds)) then begin dprint, 'Invalid time entered: ', time_string(date) ok = dialog_message(time_string(date) + ' is not a valid date. '+ $ 'Overview plots cannot be generated for future dates.', $ /center, title='Generate Overview Plot') return endif endelse if not keyword_set(device) then begin help,/device,output=plot_device plot_device=strtrim(strlowcase(strmid(plot_device(1),24)),2) if plot_device eq 'z' then device, set_resolution = [750, 800] endif else begin set_plot,device help,/device,output=plot_device plot_device=strtrim(strlowcase(strmid(plot_device(1),24)),2) if plot_device eq 'z' then device, set_resolution = [750, 800] endelse thx = 'th'+sc[0] ;need a scalar for this date=time_string(date) if ~keyword_set(dont_delete_data) then begin clear_esa_common_blocks common data_cache_com, dcache dcache = '' endif widget_control, /hourglass tn_before = tnames('*') ; clear any previously set var_labels. This prevents the SPEDAS overview plots from ; using tplot's var_label from a different plot in the same session tplot_options, var_label = '' ; set the suffix to identify different calls to overplot osuffix = ('_op' + strcompress(string(*oplot_calls + 1), /remove_all))[0] timespan,date,dur ; load gmag data ;---------------- load_position='gmag' ;----------------------------------------------------------------------- ;Comment this block of code to disable AE index generation from GMAG ;thm_load_gmag,/subtract_median ;If(is_string(tnames('thg_mag_*'))) Then Begin ; split_vec,'thg_mag_*' ; superpo_histo,'thg_mag_*_x', dif='thg_pseudoAE', res=60.0 ; tdespike_AE,-2000.0,1500.0 ; clean_spikes, 'thg_pseudoAE_despike', new_name = 'thg_pseudoAE', thresh = 5 ; options,'thg_pseudoAE',ytitle='SPEDAS AE Index' ;Endif Else Begin ; filler=fltarr(2) ; filler[*]=float('NaN') ; store_data,'thg_pseudoAE',data={x:time_double(date)+findgen(2),y:filler} ; options,'thg_pseudoAE',ytitle='SPEDAS AE Index' ;Endelse ;spd_ui_update_dlimits,'thg_pseudoAE',sc,'ae_index','none','none' ;copy_data,'thg_pseudoAE','thg_pseudoAE'+osuffix ;----------------------------------------------------------------------- ;----------------------------------------------------------------------- ;Uncomment this block of code to implement AE index loading from CDF thm_load_pseudoAE If ~(is_string(tnames('thg_idx_ae'))) then begin filler=fltarr(2) filler[*]=float('NaN') store_data,'thg_idx_ae',data={x:[time_double(date),time_double(date)+dur*60*60*24],y:filler} endif options,'thg_idx_ae',ytitle='SPEDAS_AE Index' spd_ui_update_dlimits,'thg_idx_ae',sc,'ae_index','none','none' copy_data,'thg_idx_ae','thg_pseudoAE'+osuffix ;---------------------------------------------------------------------- tplot_gui,'thg_pseudoAE'+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_GMAG_LOAD: load_position='roi' thm_load_state,probe=sc,/get_support If(is_string(tnames(thx+'_state_roi'))) Then Begin ;test for state, jmm, 21-apr-2009 roi_bar = thm_roi_bar(thx+'_state_roi') ;change roi_bar to a 16xN tplot variable to avoid the need to write a ;GUI compatible bitplot routine, jmm, 9-apr-2009 get_data,roi_bar,data=ppp,limit=l If(is_struct(ppp)) Then Begin bits2,ppp.y,new_bar ;the extracts the bit values into a bytarr of 16xntimes nbv=n_elements(new_bar[*,0]) nxv=n_elements(ppp.x) new_bar=rebin(1+bindgen(nbv), nbv, nxv)*new_bar ;this assigns each non-zero bit the value equal to its bit number new_bar=float(new_bar) ;to allow for NaN's to plot correctly zerov=where(new_bar Eq 0,nzerov) If(nzerov Gt 0) Then new_bar[zerov]=!values.f_nan str_element,ppp,'y',transpose(new_bar),/add_replace ;need transpose to make it ntimesX16 str_element,l,'panel_size',/delete ;panel sizing handled using different method in gui overview plots store_data,roi_bar,data=ppp,limit=l options,roi_bar,tplot_routine='mplot' ;reset this for testing Endif Else Goto, dummy_roi_bar Endif Else Begin dummy_roi_bar: filler=fltarr(2,16) filler[*,*]=float('NaN') store_data,thx+'_roi_bar',data={x:time_double(date)+findgen(2),y:filler} roi_bar=thx+'_roi_bar' Endelse spd_ui_update_dlimits,roi_bar,sc,'roi','none','none' copy_data,roi_bar,roi_bar+osuffix tplot_gui,roi_bar+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_ROI_LOAD: ; load ASK data and plot 3 specific ones (can be changed) ;--------------------------------------------------------- load_position='asi' thm_load_ask, site='fsmi' filler = fltarr(2, 10) ; (10 chosen arbitrarily) filler(*, *) = float('NaN') ;Harald requests using FSMI as first choice: fsmi_site = tnames('*ask*fsmi*') if (fsmi_name = check_keogram_dims(fsmi_site)) ne '' then begin copy_data, fsmi_name, 'Keogram' keo_site = 'fsmi' endif else begin thm_load_ask,site=site_list,/valid_names found_k_site = 0 for i = 0,n_elements(site_list)-1 do begin thm_load_ask,site=site_list[i] asi_sites = tnames('*ask*') if (asi_name = check_keogram_dims(asi_sites)) ne '' then begin copy_data, asi_name, 'Keogram' split_name = strsplit(asi_name,'_',/extract) ;don't want entire name of quantity just site, which should be the last part keo_site = split_name[n_elements(split_name)-1] found_k_site = 1 break endif endfor if ~keyword_set(found_k_site) then begin store_data, 'Keogram', data = {x:time_double(date)+findgen(2), y:filler, v:findgen(10)} keo_site = 'filler' endif endelse options, 'Keogram', 'ytitle', ' ' spd_ui_update_dlimits,'Keogram',keo_site,'ask','none','counts' copy_data,'Keogram','Keogram'+osuffix tplot_gui,'Keogram'+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_ASI_LOAD: ;load magnetic field fit data ;----------------------------- load_position='fgm' thm_load_state,probe=sc,/get_support thm_load_fit,lev=1,probe=sc,/get_support ;kluge to prevent missing data from crashing things index_fit=where(thx+'_fgs' eq tnames()) index_state=where(thx+'_state_spinras' eq tnames()) if (index_fit(0) eq -1 or index_state(0) eq -1) then begin filler=fltarr(2,3) filler(*,*)=float('NaN') store_data,thx+'_fgs_gse',data={x:time_double(date)+findgen(2),y:filler} ylim,thx+'_fgs_gse',-100,100,0 endif else begin thm_cotrans,thx+'_fgs',out_suf='_gse', in_c='dsl', out_c='gse' endelse ;clip data ;tclip, thx+'_fgs_gse', -100.0, 100.0, /overwrite name = thx+'_fgs_gse' options, name, 'ytitle', 'B FIT_GSE (nT)' options, name, 'labels', ['Bx', 'By', 'Bz'] options, name, 'labflag', 1 options, name, 'colors', [2, 4, 6] options, name, 'yrange', [-100,100] copy_data,name,name+osuffix tplot_gui,name+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_FGM_LOAD: ;load SST spectrograms ;---------------------- load_position='sst' thm_load_sst, probe = sc, level = 'l2' ;If Level 2 data didn't show up, check for L1 index_sst_e = where(thx+'_psef_en_eflux' eq tnames()) index_sst_i = where(thx+'_psif_en_eflux' eq tnames()) if(index_sst_e[0] eq -1 Or index_sst_i[0] Eq -1) then begin thm_part_load,probe=probe,trange=trange,datatype='psif' thm_part_load,probe=probe,trange=trange,datatype='psef' thm_part_moments, probe = sc, instrument = ['psif', 'psef'], $ moments = ['density', 'velocity', 't3'],method_clean='automatic' endif SKIP_SST_LOAD: ;kluge to prevent missing data from crashing things index_sst=where(thx+'_psif_en_eflux' eq tnames()) if index_sst eq -1 then begin filler = fltarr(2, 16) filler[*,*]=float('NaN') store_data, thx+'_psif_en_eflux', $ data = {x:time_double(date)+findgen(2)*86400., y:filler, v:findgen(16)} name = thx+'_psif_en_eflux' options, name, 'spec', 1 ylim, name, 1, 1000, 1 zlim, name, 1d1, 5d2, 1 options, name, 'ytitle', thx+'!CSST ions!CeV' options, name, 'ysubtitle', '' ; options, name, 'ztitle', 'Eflux !C eV/cm!U2!N!C-s-sr-eV' options, name, 'ztitle', 'Eflux, EFU' endif else begin ;SST ion panel name = thx+'_psif_en_eflux' tdegap, name, /overwrite, dt = 600.0 options, name, 'spec', 1 options, name, 'ytitle', thx+'!CSST ions!CeV' options, name, 'ysubtitle', '' ; options, name, 'ztitle', 'Eflux !C eV/cm!U2!N!C-s-sr-eV' options, name, 'ztitle', 'Eflux, EFU' options, name, 'y_no_interp', 1 options, name, 'x_no_interp', 1 zlim, name, 1d1, 5d2, 1 endelse index_sst = where(thx+'_psef_en_eflux' eq tnames()) if index_sst eq -1 then begin filler = fltarr(2, 16) filler[*, *] = float('NaN') store_data, thx+'_psef_en_eflux', $ data = {x:time_double(date)+findgen(2), y:filler, v:findgen(16)} name = thx+'_psef_en_eflux' options, name, 'spec', 1 ylim, name, 1, 1000, 1 zlim, name, 1d1, 5d2, 1 options, name, 'ytitle', thx+'!CSST elec!CeV' options, name, 'ysubtitle', '' ; options, name, 'ztitle', 'Eflux !C eV/cm!U2!N!C-s-sr-eV' options, name, 'ztitle', 'Eflux, EFU' endif else begin ;SST electron panel name = thx+'_psef_en_eflux' tdegap, name, /overwrite, dt = 600.0 options, name, 'spec', 1 options, name, 'ytitle', thx+'!CSST elec!CeV' options, name, 'ysubtitle', '' ; options, name, 'ztitle', 'Eflux !C eV/cm!U2!N!C-s-sr-eV' options, name, 'ztitle', 'Eflux, EFU' options, name, 'y_no_interp', 1 options, name, 'x_no_interp', 1 zlim, name, 1d1, 5d2, 1 endelse ;load ESA spectrograms and moments ;---------------------------------- load_position='esa' ;load both full and reduced data: mtyp = ['f', 'r'] ok_esai_flux = bytarr(2) ok_esae_flux = bytarr(2) ok_esai_moms = bytarr(2) ok_esae_moms = bytarr(2) For j = 0, 1 Do Begin thm_load_esa, probe = sc, datatype = 'pe?'+mtyp[j]+'*', level = 'l2' itest = thx+'_pei'+mtyp[j] etest = thx+'_pee'+mtyp[j] ;If Level 2 data didn't show up, check for L1 index_esa_e_en = where(etest+'_en_eflux' eq tnames()) index_esa_e_d = where(etest+'_density' eq tnames()) index_esa_e_v = where(etest+'_velocity_dsl' eq tnames()) index_esa_e_t = where(etest+'_t3' eq tnames()) index_esa_i_en = where(itest+'_en_eflux' eq tnames()) index_esa_i_d = where(itest+'_density' eq tnames()) index_esa_i_v = where(itest+'_velocity_dsl' eq tnames()) index_esa_i_t = where(itest+'_t3' eq tnames()) if(index_esa_e_en[0] eq -1 Or index_esa_i_en[0] Eq -1) then begin thm_load_esa_pkt, probe = sc instr_all = ['pei'+mtyp[j], 'pee'+mtyp[j]] for k = 0, 1 do begin test_index = where(thx+'_'+instr_all[k]+'_en_counts' eq tnames()) If(test_index[0] Ne -1) Then Begin thm_part_moments, probe = sc, instrument = instr_all[k], $ moments = '*' copy_data, thx+'_'+instr_all[k]+'_velocity', $ thx+'_'+instr_all[k]+'_velocity_dsl' Endif endfor index_esa_e_en = where(etest+'_en_eflux' eq tnames()) index_esa_e_d = where(etest+'_density' eq tnames()) index_esa_e_v = where(etest+'_velocity_dsl' eq tnames()) index_esa_e_t = where(etest+'_t3' eq tnames()) index_esa_i_en = where(itest+'_en_eflux' eq tnames()) index_esa_i_d = where(itest+'_density' eq tnames()) index_esa_i_v = where(itest+'_velocity_dsl' eq tnames()) index_esa_i_t = where(itest+'_t3' eq tnames()) endif if index_esa_i_en[0] eq -1 then begin filler = fltarr(2, 32) filler[*, *] = float('Nan') name1 = itest+'_en_eflux' store_data, name1, data = {x:time_double(date)+findgen(2), y:filler, v:findgen(32)} zlim, name1, 1d3, 7.5d8, 1 ylim, name1, 3., 40000., 1 ; options, name1, 'ztitle', 'Eflux !C!C eV/cm!U2!N!C-s-sr-eV' options, name1, 'ztitle', 'Eflux, EFU' ; options, name1, 'ytitle', 'ESA i+ '+thx+'!C eV' options, name1, 'ysubtitle', '' options, name1, 'spec', 1 endif else begin name1 = itest+'_en_eflux' tdegap, name1, /overwrite, dt = 600.0 zlim, name1, 1d3, 7.5d8, 1 ylim, name1, 3., 40000., 1 ; options, name1, 'ztitle', 'Eflux !C!C eV/cm!U2!N!C-s-sr-eV' options, name1, 'ztitle', 'Eflux, EFU' options, name1, 'ytitle', 'ESA i+ '+thx+'!C eV' options, name1, 'ysubtitle', '' options, name1, 'spec', 1 options, name1, 'x_no_interp', 1 options, name1, 'y_no_interp', 1 ok_esai_flux[j] = 1 endelse if index_esa_i_d[0] eq -1 then begin filler = fltarr(2) filler[*] = float('Nan') store_data, itest+'_density', data = {x:time_double(date)+findgen(2), y:filler} ; options, itest+'_density', 'ytitle', 'Ni '+thx+'!C!C1/cm!U3' options, itest+'_density', 'ytitle', 'Ni '+thx endif else begin name1 = itest+'_density' tdegap, name1, /overwrite, dt = 600.0 ylim, name1, .1, nmax, 1 options, name1, 'ytitle', 'Ni '+thx ok_esai_moms[j] = 1 endelse if index_esa_i_v[0] eq -1 then begin filler = fltarr(2, 3) filler[*, *] = float('Nan') store_data, itest+'_velocity_dsl', data = {x:time_double(date)+findgen(2), y:filler} options, itest+'_velocity_dsl', 'ytitle', 'VI '+thx+'!Ckm/s' options, itest+'_velocity_dsl', 'ysubtitle', '' endif else begin name1 = itest+'_velocity_dsl' tdegap, name1, /overwrite, dt = 600.0 itstrg=[t0,t1] get_ylimits, name1, itslimits, itstrg minmaxvals=itslimits.yrange maxvel=max(abs(minmaxvals)) maxlim=min([maxvel,2000.]) minlim=0.-maxlim if maxvel le 100. then ylim, name1, -50,50,0 else ylim, name1, minlim, maxlim, 0 options, name1, 'colors', [2, 4, 6] options, name1, 'labflag', 1 options, name1, 'ytitle', 'VI '+thx+'!Ckm/s' options, name1, 'ysubtitle', '' ;; options, name1, labels = ['Vi!dx!n', 'Vi!dy!n', 'Vi!dz!n'], constant = 0. options, name1, labels = ['VIx', 'VIy', 'VIz'], constant = 0. endelse if index_esa_i_t[0] eq -1 then begin filler = fltarr(2, 6) filler[*, *] = float('Nan') store_data, itest+'_t3', data = {x:time_double(date)+findgen(2), y:filler} options, itest+'_t3', 'ytitle', 'Ti '+thx+'!CeV' options, itest+'_t3', 'ysubtitle', '' endif else begin name1 = itest+'_t3' tdegap, name1, /overwrite, dt = 600.0 ylim, name1, 10, 10000., 1 options, name1, 'colors', [2, 4, 6, 0] options, name1, 'ytitle', 'Ti '+thx+'!C eV' options, name1, 'ysubtitle', '' endelse index_esa_e_en = where(etest+'_en_eflux' eq tnames()) if index_esa_e_en[0] eq -1 then begin filler = fltarr(2, 32) filler[*, *] = float('Nan') name1 = etest+'_en_eflux' store_data, name1, data = {x:time_double(date)+findgen(2), y:filler, v:findgen(32)} zlim, name1, 1d4, 7.5d8, 1 ylim, name1, 3., 40000., 1 ; options, name1, 'ztitle', 'Eflux !C!C eV/cm!U2!N!C-s-sr-eV' options, name1, 'ztitle', 'Eflux, EFU' options, name1, 'ytitle', 'ESA e- '+thx+'!C eV' options, name1, 'ysubtitle', '' options, name1, 'spec', 1 endif else begin name1 = etest+'_en_eflux' tdegap, name1, /overwrite, dt = 600.0 zlim, name1, 1d4, 7.5d8, 1 ylim, name1, 3., 40000., 1 ; options, name1, 'ztitle', 'Eflux !C!C eV/cm!U2!N!C-s-sr-eV' options, name1, 'ztitle', 'Eflux, EFU' options, name1, 'ytitle', 'ESA e- '+thx+'!C eV' options, name1, 'ysubtitle', '' options, name1, 'spec', 1 options, name1, 'x_no_interp', 1 options, name1, 'y_no_interp', 1 ok_esae_flux[j] = 1 endelse if index_esa_e_d[0] eq -1 then begin filler = fltarr(2) filler[*] = float('Nan') store_data, etest+'_density', data = {x:time_double(date)+findgen(2), y:filler} ; options, etest+'_density', 'ytitle', 'Ne '+thx+'!C!C1/cm!U3' options, etest+'_density', 'ytitle', 'Ne '+thx+'!C1/cc' options, etest+'_density', 'ysubtitle', '' no_npot: filler = fltarr(2) filler[*] = float('Nan') store_data, etest+'_density_npot', data = {x:time_double(date)+findgen(2), y:filler} options, etest+'_density_npot', 'ytitle', 'Ne '+thx+'!C1/cc' options, etest+'_density_npot', 'ysubtitle', '' endif else begin name1 = etest+'_density' ylim, name1, .1, nmax, 1 ; options, name1, 'ytitle', 'Ne '+thx+'!C!C1/cm!U3' options, name1, 'ytitle', 'Ne '+thx+'!C1/cc' options, name1, 'ysubtitle', '' ok_esae_moms[j] = 1 ;Npot calculation, 2009-10-12, jmm thm_scpot2dens_opt_n, probe = sc, /no_data_load, datatype_esa = 'pee'+mtyp[j] ;degap after npot calculation tdegap, name1, /overwrite, dt = 600.0 name1x = tnames(etest+'_density_npot') get_data, name1x, data = npot_test If(is_struct(temporary(npot_test)) Eq 0) Then Goto, no_npot tdegap, name1x, /overwrite, dt = 600.0 options, name1x, 'ytitle', 'Ne '+thx+'!C1/cc' options, name1x, 'ysubtitle', '' endelse if index_esa_e_v[0] eq -1 then begin filler = fltarr(2, 3) filler[*, *] = float('Nan') store_data, etest+'_velocity_dsl', data = {x:time_double(date)+findgen(2), y:filler} ; options, etest+'_velocity_dsl', 'ytitle', 'Ve '+thx+'!C!Ckm/s' options, etest+'_velocity_dsl', 'ytitle', 'VE '+thx+'!Ckm/s' options, etest+'_velocity_dsl', 'ysubtitle', '' endif else begin name1 = etest+'_velocity_dsl' tdegap, name1, /overwrite, dt = 600.0 ylim, name1, -500, 200., 0 ; options, name1, 'ytitle', 'Ve '+thx+'!C!Ckm/s' options, name1, 'ytitle', 'VE '+thx+'!Ckm/s' options, name1, 'ysubtitle', '' endelse if index_esa_e_t[0] eq -1 then begin filler = fltarr(2, 6) filler[*, *] = float('Nan') store_data, etest+'_t3', data = {x:time_double(date)+findgen(2), y:filler} options, etest+'_t3', 'ytitle', 'Te '+thx+'!CeV' options, etest+'_t3', 'ysubtitle', '' endif else begin ;options,name1,'colors',[cols.blue,cols.green,cols.red] ; options, name1, labels = ['V!dex!n', 'V!dey!n', 'V!dez!n'], constant = 0. name1 = etest+'_t3' tdegap, name1, /overwrite, dt = 600.0 options, name1, labels = ['TEx', 'TEy', 'TEz'], constant = 0. options, name1, labels = ['Ti_para', 'Ti_perp', 'Te_para', 'Te_perp'] ylim, name1, 10, 10000., 1 options, name1, 'colors', [2, 0, 4, 6] options, name1, 'ytitle', 'TE '+thx+'!CeV' options, name1, 'ysubtitle', '' endelse ; plot quantities (manipulating the plot quantities for the sake of plot aesthetics) ;kluge for labeling the density, added Npot, 2009-10-12, jmm get_data, etest+'_density', data = d get_data, etest+'_density_npot', data = d1 Ne_kluge_name = 'Ne_'+etest+'_kluge' If(n_elements(d1.x) Eq n_elements(d.x)) Then Begin dummy = fltarr(n_elements(d.y), 3) dummy[*, 0] = d1.y dummy[*, 1] = d.y dummy[*, 2] = d.y store_data, Ne_kluge_name, data = {x:d.x, y:dummy} options, Ne_kluge_name, labels = ['Npot', 'Ni', 'Ne'] options, Ne_kluge_name, colors = [2, 0, 6] options, Ne_kluge_name, 'labflag', 1 Endif Else Begin dummy = fltarr(n_elements(d.y), 2) dummy[*, 0] = d.y dummy[*, 1] = d.y store_data, Ne_kluge_name, data = {x:d.x, y:dummy} options, Ne_kluge_name, labels = ['Ni', 'Ne'] options, Ne_kluge_name, colors = [0, 6] options, Ne_kluge_name, 'labflag', 1 Endelse ;copy data to a different name so that it doesn't overwrite gui data when imported copy_data,itest+'_density', itest+'_density'+osuffix copy_data,Ne_kluge_name, Ne_kluge_name+osuffix store_data, thx+'_Nie'+mtyp[j], data = [itest+'_density', Ne_kluge_name]+osuffix ; options, thx+'_Nie'+mtyp[j], 'ytitle', 'Ni,e '+thx+'!C1/cm!U3' ; options, thx+'_Nie'+mtyp[j], 'ytitle', 'Ni,e '+thx options, thx+'_Nie'+mtyp[j], 'ytitle', 'Ni,e '+thx+'!C1/cc' options, thx+'_Nie'+mtyp[j], 'ysubtitle', '' nameti=itest+'_t3' namete=etest+'_t3' ;copy data to a different name so that it doesn't overwrite gui data when imported copy_data,nameti,nameti+osuffix copy_data,namete,namete+osuffix store_data, thx+'_Tie'+mtyp[j], data = [nameti,namete]+osuffix options, thx+'_Tie'+mtyp[j], 'ytitle', 'Ti,e '+thx+'!CeV' options, thx+'_Tie'+mtyp[j], 'ysubtitle', '' options,nameti,'labels',['Ti!9'+string(120B)+'!X','','Ti!9'+string(35B)+'!X'] options,namete,'labels',[' Te!9'+string(120B)+'!X',' ',' Te!9'+string(35B)+'!X'] options, thx+'_Tie'+mtyp[j], 'labflag', 1 options,nameti, 'colors', [2,2, 4] options,namete, 'colors', [6,6, 0] Endfor SKIP_ESA_LOAD: load_position='mode' ; make tplot variable tracking the sample rate (0=SS,1=FS,2=PB,3=WB) ;------------------------------------------------------------------- sample_rate_var = thm_sample_rate_bar(date, dur, sc, /outline) get_data,sample_rate_var,limit=l str_element,l,'panel_size',/delete ;panel arrangement handled through kluge below. Using proper keywords breaks results store_data,sample_rate_var,limit=l options, sample_rate_var,'ytitle','' ;copy data to prevent overwriting data when importing into GUI get_data,sample_rate_var,data=d for sample_loop_var = 0,n_elements(d)-1 do begin copy_data,d[sample_loop_var],d[sample_loop_var]+osuffix endfor ;run degap procedure on sample_rate_bar variables tdegap, ['wave_burst_sym_', 'particle_burst_sym_']+sc[0], dt=600.0, /overwrite spd_ui_update_dlimits,'sample_rate_'+sc,sc,'sample_rate','none','none' SKIP_SURVEY_MODE: load_position='bound' ; final tplot preparations ;-------------------------- ; plot it! thm_spec_lim4overplot, thx+'_peif_en_eflux', zlog = 1, ylog = 1, /overwrite thm_spec_lim4overplot, thx+'_peef_en_eflux', zlog = 1, ylog = 1, /overwrite thm_spec_lim4overplot, thx+'_peir_en_eflux', zlog = 1, ylog = 1, /overwrite thm_spec_lim4overplot, thx+'_peer_en_eflux', zlog = 1, ylog = 1, /overwrite ssti_name=thx+'_psif_en_eflux' sste_name=thx+'_psef_en_eflux' thm_spec_lim4overplot, ssti_name, zlog = 1, ylog = 1, /overwrite ;reset sst ylimit maxima to 3.0e6 get_data, ssti_name, data = d If(is_struct(d)) Then ylim, ssti_name, min(d.v), 3.0e6, 1 thm_spec_lim4overplot, sste_name, zlog = 1, ylog = 1, /overwrite get_data, sste_name, data = d If(is_struct(d)) Then ylim, sste_name, min(d.v), 3.0e6, 1 ;tplot_options, 'lazy_ytitle', 0 ; prevent auto formatting on ytitle (namely having carrage returns at underscores) ;!p.background=255. ;!p.color=0. ;time_stamp,/off ;loadct2,43 ;!p.charsize=0.6 scv = strcompress(strlowcase(sc[0]),/remove_all) pindex = where(vsc Eq scv) ;this is always true for one probe by the time we are here ;tplot_options,'ygap',0.0D ;For esa data we would like to plot full mode if possible, but reduced ;mode if no full mode is available esaif_flux_name = thx+'_peif_en_eflux' If(ok_esai_flux[0] Eq 0) Then Begin ;esa ion flux is not present full resolution, If(ok_esai_flux[1]) Then esaif_flux_name = thx+'_peir_en_eflux' Endif esaif_v_name = thx+'_peif_velocity_dsl' If(ok_esai_moms[0] Eq 0) Then Begin If(ok_esai_moms[1]) Then esaif_v_name = thx+'_peir_velocity_dsl' Endif esaef_flux_name = thx+'_peef_en_eflux' If(ok_esae_flux[0] Eq 0) Then Begin ;esa electron flux is not present full resolution, rename if possible If(ok_esae_flux[1]) Then esaef_flux_name = thx+'_peer_en_eflux' Endif esaef_v_name = thx+'_peef_velocity_dsl' If(ok_esae_moms[0] Eq 0) Then Begin If(ok_esae_moms[1]) Then esaef_v_name = thx+'_peer_velocity_dsl' Endif esaf_t_name = thx+'_Tief' ;T and N are done for ions, electrons together esaf_n_name = thx+'_Nief' If(ok_esai_moms[0] Eq 0) Then Begin If(ok_esai_moms[1]) Then Begin esaf_t_name = thx+'_Tier' esaf_n_name = thx+'_Nier' Endif Endif ; kludge to make y titles consistent with the summary plots from the website ;options, esaf_n_name, 'ytitle', ' ' ;options, esaif_v_name, 'ytitle', ' ' options, roi_bar, 'ytitle', ' ' options, thg_pseudoAE, 'ytitle', ' ' options, 'sample_rate_'+sc, 'ytitle', ' ' options, ssti_name, 'ysubtitle', ' ' options, esaif_flux_name, 'ysubtitle', ' ' options, sste_name, 'ysubtitle', ' ' options, esaef_flux_name, 'ytitle', ' ' options, ssti_name, 'ytitle', ' ' options, esaif_flux_name, 'ytitle', ' ' options, sste_name, 'ytitle', ' ' options, esaef_flux_name, 'ytitle', ' ' names = [esaf_n_name, esaif_v_name, esaf_t_name, 'sample_rate_'+sc, $ ssti_name, esaif_flux_name, sste_name, $ esaef_flux_name] for i = 0,n_elements(names)-1 do begin copy_data,names[i],names[i]+osuffix endfor tplot_gui,names+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_BOUNDS: ;load FBK data ;-------------- load_position='fbk' thm_load_fbk,probe=sc ;fbk_tvars=tnames(thx+'_fb_*') ;this should give us two tplot variables (but sometimes more) ;Set fbk variables to spectrograms, in dlimits fbk_tvars = [tnames(thx+'_fb_e*'), tnames(thx+'_fb_s*')] if is_string(fbk_tvars) eq 0 then fbk_tvars=['fbk_filler','fbk_filler'] ;need two blank panels if n_elements(fbk_tvars) eq 1 then fbk_tvars=[fbk_tvars,'fbk_filler'] ;need one blank panel for i=0,n_elements(fbk_tvars)-1 do begin ;kluge to prevent missing data from crashing things get_data,fbk_tvars[i],data=dd,dlimits=dl if size(dd,/type) ne 8 then begin filler = fltarr(2, 6) filler[*, *] = float('NaN') name = thx+'_fb_NaN'+strcompress(string(i+1), /remove_all) store_data, name, data = {x:time_double(date)+findgen(2), y:filler, v:findgen(6)} options, name, 'spec', 1 ylim, name, 1, 1000, 1 zlim, name, 0, 0, 1 fbk_tvars[i] = name endif else begin options, fbk_tvars[i], 'spec', 1 options, fbk_tvars[i], 'zlog', 1 ylim, fbk_tvars[i], 2.0, 2048.0, 1 thm_spec_lim4overplot, fbk_tvars[i], ylog = 1, zlog = 1, /overwrite options, fbk_tvars[i], 'ysubtitle', ' ' options, fbk_tvars[i], 'ytitle', ' ' ; options, fbk_tvars[i], 'ytitle', thx+'_FBK '+strmid(fbk_tvars[i], 7) ;for ztitle, we need to figure out which type of data is there ; for V channels, <|V|>. ; for E channels, <|mV/m|>. ; for SCM channels, <|nT|>. x1 = strpos(fbk_tvars[i], 'scm') If(x1[0] Ne -1) Then Begin options, fbk_tvars[i], 'ztitle', '<|nT|>' ;reset the upper value of zlimit to 2.0, jmm, 30-nov-2007 get_data, fbk_tvars[i], data = d If(is_struct(d)) Then zlim, fbk_tvars[i], min(d.y), 2.0, 1 Endif xv = strpos(fbk_tvars[i], 'v') If(xv[0] Ne -1) Then options, fbk_tvars[i], 'ztitle', '<|V|>' xe = strpos(fbk_tvars[i], 'e') If(xe[0] Ne -1) Then Begin options, fbk_tvars[i], 'ztitle', '<|mV/m|>' ;reset the upper value of zlimit to 2.0, jmm, 30-nov-2007 get_data, fbk_tvars[i], data = d If(is_struct(d)) Then zlim, fbk_tvars[i], min(d.y), 2.0, 1 Endif ;Explicitly set fbk variables to spectrograms, in dlimits If(is_struct(dl)) Then Begin str_element, dl, 'spec', 1, /add_replace str_element, dl, 'log', 1, /add_replace store_data, fbk_tvars[i], dlimits = dl Endif endelse endfor for i = 0,n_elements(fbk_tvars)-1 do begin copy_data,fbk_tvars[i],fbk_tvars[i]+osuffix endfor tplot_gui,fbk_tvars+osuffix,/no_verify,/no_update,/add_panel,no_draw=keyword_set(no_draw) spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data SKIP_FBK_LOAD: ; Get position info ;--------------------------------------------------------- load_position='pos' thm_load_state,probe=sc thm_cotrans,thx+'_state_pos',out_suf='_gse',in_coord='gei',out_coord='gse' get_data, thx+'_state_pos_gse',data=tmp ;creating new state variables in unnecessary now, jmm, 15-apr-2009 If(is_struct(tmp)) Then begin store_data,thx+'_state_pos_gse',data={x:tmp.x,y:tmp.y/6370.} endif else begin store_data,thx+'_state_pos_gse',data={x:time_double(date)+findgen(2),y:dblarr(2,3)+!VALUES.D_NAN} endelse SKIP_POS_LOAD: load_position='plot' vars_full = ['thg_pseudoAE', roi_bar, 'Keogram', thx+'_fgs_gse', $ esaf_n_name, esaif_v_name, esaf_t_name, 'sample_rate_'+sc, $ ssti_name, esaif_flux_name, sste_name, $ esaef_flux_name, fbk_tvars] ; add suffix to tplot vars so that multiple loads won't interfere w/each other ;nvars = n_elements(vars_full) ;;osuffix = '_test' ;vars_full_sfx = vars_full + osuffix ;for i=0, nvars-1 do begin ; ; get_data, vars_full[i], data=d, dlimits=dl, limits=l ; if size(d,/type) eq 7 then begin ; ;handle pseudovariables ; ; for j=0, n_elements(d)-1 do begin ; copy_data, d[j], d[j] + osuffix ; del_data, d[j] ; endfor ; store_data, vars_full[i], data=d+osuffix, dlimits=dl, limits=l ; endif ; copy_data, vars_full[i], vars_full_sfx[i] ;endfor ;del_data, vars_full ; ;tplot_gui, vars_full_sfx, /no_verify,no_draw=no_draw probes_title = ['P5', 'P1', 'P2', 'P3', 'P4'] ;Ok, this is the 2nd plot, what you need to do here is get the GUI ;info object and redo the plots once. if ~keyword_set(no_draw) then begin ;Only one active window at a time wo = windowstorage -> getactive() If(n_elements(wo) gt 0 && obj_valid(wo[0])) Then Begin wo[0] -> setproperty, tracking = 0 wo[0] -> getproperty, panels = panelsj panel = panelsj -> get(/all) Endif current_row = 1 for i = 0,n_elements(panel)-1 do begin panel[i]->getProperty,xAxis=xobj,yAxis=yobj,zAxis=zObj,settings=panel_settings ;fix the range of the plot to the requested range. if i eq 0 then begin xobj->setProperty,rangeOption=2,minFixedRange=time_double(date),maxFixedRange=time_double(date)+dur*24*60*60 endif ;decrease the number of ticks on the panels to avoid clutter if i eq 1 || i eq 7 then begin xObj->setProperty,majorLength=2,minorLength=1 yobj->setProperty,numMajorTicks=0,autoticks=0,numMinorTicks=0 endif else begin xObj->setProperty,majorLength=4,minorLength=2 yobj->setProperty,numMajorTicks=3,numMinorTicks=1 endelse ; ; ;autoticks = 0 guarantees that our y/z axis tick settings will not be overridden ; yObj->setProperty,autoticks=1 ; ; if obj_valid(zObj) then begin ; ; zObj->setProperty,tickNum=3,minorTickNum=1,autoticks=1 ; ; endif ;modify layout to make status bars smaller if i eq 1 || i eq 7 then begin if i eq 1 then begin panel_settings->setProperty,row=current_row,rSpan=2 current_row += 2 endif else begin panel_settings->setProperty,row=current_row,rSpan=1 current_row += 1 endelse endif else begin panel_settings->setProperty,row=current_row,rSpan=4 current_row += 4 endelse endfor ;set the vertical spacing to a smaller number, to save space on a many panel layout wo[0]->getProperty,settings=page page->setProperty,yPanelSpacing=0 ;set the total number of rows wo[0]->setProperty,nRows=current_row-1 ;get trace for the first panel and change the color to black, and set ;the title to = obj_new('spd_ui_text') to -> setproperty, value = probes_title[pindex[0]]+' (TH-'+strupcase(sc)+')' to -> setproperty, size=8. ;panel[0] -> getproperty, settings = p0settings ;p0settings -> setproperty, titleobj = to page->setproperty, title=to panel[0] -> getproperty, tracesettings = obj0 trace_obj = obj0 -> get(/all) if obj_valid(trace_obj[0]) then begin trace_obj[0] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [0b, 0b, 0b] endif ;set labels ; quick_set_panel_labels, panel[0], ['SPEDAS AE Index'] ;quick_set_panel_labels, panel[0], ['SPEDAS', 'AE Index'] ;ROI panel now, this needs 1 black label ;ROI also needs 1:1 labels:traces ; quick_set_panel_labels, panel[1], ['ROI','','','','','','','','','','','','','','',''] ; quick_set_panel_labels, panel[1], ['Region_Panel'] panel[1]->getProperty,yaxis=yobj yobj->setProperty,annotateAxis=0,lineatzero=0 yobj->getProperty, labels = ylbls if obj_valid(ylbls) then begin lobj = ylbls->get(/all) lobj0 = lobj[0] lobj0->setProperty, size=8.0 ; lobj0->setProperty, value='Region_Panel' endif ; to = obj_new('spd_ui_text') ; to->setProperty, value='Region_Panel' ; to->setProperty, size=9 ; ylbls->add, to ; yobj->setProperty, labels=ylbls ; panel[1]->setProperty, yaxis=yobj ; quick_set_panel_labels, panel[1], ['Region_Panel'] panel[1] -> getproperty, tracesettings = obj0 trace_obj = obj0 -> get(/all) ntr = 11 ; 11 is based on the total of the bit mask in thm_roi_bar.pro ; setup colors for ROI plot ;this code will break if spedas color table is not loaded. ;roi_cols = 10+indgen(ntr)*7*245/ntr mod 245 ;tvlct, r,g,b,/get ;ctbl = [[r],[g],[b]] ctbl = transpose([[7,0,5],$ [235,255,0],$ [0,97,255],$ [253,0,0],$ [90,255,0],$ [43,0,232],$ [255,103,0],$ [0,255,133],$ [83,0,117],$ [255,199,0],$ [0,235,254]]) for i = 0,ntr-1 do begin trace_obj[i]->getProperty,linestyle=linestyleObj lineStyleObj->setProperty,thickness=2, color=ctbl[i,*] endfor ;Fix the keogram panel, needs linear Z scaling & labels ;;;;;;;;;;;;;;;;;;;; eric's test stuff ;vars_full = ['thg_pseudoAE', roi_bar, 'Keogram', thx+'_fgs_gse', $ ; esaf_n_name, esaif_v_name, esaf_t_name, 'sample_rate_'+sc, $ ; ssti_name, esaif_flux_name, sste_name, $ ; esaef_flux_name, fbk_tvars] ;;;;;;;;;;;;;;;;;;;; quick_set_panel_labels, panel[2], 'Keogram' quick_set_panel_labels, panel[2], ' ', /zaxis panel[2] -> getproperty, zaxis = zobj if(obj_valid(zobj)) then zobj -> setproperty, scaling = 0 ;Reset labels for FGS data: ; quick_set_panel_labels, panel[3], ['B FIT_GSE(nT)','',''], $ ; quick_set_panel_labels, panel[3], ['Bx', 'By', 'Bz'], $ ; colors_in = [[0b, 0b, 255b], [0b, 255b, 0b], [255b, 0b, 0b]] panel[3]->getProperty, yaxis=yobj yobj->setProperty, stackLabels = 1, orientation = 0 quick_set_panel_labels, panel[3], ['Bx', 'By', 'Bz'], $ colors_in = [2, 4, 6], /zaxis ;Density and Temperature variables need log scaling, and new colors ;and labels panel[4] -> getproperty, tracesettings = obj0 trace_obj = obj0 -> get(/all) If(n_elements(trace_obj) Eq 4) Then Begin;we have an n_pot trace_obj[0] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [0b, 0b, 0b] trace_obj[1] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [0b, 0b, 255b] trace_obj[2] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [255b, 0b, 0b] panel[4] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 ; yobj -> getProperty, labels=labels ;quick_set_panel_labels, panel[4], ['ESA Ni', 'ESA Ne', 'ESA Npot', '[cm^-3]'],$ ; quick_set_panel_labels, panel[4], ['Ni,e_1/cc','', '', ''],$ ; colors_in = [[0b, 0b, 0b], [255b, 0b, 0b], [0b, 0b, 255b], [0b, 0b, 0b]] yobj->setProperty, orientation = 0 yobj->setProperty, stackLabels = 1 quick_set_panel_labels, panel[3], ['Ne', 'Ni', 'Npot'], $ colors_in = [2, 0, 6], /zaxis quick_set_panel_labels, panel[3], ['Ne', 'Ni', 'Npot'],/zaxis Endif Else Begin ;no n_pot trace_obj[0] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [0b, 0b, 0b] trace_obj[1] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [255b, 0b, 0b] panel[4] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 quick_set_panel_labels, panel[4], ['Ni,e_1/cc','', ''],$ colors_in = [[0b, 0b, 0b], [255b, 0b, 0b], [0b, 0b, 0b]] Endelse ;reset velocity labels ;quick_set_panel_labels, panel[5], ['Vx DSL', 'Vy DSL', 'Vz DSL', '[km/s]', 'ESA Ion'], $ ; quick_set_panel_labels, panel[5], ['Vi_km/s', '', '', '', ''], $ ; colors_in = [[0b, 0b, 255b], [0b, 255b, 0b], [255b, 0b, 0b], $ ; [0b, 0b, 0b], [0b, 0b, 0b]] panel[5]->getProperty, yaxis=yobj yobj->setProperty, stackLabels = 1, orientation = 0 quick_set_panel_labels, panel[5], ['Viz', 'Viy', 'Viz'], $ colors_in = [2, 4, 6], /zaxis, /zhorizontal ;t3 needs labels, and logs panel[6] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 yobj->setProperty, stackLabels = 2, orientation = 0 ; quick_set_panel_labels, panel[6], ['Ti_para', 'Te_para', 'Ti_perp', 'Te_perp'], colors_in = [2, 0, 4, 6], /zaxis ; quick_set_panel_labels, panel[6], ['Ti,e_eV', '', $ ; '', '', $ ; ''], $ ;statusbar panel needs different colors panel[7] -> getproperty, tracesettings = obj0 trace_obj = obj0 -> get(/all) ntraces = n_elements(trace_obj) If(ntraces Gt 1) Then Begin If(obj_valid(trace_obj[1])) Then Begin trace_obj[1] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [255b, 0b, 0b] Endif If(obj_valid(trace_obj[2])) Then Begin trace_obj[2] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [255b, 255b, 0b] Endif If(obj_valid(trace_obj[4])) Then Begin trace_obj[4] -> getproperty, linestyle = linestyleobj linestyleobj -> setproperty, color = [0b, 0b, 0b] Endif for i = 0,n_elements(trace_obj)-1 do begin trace_obj[i]->getProperty,linestyle=linestyleObj,symbol=symbolobj lineStyleObj->setProperty,thickness=2 lineStyleObj->getProperty,color=linecolor symbolobj->setProperty,show=0,name=symbolobj->getSymbolName(symbolid=4),id=4,color=linecolor endfor Endif quick_set_panel_labels, panel[7], ['', '','','',''] panel[7]->getProperty,yaxis=yobj ;, tracesettings = trobj yobj->setProperty,annotateAxis=0 ;;make sure SST and ESA panels have same z-range ;panel[8]->getproperty, zaxis=zobj_sst ;zobj_sst->getproperty, minRange=zmin0, maxRange=zmax0 ;panel[9]->getproperty, zaxis=zobj_esa ;zobj_esa->getproperty, minRange=zmin1, maxRange=zmax1 ;zmin=min([zmin0,zmin1]) ;if zmin le 0 then zmin=1. ;zmax=max([zmax0,zmax1]) ;zobj_sst->setproperty, minRange=zmin, maxRange=zmax, fixed=1 ;zobj_esa->setproperty, minRange=zmin, maxRange=zmax, fixed=1 ;panel[8]->setproperty, zaxis=zobj_sst ;panel[9]->setproperty, zaxis=zobj_esa ; ;panel[10]->getproperty, zaxis=zobj_sst ;zobj_sst->getproperty, minRange=zmin0, maxRange=zmax0 ;panel[11]->getproperty, zaxis=zobj_esa ;zobj_esa->getproperty, minRange=zmin1, maxRange=zmax1 ;zmin=min([zmin0,zmin1]) ;if zmin le 0 then zmin=1. ;zmax=max([zmax0,zmax1]) ;zobj_sst->setproperty, minRange=zmin, maxRange=zmax, fixed=1 ;zobj_esa->setproperty, minRange=zmin, maxRange=zmax, fixed=1 ;panel[10]->setproperty, zaxis=zobj_sst ;panel[11]->setproperty, zaxis=zobj_esa ;Keep ESA and SST spectra's z-range consistant panel[8]->getproperty, zaxis = zobj_ssti zobj_ssti->setproperty, minrange = 1d0, maxrange = 5d7, fixed=1 panel[9]->getproperty, zaxis = zobj_esai zobj_esai->setproperty, minrange = 1d3, maxrange = 7.5d8, fixed=1 panel[10]->getproperty, zaxis = zobj_sste zobj_sste->setproperty, minrange = 1d0, maxrange = 5d7, fixed=1 panel[11]->getproperty, zaxis = zobj_esae zobj_esae->setproperty, minrange = 1d4, maxrange = 7.5d8, fixed=1 ;SST needs to have y scaling set to log panel[8] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 yobj -> setproperty, rangeoption = 2 get_data, ssti_name, data = d yobj -> setproperty, maxfixedrange = 3.0e6 If(is_struct(d) && tag_exist(d, 'v')) Then yobj -> setproperty, minfixedrange = min(d.v) $ Else yobj -> setproperty, minfixedrange = 3.0e4 quick_set_panel_labels, panel[8], ['SSTi_Eflux_[eV]'] ; quick_set_panel_labels, panel[8], 'Eflux!CeV/cm!U2!N!C-s-sr-eV', /zaxis, /zhorizontal quick_set_panel_labels, panel[8], 'Eflux, EFU', /zaxis ;ESA needs logs too panel[9] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 quick_set_panel_labels, panel[9], ['ESAi_Eflux_[eV]'] ; quick_set_panel_labels, panel[9], 'Eflux!CeV/cm!U2!N!C-s-sr-eV', /zaxis, /zhorizontal quick_set_panel_labels, panel[9], 'Eflux, EFU', /zaxis ;SST electrons panel[10] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 yobj -> setproperty, rangeoption = 2 get_data, sste_name, data = d yobj -> setproperty, maxfixedrange = 3.0e6 If(is_struct(d) && tag_exist(d, 'v')) Then yobj -> setproperty, minfixedrange = min(d.v) $ Else yobj -> setproperty, minfixedrange = 3.0e4 quick_set_panel_labels, panel[10], ['SSTe_eflux_[eV]'] ; quick_set_panel_labels, panel[10], 'Eflux!CeV/cm!U2!N!C-s-sr-eV', /zaxis, /zhorizontal quick_set_panel_labels, panel[10], 'Eflux, EFU', /zaxis ;ESA electrons panel[11] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 quick_set_panel_labels, panel[11], ['ESAe_eflux_[eV]'] ; quick_set_panel_labels, panel[11], 'Eflux!CeV/cm!U2!N!C-s-sr-eV', /zaxis, /zhorizontal quick_set_panel_labels, panel[11], 'Eflux, EFU', /zaxis ;FBK panels npanels = n_elements(panel) For j = 0, n_elements(fbk_tvars)-1 Do Begin jp = j+12 If(jp Le npanels-1) Then Begin panel[jp] -> getproperty, yaxis = yobj yobj -> setproperty, scaling = 1 get_data, fbk_tvars[j], dlimits = dl, limits = al lbl0 = strupcase(strmid(fbk_tvars[j], 4)) & lbl1 = ' ' If(is_struct(al) && tag_exist(al, 'ztitle')) Then lbl1 = al.ztitle ; quick_set_panel_labels, panel[jp], [lbl0, '[Hz]'] quick_set_panel_labels, panel[jp], lbl1, /zaxis quick_set_panel_labels, panel[jp], thx+'_FBK '+strmid(fbk_tvars[j], 7) Endif quick_set_panel_labels, panel[12], '<|mV/m|>', /zaxis quick_set_panel_labels, panel[13], '<|nT|>', /zaxis Endfor endif ;for variable options, the state data needs to get into the ;loaded_data object: ; but first we need to add suffix copy_data, thx+'_state_pos_gse', thx+'_state_pos_gse'+osuffix If(is_string(tnames(thx+'_state_pos_gse'+osuffix))) Then Begin ok = loadeddata -> add(thx+'_state_pos_gse'+osuffix) if ~keyword_set(no_draw) then begin vo = objarr(3) panel[n_elements(panel)-1] -> getproperty, variables = var_container lbl = 'GSE '+['Z', 'Y', 'X']+' (Re)' vname = thx+'_state_pos_gse'+osuffix+'_'+['z', 'y', 'x'] For j = 0, 2 Do Begin to0 = obj_new('spd_ui_text') to0 -> setproperty, value = lbl[j] to0 -> setproperty, size = 8.0 vo[j] = obj_new('spd_ui_variable') vo[j] -> setproperty, text = to0 vo[j] -> setproperty, controlname = thx+'_state_pos_gse'+osuffix+'_time' vo[j] -> setproperty, fieldname = vname[j] Endfor var_container -> add, vo endif Endif ;WARNING: ; Any code that requires windowStorage or loadedData should be put in the "~keyword_set(no_draw)" block. ; spd_ui_overplot is called without valid values for these parameters during document load if ~keyword_set(no_draw) then begin ;update the draw object and redraw the panels spd_ui_orientation_update,drawObject,windowStorage drawObject -> update, windowStorage, loadedData drawObject -> draw endif SKIP_DAY: spd_ui_overplot_data_clean,tn_before,dont_delete_data=dont_delete_data message, /info, 'Returning:' error=0 Return end