Re: gc_clkwise failure

A slight variation ..... If the following is done
    lats = (/ 15, 40, 0, -30, -30 /)
    lats = lats + 30

Both return False!

On 4/15/13 6:06 PM, Dave Allured - NOAA Affiliate wrote:
> NCL team,
>
> The following demonstrates two large simple polygons over the Pacific
> ocean. Both are written in clockwise order, yet the function
> gc_clkwise returns the wrong value for the first polygon.
>
> lats = (/ 15, 40, 0, -30, -30 /)
> lons = (/ 150, 210, 260, 265, 160 /)
> func = gc_clkwise (lats, lons)
> print ("Result #1 = " + func)
>
> lons(3) = 250
> func = gc_clkwise (lats, lons)
> print ("Result #2 = " + func)
>
> mac56:~/ncl/bugs/gc_clkwise 63> ncl demo.gc_clkwise.ncl
> NCAR Command Language Version 6.1.2
> (0) Result #1 = False
> (0) Result #2 = True
>
> I believe that both polygons conform to the documented constraint that
> each must be containted in a hemisphere. We get the same behavior
> with NCL 6.0.0 and 6.1.2 on Mac OS 10.8.3, 64-bit x86 CPU. Am I
> mistaken, or is this a bug in the function? Thanks for looking into
> this.
>
> --Dave
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Apr 16 08:14:03 2013