;+ ; Batch File: THM_CRIB_EFI ; ; Purpose: Demonstrate the loading, calibration, and plotting ; of THEMIS EFI data. ; ; Calling Sequence: ; .run thm_crib_fbk, or using cut-and-paste. ; ; Arguements: ; None. ; ; Notes: ; ***WARNING: Running THM_CRIB_EFI (or likely just calling THM_LOAD_EFI) after THM_CRIB_CLEANEFW will result in the wrong plot labels. ; The problem is probably in the way that some of the LASP code handles, or does not handle, the TPLOT labelling. ; None. ; ; $LastChangedBy: mfeuerstein $ ; $LastChangedDate: 2009-08-26 17:06:41 -0700 (Wed, 26 Aug 2009) $ ; $LastChangedRevision: 6642 $ ; $URL $ ;- ; EFI waveform data example. tplot_title = 'THEMIS EFI Waveform Examples' tplot_options, 'title', tplot_title tplot_options, 'xmargin', [ 15, 10] tplot_options, 'ymargin', [ 5, 5] ; Data parameters, load data, and set color table: ; dur = 1.0 ;days ;tr_string = '2008-05-15' tr_string = '2007-07-20' timespan, tr_string, dur, /day tr = timerange() probe = 'c d' ; converting a to d below. ;; state data needed by thm_load_fgm, which is needed for dot0 quantities. thm_load_state , probe = probe, /get_support thm_load_efi, level=1, suffix='_raw', type='raw', /get_support, probe = probe;, /no_download loadct2, 39 ; the following are optional, since they should be set ; by thm_load_efi as the default plotting options (dlimits). ; you can use these as a starting point if you want to modify the labels/colors. ;options, 'th?_v??_raw', 'colors', [ 1, 2, 3, 4, 5, 6] ;options, 'th?_v??_raw', 'labels', [ 'V1', 'V2', 'V3', 'V4', 'V5', 'V6'] ;options, 'th?_v??_raw', 'labflag', 1 ;options, 'th?_e??_raw', 'colors', [ 2, 4, 6] ;options, 'th?_e??_raw', 'labels', [ 'e12', 'e34', 'e56'] ;options, 'th?_e??_raw', 'labflag', 1 tplot, [ 'th[cde]_vaf_raw', 'th[cde]_eff_raw'] print, 'The plot shows raw data for all probe for which the booms have print, 'been deployed as of 2007-06-30.' print,'Enter ".c" to print the contents of DLIMITS.DATA_ATT structure for' print,'"thd_vaf_raw", and "thd_eff_raw".' stop get_data,'thd_vaf_raw',dlimits=dl help,dl.data_att,/st get_data,'thd_eff_raw',dlimits=dl help,dl.data_att,/st print,'Enter ".c" to zoom in.' stop ;tlimit, [ tr[0] + 10d*3600 + 10d*60 , tr[0] + 13d*3600 + 56d*60 ] ;10:10:00 to 13:56:00 tlimit, [ tr[0] + 17d*3600 + 36d*60 , tr[0] + 17d*3600 + 45d*60 ] ;17:36:00 to 17:45:00 Mozer intervals. ;print, 'zoomed into a period where all satellites have some good data ' print, 'zoomed into a period of interest' print, 'Enter ".c" to continue.' stop thm_cal_efi, /verbose, in_suffix='_raw', onthefly_edc_offset = onthefly_edc_offset, gap_begin = gap_begin, gap_end = gap_end, $ out_suffix='_dsl _gei', coord = 'dsl gei' tplot, [ 'thc_vaf_raw', 'thc_vaf', 'thc_eff_raw', 'thc_eff', 'thc_eff_0', 'thc_eff_dot0_gei' ] print,'Enter ".c" to print the contents of DLIMITS.DATA_ATT structure for' print,'"thc_vaf", and "thc_eff" (note updated and additional fields).' stop get_data,'thd_vaf',dlimits=dl help,dl.data_att,/st get_data,'thd_eff',dlimits=dl help,dl.data_att,/st print, 'Enter ".c" to view the structure and substructure of the optional ONTHFLY_EDC_OFFSET keyword (GAP_BEGIN and GAP_END are similar).' stop help, onthefly_edc_offset, onthefly_edc_offset.d, /st print,'Enter ".c" to view particle burst data.' stop ;tlimit, /full timespan, tr[0], tr[1] - tr[0], /seconds tplot, [ 'thd_vap_raw', 'thd_vap', 'thd_efp_raw', 'thd_efp_dsl', 'thd_efp_gei', 'thd_efp_0', 'thd_efp_dot0_dsl', 'thd_efp_dot0_gei' ] print, 'particle burst' print, 'use the mouse to zoom in' tlimit print, 'Enter ".c" to view wave burst data.' stop ;tplot, [ 'thd_vaw_raw', 'thd_vaw', 'thd_efw_raw', 'thd_efw', 'thd_efw_0', 'thd_efw_dot0' ] tplot, [ 'thd_vaw_raw', 'thd_vaw', 'thd_efw_raw', 'thd_efw_dsl', 'thd_efw_gei', 'thd_efw_0', 'thd_efw_dot0_dsl', 'thd_efw_dot0_gei' ] print, 'wave burst occurs within particle burst' print, 'use the mouse to zoom in' tlimit end