Re: gc_clkwise failure

I have filed a JIRA ticket, NCL-1699

--Mary

On Apr 16, 2013, at 8:13 AM, Dennis Shea wrote:

> 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

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Apr 17 08:07:01 2013