timespan,'8 1 16/14:00',6,/hours
;
sc='d'
;
thm_load_state,probe=sc,/get_supp
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
thm_load_sst,level=2,probe=sc
;
; plot density from spacecraft potential
;
tinterpol_mxn,'thd_peer_t3','thd_pxxm_pot',newname='thd_peer_t3_int'
get_data,'thd_pxxm_pot',data=thd_pxxm_pot,dl=dl
get_data,'thd_peer_t3_int',data=thd_peer_t3_int
thm_pot2dens,thd_pxxm_pot.y,thd_pxxm_potdens,Te=total(thd_peer_t3_int.y,2)/3.
store_data,'thd_pxxm_potdens', $
  data={x:thd_pxxm_pot.x,y:thd_pxxm_potdens},dl=dl
;
; plot uncorrected densities
;
store_data,'thd_peer_en_eflux_pot',data='thd_peer_en_eflux thd_pxxm_pot'
options,'thd_fgs_gsm',yrange=[-150,150]
options,'thd_peer_density',colors=['r']
options,'thd_peir_density',colors=['b']
options,'thd_pxxm_potdens',colors=['g']
options,'thd_pxxm_potdens',ylog=1
options,'thd_peer_t3',ylog=0
options,'thd_pxxm_pot',ylog=0
options,'thd_pe?r_en_eflux*',yrange=[7.,25000.]
store_data,'thd_densities',data='thd_peir_density thd_peer_density thd_pxxm_potdens'
tplot,'thd_fgs_gsm thd_peer_t3 thd_pxxm_pot thd_densities thd_psef_en_eflux thd_peer_en_eflux_pot thd_peir_en_eflux'
;
; CORRECT DENSITIES
;
; load L0 omni spectra, all ESA data in memory
thm_load_esa_pkt,probe=sc
;
; PEIR MOMS/SPECTRA
; Remove constant background ; this removes also scattered electrons
thm_part_moments, probe = sc, instrum = 'peir', scpot_suffix = '_pxxm_pot', $
   trange=['8 1 16/14:00','8 1 16/20:00'], erange=[0,31], $
   mag_suffix = '_fgs_dsl', tplotnames = tn, verbose = 2, $
   /bgnd_remove ; names are output into variable tn
;
; PEER MOMS/SPECTRA
; Remove constant background and integrate above 40eV to remove scattered electrons
thm_part_moments, probe = sc, instrum = 'peer', scpot_suffix = '_esa_pot', $
   trange=['8 1 16/14:00','8 1 16/20:00'], erange=[0,24], $
   mag_suffix = '_fgs_dsl', tplotnames = tn, verbose = 2, $
   /bgnd_remove ; names are output into variable tn
;
; scpot determination of density, with (now/see above) better temperature
;
tinterpol_mxn,'thd_peer_t3','thd_pxxm_pot',newname='thd_peer_t3_int'
get_data,'thd_pxxm_pot',data=thd_pxxm_pot,dl=dl
get_data,'thd_peer_t3_int',data=thd_peer_t3_int
thm_pot2dens,thd_pxxm_pot.y,thd_pxxm_potdens, $
  Te=total(thd_peer_t3_int.y,2)/3.
store_data,'thd_pxxm_potdens', $
  data={x:thd_pxxm_pot.x,y:thd_pxxm_potdens},dl=dl

tplot,'thd_fgs_gsm thd_peer_t3 thd_pxxm_pot thd_densities ' + $
  'thd_psef_en_eflux thd_peer_en_eflux_pot thd_peir_en_eflux'
;
;