求助！！用gsn_csm_vector_map画WRF风场

Rita

求教各位大侠！不用wrf_vector， 使用gsn_csm_vector_map画WRF的10米风场，风向杆画不出来，这是咋回事啊？

脚本如下：

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/wrf/WRFUserARW.ncl"
begin
  a = addfile("./wrfout_d01_2011-07-22_10:15:00.nc","r")
  figname = "wrf_spatial_uv10_barb_urban_1815LST_4"
  type = "pdf"
  u10Urban = a->U10(0,:,:)
  v10Urban = a->V10(0,:,:)
  hgt = a->HGT(0,:,:)
  u10Urban@lat2d = a->XLAT(0,:,:)
  u10Urban@lon2d = a->XLONG(0,:,:)
  u10Urban@units = "m/s"
  u10Urban@description = "10m U"
  copy_VarMeta(u10Urban,v10Urban)
  v10Urban@description = "10m V"
  copy_VarMeta(v10Urban,hgt)
  hgt@units = "m"
  hgt@description = "topography"
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  wks = gsn_open_wks(type,figname)
; Set some basic resources
  res = True
  res@mpLimitMode = "LatLon"
  res@mpMinLatF = 39.6
  res@mpMaxLatF = 40.21
  res@mpMinLonF = 116.0
  res@mpMaxLonF = 116.8
;    res@tmXTOn = True
;    res@tmYROn = True
;    res@tmXTLabelsOn = True
;    res@tmYRLabelsOn = True
  res@gsnAddCyclic            = False   ; regional data
  res = res
  res@vcGlyphStyle = "WindBarb"         ; choose wind barbs
  res@vcMonoWindBarbColor = True       ; color barbs by scalar
  res@vcWindBarbScaleFactorF=10
; res@vcMinDistanceF = 0.025            ; thin out windbarbs
  res@vcWindBarbTickLengthF = 0.4       ; length of the barb tick
  res@vcWindBarbTickSpacingF = 0.2      ; distance between ticks
  res@vcRefMagnitudeF         = 3.                ; make vectors larger
  res@vcRefAnnoOn = False               ; turn off ref wind barb
  plot = gsn_csm_vector_map(wks,u10Urban,v10Urban,res)
  draw(plot)
  frame(wks)
end

Rita

数据print出来没有问题，用wrf_vector可以正确出图，为什么gsn_csm_vector_map不行呢？

zly4814624

gsn_csm_vector

zly4814624

res@gsnDraw               = False
res@gsnFrame             = False
你把这两句加上去试试看。