Wind vetcor plotting

[Matt Georgescu]

Vijay,

there is no attached figure in your e-mail.

As regards your question, Eric L. Altshuler provided me with an answer to
a similar query I had some time ago.

If you'd like to take a time average over separate wind components, try
this:

*Define "skipped" version of u and v-wind separately
'define x = skip(ave(uwnd.1, t=1, t=30),5,5)'
'define y = skip(ave(vwnd.2, t=1, t=30),5,5)'
*Resulting vector of "skipped" u and v-wind.
'd x;y'

Here is more on the "skip" function from the GrADS Documentation:
http://grads.iges.org/grads/gadoc/gadocindex.html

matt

On Thu, 25 Oct 2007, Kanawade, V.P. wrote:

> Dear All,
> Greetings from Leicester University.
> I am trying to plot mean wind vector for July 2004. Could you please
> have a look at the GrADS code to generate wind vector plot and please
> correct and suggest to produce nice and informative plot.
>
> 'reinit'
> 'set display color white'
> 'c'
> 'set mpdraw on'
> 'set arrowhead 0.15
> 'set arrlab on'
> 'set map 1 1 7'
> 'set gxout vector'
> 'sdfopen F:/uwnd05_1000mb.nc   (uwind file)'
> 'sdfopen F:/vwnd05_1000mb.nc   (vwind file)'
> 'set lat -90 90'
> 'set lon -180 180'
> 'set lev 1000'
> 'd uwnd.1;vwnd.2'
> 'draw title Wind Vector at 1000 hPa   July 2004'
> 'printim windjul04_1000mb.gif gif'
>
> please see attached wind vector plot.
>
> Then I tired to average u and v component of windfields using following
> command but it doesn't work.
>
> 'define x =ave(uwnd.1, t=1, t=30)'
> 'define y=ave(vwnd.2, t=1, t=30)'
> 'd x;y'
>
> Any suggestions are very much appreciated.
> Thank you.
>
>
> Best regards,
> Vijay Kanawade
>
> Ph D Student
> Earth Observation Science Group
> Space Research Centre
> University of Leicester
> Leicester, LE1 7RH
> Tel: +44 (0) 116 252 2771
> Email: vpk1 at le.ac.uk
>