;Add fce, fce/2, flh and fci lines to spectral plots ;Requirements: ; Can load some EFW spectral data with ; rbsp_load_efw_spec.pro or rbsp_load_efw_xspec.pro ;speclist -> list of tplot names of spectral data ;sc -> Van Allen probe 'a' or 'b' ;period -> Downsample the EMFISIS data to this time period. ; Defaults to 11 seconds (~ Van Allen Probe spin period) ;Written by: ;Aaron W Breneman, Dec 2012 ;Modified to use EMFISIS L3 data 2013-01-25 AWB pro rbsp_add_fce2spec,speclist,sc,period=period if ~keyword_set(period) then period = 11. ;sc rotation period ;First find out if EMFISIS data is loaded ;If not, then load it ;get_data,'rbsp'+sc+'_emfisis_quicklook_Magnitude',data=mag get_data,'rbsp'+sc+'_emfisis_l3_4sec_gse_Magnitude',data=mag if ~is_struct(mag) then begin ; rbsp_load_emfisis,probe=sc,/quicklook ; get_data,'rbsp'+sc+'_emfisis_quicklook_Magnitude',data=mag rbsp_load_emfisis,probe=sc,coord='gse',cadence='4sec',level='l3' get_data,'rbsp'+sc+'_emfisis_l3_4sec_gse_Magnitude',data=mag endif if keyword_set(period) then begin ; rbsp_downsample,['rbsp'+sc+'_emfisis_quicklook_Magnitude'],1/period,suffix='_DS_tmp' ; get_data,'rbsp'+sc+'_emfisis_quicklook_Magnitude_DS_tmp',data=mag rbsp_downsample,'rbsp'+sc+'_emfisis_l3_4sec_gse_Magnitude',1/period,suffix='_DS_tmp' get_data,'rbsp'+sc+'_emfisis_l3_4sec_gse_Magnitude_DS_tmp',data=mag endif ;Loop through each spec variable for i=0,n_elements(speclist)-1 do begin get_data,speclist[i],data=dat,dlimits=dlim,limits=lim ;Check to see if tplot variable is spectral data if is_struct(dat) then begin if dlim.spec eq 1 then begin print,'Adding fce lines to ' + speclist[i] fce = 28.*mag.y fce = interpol(fce,mag.x,dat.x) store_data,'fce',data={x:dat.x,y:fce} store_data,'fce_2',data={x:dat.x,y:fce/2.} store_data,'fci',data={x:dat.x,y:fce/1836.} store_data,'flh',data={x:dat.x,y:sqrt(fce*fce/1836.)} store_data,speclist[i]+'_fce',data=[speclist[i],'fce','fce_2','fci','flh'] ylim,speclist[i]+'_fce',3,10000,1 ; if keyword_set(nochange) then begin ; ; get_data,speclist[i]+'_fce',data=dat ; store_data,speclist[i],data=dat ; ylim,speclist[i],3,10000,1 ; ; ; tplot,[speclist[i]+'_fce',speclist[i]] ; get_data,speclist[i],data=dd1 ; get_data,speclist[i]+'_fce',data=dd2 ; ; endif endif endif endfor end