<div dir="ltr"><div><div><div><div>Stuart,<br><br></div>This is part of my script that I use for plotting temperature at specific latitude/longitude points:<br><br>-----------------------------------------------------------------------------------------------------------------------<br><br>'set strsiz 0.08'<br>'set string 1 c 6'<br>#'set gxout print'<br>#'set<br># Hamilton<br>***First query to get x,y from lat/lon:<br>lon = -79.93<br>lat = 43.17<br>'q w2xy 'lon' 'lat<br>x = subwrd(result,3)<br>y = subwrd(result,6)<br>***Second query to get temp at i,j rather than x,y:<br>'q w2gr 'lon' 'lat<br>x2 = subwrd(result,3)<br>y2 = subwrd(result,6)<br>'define t2 = (TMP2m(t='dis_t')-273.15)'<br>***Query (defval) value of temp at i,j<br>'q defval t2 'x2' 'y2<br>t2m = subwrd(result,3)<br>**Round temperature value:<br>t2m = math_nint(t2m)<br>'draw string 'x' 'y' 't2m<br>#<br># London<br>***First query to get x,y from lat/lon:<br>lon = -81.15<br>lat = 43.03<br>'q w2xy 'lon' 'lat<br>x = subwrd(result,3)<br>y = subwrd(result,6)<br>***Second query to get temp at i,j rather than x,y:<br>'q w2gr 'lon' 'lat<br>x2 = subwrd(result,3)<br>y2 = subwrd(result,6)<br>'define t2 = (TMP2m(t='dis_t')-273.15)'<br>***Query (defval) value of temp at i,j<br>'q defval t2 'x2' 'y2<br>t2m = subwrd(result,3)<br>**Round temperature value:<br>t2m = math_nint(t2m)<br>'draw string 'x' 'y' 't2m<br><br>---------------------------------------------------------------------------------------<br><br></div>The point plots should look like the enclosed image.<br><br></div>- Brian Bernard<br><br></div><a href="http://www.goldenhorseshoewx.ca">www.goldenhorseshoewx.ca</a><br><div><div><div><br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 29, 2017 at 9:31 AM, <span dir="ltr"><<a href="mailto:gradsusr-request@gradsusr.org" target="_blank">gradsusr-request@gradsusr.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Send gradsusr mailing list submissions to<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
or, via email, send a message with subject or body 'help' to<br>
<a href="mailto:gradsusr-request@gradsusr.org">gradsusr-request@gradsusr.org</a><br>
<br>
You can reach the person managing the list at<br>
<a href="mailto:gradsusr-owner@gradsusr.org">gradsusr-owner@gradsusr.org</a><br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of gradsusr digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
1. Plotting Values via LAT/LON (CW Weather)<br>
<br>
<br>
------------------------------<wbr>------------------------------<wbr>----------<br>
<br>
Message: 1<br>
Date: Thu, 29 Jun 2017 14:24:25 +0100<br>
From: CW Weather <<a href="mailto:weatherstu@chorleyweather.com">weatherstu@chorleyweather.com</a><wbr>><br>
Subject: [gradsusr] Plotting Values via LAT/LON<br>
To: GrADS Users Forum <<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a>><br>
Message-ID: <<a href="mailto:af14690d19f155b98cc392beb21798a2@chorleyweather.com">af14690d19f155b98cc392beb2179<wbr>8a2@chorleyweather.com</a>><br>
Content-Type: text/plain; charset="us-ascii"<br>
<br>
Hi there,<br>
<br>
Im trying to plot t2m temperatures at certain LAT-LON but with no<br>
success. I have found some code that is used with WRF but isnt working<br>
on my grads when running GFS data. Can anybody shed more light on this<br>
or share some code for doing this.<br>
<br>
Code Example:<br>
<br>
***First query to get x,y from lat/lon:<br>
lon = -0.118092'<br>
lat = 51.509865'<br>
'q w2xy 'lon' 'lat<br>
x = subwrd(result,3)<br>
y = subwrd(result,6)<br>
***Second query to get temp at i,j rather than x,y:<br>
'q w2gr 'lon' 'lat<br>
x2 = subwrd(result,3)<br>
y2 = subwrd(result,6)<br>
'define t2m = const((tmp2m-273.16),0,-u)'<br>
***Query (defval) value of temp at i,j<br>
'q defval t2m 'x2' 'y2<br>
tmp2m = subwrd(result,3)<br>
**Round temperature value:<br>
t2m = math_nint(t2m)<br>
'draw string 'x' 'y' 't2m<br>
<br>
--<br>
Kind regards,<br>
<br>
STUART MARKHAM<br>
<br>
_Founder of CW Weather_<br>
-------------- next part --------------<br>
An HTML attachment was scrubbed...<br>
URL: <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20170629/ab18760e/attachment.html" rel="noreferrer" target="_blank">http://gradsusr.org/pipermail/<wbr>gradsusr/attachments/20170629/<wbr>ab18760e/attachment.html</a><br>
-------------- next part --------------<br>
A non-text attachment was scrubbed...<br>
Name: 9e849405.png<br>
Type: image/png<br>
Size: 22440 bytes<br>
Desc: not available<br>
Url : <a href="http://gradsusr.org/pipermail/gradsusr/attachments/20170629/ab18760e/attachment.png" rel="noreferrer" target="_blank">http://gradsusr.org/pipermail/<wbr>gradsusr/attachments/20170629/<wbr>ab18760e/attachment.png</a><br>
<br>
------------------------------<br>
<br>
______________________________<wbr>_________________<br>
gradsusr mailing list<br>
<a href="mailto:gradsusr@gradsusr.org">gradsusr@gradsusr.org</a><br>
<a href="http://gradsusr.org/mailman/listinfo/gradsusr" rel="noreferrer" target="_blank">http://gradsusr.org/mailman/<wbr>listinfo/gradsusr</a><br>
<br>
<br>
End of gradsusr Digest, Vol 88, Issue 30<br>
******************************<wbr>**********<br>
</blockquote></div><br></div>