;+ ; 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 number of values ; the local mean is away from the local standard deviation. 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 100. ;; 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: aaronbreneman $ ;$LastChangedDate: 2014-04-21 13:19:06 -0700 (Mon, 21 Apr 2014) $ ;$LastChangedRevision: 14901 $ ;$URL: svn+ssh://thmsvn@ambrosia.ssl.berkeley.edu/repos/spdsoft/tags/spedas_4_1/general/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.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. ;the ratio of the local value of the peak FBK value divided by the standard deviation if ~keyword_set(maxrat) then maxrat = 30 ;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 based on keeping peak values within a certain number ;of standard deviations. Avoids having the scaling set by spiky noise ;----------------------------------------------- ;--------- ;FBK1 - 13 ;--------- 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 nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,13) maxv = fltarr(nchunks,13) ;peak values varv = fltarr(nchunks,13) ;variances 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 maxv[i,0] = max(goo0.y[i*sz:i*sz+sz]) maxv[i,1] = max(goo1.y[i*sz:i*sz+sz]) maxv[i,2] = max(goo2.y[i*sz:i*sz+sz]) maxv[i,3] = max(goo3.y[i*sz:i*sz+sz]) maxv[i,4] = max(goo4.y[i*sz:i*sz+sz]) maxv[i,5] = max(goo5.y[i*sz:i*sz+sz]) maxv[i,6] = max(goo6.y[i*sz:i*sz+sz]) maxv[i,7] = max(goo7.y[i*sz:i*sz+sz]) maxv[i,8] = max(goo8.y[i*sz:i*sz+sz]) maxv[i,9] = max(goo9.y[i*sz:i*sz+sz]) maxv[i,10] = max(goo10.y[i*sz:i*sz+sz]) maxv[i,11] = max(goo11.y[i*sz:i*sz+sz]) maxv[i,12] = max(goo12.y[i*sz:i*sz+sz]) tst = moment(goo0.y[i*sz:i*sz+sz]) varv[i,0] = sqrt(tst[1]) tst = moment(goo1.y[i*sz:i*sz+sz]) varv[i,1] = sqrt(tst[1]) tst = moment(goo2.y[i*sz:i*sz+sz]) varv[i,2] = sqrt(tst[1]) tst = moment(goo3.y[i*sz:i*sz+sz]) varv[i,3] = sqrt(tst[1]) tst = moment(goo4.y[i*sz:i*sz+sz]) varv[i,4] = sqrt(tst[1]) tst = moment(goo5.y[i*sz:i*sz+sz]) varv[i,5] = sqrt(tst[1]) tst = moment(goo6.y[i*sz:i*sz+sz]) varv[i,6] = sqrt(tst[1]) tst = moment(goo7.y[i*sz:i*sz+sz]) varv[i,7] = sqrt(tst[1]) tst = moment(goo8.y[i*sz:i*sz+sz]) varv[i,8] = sqrt(tst[1]) tst = moment(goo9.y[i*sz:i*sz+sz]) varv[i,9] = sqrt(tst[1]) tst = moment(goo10.y[i*sz:i*sz+sz]) varv[i,10] = sqrt(tst[1]) tst = moment(goo11.y[i*sz:i*sz+sz]) varv[i,11] = sqrt(tst[1]) tst = moment(goo12.y[i*sz:i*sz+sz]) varv[i,12] = sqrt(tst[1]) 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 stdevv = sqrt(varv) ratio = maxv/stdevv ;q=11 ;!p.multi = [0,0,4] ;plot,maxv[*,q] ;plot,mom[*,q] ;plot,stdevv[*,q] ;plot,maxv[*,q]/stdevv[*,q] ;Reject values (for the y-scaling) where the max value is >> than the standard deviation good = bytarr(nchunks,13) good[*] = 1b goo = where(ratio[*,0] ge maxrat) if goo[0] ne -1 then good[goo,0] = 0b goo = where(ratio[*,1] ge maxrat) if goo[0] ne -1 then good[goo,1] = 0b goo = where(ratio[*,2] ge maxrat) if goo[0] ne -1 then good[goo,2] = 0b goo = where(ratio[*,3] ge maxrat) if goo[0] ne -1 then good[goo,3] = 0b goo = where(ratio[*,4] ge maxrat) if goo[0] ne -1 then good[goo,4] = 0b goo = where(ratio[*,5] ge maxrat) if goo[0] ne -1 then good[goo,5] = 0b goo = where(ratio[*,6] ge maxrat) if goo[0] ne -1 then good[goo,6] = 0b goo = where(ratio[*,7] ge maxrat) if goo[0] ne -1 then good[goo,7] = 0b goo = where(ratio[*,8] ge maxrat) if goo[0] ne -1 then good[goo,8] = 0b goo = where(ratio[*,9] ge maxrat) if goo[0] ne -1 then good[goo,9] = 0b goo = where(ratio[*,10] ge maxrat) if goo[0] ne -1 then good[goo,10] = 0b goo = where(ratio[*,11] ge maxrat) if goo[0] ne -1 then good[goo,11] = 0b goo = where(ratio[*,12] ge maxrat) if goo[0] ne -1 then good[goo,12] = 0b good = float(good) ylim,'rbsp'+probe +'_fbk1_13pk_0',0.,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk1_13pk_1',0.,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk1_13pk_2',0.,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk1_13pk_3',0.,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk1_13pk_4',0.,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk1_13pk_5',0.,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk1_13pk_6',0.,max(maxv[*,6]*good[*,6]),0 ylim,'rbsp'+probe +'_fbk1_13pk_7',0.,max(maxv[*,7]*good[*,7]),0 ylim,'rbsp'+probe +'_fbk1_13pk_8',0.,max(maxv[*,8]*good[*,8]),0 ylim,'rbsp'+probe +'_fbk1_13pk_9',0.,max(maxv[*,9]*good[*,9]),0 ylim,'rbsp'+probe +'_fbk1_13pk_10',0.,max(maxv[*,10]*good[*,10]),0 ylim,'rbsp'+probe +'_fbk1_13pk_11',0.,max(maxv[*,11]*good[*,11]),0 ylim,'rbsp'+probe +'_fbk1_13pk_12',0.,max(maxv[*,12]*good[*,12]),0 ; 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,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk1_13comb_1',0,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk1_13comb_2',0,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk1_13comb_3',0,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk1_13comb_4',0,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk1_13comb_5',0,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk1_13comb_6',0,max(maxv[*,6]*good[*,6]),0 ylim,'rbsp'+probe +'_fbk1_13comb_7',0,max(maxv[*,7]*good[*,7]),0 ylim,'rbsp'+probe +'_fbk1_13comb_8',0,max(maxv[*,8]*good[*,8]),0 ylim,'rbsp'+probe +'_fbk1_13comb_9',0,max(maxv[*,9]*good[*,9]),0 ylim,'rbsp'+probe +'_fbk1_13comb_10',0,max(maxv[*,10]*good[*,10]),0 ylim,'rbsp'+probe +'_fbk1_13comb_11',0,max(maxv[*,11]*good[*,11]),0 ylim,'rbsp'+probe +'_fbk1_13comb_12',0,max(maxv[*,12]*good[*,12]),0 ; 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 ;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 ;--------- ;FBK2 - 13 ;--------- 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 nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,13) maxv = fltarr(nchunks,13) ;peak values varv = fltarr(nchunks,13) ;variances 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 maxv[i,0] = max(goo0.y[i*sz:i*sz+sz]) maxv[i,1] = max(goo1.y[i*sz:i*sz+sz]) maxv[i,2] = max(goo2.y[i*sz:i*sz+sz]) maxv[i,3] = max(goo3.y[i*sz:i*sz+sz]) maxv[i,4] = max(goo4.y[i*sz:i*sz+sz]) maxv[i,5] = max(goo5.y[i*sz:i*sz+sz]) maxv[i,6] = max(goo6.y[i*sz:i*sz+sz]) maxv[i,7] = max(goo7.y[i*sz:i*sz+sz]) maxv[i,8] = max(goo8.y[i*sz:i*sz+sz]) maxv[i,9] = max(goo9.y[i*sz:i*sz+sz]) maxv[i,10] = max(goo10.y[i*sz:i*sz+sz]) maxv[i,11] = max(goo11.y[i*sz:i*sz+sz]) maxv[i,12] = max(goo12.y[i*sz:i*sz+sz]) tst = moment(goo0.y[i*sz:i*sz+sz]) varv[i,0] = sqrt(tst[1]) tst = moment(goo1.y[i*sz:i*sz+sz]) varv[i,1] = sqrt(tst[1]) tst = moment(goo2.y[i*sz:i*sz+sz]) varv[i,2] = sqrt(tst[1]) tst = moment(goo3.y[i*sz:i*sz+sz]) varv[i,3] = sqrt(tst[1]) tst = moment(goo4.y[i*sz:i*sz+sz]) varv[i,4] = sqrt(tst[1]) tst = moment(goo5.y[i*sz:i*sz+sz]) varv[i,5] = sqrt(tst[1]) tst = moment(goo6.y[i*sz:i*sz+sz]) varv[i,6] = sqrt(tst[1]) tst = moment(goo7.y[i*sz:i*sz+sz]) varv[i,7] = sqrt(tst[1]) tst = moment(goo8.y[i*sz:i*sz+sz]) varv[i,8] = sqrt(tst[1]) tst = moment(goo9.y[i*sz:i*sz+sz]) varv[i,9] = sqrt(tst[1]) tst = moment(goo10.y[i*sz:i*sz+sz]) varv[i,10] = sqrt(tst[1]) tst = moment(goo11.y[i*sz:i*sz+sz]) varv[i,11] = sqrt(tst[1]) tst = moment(goo12.y[i*sz:i*sz+sz]) varv[i,12] = sqrt(tst[1]) 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 stdevv = sqrt(varv) ratio = maxv/stdevv ;Reject values (for the y-scaling) where the max value is >> than the standard deviation good = bytarr(nchunks,13) good[*] = 1b goo = where(ratio[*,0] ge maxrat) if goo[0] ne -1 then good[goo,0] = 0b goo = where(ratio[*,1] ge maxrat) if goo[0] ne -1 then good[goo,1] = 0b goo = where(ratio[*,2] ge maxrat) if goo[0] ne -1 then good[goo,2] = 0b goo = where(ratio[*,3] ge maxrat) if goo[0] ne -1 then good[goo,3] = 0b goo = where(ratio[*,4] ge maxrat) if goo[0] ne -1 then good[goo,4] = 0b goo = where(ratio[*,5] ge maxrat) if goo[0] ne -1 then good[goo,5] = 0b goo = where(ratio[*,6] ge maxrat) if goo[0] ne -1 then good[goo,6] = 0b goo = where(ratio[*,7] ge maxrat) if goo[0] ne -1 then good[goo,7] = 0b goo = where(ratio[*,8] ge maxrat) if goo[0] ne -1 then good[goo,8] = 0b goo = where(ratio[*,9] ge maxrat) if goo[0] ne -1 then good[goo,9] = 0b goo = where(ratio[*,10] ge maxrat) if goo[0] ne -1 then good[goo,10] = 0b goo = where(ratio[*,11] ge maxrat) if goo[0] ne -1 then good[goo,11] = 0b goo = where(ratio[*,12] ge maxrat) if goo[0] ne -1 then good[goo,12] = 0b good = float(good) ylim,'rbsp'+probe +'_fbk2_13pk_0',0.,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk2_13pk_1',0.,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk2_13pk_2',0.,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk2_13pk_3',0.,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk2_13pk_4',0.,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk2_13pk_5',0.,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk2_13pk_6',0.,max(maxv[*,6]*good[*,6]),0 ylim,'rbsp'+probe +'_fbk2_13pk_7',0.,max(maxv[*,7]*good[*,7]),0 ylim,'rbsp'+probe +'_fbk2_13pk_8',0.,max(maxv[*,8]*good[*,8]),0 ylim,'rbsp'+probe +'_fbk2_13pk_9',0.,max(maxv[*,9]*good[*,9]),0 ylim,'rbsp'+probe +'_fbk2_13pk_10',0.,max(maxv[*,10]*good[*,10]),0 ylim,'rbsp'+probe +'_fbk2_13pk_11',0.,max(maxv[*,11]*good[*,11]),0 ylim,'rbsp'+probe +'_fbk2_13pk_12',0.,max(maxv[*,12]*good[*,12]),0 if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk2_13comb_0',0,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk2_13comb_1',0,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk2_13comb_2',0,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk2_13comb_3',0,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk2_13comb_4',0,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk2_13comb_5',0,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk2_13comb_6',0,max(maxv[*,6]*good[*,6]),0 ylim,'rbsp'+probe +'_fbk2_13comb_7',0,max(maxv[*,7]*good[*,7]),0 ylim,'rbsp'+probe +'_fbk2_13comb_8',0,max(maxv[*,8]*good[*,8]),0 ylim,'rbsp'+probe +'_fbk2_13comb_9',0,max(maxv[*,9]*good[*,9]),0 ylim,'rbsp'+probe +'_fbk2_13comb_10',0,max(maxv[*,10]*good[*,10]),0 ylim,'rbsp'+probe +'_fbk2_13comb_11',0,max(maxv[*,11]*good[*,11]),0 ylim,'rbsp'+probe +'_fbk2_13comb_12',0,max(maxv[*,12]*good[*,12]),0 endif dmessage='ylim fbk2_13pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;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 ;--------- ;FBK1 - 7 ;--------- 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 nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,7) maxv = fltarr(nchunks,7) ;peak values varv = fltarr(nchunks,7) ;variances 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 maxv[i,0] = max(goo0.y[i*sz:i*sz+sz]) maxv[i,1] = max(goo1.y[i*sz:i*sz+sz]) maxv[i,2] = max(goo2.y[i*sz:i*sz+sz]) maxv[i,3] = max(goo3.y[i*sz:i*sz+sz]) maxv[i,4] = max(goo4.y[i*sz:i*sz+sz]) maxv[i,5] = max(goo5.y[i*sz:i*sz+sz]) maxv[i,6] = max(goo6.y[i*sz:i*sz+sz]) tst = moment(goo0.y[i*sz:i*sz+sz]) varv[i,0] = sqrt(tst[1]) tst = moment(goo1.y[i*sz:i*sz+sz]) varv[i,1] = sqrt(tst[1]) tst = moment(goo2.y[i*sz:i*sz+sz]) varv[i,2] = sqrt(tst[1]) tst = moment(goo3.y[i*sz:i*sz+sz]) varv[i,3] = sqrt(tst[1]) tst = moment(goo4.y[i*sz:i*sz+sz]) varv[i,4] = sqrt(tst[1]) tst = moment(goo5.y[i*sz:i*sz+sz]) varv[i,5] = sqrt(tst[1]) tst = moment(goo6.y[i*sz:i*sz+sz]) varv[i,6] = sqrt(tst[1]) 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 stdevv = sqrt(varv) ratio = maxv/stdevv ;Reject values (for the y-scaling) where the max value is >> than the standard deviation good = bytarr(nchunks,7) good[*] = 1b goo = where(ratio[*,0] ge maxrat) if goo[0] ne -1 then good[goo,0] = 0b goo = where(ratio[*,1] ge maxrat) if goo[0] ne -1 then good[goo,1] = 0b goo = where(ratio[*,2] ge maxrat) if goo[0] ne -1 then good[goo,2] = 0b goo = where(ratio[*,3] ge maxrat) if goo[0] ne -1 then good[goo,3] = 0b goo = where(ratio[*,4] ge maxrat) if goo[0] ne -1 then good[goo,4] = 0b goo = where(ratio[*,5] ge maxrat) if goo[0] ne -1 then good[goo,5] = 0b goo = where(ratio[*,6] ge maxrat) if goo[0] ne -1 then good[goo,6] = 0b good = float(good) ylim,'rbsp'+probe +'_fbk1_7pk_0',0.,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk1_7pk_1',0.,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk1_7pk_2',0.,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk1_7pk_3',0.,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk1_7pk_4',0.,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk1_7pk_5',0.,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk1_7pk_6',0.,max(maxv[*,6]*good[*,6]),0 if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk1_7comb_0',0,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk1_7comb_1',0,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk1_7comb_2',0,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk1_7comb_3',0,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk1_7comb_4',0,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk1_7comb_5',0,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk1_7comb_6',0,max(maxv[*,6]*good[*,6]),0 endif dmessage='ylim fbk1_7pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;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 ;--------- ;FBK2 - 7 ;--------- 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 nchunks = floor(n_elements(goo0.y)/sz) mom = fltarr(nchunks,7) maxv = fltarr(nchunks,7) ;peak values varv = fltarr(nchunks,7) ;variances 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 maxv[i,0] = max(goo0.y[i*sz:i*sz+sz]) maxv[i,1] = max(goo1.y[i*sz:i*sz+sz]) maxv[i,2] = max(goo2.y[i*sz:i*sz+sz]) maxv[i,3] = max(goo3.y[i*sz:i*sz+sz]) maxv[i,4] = max(goo4.y[i*sz:i*sz+sz]) maxv[i,5] = max(goo5.y[i*sz:i*sz+sz]) maxv[i,6] = max(goo6.y[i*sz:i*sz+sz]) tst = moment(goo0.y[i*sz:i*sz+sz]) varv[i,0] = sqrt(tst[1]) tst = moment(goo1.y[i*sz:i*sz+sz]) varv[i,1] = sqrt(tst[1]) tst = moment(goo2.y[i*sz:i*sz+sz]) varv[i,2] = sqrt(tst[1]) tst = moment(goo3.y[i*sz:i*sz+sz]) varv[i,3] = sqrt(tst[1]) tst = moment(goo4.y[i*sz:i*sz+sz]) varv[i,4] = sqrt(tst[1]) tst = moment(goo5.y[i*sz:i*sz+sz]) varv[i,5] = sqrt(tst[1]) tst = moment(goo6.y[i*sz:i*sz+sz]) varv[i,6] = sqrt(tst[1]) 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 stdevv = sqrt(varv) ratio = maxv/stdevv ;Reject values (for the y-scaling) where the max value is >> than the standard deviation good = bytarr(nchunks,7) good[*] = 1b goo = where(ratio[*,0] ge maxrat) if goo[0] ne -1 then good[goo,0] = 0b goo = where(ratio[*,1] ge maxrat) if goo[0] ne -1 then good[goo,1] = 0b goo = where(ratio[*,2] ge maxrat) if goo[0] ne -1 then good[goo,2] = 0b goo = where(ratio[*,3] ge maxrat) if goo[0] ne -1 then good[goo,3] = 0b goo = where(ratio[*,4] ge maxrat) if goo[0] ne -1 then good[goo,4] = 0b goo = where(ratio[*,5] ge maxrat) if goo[0] ne -1 then good[goo,5] = 0b goo = where(ratio[*,6] ge maxrat) if goo[0] ne -1 then good[goo,6] = 0b good = float(good) ylim,'rbsp'+probe +'_fbk2_7pk_0',0.,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk2_7pk_1',0.,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk2_7pk_2',0.,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk2_7pk_3',0.,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk2_7pk_4',0.,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk2_7pk_5',0.,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk2_7pk_6',0.,max(maxv[*,6]*good[*,6]),0 if keyword_set(combine) then begin ylim,'rbsp'+probe +'_fbk2_7comb_0',0,max(maxv[*,0]*good[*,0]),0 ylim,'rbsp'+probe +'_fbk2_7comb_1',0,max(maxv[*,1]*good[*,1]),0 ylim,'rbsp'+probe +'_fbk2_7comb_2',0,max(maxv[*,2]*good[*,2]),0 ylim,'rbsp'+probe +'_fbk2_7comb_3',0,max(maxv[*,3]*good[*,3]),0 ylim,'rbsp'+probe +'_fbk2_7comb_4',0,max(maxv[*,4]*good[*,4]),0 ylim,'rbsp'+probe +'_fbk2_7comb_5',0,max(maxv[*,5]*good[*,5]),0 ylim,'rbsp'+probe +'_fbk2_7comb_6',0,max(maxv[*,6]*good[*,6]),0 endif dmessage='ylim fbk2_7pk runtime (s): '+string(systime(1)-ylim_start_time) dprint,verbose=verbose,dlevel=2,dmessage ;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