PRO STEREO_plot, doy1, doy2, DT=dt

; DT = Decimal day shift of STEREO B data

@'D:\idl\pc_colors.pic'

if not keyword_set(DT) then dt = 0.0

STEREO_read, sta, stb

stb.doy = stb.doy + dt

ix = where(sta.doy ge doy1 and sta.doy le doy2)
sta = sta[ix]

ix = where(stb.doy ge doy1 and stb.doy le doy2)
stb = stb[ix]

; Read in model data

Xycols2,'D:\Meetings\CCMC Arecibo 2007\ENLIL_Line_20070821_BBrV_Lon-11_Lat7.txt', $
  6, 1, br, 4, bbr, 5, bb, skip=7
Xycols2,'D:\Meetings\CCMC Arecibo 2007\ENLIL_Line_20070821_BBrV_Lon15_Lat7.txt', $
  6, 1, ar, 4, abr, 5, ab, skip=7

; Using the Parker formula normalize B to 1 AU
; We take 500 km/s constant solar wind speed

abmod = 1.305d0*ar*ar*ab/sqrt(1.d0+0.704d0*ar*ar)
abrmod = ar*ar*abr
bbmod = 1.305d0*br*br*bb/sqrt(1.d0+0.704d0*br*br)
bbrmod = br*br*bbr

; Convert the distance scale to time line assuming a constant 500 km/s speed

adt = (1.d0-ar)*1.495979d8/500.d0/86400.d0
bdt = (1.d0-br)*1.495979d8/500.d0/86400.d0

; The simulation is at 2007/08/21 0600 UT = 233.25 doy

a_enlil_time = 233.25 + adt
b_enlil_time = 233.25 + bdt

pos = get_pos(4, sp=[0.01,0.01,0.01])
defplotn, sta.doy, sta.b, pos=pos[0,*], xr=[doy1, doy2], xtitle='!17Days of 2007', $
 ytitle='!9!!!17B!9!!!17 (nT)', yr=[0,12], ytitpos=0.03
oplot, sta.doy, sta.b, color=red, thick=1
oplot, stb.doy, stb.b, color=blue, thick=1
;oplot, a_enlil_time, abmod, color=green, thick=3
;oplot, b_enlil_time, bbmod, color=purple, thick=3

defplotn, sta.doy, sta.bn, pos=pos[1,*], xr=[doy1, doy2], yr=[-10,10], ytitpos=0.03, $
 ytitle='!17B!dN!n (nT)', /noxtick
oplot, sta.doy, sta.bn, color=red, thick=1
oplot, stb.doy, stb.bn, color=blue, thick=1

defplotn, sta.doy, sta.bt, pos=pos[2,*], xr=[doy1, doy2], yr=[-10,10], ytitpos=0.03, $
 ytitle='!17B!dT!n (nT)', /noxtick
oplot, sta.doy, sta.bt, color=red, thick=1
oplot, stb.doy, stb.bt, color=blue, thick=1

defplotn, sta.doy, sta.br, pos=pos[3,*], xr=[doy1, doy2], yr=[-10,10], ytitpos=0.03, $
 ytitle='!17B!dR!n (nT)', /noxtick
oplot, sta.doy, sta.br, color=red, thick=1
oplot, stb.doy, stb.br, color=blue, thick=1
;oplot, a_enlil_time, -abrmod, color=green, thick=3
;oplot, b_enlil_time, -bbrmod, color=purple, thick=3

end