;; ID/password of the ERG working group ;;uname = '?????' ;;pass = '?????' ;; Slide #9 timespan, '2017-05-01' erg_load_lepe_pa, uname=uname, pass=pass tplot, 'erg_lepe_l3_prov_PA_' + ['energy_15', 'pitchangle_09' ] ;; Slide #12 timespan, '2017-04-12' erg_load_lepe, datatype='3dflux', uname=uname, pass=pass, /no_sort_enebin erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy' tplot, 'erg_lepe_l2_3dflux_FEDU_energy' ;; Slide #13 timespan, '2017-04-12/04:00', 8, /hour & get_timespan, tr erg_load_mgf & set_erg_var_label magvn = 'erg_mgf_l2_mag_8sec_dsi' & posvn = 'erg_orb_l2_pos_gse' erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='pa', $ energy=[1000., 1500.], trange=tr, mag=magvn, pos=posvn tplot, 'erg_lepe_l2_3dflux_FEDU_*' ;; Slide #14 erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy', $ trange=tr, mag=magvn, pos=posvn, pitch=[80, 100] tplot, 'erg_lepe_l2_3dflux_FEDU_*' ;; Slide #15 erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='moments', mag=magvn, pos=posvn, trange=tr tplot, 'erg_lepe_l2_3dflux_FEDU_' + ['energy','pa','ptens', 'density'] ;; From here to bottom, uname and pass should be the ERG working group ID/pass below ;; Slide #17 ;; For Figure 1 of Kazama+2018 timespan, '2017-04-12/16:00', 2, /hour get_timespan, tr erg_load_lepe, datatype='3dflux', uname=uname, pass=pass, /no_sort_enebin erg_load_mgf set_erg_var_label erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy' erg_load_pwe_hfa, datatype='spec' erg_load_pwe_ofa, datatype='spec' options, 'erg_pwe_hfa_l2_lm_spectra_e_mix', zrange=[1e-8, 1e-2], zlog=1, ytitle='PWE/HFA!C(E_MIX)', ysubtitle='[Hz]', auto_downsample=1 options, 'erg_pwe_ofa_spec_l2_E_spectra_merged', zrange=[1e-8, 1e-2], zlog=1, ytitle='PWE/OFA!Cspec', ysubtitel='[Hz]', auto_downsample=1 tplot, ['erg_pwe_hfa_l2_lm_spectra_e_mix', 'erg_pwe_ofa_spec_l2_E_spectra_merged', 'erg_lepe_l2_3dflux_FEDU_energy' ] ;; Slide #18 erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy', $ trange=tr, mag=magvn, pos=posvn, pitch=[0, 15], suffix='_pa00-15' erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy', $ trange=tr, mag=magvn, pos=posvn, pitch=[85, 95], suffix='_pa85-95' zlim, ['erg_lepe_l2_3dflux_FEDU_energy_pa*'], 1e+2, 1e+6, 1 tplot, ['erg_pwe_hfa_l2_lm_spectra_e_mix', 'erg_pwe_ofa_spec_l2_E_spectra_merged', 'erg_lepe_l2_3dflux_FEDU_energy_pa*'] ;; Slide #19 ;; For Figure 2 of Kazama+2018 dists = erg_lepe_get_dist( 'erg_lepe_l2_3dflux_FEDU', trange=tr ) erg_part_en_pa_spec_plot, dists[450], units='flux' erg_part_en_pa_spec_plot, dists, time='2017-04-12/17:00:00', units='eflux' erg_part_en_pa_spec_plot, dists, time=['2017-04-12/17:05:35','2017-04-12/17:06:08'], /with_contour, units='df_km' ;; Slide #20 erg_part_en_pa_spec_plot, dists, time=['2017-04-12/17:05:35','2017-04-12/17:06:08'], /noplot, rslt=rslt, units='df_km' help, rslt yid = nn( rslt.y_pitchangle, 30. ) ;; for PA=30 deg plot, rslt.x_energy, rslt.z_hist[*, yid], /xlog, /ylog, psym=-6, xtitle='Energy [eV]', ytitle='PSD [s!U3!N/km!U6!N]' yid2 = nn( rslt.y_pitchangle, 90 ) ;; for PA=90 deg oplot, rslt.x_energy, rslt.z_hist[*, yid2], psym=-6, color=spd_get_color('purple') ;; Slide #21 erg_part_maxwellian_fit, rslt.x_energy[2:14], rslt.z_hist[2:14, yid2], n90, kbt90 df_km = erg_part_get_maxwellian( rslt.x_energy, n90, kbt90 ) oplot, rslt.x_energy, df_km, color=spd_get_color('purple'), linestyle=2, thick=2.5 print, 'n = ', n90, ' [/cc]', ' kBT = ', kbt90, ' [eV]' erg_part_maxwellian_fit, rslt.x_energy[2:14], rslt.z_hist[2:14, yid], n30, kbt30 df_km = erg_part_get_maxwellian( rslt.x_energy, n30, kbt30 ) oplot, rslt.x_energy, df_km, color=spd_get_color('black'), linestyle=2, thick=2.5 print, 'n = ', n30, ' [/cc]', ' kBT = ', kbt30, ' [eV]' ;; Slide #22 ;; For Figure 3 of Kazama+2018 timespan, '2017-04-12/16:50', 40, /min & get_timespan, tr erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='pa', units='df_km', $ energy=[90., 92.], suffix='_91ev', trange=tr, mag=magvn, pos=posvn erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='pa', units='df_km', $ energy=[341., 343.], suffix='_342ev', trange=tr, mag=magvn, pos=posvn erg_load_pwe_efd, datatype='pot', uname=uname, pass=pass tplot, ['erg_pwe_efd_l2_Vave', 'erg_lepe_l2_3dflux_FEDU_pa_342ev', 'erg_lepe_l2_3dflux_FEDU_pa_91ev' ] ;; Slide #23 erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy', units='eflux', $ pitch=[85., 95.], suffix='_pa85-95', trange=tr, mag=magvn, pos=posvn erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='energy', units='eflux', $ pitch=[0., 15.], suffix='_pa00-15', trange=tr, mag=magvn, pos=posvn get_data, 'erg_lepe_l2_3dflux_FEDU_energy_pa85-95', data=dperp, dl=dl, lim=lim get_data, 'erg_lepe_l2_3dflux_FEDU_energy_pa00-15', data=dpara, dl=dl, lim=lim !null = min( abs( dperp.v[0, *] - 342. ), id_342ev,/nan ) !null = min( abs( dperp.v[0, *] - 91 ), id_91ev, /nan ) store_data, 'erg_lepe_eflux_342ev', data={x:dperp.x, y:[ [dperp.y[*, id_342ev]], [dpara.y[*, id_342ev]] ] } store_data, 'erg_lepe_eflux_91ev', data={x:dperp.x, y:[ [dperp.y[*, id_91ev]], [dpara.y[*, id_91ev]] ] } ylim, 'erg_lepe_eflux_*ev', 0, 0, 0 options, 'erg_lepe_eflux_*ev', labels=['perp', 'para'], colors=[2, 187], ysubtitle=dl.ztitle tplot, ['erg_pwe_efd_l2_Vave', 'erg_lepe_eflux_342ev', 'erg_lepe_l2_3dflux_FEDU_pa_342ev', 'erg_lepe_eflux_91ev', 'erg_lepe_l2_3dflux_FEDU_pa_91ev' ] ;; Slide #24 ;; For Figure 4 of Kazama+2018 erg_load_pwe_hfa, level='l3', uname=uname, pass=pass ylim, 'erg_pwe_hfa_l2_lm_spectra_e_mix', 10., 100., 1 zlim, 'erg_pwe_hfa_l2_lm_spectra_e_mix', 1e-9, 1e-2, 1 calc, '"erg_pwe_efd_ne_vsc" = 26.39 * exp( -0.6546 * "erg_pwe_efd_l2_Vave" )' erg_lep_part_products, 'erg_lepe_l2_3dflux_FEDU', output='moments', energy=[40., 750.], mag=magvn, pos=posvn, trange=tr ylim, 'erg_lepe_l2_3dflux_FEDU_magt3', 1e+1, 1e+3, 1 store_data, 'erg_lepe_efd_density', data=['erg_pwe_efd_ne_vsc', 'erg_lepe_l2_3dflux_FEDU_density' ] ylim, 'erg_lepe_efd_density', 1e-2, 1e+2, 1 options, 'erg_lepe_efd_density', ytitle='Density', ysubtitle='[/cc]', colors=[2, 6], labels=['total', 'hot'], labflag=-1, constant=[0.1, 1., 10.] tplot, [ 'erg_pwe_hfa_l3_ne_mgf', 'erg_pwe_hfa_l2_lm_spectra_e_mix', 'erg_lepe_efd_density', 'erg_lepe_l2_3dflux_FEDU_magt3' ]