<div dir="ltr"><div class="gmail_default" style="font-size:small">Andries,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Sorry for the delay in response. We are working on getting a new release out and have a tight deadline.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I always thought gc_inout handled points across the equator just fine, but I haven't specifically tried this out. If your longitudes go from -180 to 180, then one thing you could try is converting them to go from 0 to 360, which may take care of that big jump:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">lon = where(lon.lt.0,lon+360,lon)</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">If this doesn't help, can you send us a clean script that shows the issue?</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Dec 12, 2016 at 2:41 PM, 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 users/developers,<br>
<br>
I ran into a problem with ncl (v.6.3.0) using the function gc_inout. I used the function to identify grid points on a gaussian grid (512x256) within a polygon. When the polygon has a zonal extent larger than 180 degrees and extends over the equator (covers
partly both northern and southern hemispheres), a problem occurs; on the edges gridpoints are missing that should have been identified.
<div><br>
</div>
<div>I attached 2 simple idealized plots that demonstrate the issue; all gridpoints should be identified within the boxes:<br>
plot1) box lon-100:100, lat 10:30; everything fine<br>
plot2) box lon-100:100, lat -10:10; gridpoints on both edges are missing<br>
<br>
Is this a bug or a known issue as the description of gc_inout says: "and that the vertices can be contained in a hemisphere"?<br>
<br>
<div>Any suggestions how can I work around this issue?<br>
<br>
Looking forward to suggestions!<br>
Cheers,<br>
Andries</div>
</div>
</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>