how to make maximum wind plot along with level showing maximum wind for each grid

[Phillip Lin]

Hi Chuck,
I also met this problem. And there are something unclear to me (only to
me!)in your last email, because:
(1) the level that has minimum temp isn't the same as the level where
there is maximum zonal wind speed;
(2) what does the meaning of the step (3) you said; sorry how can plot the
values (xloc,j,yloc,j) on a (y,z)plot (this is the cross section in the
height-latitude plane?
(3)could you be kind to tell me the basic ideas to get the maximum jet
stream, or... this means the contours of zonal wind speed in a specific
level as 250mb or ...?
(4)could you be kind to send me the script to do this plotting, and data
file too for a test.
I am quite confused about this! Thank you so much in advance!
Phillip


> Sushant,
>
> Here's some code that:
> 1) finds the height level "zmin" of the minimum temperature (represented
> by variable ztav = zonal average of weighted time average temperature)
> at each horizontal location "j" in a (y,z) physical domain
> 2) uses GrADS function "gr2xy"
> (http://grads.iges.org/grads/gadoc/script.html#commands) to transform
> the grid point location (j,zmin) to plot location (xloc.j,yloc.j)
> 3) plots the values (xloc.j,yloc.j) on a (y,z) plot
>
>     say '   find and save location of minimum temperature...'
>
>     j=1
>     while ( j <= ny )
>       'set y 'j
>       'set gxout print'
>       'd minloc(ztav,z=1,z='nz')' ; rec3 = sublin(result,3) ; zmin =
> subwrd(rec3,1)
> *      say result
>       say '   ...j, zmin = 'j', 'zmin
>       'query gr2xy 'j' 'zmin
>       rec = sublin(result,1) ; xloc.j = subwrd(rec,3) ; yloc.j =
> subwrd(rec,6)
> *      say result
>       say '      xloc, yloc = 'xloc.j', 'yloc.j
>       j=j+1
>     endwhile
> ...
>   if( plot_min_temp = 'yes' )
>     say '   plot level of minimum temperature... '
>     j=1
>     while ( j <= ny )
>       'set y 'j
> *       say '      xloc, yloc = 'xloc.j', 'yloc.j
>       'set line 1'
>       'draw mark 3 'xloc.j' 'yloc.j' 0.03'
>       j=j+1
>     endwhile
>   endif
>
> Please contact me if you have any questions about this code, or would
> like a script that does a plot.
>
> Hope this helps,
> Chuck
>
>
> sushant puranik wrote:
>> Hello Sir
>> I want to study the structure of a jet stream for that purpose i wish
>> to make a plot showing maximum wind speed for each grid box alongwith
>> the level in which the speed is maximum. Means it should print level
>> with maximum wind speed value. I am trying to solve this problem since
>> last 15 days but i am not able to solve it. Does anyone has any idea
>> how i can make such plot using GrADS. I am using ECMWF interim
>> reanalysis data for this purpose.
>>
>> Any suggestion is appreciated
>>
>> Thanking you.
>>
>> --
>> Sushant Puranik
>> Junior Research Fellow
>> Dept. of Atmospheric & Space Sciences,
>> University of Pune,
>> Pune-07,
>> India.
>
> --
>
> Please note that Charles.Seman at noaa.gov should be considered my NOAA
> email address, not cjs at gfdl.noaa.gov.
>
> ********************************************************************
>  Charles Seman                                Charles.Seman at noaa.gov
>  U.S. Department of Commerce / NOAA / OAR
>  Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
>  201 Forrestal Road                              fax: (609) 987-5063
>  Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
> ********************************************************************
>
> "The contents of this message are mine personally and do not necessarily
> reflect any position of the Government or NOAA."
>