timespan,'8 2 22/04:00',1,/hours
sc='b'
thm_load_state,probe=sc,/get_supp
thm_load_fit,probe=sc,coord='dsl',suff='_dsl'
thm_load_fit,probe=sc,coord='gsm',suff='_gsm'
;
;Exoffset=1.125
;Eyoffset=0.223
Exoffset=1. ; <---- modify this from a time of zero field
Eyoffset=0.1; <---- modify this from a time of zero field
angle=10. ; degrees
tanangle=tan(angle*!PI/180.)
get_data,'th'+sc+'_efs_0_dsl',data=thx_efs_dsl
get_data,'th'+sc+'_fgs_dsl',data=thx_fgs_dsl
igood=where(abs(thx_fgs_dsl.y(*,2)/sqrt(thx_fgs_dsl.y(*,0)^2+thx_fgs_dsl.y(*,1)^2)) ge tanangle,janygood)
ibad=where(abs(thx_fgs_dsl.y(*,2)/sqrt(thx_fgs_dsl.y(*,0)^2+thx_fgs_dsl.y(*,1)^2)) lt tanangle,janybad)
;
thx_efs_dsl.y(*,0)=thx_efs_dsl.y(*,0)-Exoffset
thx_efs_dsl.y(*,1)=thx_efs_dsl.y(*,1)-Eyoffset
thx_efs_dot0_dsl=thx_efs_dsl
if (janybad ge 1) then thx_efs_dot0_dsl.y(ibad,*)=!VALUES.F_NAN
if (janygood lt 1) then print,'*****WARNING: NO GOOD 3D ExB data'
if (janygood ge 1) then thx_efs_dot0_dsl.y(igood,2)=$
  -(thx_efs_dsl.y(igood,0)*thx_fgs_dsl.y(igood,0)+$
    thx_efs_dsl.y(igood,1)*thx_fgs_dsl.y(igood,1)+$
    thx_efs_dsl.y(igood,2)*thx_fgs_dsl.y(igood,2))/ $
    thx_fgs_dsl.y(igood,2)
;
thx_exb_dot0_dsl=thx_efs_dot0_dsl
store_data,'th'+sc+'_efs_dot0_dsl',data={x:thx_efs_dot0_dsl.x,y:thx_efs_dot0_dsl.y}
options,'th'+sc+'_efs_dot0_dsl','colors',[2,4,6]
tcrossp,'th'+sc+'_efs_dot0_dsl','th'+sc+'_fgs',newname='thx_exb_temp'
tdotp,'th'+sc+'_fgs','th'+sc+'_fgs',newname='thx_bdotb_temp'
div_data,'thx_exb_temp','thx_bdotb_temp',newname='thx_exbob2_temp'
get_data,'thx_exbob2_temp',data=thx_exbob2_temp
store_data,'th'+sc+'_Vexb_dot0_dsl',data={x:thx_exbob2_temp.x,y:thx_exbob2_temp.y*1000.}
options,'th'+sc+'_Vexb_dot0_dsl','colors',[2,4,6]
tnormalize,'th'+sc+'_fgs_gsm',newname='th'+sc+'_fgs_unit'
;
thm_cotrans,'th'+sc+'_Vexb_dot0',in_coord='dsl',out_coord='gsm',in_suff='_dsl',out_suff='_gsm'
thm_cotrans,'th'+sc+'_efs_dot0',in_coord='dsl',out_coord='gsm',in_suff='_dsl',out_suff='_gsm'
;
tvectot,'thb_fgs_gsm',newname='thb_fgs_gsmt'
tplot,'thb_fgs_gsm th'+sc+'_efs_0_dsl th'+sc+'_efs_dot0_gsm th'+sc+'_Vexb_dot0_gsm'
;