timespan,'1998-11-14',2,/days             ; Set the time span (2 days)

wi_or_load                                ; Load orbit data
wi_mfi_load                               ; Load 3 second mag data (creates variable 'wi_h0_mfi_B3GSE')

xyz_to_polar,'wi_h0_mfi_B3GSE'
store_data,'wi_h0_mfi',data='wi_h0_mfi_B3GSE wi_h0_mfi_B3GSE_mag',dlimits={panel_size:2}
ylim,'wi_h0_mfi',-26,26                    ; Change y- limits
tplot_names                               ; Display all loaded names
tplot,'wi_h0_mfi wi_h0_mfi_B3GSE_phi wi_h0_mfi_B3GSE_th'

dprint,'Use Cursor to select begin and end times. Use right button to end'
tlimit

wi_3dp_load,datatype = 'pm'     ; 3 second proton moments
wi_3dp_load,datatype = 'elpd'  ; Electron reduced values

tplot,'wi_3dp_elpd_FLUX wi_3dp_pm_P_VELS wi_3dp_pm_P_DENS',/add

; Some code examples:

reduce_pads,'wi_3dp_elpd_FLUX',1,3,3   ;  Create PAD spectrogram for energy step 3  (integral from 3 to 3)

reduce_pads,'wi_3dp_elpd_FLUX',2,0,1   ;  Create Energy spectrogram for pitch angles steps 0 to 1  (Sum of PA steps 0 and 1)


; Display the names of data that have been loaded:
tplot_names


; Display the names of subset of data variables (variables that match '*elpd*')
tplot_names,'*elpd*'

;Display info on all the data stored for one variable:
tplot_names,'wi_3dp_elpd_FLUX',/verbose


; There are multiple ways to get the data (within IDL). I will show only one method:
get_data,'wi_3dp_elpd_FLUX',data=dat    ; all data is stored in the structure "dat"


help,/structure,dat

; or alternatively use:
printdat,dat


;Explanation of the "dat" structure:

; DAT.X   : Time ( seconds since 1970 )
; DAT.Y   : flux ( #/cm2/sec/ster/eV)
; DAT.V1  : Energy steps
; DAT.V2  : Pitch angles







end