timespan,'7 6 8/21:00',3,/hours
;
sc='c'
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
;
; plot density from spacecraft potential
;
tinterpol_mxn,'thc_peer_t3','thc_pxxm_pot',newname='thc_peer_t3_int'
get_data,'thc_pxxm_pot',data=thc_pxxm_pot,dl=dl
get_data,'thc_peer_t3_int',data=thc_peer_t3_int
thm_pot2dens,thc_pxxm_pot.y,thc_pxxm_potdens,Te=total(thc_peer_t3_int.y,2)/3.
store_data,'thc_pxxm_potdens', $
  data={x:thc_pxxm_pot.x,y:thc_pxxm_potdens},dl=dl
;
; plot uncorrected densities
;
store_data,'thc_peer_en_eflux_pot',data='thc_peer_en_eflux thc_pxxm_pot'
options,'thc_fgs_gsm',yrange=[-70,100]
options,'thc_peer_density',colors=['r']
options,'thc_peir_density',colors=['b']
options,'thc_pxxm_potdens',colors=['g']
options,'thc_pxxm_potdens',ylog=1
options,'thc_peer_t3',ylog=0
options,'thc_pxxm_pot',ylog=0
options,'thc_pe?r_en_eflux*',yrange=[7.,25000.]
store_data,'thc_densities',data='thc_peir_density thc_peer_density thc_pxxm_potdens'
tplot,'thc_fgs_gsm thc_peir_velocity_gsm thc_densities thc_psef_en_eflux thc_peer_en_eflux_pot thc_peir_en_eflux'
tlimit,['7 6 8/21:00','7 6 8/24:00']
;
tlimit,['7 6 8/22:00','7 6 8/22:30']
;
tvectot,'thc_peir_velocity_gsm',newname='thc_peir_velocity_gsmt'
tvectot,'thc_peir_velocity_gsm',tot='thc_peir_velocity_t'
tinterpol_mxn,'thc_peir_velocity_t','thc_pxxm_pot',newname='thc_peir_velocity_tint'
get_data,'thc_peir_velocity_tint',data=thc_peir_velocity_tint
get_data,'thc_pxxm_pot',data=thc_pxxm_pot
eflow=1000.*(thc_peir_velocity_tint.y/310.)^2 - thc_pxxm_pot.y; in eV
store_data,'thc_eflow',data={x:thc_peir_velocity_tint.x,y:eflow}
store_data,'thc_peir_en_eflux-n-flow',data='thc_peir_en_eflux thc_eflow'
options,'thc_peir_en_eflux*',yrange=[7.,25000.]
;
tplot,'thc_fgs_gsm thc_peir_velocity_gsmt thc_densities thc_psef_en_eflux thc_peer_en_eflux_pot thc_peir_en_eflux-n-flow'
tlimit,['7 6 8/22:00','7 6 8/22:30']