pro erg_crib_training_beginner_20190306
  
; *****************
; settings for the training (P11)
; *****************
  erg_init
  
; *****************
; set time span (P11)
; *****************
  timespan, '2017-05-27', 07
  
; *****************
; load ERG MGF CDF files (P11) 
; *****************
  erg_load_mgf
  
; *****************
; plot MGF time series data (P11) 
; *****************
  tplot, ['erg_mgf_l2_mag_8sec_dsi','erg_orb_l2_magt_8sec'] 
  
   stop
  
; *****************
; Set a date/time range (P12) 
; *****************  
 timespan, '2017-05-27'
 tplot
 stop
 timespan, '2017-05-29/03:25:30',90,/min
 tplot
 stop
  
; *****************
; Listing tplot data & viewing the content (P14) 
; *****************  
  tplot_names
  print_tinfo,'erg_mgf_l2_mag_8sec_dsi'
  stop
  
; *****************
; Listing tplot data & viewing the content (P15) 
; *****************  
  tplot_names,'erg_mgf_l2_mag_8sec_dsi',/verbose
  stop
  tplot_names,'erg_mgf_l2_mag_8sec_dsi',/v
  stop
  
  
; *****************
; Plotting a tplot data (P16) 
; *****************  
  timespan, '2017-05-27', 07
  tplot,'erg_mgf_l2_mag_8sec_dsi'
  stop
  tplot,['erg_mgf_l2_mag_8sec_dsi','erg_mgf_l2_mag_8sec_gse']
  stop
  tplot,1
  stop
  tplot,[2,1]
  stop
  tplot,'erg_mgf_l2_mag_8sec_*'
  stop

; *****************
; Decorate the plot panel for each tplot variable (P17)
; *****************
  options, 'erg_mgf_l2_mag_8sec_gse', ytitle= 'Arase mag' , ysubtitle='[nT]'
  options, 'erg_mgf_l2_mag_8sec_*', labels=['Bx-comp', 'By-comp', 'Bz-comp' ] 
  tplot, 'erg_mgf_l2_mag_8sec_gse'
  stop
  
; *****************
; Separate a tplot variable with vector data (P18)
; *****************
  split_vec, 'erg_mgf_l2_mag_8sec_gse'
  tplot,'erg_mgf_l2_mag_8sec_gse_?' 
  
; *****************
; Change the time range of a plot (P19)
; *****************
  tlimit
  tlimit, '2017-05-29/03:30','2017-05-29/08:30'
  stop
  tlimit,/last
  stop
  tlimit,/full 
  stop
 
; *****************
; Change the  vertical scale of a plot (P20)
; *****************
tlimit, '2017-05-29/03:30' , '2017-05-29/08:30'
ylim, 'erg_mgf_l2_mag_8sec_gse_x',180,250,0
tplot
stop
 
; *****************
; Dump to png, postscript, and Ascii files (P21)
; *****************
cwd
tplot,'erg_mgf_l2_mag_8sec_gse_?'
makepng, 'erg_mgf_plot'
	  
popen, 'erg_mgf_plot'
tplot 
pclose

tplot_ascii, 'erg_mgf_l2_mag_8sec_gse'

stop


; *****************
; baxcar-average data (P24)
; *****************
del_data,'*'
timespan,'2017-05-29-03:30',70,/min
erg_load_mgf
split_vec,'erg_mgf_l2_mag_8sec_gse'
avg_data,'erg_mgf_l2_mag_8sec_gse_x',600.
tplot,['erg_mgf_l2_mag_8sec_gse_x','erg_mgf_l2_mag_8sec_gse_x_avg']

stop


; *****************
; smoothing  data (P25)
; *****************
tsmooth_in_time,'erg_mgf_l2_mag_8sec_gse_x',600.
tplot_names
tplot,[19,23]
stop


; *****************
; High-pass fileter in time (P26)
; *****************
thigh_pass_filter,'erg_mgf_l2_mag_8sec_gse',600.
tplot_names
tplot,[2,24]
stop


; *****************
; Dynamics spectra in time (P27)
; *****************
timespan,'2017-05-29-03:00',180,/min
tdpwrspc,'erg_mgf_l2_mag_8sec_gse_hpfilt'
tplot_names
tplot,[25,28]
stop

; *****************
; 3-day plot of HEP data and geomag. indices (P32)
; *****************
timespan,'2017-06-16',3,/day
set_erg_var_label
erg_load_hep,datatype='omniflux'
omni_hro_load,/res5min
tplot,['OMNI_HRO_5min_SYM_H','OMNI_HRO_5min_AL_INDEX','erg_hep_l2_FEDO_*']

stop

; *****************
; Coordinate transformation of MGF data (P33)
; *****************
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_sm', in_coord='j2000', out_coord='sm'
tplot, ['mgf_sm', 'erg_orb_l2_pos_blocal']
stop

; *****************
; Frequency analysis of MGF data (P34)
; *****************
timespan, '2017-03-28', 1, /day
tdpwrspc, 'mgf_sm'
wav_data, 'mgf_sm_z'
tplot,'mgf_sm_z*'

stop

; *****************
;Level-2 particle omni-flux data  (P35)
; *****************

;; 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'    ]
stop

; *****************
;Change the time range of a plot  (P36)
; *****************
tlimit
stop
tlimit,'2017-04-19/00:00','2017-04-19/08:00'
stop
tlimit,/last
stop
tlimit,/full
stop

; *****************
;Add MGF data with /add keyword (P37)
; *****************
erg_load_mgf
tplot,'erg_mgf_l2_mag_8sec_sm',/add
stop

; *****************
;Set time range and load ERG orbit data (P38-P39)
; *****************
timespan,'2017-03-28/00:00:00',3,/day
erg_load_orb 
tplot_names
tplot, 'erg_orb_l2_pos_gsm'
tplotxy, 'erg_orb_l2_pos_gsm' 
split_vec, 'erg_orb_l2_pos_rmlatmlt'
tplot, 'erg_orb_l2_pos_rmlatmlt_?'
stop

; *****************
;Separate a tplot variable with vector data (P40)
; *****************
 split_vec, 'erg_orb_l2_pos_rmlatmlt'
 tplot, 'erg_orb_l2_pos_rmlatmlt_?'
stop

; *****************
; Insert orbit values below a time-series plot (P41)
; *****************
set_erg_var_label
tplot 
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


end