<div dir="ltr">your code looks good, I think the problem is that Grads defines longitude in values from 0 - 360 so "if lon < 0 then lon=lon + 180"</div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jun 29, 2017 at 9:24 AM, CW Weather <span dir="ltr"><<a href="mailto:weatherstu@chorleyweather.com" target="_blank">weatherstu@chorleyweather.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-size:10pt;font-family:Tahoma,Arial,Helvetica,sans-serif">
<p>Hi there, </p>
<p>Im trying to plot t2m temperatures at certain LAT-LON but with no success. I have found some code that is used with WRF but isnt working on my grads when running GFS data. Can anybody shed more light on this or share some code for doing this. </p>
<p> </p>
<p>Code Example:</p>
<p>***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</p><span class="HOEnZb"><font color="#888888">
<div>-- <br>
<p>Kind regards,</p>
<p><strong>Stuart Markham</strong></p>
<p><em>Founder of CW Weather</em></p>
<p><img src="cid:9e849405edb0a79e162c211a585f3245@chorleyweather.com" alt="" width="147" height="65"></p>
</div>
</font></span></div>
<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></blockquote></div><br></div>