NCL 风玫瑰图warning

Gabriella

画风玫瑰图，但是总是会有warning:TransformPostDraw: tfPolyDrawList element 5 is invalid,然后导致风玫瑰图中有些线不显示，比如鞋画了一半就没有了。
比如图一正北方向画了一半就没有了。
脚本如下，请大家帮忙看一下是什么问题，谢谢。


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/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"

begin
wks_type="png"
;wks_type@wkWidth=4500
;wks_type@wkHeight=4500
res                         = True
res@gsnDraw                 = True
res@gsnFrame                = True
;8aws的风向panel
;文件名循环打开 先把8站文件名给name
fname ="D:/WR/winterave/id/8id.txt"
name  = asciiread(fname,-1,"string")

do ii=0,7
wks     =gsn_open_wks(wks_type,"D:/WR/wind_rose/winter3months/8aws/winter_"+name(ii))
wrname="D:/WR/awsnew/wind/dd2ff28aws/"+name(ii)+".txt"
wrData=asciiread(wrname,-1,"string")
wspd=tofloat(str_get_field(wrData,6," "))
wdir=tofloat(str_get_field(wrData,5," "))
wspd@longname="wind speed"
wspd@units   ="m/s"
wdir@longname="wind direction"
numPetals  =16
circFr     =5   ;频数圈间隔
spdBounds  =(/0.5,1.5,3.,5.,8.,15./)*1.0
colorBounds=(/"black","blue","green","red","yellow","orange"/)
res@tiMainString ="wind rose of "+name(ii)

lbres                    = True
lbres@lbAutoManage       = False         ; Necessary to control sizes
lbres@vpWidthF           = 0.05          ; labelbar width
lbres@vpHeightF          = 0.15          ; labelbar height
lbres@vpXF               = 0.78          ; labelbar position in x direction
lbres@vpYF               = 0.58          ; labelbar position in y direction
lbres@lbBoxMajorExtentF  = 0.80          ; puts space between color boxes  样条间隔
lbres@lbFillColors       = colorBounds   ; labelbar colors        
lbres@lbMonoFillPattern  = True          ; Solid fill pattern
lbres@lbLabelFontHeightF = 0.01         ; font height. default is small
lbres@lbPerimOn          = False         ;关闭样条边框
lbres@lbTitleString      ="Wind Speed(m/s)"
lbres@lbTitleFontHeightF = 0.01
;绘制标签
nbar   = dimsizes(spdBounds)
labels = new(nbar,string)
labels(0) = 0 + "-" + spdBounds(0)
do i=1,nbar-1
     labels(i) = spdBounds(i-1) + "-" +spdBounds(i)
end do
lbid = gsn_create_labelbar(wks,nbar,labels,lbres)

draw(lbid)
plot = WindRoseColor(wks,wspd,wdir,numPetals,\
                        circFr,spdBounds,colorBounds,res)
delete(wrData)  
delete(wspd)
delete(wdir)                       

end do
frame(wks)
end