Re: polygon problem

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
Received on Wed Jan 21 2009 - 13:06:05 MST