timespan,'7 11 07/10',2,/hours
sc='a'
;
thm_load_state,probe=sc,/get_supp
thm_load_fit,probe=sc,data='fgs',coord='gsm',suff='_gsm'
thm_load_fit,probe=sc,data='fgs',coord='dsl',suff='_dsl'
thm_load_mom,probe=sc ; L2: onboard processed moms
thm_load_esa,probe=sc ; L2: ground processed gmoms, omni spectra
;
; Introduce and correct spacecraft potential
;
thm_load_esa_pkt,probe=sc
get_data,'tha_pxxm_pot',data=tha_pxxm_pot,dlimits=dlimits
tha_pxxm_pot.y(*)=10. ; eV
store_data,'tha_pxxm_pot_corr',data={x:tha_pxxm_pot.x,y:tha_pxxm_pot.y},dlimits=dlimits
;
thm_part_moments, probe = sc, instrum = 'peer', scpot_suffix = '_pxxm_pot_corr', mag_suffix = '_fgs_dsl', tplotnames = tn, verbose = 2 ; names are output into variable tn
;
options,'tha_peer_density','colors',['b']
options,'tha_peim_density','colors',['r']
store_data,'tha_pexm_density',data='tha_peer_density tha_peim_density'
options,'tha_pexm_density','colors',['b','r']
options,'tha_pexm_density',yrange=[0,2]
options,'tha_pexm_density',ylog=0
;
tplot,'tha_fgs_gsm tha_pexm_density tha_pe?r_en_eflux'