[GrADS] How to mark the cities in the contour map?

[Charles Seman]

Qingfang Dai,

Please try something like this:

'query w2xy 53.32 -113.58'
xpos = subwrd(result,3)
ypos = subwrd(result,6)
say '    (xpos,ypos) = ('xpos','ypos')'
'draw mark 3 'xpos' 'ypos' 2'


I hope this helps,
Chuck

Qingfang Dai wrote:
> Hi,
>
> I revised my gs file as
> 'reinit'
> 'open ../Climatology.ctl'
> 'set grads off'
> 'set mproj nps'
> 'set grid off'
> 'set frame off'
> 'set cstyle 1'
> 'set csmooth off'
> 'set cmin 0.0'
> 'set clab forced'
> 'set t 1 1'
> 'set mpdset ../province'
> 'd prep'
> 'set font 1'
> 'set string 1 tc 3'
> 'set strsiz 0.15 0.15'
> 'query w2xy 53.32 -113.58'
> 'draw mark 3 x y 2'
> 'printim ../climatology.gif gif white'
>
> But it's still not working.
>
> Thanks.
>
> Qingfang
>
>
> On Thu, 17 Jan 2008 23:31:14 +0530, M.Mahakur wrote
>
>> Hi,
>>
>> Convert the lon, lat first to vpage co-ordinate using 'query w2xy lon
>> lat' and later use the results (x and y) in the 'draw mark mark_type
>> x y size' to mark your desired point(s).
>>
>> Regards,
>>
>> Mahakur
>>
>> Qingfang Dai wrote:
>>
>>
>>> Hi All,
>>>
>>> I want to mark 2 cities in a contour map. The lattitude and longitude of
>>> these two cities are (53.32,-113.58 ) and (51.12,-114.02). So I used the
>>> command 'draw mark 3 53.32 -113.58 1' and 'draw mark 3 51.12 -114.02 1' in
>>> the .gs file but I they didn't work. Please tell me how I can mark the two
>>> cities in the contour map with closed circle and cities name. Thanks.
>>>
>>> Qingfang
>>>
>>>
>>>
>>>
>> ----------------------------------------------------------
>> M. Mahakur
>> FORECASTING RESEARCH DIVISION
>> INDIAN INSTITUTE OF TROPICAL METEOROLOGY
>> (Min. of Earth Sciences, Govt. of India)
>> DR. HOMI BHABHA ROAD
>> NCL POST, PASHAN
>> PUNE - 411 008
>> INDIA.
>> -----------------------------------------------------------
>>
>
>
> --
> Department of Math and Stat Sciences (http://www.math.ualberta.ca)
>

--

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 reflect
any position of the Government or NOAA."