;***********************************************************************************
; ESA data - all s/c


startdate = '2007-03-23/0:00'

ndays=2
thm_init
loadct2,43
timespan,startdate,ndays
probe=['a','b','c','d','e']

thm_load_state,/get_support_data
thm_load_esa_pkt

.r
for i=0,4 do begin
	sc=probe(i)
	cotrans,'th'+sc+'_state_pos','th'+sc+'_state_pos'+'_gse',/gei2gse
	get_data,'th'+sc+'_state_pos_gse',data=tmp
	store_data,'th'+sc+'_state_pos_gse_x',data={x:tmp.x,y:tmp.y(*,0)/6370.}
		options,'th'+sc+'_state_pos_gse_x','ytitle','th'+sc+'_X-GSE'
	store_data,'th'+sc+'_state_pos_gse_y',data={x:tmp.x,y:tmp.y(*,1)/6370.}
		options,'th'+sc+'_state_pos_gse_y','ytitle','th'+sc+'_Y-GSE'
	store_data,'th'+sc+'_state_pos_gse_z',data={x:tmp.x,y:tmp.y(*,2)/6370.}
		options,'th'+sc+'_state_pos_gse_z','ytitle','th'+sc+'_Z-GSE'
	name1='th'+sc+'_peif_en_counts'
		zlim,name1,1.,1.e4,1
	name1='th'+sc+'_peir_en_counts'
		zlim,name1,1.,1.e4,1
	name1='th'+sc+'_peef_en_counts'
		zlim,name1,10.,1.e5,1
	name1='th'+sc+'_peer_en_counts'
		zlim,name1,10.,1.e5,1
endfor
end

; plot reduced distribution

sc='b'
tplot,['th?_peir*counts','th?_peer*counts'],title='Themis',var_label=['th'+sc+'_state_pos_gse_x','th'+sc+'_state_pos_gse_y','th'+sc+'_state_pos_gse_z']

; plot full distribution

sc='b'
tplot,['th?_peif*counts','th?_peef*counts'],title='Themis',var_label=['th'+sc+'_state_pos_gse_x','th'+sc+'_state_pos_gse_y','th'+sc+'_state_pos_gse_z']

; plot burst distribution

sc='b'
tplot,['th?_peib*counts','th?_peeb*counts'],title='Themis',var_label=['th'+sc+'_state_pos_gse_x','th'+sc+'_state_pos_gse_y','th'+sc+'_state_pos_gse_z']