pro training_20180328_beginner
  
; *****************
; settings for the training
; *****************
  erg_init, remote_data_dir='https://ergsc.isee.nagoya-u.ac.jp/data/ergsc_training/nagoya_201803/'
  uname = '***************'
  pass = '****************'
  
; *****************
; set time span
; *****************
  timespan, '2017-03-28', 3, /day
  
; *****************
; load orbit CDF 
; *****************
  erg_load_orb
  tplot_names
  
; *****************
; plot orbit time series data 
; *****************
  tplot, 'erg_orb_l2_pos_gsm'
  
  
; *****************
; plot orbit data in the X-Y plane
; *****************
  tplotxy, 'erg_orb_l2_pos_gsm' 
  
; *****************
; load orbit CDF & set var label
; *****************
  set_erg_var_label
  
; *****************
; Change the time span
; *****************
  timespan, '2017-03-28'
  timespan, '2017-03-28/01:31:41', 10, /min
  
; *****************
; Separate a tplot variable with vector data
; *****************
  split_vec, 'erg_orb_l2_pos_rmlatmlt'
  tplot, 'erg_orb_l2_pos_rmlatmlt_?'
  
  
; *****************
; Insert orbit values below a time-series plot
; *****************
  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
  
; *****************
; load MGF L2 CDF
; *****************
  timespan, '2017-03-28'
  erg_load_mgf, datatype='8sec', uname=uname, pass=pass
  
; *****************
; Change the vertical scale of a plot
; *****************
  ylim, 'erg_mgf_l2_mag_8sec_gse', -300,300,  0
  tplot, 'erg_mgf_l2_mag_8sec_gse'
  
; *****************
; load PWE Provisional CDF and plot
; *****************
  erg_load_pwe_pre, uname=uname, pass=pass
  tplot, [ 'erg_mgf_l2_mag_8sec_gse' , 'erg_pwe_pre_HFA-merged' , $
           'erg_pwe_pre_OFA_E_spectra_132' , 'erg_pwe_pre_OFA_B_spectra_132' , $
           'erg_pwe_pre_EFD_DPB_spectra' ]
  
; *****************
; Change the time range
; *****************
  tlimit
  tlimit, '2017-03-28/00:00' , '2017-03-28/08:00'
  tlimit, /last
  tlimit, /full
  
  
; *****************
; load PWE/OFA L2 CDF and plot
; *****************
  erg_load_pwe_ofa, uname=uname, pass=pass
  tplot, [ 'erg_mgf_l2_mag_8sec_gse' , 'erg_pwe_pre_HFA-merged' , $
           'erg_pwe_ofa_spec_l2_E_spectra_132' , 'erg_pwe_ofa_spec_l2_B_spectra_132' , $
           'erg_pwe_pre_EFD_DPB_spectra' ]
  
; *****************
; checking CDF attributes
; *****************
  show_cdf_att,'erg_pwe_ofa_spec_l2_E_spectra_132'
  
; *****************
; overplot fce
; *****************
; calculate fce
  get_data, 'erg_mgf_l2_magt_8sec', data=data
  fce = data.y/ 10^(9.) * 1.6 * 10^(-19.) / (9.1093D * 10^(-31.)) / 2. / !pi / 1000.
  fce_half = fce / 2.
  
  store_data, 'fce', data={x:data.x, y:fce}, dlim={colors:fsc_color('yellow')}
  store_data, 'fce_half', data={x:data.x, y:fce_half}, dlim={colors:fsc_color('magenta')}
  
  store_data, 'erg_pwe_pre_HFA_merged_gyro', data = ['erg_pwe_pre_HFA_L_spectra_ey','erg_pwe_pre_HFA_H_spectra_ey','fce']
  store_data, 'erg_pwe_ofa_spec_l2_E_gyro', data = ['erg_pwe_ofa_spec_l2_E_spectra_132',  'fce','fce_half']
  store_data, 'erg_pwe_ofa_spec_l2_B_gyro', data = ['erg_pwe_ofa_spec_l2_B_spectra_132',  'fce','fce_half']
  
  ylim, 'erg_pwe_pre_HFA_merged_gyro', 1, 1e4, 1 ; kHz
  ylim, 'erg_pwe_ofa_spec_l2_E_gyro', 0.032, 20, 1 ; kHz
  ylim, 'erg_pwe_ofa_spec_l2_B_gyro', 0.032, 20, 1 ; kHz
  
  tplot, ['erg_pwe_pre_HFA_merged_gyro', 'erg_pwe_ofa_spec_l2_E_gyro','erg_pwe_ofa_spec_l2_B_gyro']
  
; *****************
; Dump to png and postscript files
; *****************
  cwd
  tplot
  makepng, ’erg_pwe_plot'
  popen, ’erg_pwe_plot'
  tplot
  pclose
  
; *****************
; Power spectrum at a particular time
; *****************

  t_index=nn('erg_pwe_ofa_spec_l2_E_spectra_132', '2017-03-28/18:53')
  get_data,'erg_pwe_ofa_spec_l2_E_spectra_132',data=ofa_e_spec
  tf_index=nn('fce_half', '2017-03-28/18:53')
  get_data,'fce_half',data=data_fce_half
  
  plot, ofa_e_spec.v, ofa_e_spec.y[t_index,*],$
        /ylog,/xlog,xstyle=1,xrange=[0.032,20],$
        xtitle='Frequency [kHz]', ytitle='Power [mV^2/m^2/Hz]'
  
  oplot,[data_fce_half.y[tf_index],data_fce_half.y[tf_index]],[1e-8,1],color=fsc_color('green')

end