Re: polygon problem

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Mon, 26 Jan 2009 07:56:30 -0700 (MST)

This problem has been fixed and will be in the 5.1.0 release. It had
to do with floating point precision in some low-level map code.

Michael, if you need the fix now, email me separately with your
system information.

--Mary

On Wed, 21 Jan 2009, Mary Haley wrote:

> Hi Michael,
>
> I was able to use this script and narrow it down a little bit. Fred
> Clare is going to look into it. We think maybe the ordering of the
> points may be the problem.
>
> One of us will post back here when we know more.
>
> --Mary
>
> On Wed, 21 Jan 2009, Michael Notaro wrote:
>
>> This is a simpler version of the code. It just gives a map fully covered
>> in black.
>> For some reason, that set of lats/lons results in a problem.
>>
>> Mike
>>
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
>> begin
>>
>> wks = gsn_open_wks ("ps", "test")
>> plot=new(1,graphic)
>> setvalues wks
>> "wkColorMap" : "BlWhRe"
>> end setvalues
>> res = True
>> res_at_gsnFrame=False
>> res_at_tiMainFontHeightF = 0.033
>> res_at_tiMainString = "Plant Distribution"
>> res_at_tmYROn = False
>> res_at_tmXTOn = False
>> res_at_mpOceanFillColor = -1
>> res_at_mpLandFillColor = -1
>> res_at_mpMinLatF = 10.
>> res_at_mpMaxLatF = 75.
>> res_at_mpMinLonF = -180.00
>> res_at_mpMaxLonF = -50.00
>> res_at_mpCenterLonF = -115.
>> res_at_mpOutlineBoundarySets = "AllBoundaries"
>> res_at_mpGeophysicalLineThicknessF = 2.5
>> res_at_mpOutlineOn = True
>> plot = gsn_csm_map(wks,res)
>>
>> lat=new((/8/),double)
>> lon=new((/8/),double)
>> lat=(/
>> 54.4126815796,54.4171485901,54.4163818359,54.4119338989,54.4063796997,54
>> .4074935913,54.4108200073,54.4126815796/)
>> lon=
>> (/-130.6879730225,-130.674041748,-130.6725158691,-130.6683349609,-130.66
>> 88842773,-130.6794433594,-130.686126709,-130.6879730225/)
>>
>> lat!0="lat"
>> lat_at_units="degrees_north"
>> lat&lat=lat
>> lon!0="lon"
>> lon_at_units="degrees_east"
>> lon&lon=lon
>>
>> res1=True
>> res1_at_gsFillScaleF = 1.0
>> res1_at_gsFillColor = "forestgreen"
>> gsn_polygon(wks,plot,lon,lat,res1)
>>
>> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 26 2009 - 07:56:30 MST

This archive was generated by hypermail 2.2.0 : Wed Jan 28 2009 - 09:24:27 MST