;+
;
; erg_crib_bgn_training_201809.pro
;
; The crib sheet used for the beginner's cource of the SPEDAS training
; held in Tohoku Univ. on September 20, 2018.
;
;
;-

;;; Slide #11

idl
erg_init 
timespan, '2017-06-16', 3
thm_load_gmag, site='atha pokr'

tplot, [ 'thg_mag_atha', 'thg_mag_pokr' ] 


;;; Slide #12

timespan, '2017-06-16'
timespan, '2017-06-17/01:31:41', 10, /min


;;; Slide #14

tplot_names
print_tinfo, 'thg_mag_pokr'


;;; Slide #15
tplot_names, 'thg_mag_atha', /verbose
tplot_names, 1, /v


;;; Slide #16

tplot, 'thg_mag_atha'
tplot, [ 'thg_mag_atha' , 'thg_mag_pokr' ]

tplot, 1
tplot, [2, 1]

tplot, 'thg_mag_*'


;;; Slide #17

options, 'thg_mag_atha', ytitle='ATHA geomag' , ysubtitle='[nT]'
options, 'thg_mag_????', labels=[ 'H-comp', 'D-comp', 'Z-comp' ] 
tplot, 'thg_mag_atha'


;;; SLide #18

split_vec, 'thg_mag_atha'
tplot, 'thg_mag_atha_?'


;;; Slide #19

tlimit
tlimit, '2017-06-16/00:00', '2017-06-17/00:00'
tlimit, /last
tlimit, /full


;;; Slide #20

ylim, 'thg_mag_atha_x', 13950,14050,  0
tplot


;;; Slide #21

cwd         ;Display the current directory

tplot, 'thg_mag_atha_?'
makepng, 'atha_plot'

popen, 'atha_plot'
tplot         ;Redo the last plot
pclose

tplot_ascii, 'thg_mag_atha'


;;; Slide #24

del_data, '*'
timespan, '2017-06-16'
thm_load_gmag, site='atha' 
tsub_average, 'thg_mag_atha' 
tplot, ['thg_mag_atha', 'thg_mag_atha-d' ]


;;; Slide #25

avg_data, 'thg_mag_atha-d', 1000.
tplot, 'thg_mag_atha-d*'


;;; Slide #26

tsmooth_in_time, 'thg_mag_atha-d', 1000. 
tplot_names

tplot, [ 2, 4 ]


;;; Slide #27

thigh_pass_filter, 'thg_mag_atha-d', 1000.
tplot_names

tplot, [ 2, 5 ]


;;; Slide #29

tdpwrspc, 'thg_mag_atha-d' 
tplot_names

tplot, [ 5, 9 ]


;;; Slide #30

split_vec, 'thg_mag_atha-d'
avg_data, 'thg_mag_atha-d_x', 5.
wav_data, 'thg_mag_atha-d_x_avg'

tplot, 'thg_mag_atha-d_x_avg*' 

zlim, 13, 100, 1000000, 1 
tplot, 'thg_mag_atha-d_x_avg*'


;;; Slide #33

timespan, '2017-06-16', 3, /day

erg_load_hep, datatype='omniflux'

omni_hro_load, /res5min
tplot,[ 'OMNI_HRO_5min_SYM_H' , 'OMNI_HRO_5min_AL_INDEX' , 'erg_hep_l2_FEDO_*' ]


;;; Slide #34

timespan, '2017-03-28' & erg_load_orb & erg_load_mgf
erg_cotrans, 'erg_mgf_l2_mag_8sec_dsi', 'mgf_j2000', in_coord='dsi', out_coord='j2000'

spd_cotrans, 'mgf_j2000', 'mgf_gsm', in_coord='j2000', out_coord='gsm'
tplot, ['mgf_gsm', 'erg_orb_l2_pos_blocal']


;;; Slide #35

timespan, '2017-03-28', 1, /day
tdpwrspc, 'mgf_sm'
wav_data, 'mgf_sm_z' 
tplot,'mgf_sm_z*'


;;; Slide #36

;; Set the time span
timespan, '2017-04-19'
;; Load data
erg_load_xep, datatype='omniflux'
erg_load_hep, datatype='omniflux'
erg_load_mepe, datatype='omniflux'
erg_load_lepe, datatype='omniflux'
erg_load_mepi_nml, datatype='omniflux'
erg_load_lepi_nml, datatype='omniflux'
tplot_names

tplot,[ 'erg_xep_l2_FEDO_SSD', $
  'erg_hep_l2_FEDO_H','erg_hep_l2_FEDO_L', $
  'erg_mepe_l2_omniflux_FEDO', $
  'erg_lepe_l2_omniflux_FEDO', $
  'erg_mepi_l2_omniflux_FPDO', $
  'erg_lepi_l2_omniflux_FPDO'    ]


;;; Slide #37

tlimit
tlimit, '2017-04-19/00:00', '2017-04-19/08:00'
tlimit, /last
tlimit, /full


;;; SLide #38

erg_load_mgf, uname=uname, pass=pass
tplot, 'erg_mgf_l2_mag_8sec_sm' , /add


;;; Slide #39

timespan,'2017-03-28/00:00:00',3,/day 

erg_load_orb 
tplot_names


;;; Slide #40

tplot, 'erg_orb_l2_pos_gsm'

tplotxy, 'erg_orb_l2_pos_gsm' 


;;; Slide #41

split_vec, 'erg_orb_l2_pos_rmlatmlt'
tplot, 'erg_orb_l2_pos_rmlatmlt_?'


;;; Slide #42

set_erg_var_label
tplot           ;just type "tplot" to replot the previous panels

split_vec, 'erg_orb_l2_pos_rmlatmlt'
options, 'erg_orb_l2_pos_rmlatmlt_x', ytitle='R'
options, 'erg_orb_l2_pos_rmlatmlt_y', ytitle='MLAT'
options, 'erg_orb_l2_pos_rmlatmlt_z', ytitle='MLT'
tplot_options, var_label=['erg_orb_l2_pos_rmlatmlt_z', 'erg_orb_l2_pos_rmlatmlt_y', 'erg_orb_l2_pos_rmlatmlt_x']
tplot


;;; Slide #44

get_data, 'thg_mag_atha', data=data, dlimits=dlimits, lim=lim

help, data.x , data.y


;;; Slide #45

timestr = '2017-06-16/12:30:00'
spedastime = time_double( timestr )
print, spedastime

print, time_string( spedastime )