;+ ; NAME: rbsp_split_fbk ; ; ; PURPOSE: Split the filterbank data into separate tplot variables for ; each channel. Can also combine peak and average on a single panel ; ; KEYWORDS: ; probe = 'a' or 'b' NOTE: single spacecraft only, does not accept ['a b'] ; combine -> set to combine peak and average onto a single plot ; meansz -> number of data points over which to calculate the mean. This is used ; to set y-scaling. If not done then the yscale of the FBK plots ; is often dominated by single large amplitude spikes, masking the majority ; of the data. Default is 20. ; ysc -> Scale factor for y-scaling. Default is 1. This should be set to ; values greater than 1 for meansz larger than the default. For ; example, ysc=3. works well for meansz=1000. ; ; CREATED: Aaron Breneman 11/07/2012 ; ; MODIFIED: changed y-scaling based on a running average rather than the max value ; for each FBK channel. This avoids having a ridiculous yscaling based on a few ; very large amplitude spiky events. ; ; VERSION: ;$LastChangedBy: kersten $ ;$LastChangedDate: 2013-11-06 08:06:42 -0800 (Wed, 06 Nov 2013) $ ;$LastChangedRevision: 13498 $ ;$URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/ssl_general/trunk/missions/rbsp/efw/rbsp_split_fbk.pro $ ; ;- pro rbsp_split_fbk,probe,combine=combine,meansz=sz,verbose=verbose,ysc=ysc vb = keyword_set(verbose) ? verbose : 0 vb = vb > !rbsp_efw.verbose start_time=systime(1) ; fbk13_bins=['0.8-1.5Hz', '1.5-3Hz', '3-6Hz', '6-12Hz', '12-25Hz', '25-50Hz', $ ; '50-100Hz', '100-200Hz', '200-400Hz', '400-800Hz', $ ; '0.8-1.6kHz', '1.6-3.2kHz', '3.2-6.5kHz'] fbk13_bins=['0.8-1.5', '1.5-3', '3-6', '6-12', '12-25', '25-50', $ '50-100', '100-200', '200-400', '400-800', $ '0.8-1.6', '1.6-3.2', '3.2-6.5'] fbk7_bins=fbk13_bins[lindgen(7)*2] ;Determine ylim based on the mean value. Doing this avoids large y-scalings ;based on very spiky FBK data if ~keyword_set(sz) then sz = 100. if ~keyword_set(ysc) then ysc=1. ;grab dlimits structure for each of the channels to determine the source get_data,'rbsp'+probe+'_efw_fbk_7_fb1_pk',dlimits=fb7_fb1 get_data,'rbsp'+probe+'_efw_fbk_7_fb2_pk',dlimits=fb7_fb2 get_data,'rbsp'+probe+'_efw_fbk_13_fb1_pk',dlimits=fb13_fb1 get_data,'rbsp'+probe+'_efw_fbk_13_fb2_pk',dlimits=fb13_fb2 get_data,'rbsp'+probe +'_efw_fbk_13_fb1_pk',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk1_13pk_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk1_13pk_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk1_13pk_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk1_13pk_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk1_13pk_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk1_13pk_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk1_13pk_6',data={x:goo.x,y:goo.y[*,6]} store_data,'rbsp'+probe +'_fbk1_13pk_7',data={x:goo.x,y:goo.y[*,7]} store_data,'rbsp'+probe +'_fbk1_13pk_8',data={x:goo.x,y:goo.y[*,8]} store_data,'rbsp'+probe +'_fbk1_13pk_9',data={x:goo.x,y:goo.y[*,9]} store_data,'rbsp'+probe +'_fbk1_13pk_10',data={x:goo.x,y:goo.y[*,10]} store_data,'rbsp'+probe +'_fbk1_13pk_11',data={x:goo.x,y:goo.y[*,11]} store_data,'rbsp'+probe +'_fbk1_13pk_12',data={x:goo.x,y:goo.y[*,12]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_13pk_0','ytitle','pk!C'+fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_1','ytitle','pk!C'+fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_2','ytitle','pk!C'+fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_3','ytitle','pk!C'+fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_4','ytitle','pk!C'+fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_5','ytitle','pk!C'+fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_6','ytitle','pk!C'+fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_7','ytitle','pk!C'+fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_8','ytitle','pk!C'+fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_9','ytitle','pk!C'+fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_10','ytitle','pk!C'+fbk13_bins[10]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_11','ytitle','pk!C'+fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13pk_12','ytitle','pk!C'+fbk13_bins[12]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13??_*','labels','RBSP'+probe+'!C '+fb13_fb1.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_13_fb2_pk',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk2_13pk_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk2_13pk_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk2_13pk_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk2_13pk_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk2_13pk_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk2_13pk_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk2_13pk_6',data={x:goo.x,y:goo.y[*,6]} store_data,'rbsp'+probe +'_fbk2_13pk_7',data={x:goo.x,y:goo.y[*,7]} store_data,'rbsp'+probe +'_fbk2_13pk_8',data={x:goo.x,y:goo.y[*,8]} store_data,'rbsp'+probe +'_fbk2_13pk_9',data={x:goo.x,y:goo.y[*,9]} store_data,'rbsp'+probe +'_fbk2_13pk_10',data={x:goo.x,y:goo.y[*,10]} store_data,'rbsp'+probe +'_fbk2_13pk_11',data={x:goo.x,y:goo.y[*,11]} store_data,'rbsp'+probe +'_fbk2_13pk_12',data={x:goo.x,y:goo.y[*,12]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_13pk_0','ytitle','pk!C'+fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_1','ytitle','pk!C'+fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_2','ytitle','pk!C'+fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_3','ytitle','pk!C'+fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_4','ytitle','pk!C'+fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_5','ytitle','pk!C'+fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_6','ytitle','pk!C'+fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_7','ytitle','pk!C'+fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_8','ytitle','pk!C'+fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_9','ytitle','pk!C'+fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_10','ytitle','pk!C'+fbk13_bins[10]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_11','ytitle','pk!C'+fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13pk_12','ytitle','pk!C'+fbk13_bins[12]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13??_*','labels','RBSP'+probe+'!C '+fb13_fb2.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_7_fb1_pk',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk1_7pk_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk1_7pk_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk1_7pk_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk1_7pk_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk1_7pk_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk1_7pk_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk1_7pk_6',data={x:goo.x,y:goo.y[*,6]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_7pk_0','ytitle','pk!C'+fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_1','ytitle','pk!C'+fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_2','ytitle','pk!C'+fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_3','ytitle','pk!C'+fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_4','ytitle','pk!C'+fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_5','ytitle','pk!C'+fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7pk_6','ytitle','pk!C'+fbk7_bins[6]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7??_*','labels','RBSP'+probe+'!C '+fb7_fb1.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_7_fb2_pk',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk2_7pk_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk2_7pk_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk2_7pk_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk2_7pk_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk2_7pk_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk2_7pk_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk2_7pk_6',data={x:goo.x,y:goo.y[*,6]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_7pk_0','ytitle','pk!C'+fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_1','ytitle','pk!C'+fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_2','ytitle','pk!C'+fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_3','ytitle','pk!C'+fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_4','ytitle','pk!C'+fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_5','ytitle','pk!C'+fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7pk_6','ytitle','pk!C'+fbk7_bins[6]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7??_*','labels','RBSP'+probe+'!C '+fb7_fb2.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_13_fb1_av',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk1_13av_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk1_13av_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk1_13av_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk1_13av_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk1_13av_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk1_13av_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk1_13av_6',data={x:goo.x,y:goo.y[*,6]} store_data,'rbsp'+probe +'_fbk1_13av_7',data={x:goo.x,y:goo.y[*,7]} store_data,'rbsp'+probe +'_fbk1_13av_8',data={x:goo.x,y:goo.y[*,8]} store_data,'rbsp'+probe +'_fbk1_13av_9',data={x:goo.x,y:goo.y[*,9]} store_data,'rbsp'+probe +'_fbk1_13av_10',data={x:goo.x,y:goo.y[*,10]} store_data,'rbsp'+probe +'_fbk1_13av_11',data={x:goo.x,y:goo.y[*,11]} store_data,'rbsp'+probe +'_fbk1_13av_12',data={x:goo.x,y:goo.y[*,12]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_13av_0','ytitle','av!C'+fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_1','ytitle','av!C'+fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_2','ytitle','av!C'+fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_3','ytitle','av!C'+fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_4','ytitle','av!C'+fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_5','ytitle','av!C'+fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_6','ytitle','av!C'+fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_7','ytitle','av!C'+fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_8','ytitle','av!C'+fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_9','ytitle','av!C'+fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_10','ytitle','av!C'+fbk13_bins[10]+'!CHkz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_11','ytitle','av!C'+fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13av_12','ytitle','av!C'+fbk13_bins[12]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13??_*','labels','RBSP'+probe+'!C '+fb13_fb1.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_13_fb2_av',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk2_13av_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk2_13av_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk2_13av_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk2_13av_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk2_13av_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk2_13av_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk2_13av_6',data={x:goo.x,y:goo.y[*,6]} store_data,'rbsp'+probe +'_fbk2_13av_7',data={x:goo.x,y:goo.y[*,7]} store_data,'rbsp'+probe +'_fbk2_13av_8',data={x:goo.x,y:goo.y[*,8]} store_data,'rbsp'+probe +'_fbk2_13av_9',data={x:goo.x,y:goo.y[*,9]} store_data,'rbsp'+probe +'_fbk2_13av_10',data={x:goo.x,y:goo.y[*,10]} store_data,'rbsp'+probe +'_fbk2_13av_11',data={x:goo.x,y:goo.y[*,11]} store_data,'rbsp'+probe +'_fbk2_13av_12',data={x:goo.x,y:goo.y[*,12]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_13av_0','ytitle','av!C'+fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_1','ytitle','av!C'+fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_2','ytitle','av!C'+fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_3','ytitle','av!C'+fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_4','ytitle','av!C'+fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_5','ytitle','av!C'+fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_6','ytitle','av!C'+fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_7','ytitle','av!C'+fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_8','ytitle','av!C'+fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_9','ytitle','av!C'+fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_10','ytitle','av!C'+fbk13_bins[10]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_11','ytitle','av!C'+fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13av_12','ytitle','av!C'+fbk13_bins[12]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13??_*','labels','RBSP'+probe+'!C '+fb13_fb2.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_7_fb1_av',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk1_7av_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk1_7av_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk1_7av_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk1_7av_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk1_7av_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk1_7av_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk1_7av_6',data={x:goo.x,y:goo.y[*,6]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_7av_0','ytitle','av!C'+fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_1','ytitle','av!C'+fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_2','ytitle','av!C'+fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_3','ytitle','av!C'+fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_4','ytitle','av!C'+fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_5','ytitle','av!C'+fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7av_6','ytitle','av!C'+fbk7_bins[6]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7??_*','labels','RBSP'+probe+'!C '+fb7_fb1.data_att.channel endif get_data,'rbsp'+probe +'_efw_fbk_7_fb2_av',data=goo,dlimits=dlim tst = size(goo,/dimensions) if tst ne 0 then begin store_data,'rbsp'+probe +'_fbk2_7av_0',data={x:goo.x,y:goo.y[*,0]} store_data,'rbsp'+probe +'_fbk2_7av_1',data={x:goo.x,y:goo.y[*,1]} store_data,'rbsp'+probe +'_fbk2_7av_2',data={x:goo.x,y:goo.y[*,2]} store_data,'rbsp'+probe +'_fbk2_7av_3',data={x:goo.x,y:goo.y[*,3]} store_data,'rbsp'+probe +'_fbk2_7av_4',data={x:goo.x,y:goo.y[*,4]} store_data,'rbsp'+probe +'_fbk2_7av_5',data={x:goo.x,y:goo.y[*,5]} store_data,'rbsp'+probe +'_fbk2_7av_6',data={x:goo.x,y:goo.y[*,6]} units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_7av_0','ytitle','av!C'+fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_1','ytitle','av!C'+fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_2','ytitle','av!C'+fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_3','ytitle','av!C'+fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_4','ytitle','av!C'+fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_5','ytitle','av!C'+fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7av_6','ytitle','av!C'+fbk7_bins[6]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7??_*','labels','RBSP'+probe+'!C '+fb7_fb2.data_att.channel endif if keyword_set(combine) then begin ;Combine FBK pk and av quantities onto a single plot tst = tnames('rbsp'+probe +'_fbk1_13pk_0') if tst ne '' then begin store_data,'rbsp'+probe +'_fbk1_13comb_0',data=['rbsp'+probe +'_fbk1_13pk_0','rbsp'+probe +'_fbk1_13av_0'] store_data,'rbsp'+probe +'_fbk1_13comb_1',data=['rbsp'+probe +'_fbk1_13pk_1','rbsp'+probe +'_fbk1_13av_1'] store_data,'rbsp'+probe +'_fbk1_13comb_2',data=['rbsp'+probe +'_fbk1_13pk_2','rbsp'+probe +'_fbk1_13av_2'] store_data,'rbsp'+probe +'_fbk1_13comb_3',data=['rbsp'+probe +'_fbk1_13pk_3','rbsp'+probe +'_fbk1_13av_3'] store_data,'rbsp'+probe +'_fbk1_13comb_4',data=['rbsp'+probe +'_fbk1_13pk_4','rbsp'+probe +'_fbk1_13av_4'] store_data,'rbsp'+probe +'_fbk1_13comb_5',data=['rbsp'+probe +'_fbk1_13pk_5','rbsp'+probe +'_fbk1_13av_5'] store_data,'rbsp'+probe +'_fbk1_13comb_6',data=['rbsp'+probe +'_fbk1_13pk_6','rbsp'+probe +'_fbk1_13av_6'] store_data,'rbsp'+probe +'_fbk1_13comb_7',data=['rbsp'+probe +'_fbk1_13pk_7','rbsp'+probe +'_fbk1_13av_7'] store_data,'rbsp'+probe +'_fbk1_13comb_8',data=['rbsp'+probe +'_fbk1_13pk_8','rbsp'+probe +'_fbk1_13av_8'] store_data,'rbsp'+probe +'_fbk1_13comb_9',data=['rbsp'+probe +'_fbk1_13pk_9','rbsp'+probe +'_fbk1_13av_9'] store_data,'rbsp'+probe +'_fbk1_13comb_10',data=['rbsp'+probe +'_fbk1_13pk_10','rbsp'+probe +'_fbk1_13av_10'] store_data,'rbsp'+probe +'_fbk1_13comb_11',data=['rbsp'+probe +'_fbk1_13pk_11','rbsp'+probe +'_fbk1_13av_11'] store_data,'rbsp'+probe +'_fbk1_13comb_12',data=['rbsp'+probe +'_fbk1_13pk_12','rbsp'+probe +'_fbk1_13av_12'] options,'rbsp'+probe +'_fbk1_13comb_*','labels','RBSP'+probe+'!C '+fb13_fb1.data_att.channel endif tst = tnames('rbsp'+probe +'_fbk2_13pk_0') if tst ne '' then begin store_data,'rbsp'+probe +'_fbk2_13comb_0',data=['rbsp'+probe +'_fbk2_13pk_0','rbsp'+probe +'_fbk2_13av_0'] store_data,'rbsp'+probe +'_fbk2_13comb_1',data=['rbsp'+probe +'_fbk2_13pk_1','rbsp'+probe +'_fbk2_13av_1'] store_data,'rbsp'+probe +'_fbk2_13comb_2',data=['rbsp'+probe +'_fbk2_13pk_2','rbsp'+probe +'_fbk2_13av_2'] store_data,'rbsp'+probe +'_fbk2_13comb_3',data=['rbsp'+probe +'_fbk2_13pk_3','rbsp'+probe +'_fbk2_13av_3'] store_data,'rbsp'+probe +'_fbk2_13comb_4',data=['rbsp'+probe +'_fbk2_13pk_4','rbsp'+probe +'_fbk2_13av_4'] store_data,'rbsp'+probe +'_fbk2_13comb_5',data=['rbsp'+probe +'_fbk2_13pk_5','rbsp'+probe +'_fbk2_13av_5'] store_data,'rbsp'+probe +'_fbk2_13comb_6',data=['rbsp'+probe +'_fbk2_13pk_6','rbsp'+probe +'_fbk2_13av_6'] store_data,'rbsp'+probe +'_fbk2_13comb_7',data=['rbsp'+probe +'_fbk2_13pk_7','rbsp'+probe +'_fbk2_13av_7'] store_data,'rbsp'+probe +'_fbk2_13comb_8',data=['rbsp'+probe +'_fbk2_13pk_8','rbsp'+probe +'_fbk2_13av_8'] store_data,'rbsp'+probe +'_fbk2_13comb_9',data=['rbsp'+probe +'_fbk2_13pk_9','rbsp'+probe +'_fbk2_13av_9'] store_data,'rbsp'+probe +'_fbk2_13comb_10',data=['rbsp'+probe +'_fbk2_13pk_10','rbsp'+probe +'_fbk2_13av_10'] store_data,'rbsp'+probe +'_fbk2_13comb_11',data=['rbsp'+probe +'_fbk2_13pk_11','rbsp'+probe +'_fbk2_13av_11'] store_data,'rbsp'+probe +'_fbk2_13comb_12',data=['rbsp'+probe +'_fbk2_13pk_12','rbsp'+probe +'_fbk2_13av_12'] options,'rbsp'+probe +'_fbk2_13comb_*','labels','RBSP'+probe+'!C '+fb13_fb2.data_att.channel endif tst = tnames('rbsp'+probe +'_fbk1_7pk_0') if tst ne '' then begin store_data,'rbsp'+probe +'_fbk1_7comb_0',data=['rbsp'+probe +'_fbk1_7pk_0','rbsp'+probe +'_fbk1_7av_0'] store_data,'rbsp'+probe +'_fbk1_7comb_1',data=['rbsp'+probe +'_fbk1_7pk_1','rbsp'+probe +'_fbk1_7av_1'] store_data,'rbsp'+probe +'_fbk1_7comb_2',data=['rbsp'+probe +'_fbk1_7pk_2','rbsp'+probe +'_fbk1_7av_2'] store_data,'rbsp'+probe +'_fbk1_7comb_3',data=['rbsp'+probe +'_fbk1_7pk_3','rbsp'+probe +'_fbk1_7av_3'] store_data,'rbsp'+probe +'_fbk1_7comb_4',data=['rbsp'+probe +'_fbk1_7pk_4','rbsp'+probe +'_fbk1_7av_4'] store_data,'rbsp'+probe +'_fbk1_7comb_5',data=['rbsp'+probe +'_fbk1_7pk_5','rbsp'+probe +'_fbk1_7av_5'] store_data,'rbsp'+probe +'_fbk1_7comb_6',data=['rbsp'+probe +'_fbk1_7pk_6','rbsp'+probe +'_fbk1_7av_6'] options,'rbsp'+probe +'_fbk1_7comb_*','labels','RBSP'+probe+'!C '+fb7_fb1.data_att.channel endif tst = tnames('rbsp'+probe +'_fbk2_7pk_0') if tst ne '' then begin store_data,'rbsp'+probe +'_fbk2_7comb_0',data=['rbsp'+probe +'_fbk2_7pk_0','rbsp'+probe +'_fbk2_7av_0'] store_data,'rbsp'+probe +'_fbk2_7comb_1',data=['rbsp'+probe +'_fbk2_7pk_1','rbsp'+probe +'_fbk2_7av_1'] store_data,'rbsp'+probe +'_fbk2_7comb_2',data=['rbsp'+probe +'_fbk2_7pk_2','rbsp'+probe +'_fbk2_7av_2'] store_data,'rbsp'+probe +'_fbk2_7comb_3',data=['rbsp'+probe +'_fbk2_7pk_3','rbsp'+probe +'_fbk2_7av_3'] store_data,'rbsp'+probe +'_fbk2_7comb_4',data=['rbsp'+probe +'_fbk2_7pk_4','rbsp'+probe +'_fbk2_7av_4'] store_data,'rbsp'+probe +'_fbk2_7comb_5',data=['rbsp'+probe +'_fbk2_7pk_5','rbsp'+probe +'_fbk2_7av_5'] store_data,'rbsp'+probe +'_fbk2_7comb_6',data=['rbsp'+probe +'_fbk2_7pk_6','rbsp'+probe +'_fbk2_7av_6'] options,'rbsp'+probe +'_fbk2_7comb_*','labels','RBSP'+probe+'!C '+fb7_fb2.data_att.channel endif ;Change colors to black and red tst = tnames('rbsp'+probe +'_fbk1_13comb_0') if tst ne '' then begin options,['rbsp'+probe+'_fbk1_13comb_0',$ 'rbsp'+probe+'_fbk1_13comb_1',$ 'rbsp'+probe+'_fbk1_13comb_2',$ 'rbsp'+probe+'_fbk1_13comb_3',$ 'rbsp'+probe+'_fbk1_13comb_4',$ 'rbsp'+probe+'_fbk1_13comb_5',$ 'rbsp'+probe+'_fbk1_13comb_6',$ 'rbsp'+probe+'_fbk1_13comb_7',$ 'rbsp'+probe+'_fbk1_13comb_8',$ 'rbsp'+probe+'_fbk1_13comb_9',$ 'rbsp'+probe+'_fbk1_13comb_10',$ 'rbsp'+probe+'_fbk1_13comb_11',$ 'rbsp'+probe+'_fbk1_13comb_12'],colors=[0,6] endif tst = tnames('rbsp'+probe +'_fbk2_13comb_0') if tst ne '' then begin options,['rbsp'+probe+'_fbk2_13comb_0',$ 'rbsp'+probe+'_fbk2_13comb_1',$ 'rbsp'+probe+'_fbk2_13comb_2',$ 'rbsp'+probe+'_fbk2_13comb_3',$ 'rbsp'+probe+'_fbk2_13comb_4',$ 'rbsp'+probe+'_fbk2_13comb_5',$ 'rbsp'+probe+'_fbk2_13comb_6',$ 'rbsp'+probe+'_fbk2_13comb_7',$ 'rbsp'+probe+'_fbk2_13comb_8',$ 'rbsp'+probe+'_fbk2_13comb_9',$ 'rbsp'+probe+'_fbk2_13comb_10',$ 'rbsp'+probe+'_fbk2_13comb_11',$ 'rbsp'+probe+'_fbk2_13comb_12'],colors=[0,6] endif tst = tnames('rbsp'+probe +'_fbk1_7comb_0') if tst ne '' then begin options,['rbsp'+probe+'_fbk1_7comb_0',$ 'rbsp'+probe+'_fbk1_7comb_1',$ 'rbsp'+probe+'_fbk1_7comb_2',$ 'rbsp'+probe+'_fbk1_7comb_3',$ 'rbsp'+probe+'_fbk1_7comb_4',$ 'rbsp'+probe+'_fbk1_7comb_5',$ 'rbsp'+probe+'_fbk1_7comb_6'],colors=[0,6] endif tst = tnames('rbsp'+probe +'_fbk2_7comb_0') if tst ne '' then begin options,['rbsp'+probe+'_fbk2_7comb_0',$ 'rbsp'+probe+'_fbk2_7comb_1',$ 'rbsp'+probe+'_fbk2_7comb_2',$ 'rbsp'+probe+'_fbk2_7comb_3',$ 'rbsp'+probe+'_fbk2_7comb_4',$ 'rbsp'+probe+'_fbk2_7comb_5',$ 'rbsp'+probe+'_fbk2_7comb_6'],colors=[0,6] endif tst = tnames('rbsp'+probe +'_fbk1_13comb_0') get_data,'rbsp'+probe+'_efw_fbk_13_fb1_pk',dlimits=dlim,data=goo if tst ne '' then begin units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_13comb_0','ytitle',fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_1','ytitle',fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_2','ytitle',fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_3','ytitle',fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_4','ytitle',fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_5','ytitle',fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_6','ytitle',fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_7','ytitle',fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_8','ytitle',fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_9','ytitle',fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_10','ytitle',fbk13_bins[10]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_11','ytitle',fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_13comb_12','ytitle',fbk13_bins[12]+'!CkHz!C['+units+']' endif tst = tnames('rbsp'+probe +'_fbk2_13comb_0') get_data,'rbsp'+probe+'_efw_fbk_13_fb2_pk',dlimits=dlim,data=goo if tst ne '' then begin units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_13comb_0','ytitle',fbk13_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_1','ytitle',fbk13_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_2','ytitle',fbk13_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_3','ytitle',fbk13_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_4','ytitle',fbk13_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_5','ytitle',fbk13_bins[5]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_6','ytitle',fbk13_bins[6]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_7','ytitle',fbk13_bins[7]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_8','ytitle',fbk13_bins[8]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_9','ytitle',fbk13_bins[9]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_10','ytitle',fbk13_bins[10]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_11','ytitle',fbk13_bins[11]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_13comb_12','ytitle',fbk13_bins[12]+'!CkHz!C['+units+']' endif tst = tnames('rbsp'+probe +'_fbk1_7comb_0') get_data,'rbsp'+probe+'_efw_fbk_7_fb1_pk',dlimits=dlim,data=goo if tst ne '' then begin units = dlim.data_att.units options,'rbsp'+probe +'_fbk1_7comb_0','ytitle',fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_1','ytitle',fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_2','ytitle',fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_3','ytitle',fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_4','ytitle',fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_5','ytitle',fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk1_7comb_6','ytitle',fbk7_bins[6]+'!CkHz!C['+units+']' endif tst = tnames('rbsp'+probe +'_fbk2_7comb_0') get_data,'rbsp'+probe+'_efw_fbk_7_fb2_pk',dlimits=dlim,data=goo if tst ne '' then begin units = dlim.data_att.units options,'rbsp'+probe +'_fbk2_7comb_0','ytitle',fbk7_bins[0]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_1','ytitle',fbk7_bins[1]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_2','ytitle',fbk7_bins[2]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_3','ytitle',fbk7_bins[3]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_4','ytitle',fbk7_bins[4]+'!CHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_5','ytitle',fbk7_bins[5]+'!CkHz!C['+units+']' options,'rbsp'+probe +'_fbk2_7comb_6','ytitle',fbk7_bins[6]+'!CkHz!C['+units+']' endif endif ;The number of y-ticks is usually too large when 7 or 13 plots are displayed at once. ;Unfortunately reducing the tick numbers means that negative numbers show up. So, I have ;to manually set the yrange options,'rbsp'+probe+'_fbk?_13pk_*','yticks',2 options,'rbsp'+probe+'_fbk?_7pk_*','yticks',2 ;If yticks are manually set, yrange isn't set or is set to a range of zero (max=min), ;ystyle is set to 1, and the data is a horizontal line (max(y)=min(y)), then the IDL ;PLOT routine crashes. This next bit checks to see if any of the pk vars are horizontal ;lines, and if so, sets yticks to 0 (automatic) for that variable. ;-------------- pk13vars=tnames('rbsp'+probe+'_fbk1_13pk_*',cnt) if cnt gt 0 then begin for i=0, n_elements(pk13vars)-1 do begin var=pk13vars[i] get_data,var,data=vardat if min(vardat.y) eq max(vardat.y) then $ options,var,'yticks',0 endfor endif pk13vars=tnames('rbsp'+probe+'_fbk2_13pk_*',cnt) if cnt gt 0 then begin for i=0, n_elements(pk13vars)-1 do begin var=pk13vars[i] get_data,var,data=vardat if min(vardat.y) eq max(vardat.y) then $ options,var,'yticks',0 endfor endif pk7vars=tnames('rbsp'+probe+'_fbk1_7pk_*',cnt) if cnt gt 0 then begin for i=0, n_elements(pk7vars)-1 do begin var=pk7vars[i] get_data,var,data=vardat if min(vardat.y) eq max(vardat.y) then $ options,var,'yticks',0 endfor endif pk7vars=tnames('rbsp'+probe+'_fbk2_7pk_*',cnt) if cnt gt 0 then begin for i=0, n_elements(pk7vars)-1 do begin var=pk7vars[i] get_data,var,data=vardat if min(vardat.y) eq max(vardat.y) then $ options,var,'yticks',0 endfor endif ;----------------------------------------------- ;Set ylimit and tick format ;----------------------------------------------- tst = tnames('rbsp'+probe +'_fbk1_13pk_0') if tst ne '' then begin get_data,'rbsp'+probe +'_fbk1_13pk_0',data=goo0 get_data,'rbsp'+probe +'_fbk1_13pk_1',data=goo1 get_data,'rbsp'+probe +'_fbk1_13pk_2',data=goo2 get_data,'rbsp'+probe +'_fbk1_13pk_3',data=goo3 get_data,'rbsp'+probe +'_fbk1_13pk_4',data=goo4 get_data,'rbsp'+probe +'_fbk1_13pk_5',data=goo5 get_data,'rbsp'+probe +'_fbk1_13pk_6',data=goo6 get_data,'rbsp'+probe +'_fbk1_13pk_7',data=goo7 get_data,'rbsp'+probe +'_fbk1_13pk_8',data=goo8 get_data,'rbsp'+probe +'_fbk1_13pk_9',data=goo9 get_data,'rbsp'+probe +'_fbk1_13pk_10',data=goo10 get_data,'rbsp'+probe +'_fbk1_13pk_11',data=goo11 get_data,'rbsp'+probe +'_fbk1_13pk_12',data=goo12 ;Determine ylim based on the mean value. Doing this avoids large y-scalings ;based on very spiky FBK data nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,13) ylim_start_time=systime(1) dmessage="using nchunks in fbk1_13pk: "+string(nchunks) dprint,verbose=verbose,dlevel=2,dmessage for i=0L,nchunks-1*sz do begin mom[i,0]=total(goo0.y[i*sz:i*sz+sz])/sz mom[i,1]=total(goo1.y[i*sz:i*sz+sz])/sz mom[i,2]=total(goo2.y[i*sz:i*sz+sz])/sz mom[i,3]=total(goo3.y[i*sz:i*sz+sz])/sz mom[i,4]=total(goo4.y[i*sz:i*sz+sz])/sz mom[i,5]=total(goo5.y[i*sz:i*sz+sz])/sz mom[i,6]=total(goo6.y[i*sz:i*sz+sz])/sz mom[i,7]=total(goo7.y[i*sz:i*sz+sz])/sz mom[i,8]=total(goo8.y[i*sz:i*sz+sz])/sz mom[i,9]=total(goo9.y[i*sz:i*sz+sz])/sz mom[i,10]=total(goo10.y[i*sz:i*sz+sz])/sz mom[i,11]=total(goo11.y[i*sz:i*sz+sz])/sz mom[i,12]=total(goo12.y[i*sz:i*sz+sz])/sz endfor ylim,'rbsp'+probe +'_fbk1_13pk_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk1_13pk_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk1_13pk_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk1_13pk_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk1_13pk_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk1_13pk_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk1_13pk_6',0,ysc*max(mom[*,6]) ylim,'rbsp'+probe +'_fbk1_13pk_7',0,ysc*max(mom[*,7]) ylim,'rbsp'+probe +'_fbk1_13pk_8',0,ysc*max(mom[*,8]) ylim,'rbsp'+probe +'_fbk1_13pk_9',0,ysc*max(mom[*,9]) ylim,'rbsp'+probe +'_fbk1_13pk_10',0,ysc*max(mom[*,10]) ylim,'rbsp'+probe +'_fbk1_13pk_11',0,ysc*max(mom[*,11]) ylim,'rbsp'+probe +'_fbk1_13pk_12',0,ysc*max(mom[*,12]) if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk1_13comb_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk1_13comb_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk1_13comb_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk1_13comb_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk1_13comb_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk1_13comb_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk1_13comb_6',0,ysc*max(mom[*,6]) ylim,'rbsp'+probe +'_fbk1_13comb_7',0,ysc*max(mom[*,7]) ylim,'rbsp'+probe +'_fbk1_13comb_8',0,ysc*max(mom[*,8]) ylim,'rbsp'+probe +'_fbk1_13comb_9',0,ysc*max(mom[*,9]) ylim,'rbsp'+probe +'_fbk1_13comb_10',0,ysc*max(mom[*,10]) ylim,'rbsp'+probe +'_fbk1_13comb_11',0,ysc*max(mom[*,11]) ylim,'rbsp'+probe +'_fbk1_13comb_12',0,ysc*max(mom[*,12]) endif dmessage='ylim fbk1_13pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;Set ylim based on maximum value ; ylim,'rbsp'+probe +'_fbk1_13pk_0',0,max(goo0.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_1',0,max(goo1.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_2',0,max(goo2.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_3',0,max(goo3.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_4',0,max(goo4.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_5',0,max(goo5.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_6',0,max(goo6.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_7',0,max(goo7.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_8',0,max(goo8.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_9',0,max(goo9.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_10',0,max(goo10.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_11',0,max(goo11.y) ; ylim,'rbsp'+probe +'_fbk1_13pk_12',0,max(goo12.y) ;Keep only 1 number after decimal place. More just clutters up plot options,'rbsp'+probe +'_fbk1_13pk_0','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_1','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_2','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_3','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_4','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_5','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_6','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_7','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_8','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_9','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_10','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_11','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_13pk_12','ytickformat','(f6.1)' endif tst = tnames('rbsp'+probe +'_fbk2_13pk_0') if tst ne '' then begin get_data,'rbsp'+probe +'_fbk2_13pk_0',data=goo0 get_data,'rbsp'+probe +'_fbk2_13pk_1',data=goo1 get_data,'rbsp'+probe +'_fbk2_13pk_2',data=goo2 get_data,'rbsp'+probe +'_fbk2_13pk_3',data=goo3 get_data,'rbsp'+probe +'_fbk2_13pk_4',data=goo4 get_data,'rbsp'+probe +'_fbk2_13pk_5',data=goo5 get_data,'rbsp'+probe +'_fbk2_13pk_6',data=goo6 get_data,'rbsp'+probe +'_fbk2_13pk_7',data=goo7 get_data,'rbsp'+probe +'_fbk2_13pk_8',data=goo8 get_data,'rbsp'+probe +'_fbk2_13pk_9',data=goo9 get_data,'rbsp'+probe +'_fbk2_13pk_10',data=goo10 get_data,'rbsp'+probe +'_fbk2_13pk_11',data=goo11 get_data,'rbsp'+probe +'_fbk2_13pk_12',data=goo12 ;Determine ylim based on the mean value. Doing this avoids large y-scalings ;based on very spiky FBK data nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,13) ylim_start_time=systime(1) dmessage="using nchunks in fbk2_13pk: "+string(nchunks) dprint,verbose=verbose,dlevel=2,dmessage for i=0L,nchunks-1*sz do begin mom[i,0]=total(goo0.y[i*sz:i*sz+sz])/sz mom[i,1]=total(goo1.y[i*sz:i*sz+sz])/sz mom[i,2]=total(goo2.y[i*sz:i*sz+sz])/sz mom[i,3]=total(goo3.y[i*sz:i*sz+sz])/sz mom[i,4]=total(goo4.y[i*sz:i*sz+sz])/sz mom[i,5]=total(goo5.y[i*sz:i*sz+sz])/sz mom[i,6]=total(goo6.y[i*sz:i*sz+sz])/sz mom[i,7]=total(goo7.y[i*sz:i*sz+sz])/sz mom[i,8]=total(goo8.y[i*sz:i*sz+sz])/sz mom[i,9]=total(goo9.y[i*sz:i*sz+sz])/sz mom[i,10]=total(goo10.y[i*sz:i*sz+sz])/sz mom[i,11]=total(goo11.y[i*sz:i*sz+sz])/sz mom[i,12]=total(goo12.y[i*sz:i*sz+sz])/sz endfor ylim,'rbsp'+probe +'_fbk2_13pk_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk2_13pk_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk2_13pk_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk2_13pk_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk2_13pk_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk2_13pk_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk2_13pk_6',0,ysc*max(mom[*,6]) ylim,'rbsp'+probe +'_fbk2_13pk_7',0,ysc*max(mom[*,7]) ylim,'rbsp'+probe +'_fbk2_13pk_8',0,ysc*max(mom[*,8]) ylim,'rbsp'+probe +'_fbk2_13pk_9',0,ysc*max(mom[*,9]) ylim,'rbsp'+probe +'_fbk2_13pk_10',0,ysc*max(mom[*,10]) ylim,'rbsp'+probe +'_fbk2_13pk_11',0,ysc*max(mom[*,11]) ylim,'rbsp'+probe +'_fbk2_13pk_12',0,ysc*max(mom[*,12]) if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk2_13comb_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk2_13comb_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk2_13comb_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk2_13comb_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk2_13comb_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk2_13comb_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk2_13comb_6',0,ysc*max(mom[*,6]) ylim,'rbsp'+probe +'_fbk2_13comb_7',0,ysc*max(mom[*,7]) ylim,'rbsp'+probe +'_fbk2_13comb_8',0,ysc*max(mom[*,8]) ylim,'rbsp'+probe +'_fbk2_13comb_9',0,ysc*max(mom[*,9]) ylim,'rbsp'+probe +'_fbk2_13comb_10',0,ysc*max(mom[*,10]) ylim,'rbsp'+probe +'_fbk2_13comb_11',0,ysc*max(mom[*,11]) ylim,'rbsp'+probe +'_fbk2_13comb_12',0,ysc*max(mom[*,12]) endif dmessage='ylim fbk2_13pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;Set ylim based on maximum value ; ylim,'rbsp'+probe +'_fbk2_13pk_0',0,max(goo0.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_1',0,max(goo1.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_2',0,max(goo2.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_3',0,max(goo3.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_4',0,max(goo4.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_5',0,max(goo5.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_6',0,max(goo6.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_7',0,max(goo7.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_8',0,max(goo8.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_9',0,max(goo9.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_10',0,max(goo10.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_11',0,max(goo11.y) ; ylim,'rbsp'+probe +'_fbk2_13pk_12',0,max(goo12.y) ;Keep only 1 number after decimal place. More just clutters up plot options,'rbsp'+probe +'_fbk2_13pk_0','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_1','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_2','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_3','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_4','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_5','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_6','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_7','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_8','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_9','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_10','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_11','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_13pk_12','ytickformat','(f6.1)' endif tst = tnames('rbsp'+probe +'_fbk1_7pk_0') if tst ne '' then begin get_data,'rbsp'+probe +'_fbk1_7pk_0',data=goo0 get_data,'rbsp'+probe +'_fbk1_7pk_1',data=goo1 get_data,'rbsp'+probe +'_fbk1_7pk_2',data=goo2 get_data,'rbsp'+probe +'_fbk1_7pk_3',data=goo3 get_data,'rbsp'+probe +'_fbk1_7pk_4',data=goo4 get_data,'rbsp'+probe +'_fbk1_7pk_5',data=goo5 get_data,'rbsp'+probe +'_fbk1_7pk_6',data=goo6 ;Determine ylim based on the mean value. Doing this avoids large y-scalings ;based on very spiky FBK data nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,7) ylim_start_time=systime(1) dmessage="using nchunks in fbk1_7pk: "+string(nchunks) dprint,verbose=verbose,dlevel=2,dmessage for i=0L,nchunks-1*sz do begin mom[i,0]=total(goo0.y[i*sz:i*sz+sz])/sz mom[i,1]=total(goo1.y[i*sz:i*sz+sz])/sz mom[i,2]=total(goo2.y[i*sz:i*sz+sz])/sz mom[i,3]=total(goo3.y[i*sz:i*sz+sz])/sz mom[i,4]=total(goo4.y[i*sz:i*sz+sz])/sz mom[i,5]=total(goo5.y[i*sz:i*sz+sz])/sz mom[i,6]=total(goo6.y[i*sz:i*sz+sz])/sz endfor ylim,'rbsp'+probe +'_fbk1_7pk_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk1_7pk_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk1_7pk_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk1_7pk_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk1_7pk_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk1_7pk_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk1_7pk_6',0,ysc*max(mom[*,6]) if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk1_7comb_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk1_7comb_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk1_7comb_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk1_7comb_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk1_7comb_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk1_7comb_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk1_7comb_6',0,ysc*max(mom[*,6]) endif dmessage='ylim fbk1_7pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;set ylim based on maximum value ; ylim,'rbsp'+probe +'_fbk1_7pk_0',0,max(goo0.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_1',0,max(goo1.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_2',0,max(goo2.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_3',0,max(goo3.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_4',0,max(goo4.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_5',0,max(goo5.y) ; ylim,'rbsp'+probe +'_fbk1_7pk_6',0,max(goo6.y) ;Keep only 1 number after decimal place. More just clutters up plot options,'rbsp'+probe +'_fbk1_7pk_0','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_1','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_2','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_3','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_4','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_5','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk1_7pk_6','ytickformat','(f6.1)' endif tst = tnames('rbsp'+probe +'_fbk2_7pk_0') if tst ne '' then begin get_data,'rbsp'+probe +'_fbk2_7pk_0',data=goo0 get_data,'rbsp'+probe +'_fbk2_7pk_1',data=goo1 get_data,'rbsp'+probe +'_fbk2_7pk_2',data=goo2 get_data,'rbsp'+probe +'_fbk2_7pk_3',data=goo3 get_data,'rbsp'+probe +'_fbk2_7pk_4',data=goo4 get_data,'rbsp'+probe +'_fbk2_7pk_5',data=goo5 get_data,'rbsp'+probe +'_fbk2_7pk_6',data=goo6 ;Determine ylim based on the mean value. Doing this avoids large y-scalings ;based on very spiky FBK data nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,7) ylim_start_time=systime(1) dmessage="using nchunks in fbk2_7pk: "+string(nchunks) dprint,verbose=verbose,dlevel=2,dmessage for i=0L,nchunks-1*sz do begin mom[i,0]=total(goo0.y[i*sz:i*sz+sz])/sz mom[i,1]=total(goo1.y[i*sz:i*sz+sz])/sz mom[i,2]=total(goo2.y[i*sz:i*sz+sz])/sz mom[i,3]=total(goo3.y[i*sz:i*sz+sz])/sz mom[i,4]=total(goo4.y[i*sz:i*sz+sz])/sz mom[i,5]=total(goo5.y[i*sz:i*sz+sz])/sz mom[i,6]=total(goo6.y[i*sz:i*sz+sz])/sz endfor ylim,'rbsp'+probe +'_fbk2_7pk_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk2_7pk_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk2_7pk_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk2_7pk_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk2_7pk_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk2_7pk_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk2_7pk_6',0,ysc*max(mom[*,6]) if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk2_7comb_0',0,ysc*max(mom[*,0]) ylim,'rbsp'+probe +'_fbk2_7comb_1',0,ysc*max(mom[*,1]) ylim,'rbsp'+probe +'_fbk2_7comb_2',0,ysc*max(mom[*,2]) ylim,'rbsp'+probe +'_fbk2_7comb_3',0,ysc*max(mom[*,3]) ylim,'rbsp'+probe +'_fbk2_7comb_4',0,ysc*max(mom[*,4]) ylim,'rbsp'+probe +'_fbk2_7comb_5',0,ysc*max(mom[*,5]) ylim,'rbsp'+probe +'_fbk2_7comb_6',0,ysc*max(mom[*,6]) endif dmessage='ylim fbk2_7pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;Set ylim based on maximum value ; ylim,'rbsp'+probe +'_fbk2_7pk_0',0,max(goo0.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_1',0,max(goo1.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_2',0,max(goo2.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_3',0,max(goo3.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_4',0,max(goo4.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_5',0,max(goo5.y) ; ylim,'rbsp'+probe +'_fbk2_7pk_6',0,max(goo6.y) ;Keep only 1 number after decimal place. More just clutters up plot options,'rbsp'+probe +'_fbk2_7pk_0','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_1','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_2','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_3','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_4','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_5','ytickformat','(f6.1)' options,'rbsp'+probe +'_fbk2_7pk_6','ytickformat','(f6.1)' endif dmessage='runtime (s): '+string(systime(1)-start_time) dprint,verbose=verbose,dlevel=2,dmessage end