how to make maximum wind plot along with level showing maximum wind for each grid
Charles Seman
Charles.Seman at NOAA.GOV
Tue Nov 24 14:31:54 EST 2009
Hi Phillip,
Please find attached a script "fig3ab.gs" and called script functions
"pause.gsf" and "scale_zdef_zlevs.gsf". "scale_zdef_zlevs.gsf" makes a
local (temporary) ctl file with scaled vertical levels for plotting
purposes: here specifically to convert the original zdef zlevs values
from meters to kilometers (there may be a better way to do this; maybe
using commands "set ylab off", followed by "set yaxis ..."). I hope all
of the scripts needed for the fig3ab.gs script are attached; if not,
please let me know. You may have questions about these scripts; please
contact me if you do...
Hope this helps,
Chuck
Phillip Lin wrote:
> Hi Chuck,
> Thanks for your dedicated input to this problem we have.
> And could you be kind to send me a script you did before for the sane or
> similar question? Then, maybe it is much helpful for me to solve this
> problem!
> Thanks much!
> Phillip
>
>
>> Hi Phillip,
>>
>> Sorry I didn't explain further: the code was meant as an example of how
>> to find and plot extrema in a field of data. To do the max wind, the
>> code would have to be modified to search in a wind speed field using
>> "maxloc"... and the exact code would also depend on the domain...
>>
>> Hope this helps,
>> Chuck
>>
>> Phillip Lin wrote:
>>
>>> 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."
>>>>
>>>>
>>>>
>> --
>>
>> 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."
>>
>>
--
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."
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: fig3ab.gs
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20091124/043a4bbe/attachment.pl
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pause.gsf
Type: application/x-font-type1
Size: 355 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20091124/043a4bbe/attachment.bin
-------------- next part --------------
A non-text attachment was scrubbed...
Name: scale_zdef_zlevs.gsf
Type: application/x-font-type1
Size: 3604 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20091124/043a4bbe/attachment-0001.bin
More information about the gradsusr
mailing list