* This one is tricky as it creates daily climatologies and then produces a hovmuller.
* Hovmuller of shaded 500hPa mean wind.
* Hovmuller of vectors is also plotted alongside the shaded wind magnitude. 
* Created Wednesday July 1, 2010

*************************************************************
*************************************************************

'reinit'
'set grid off'
'set grads off'

* Below I open the daily zonal wind file, set the level to the 500hPa level, and create daily climatologies from Jan1-Dec31
'sdfopen F:\SampeXie_files\uwnd.daily.nc'
'set lev 500'
'set t 1 365'
'define MAYdayu = ave(u,t+0,t=10593,365)'

* I do the same as above, but for the meridional wind
'sdfopen F:\SampeXie_files\vwnd.daily.nc'
'set lev 500'
'set t 1 365'
'define MAYdayv = ave(v.2,t+0,t=10593,365)'

* create a climatology of the wind magnitude from Jan1-Dec31
'set t 1 365'
'define wind = mag(MAYdayu,Maydayv)'

* Begin Hovmuller set up specifying the latitude range (12N to 30N) averaged at 80W
'set gxout shaded'
'set x 1'
'set lat 12 30'
'set t 1 365'
'define DAYwind = ave(wind,lon=-80,lon=-80)'

* Below creates the plot from May 1 (day 121) to Aug 1 (day 213)
'set t 121 213' 
'set xyrev on'
'set tlsupp year'
'set csmooth on'
'd DAYwind'
'cbar'

* over-plots the wind vectors (direction), same theory as above for hovmuller
'set gxout vector'
'set x 1'
'set lat 12 30'
'set t 121 213' 
'set clab masked'
'set cthick 7'
'set xyrev on'
'set tlsupp year'
'set cint 4'
'set arrscl 0.4 5'
'set arrowhead 0.1'
'set lon -80'
* the below line plots only every 4th vector for better viewing, or else the vectors would be too crowded!
'd skip (Maydayu,0,4);maydayv'

'draw title 500hPa mean wind 80W 1980-2008'
* output the plot to below specified location
'printim F:\GrADSex\teddyFIGURE.gif gif x800 x600 white'