!p.thick=3
!x.thick=3
!y.thick=3
!p.charthick=3
!p.charsize=1.3

plot,findgen(24),wiav(72,40,0:23)/100.,$
xst=1,yrange=[-60,40],yst=1,$
title='Vertical drift at 0Lon 11.25N',xtitle='Local time (hr)',ytitle='m/s'
oplot,findgen(24),wisdv(72,40,0:23)*4/100.,linest=2

for i=0,51 do begin
   oplot,findgen(24),wi_zp2_d(72,40,*,i)/100.,color=100,thick=.8
endfor
plot,findgen(52)+1,wi_zp2_d(72,40,10,*)/100,yrange=[-60,60],yst=1,xst=1,$
xtitle='Day of year',ytitle='m/s',title='Vertical drift: 0Lon 11.25N'
oplot,findgen(52)+1,wi_zp2_d(72,40,5,*)/100,linest=1
oplot,findgen(52)+1,wi_zp2_d(72,40,19,*)/100,linest=2

plot,findgen(24),ui_zp2(72,40,0:23)/100.,$
xst=1,yrange=[-100,100],yst=1,$
title='Zonal drift at 0Lon 11.25N',xtitle='Local time (hr)',ytitle='m/s'
for i=1,51 do begin
   oplot,findgen(24),ui_zp2(72,40,i*24:(i+1)*24-1)/100.
endfor

plot,findgen(24),shift(reform(wi_zp2(42,31,0:23)),lon(42)/15)/100.,$
xst=1,yrange=[-100,40],yst=1,$
title='Vertical drift at 75W 11.25S',xtitle='Local time (hr)',ytitle='m/s'
for i=1,51 do begin
   oplot,findgen(24),$
         shift(reform(wi_zp2(42,31,i*24:(i+1)*24-1)),lon(42)/15)/100.
endfor

plot,findgen(24),shift(reform(ui_zp2(42,31,0:23)),lon(42)/15)/100.,$
xst=1,yrange=[-100,100],yst=1,$
title='Zonal drift at 75W 11.25S',xtitle='Local time (hr)',ytitle='m/s'
for i=1,51 do begin
   oplot,findgen(24),$
         shift(reform(ui_zp2(42,31,i*24:(i+1)*24-1)),lon(42)/15)/100.
endfor

plot,findgen(24),nmf2(72,56,0:23)*1.e-5,xst=1,yrange=[0,5],$
title='NmF2 at 0Lon 51.25N',xtitle='Local time (hr)',ytitle='x10!E11!Nm!E-3!N'
for i=1,51 do begin
   oplot,findgen(24),nmf2(72,56,i*24:(i+1)*24-1)*1.e-5
endfor

plot,findgen(24),hmf2(72,56,0:23),xst=1,yrange=[200,350],$
title='hmF2 at 0Lon 51.25N',xtitle='Local time (hr)',ytitle='km'
for i=1,51 do begin
   oplot,findgen(24),hmf2(72,56,i*24:(i+1)*24-1)
endfor

end