<div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Andries,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Thanks for providing a concise description of the problem and a related link.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I've created a ticket for this, NCL-2501. Hopefully we'll get a fix in for 6.4.0, but if not, then 6.4.1.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Sep 27, 2016 at 9:26 AM, Devries, Andries <span dir="ltr"><<a href="mailto:a.devries@mpic.de" target="_blank">a.devries@mpic.de</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div>
<div style="direction:ltr;font-family:Tahoma;color:#000000;font-size:10pt">Dear NCL team,<br>
<br>
I just came across a possible bug in the gc_latlon function. Two points that are located very nearby can result in a "nan", see below. In a previous post a similar bug was reported (see link below) for two points with the exact same location. A solution was
found for this issue, but possibly two points that are located very nearby each other may still suffer the same issue. Further I noticed that a slight change in any of the four coordinate values can result in a 0 value instead of a nan. A possible way to work
around this could be using the "replace_ieeenan" function, but fixing the bug might be preferable.<br>
<br>
Looking forward to your response. Thanks!<br>
Cheers,<br>
Andries<br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<br>
<br>
Copyright (C) 1995-2015 - All Rights Reserved<br>
University Corporation for Atmospheric Research<br>
NCAR Command Language Version 6.3.0<br>
The use of this software is governed by a License Agreement.<br>
See <a href="http://www.ncl.ucar.edu/" target="_blank">http://www.ncl.ucar.edu/</a> for more details.<br>
ncl 0> print(gc_latlon(88.<wbr>67137908935547,-179.<wbr>9999847412109,88.<wbr>67137908935547,180,0,4))<br>
<br>
<br>
Variable: unnamed (return)<br>
Type: float<br>
Total Size: 4 bytes<br>
1 values<br>
Number of Dimensions: 1<br>
Dimensions and sizes: [1]<br>
Coordinates: <br>
Number Of Attributes: 4<br>
units : kilometers<br>
gclon : ( 180, 180 )<br>
gclat : ( 88.67138, 88.67138 )<br>
spacing : nan<br>
(0) nan<br>
<br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<br>
link to previously reported bug:<br>
<a href="https://www.ncl.ucar.edu/Support/talk_archives/2011/0317.html" target="_blank">https://www.ncl.ucar.edu/<wbr>Support/talk_archives/2011/<wbr>0317.html</a><br>
------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>--------------<br>
<br>
</div>
</div>
<br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>