Re: possible bug in gsn_polygon

Hi Michael,

This is definitely strange.

I was also able to get the script to work by decreasing the number of
significant digits:

   cx(0) = -2.35
   cx(1) = -4.62
   cx(2) = -2.27
   cx(3) = -2.35

   cy(0) = 65.9133
   cy(1) = 65.4328
   cy(2) = 64.9149
   cy(3) = 65.9133

I'll get this on the NCL developers' list to see if they have any
ideas on this.

--Mary

On Thu, 30 Apr 2009, Michael Duda wrote:

> Hi, Folks,
>
> I'm using NCL 5.1.0 on NCAR's bluefire (although the problem appears on other
> machines, too), and I'm running into a error when using gsn_polygon(...) to add
> a polygon to a map. The error message I see is:
>
> Copyright (C) 1995-2009 - All Rights Reserved
> University Corporation for Atmospheric Research
> NCAR Command Language Version 5.1.0
> The use of this software is governed by a License Agreement.
> See http://www.ncl.ucar.edu/ for more details.
> fatal:MapDataPolygon: ARPRAM - ALGORITHM FAILURE
>
> However, if I just changing the center lat or lon of the map (mpCenterLatF or
> mpCenterLonF) by even one degree, the problem disappears and the polygon appears
> on the plot as expected. Appended below is the script I'm using. Could anyone,
> especially an NCL developer, comment on whether this is a bug? Any other
> suggestions on a workaround would also be great.
>
> Best regards,
> Michael
>
> --------------------------------------------------------
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>
> begin
>
> wks = gsn_open_wks("pdf","triangle")
> gsn_define_colormap(wks,"gui_default")
>
> res = True
> res_at_gsnMaximize = True
>
> res_at_mpProjection = "Satellite"
> res_at_mpDataBaseVersion = "MediumRes"
> res_at_mpCenterLatF = 45.0
> res_at_mpCenterLonF = 0.0
> res_at_gsnFrame = False
>
> map = gsn_csm_map(wks, res)
>
> lres = False
>
> cx = new((/4/), float)
> cy = new((/4/), float)
>
> cx(0) = -2.356872735192781
> cx(1) = -4.627934296211549
> cx(2) = -2.270649143374328
> cx(3) = -2.356872735192781
>
> cy(0) = 65.91331700014982
> cy(1) = 65.43284608178955
> cy(2) = 64.91495523013631
> cy(3) = 65.91331700014982
>
> lres = True
> lres_at_gsFillColor = 2
> gsn_polygon(wks, map, cx, cy, lres)
>
> frame(wks)
> end
>
> _______________________________________________
> 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 Thu Apr 30 2009 - 15:29:49 MDT