ncl画风矢量图出现的问题

Cupid.

画风矢量图的时候，出现这个警告，不知道问题出在哪里，各位大神帮忙看一下吧

Cupid.

顶顶

Cupid.

顶顶

supersupergirl

只要图没出错 warning可以不理会

Cupid.

supersupergirl 发表于 2020-4-13 12:43
只要图没出错 warning可以不理会

就是感觉出错了，箭头的长度和阴影所表示的风速，不匹配

Cupid.

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"

begin
fu = addfile("2003082122.nc","r")
fv = addfile("2003082122.nc","r")
U = fu->u(10,6,:,:)
V = fv->v(10,6,:,:)
U1 =fu->u(10,6,:,:)
V1 =fv->v(10,6,:,:)
U2=short2flt(U1)
V2=short2flt(V1)
Wind = sqrt(U2*U2+V2*V2)
wks = gsn_open_wks("png","vector20")
res = True
res@gsnDraw = False
res@gsnFrame = False
res@gsnLeftString = ""
res@gsnRightString = ""

mpres = res
;mpres@gsnAddCyclic = False
mpres@mpFillOn                = True
mpres@mpOutlineOn             =True
mpres@pmTickMarkDisplayMode = "Always"
mpres@mpDataBaseVersion      ="MediumRes"
mpres@mpMinLatF = 37
mpres@mpMaxLatF= 47
mpres@mpMinLonF = -75
mpres@mpMaxLonF = -65

vcres                           = res
;vcres@gsnAddCyclic = False
vcres@vcGlyphStyle               = "CurlyVector"
vcres@vcRefLengthF               = 0.045
vcres@vcRefMagintudeF         = 1
vcres@vcRefAnnoOn               = False
;vcres@vcWindBarbLineThicknessF  = 2
;vcres@vcWindBarbColor = "black"
vcres@vcMinMagnitudeF = 0
vcres@vcMinDistanceF = 0.05
;vcres@vcWindBarbTickLengthF = 0.05

res@tiMainString = "wind"
res1=res
res1@cnFillOn               = True
res1@cnLinesOn               = False
res1@cnFillPalette               ="BlueRed"
res1@cnLevelSelectionMode    ="ManualLevels"
res1@cnMinLevelValF          ="0"
res1@cnMaxLevelValF          ="10"
res1@cnLevelSpacingF          ="0.25"
res1@gsnStringFontHeightF      =0.015
res1@gsnLeftString = "(a)2003-08-22 06:00:00"
res1@lbOrientation           ="vertical"
res1@pmLabelBarOrthogonalPosF =-0.03
res1@lbBoxMinorExtentF        =0.15
res1@lbLabelFontHeightF=0.01
res1@lbLabelFontColor  ="black"
res1@lbLabelFontThicknessF  =1.5
res1@lbLabelStride  =1
;res1@cnFillOnDrawOrder     ="PreDraw"
mapplot = gsn_csm_map(wks,mpres)
plot1 = gsn_csm_contour(wks, Wind, res1)
plot = gsn_csm_vector(wks,U,V,vcres)
overlay(mapplot,plot1)
overlay(mapplot,plot)
draw(mapplot)
frame(wks)
print("All Done.")
end

Cupid.

顶，大神们帮帮忙看一下吧

蹦蹦

37行vcres@vcRefMagintudeF应为vcres@vcRefMagnitudeF

lc2016

可不可以来张图瞧瞧啊，最近也在做类似的图。我有一点问题就是那个矢量图看你用的是plot = gsn_csm_vector(wks,U,V,vcres)，UV是short字型（有偏移）的吧，这个地方不应该用浮点字型吗，看看楼主怎么想的？我也是在学习ing。