NCL画水汽通量或风矢量时的箭头问题

timevigor

NCL画水汽通量或风矢量时的箭头问题。请问在设置箭头图例时是自动生成还是需要设置，我在运用相关命令后图中还是没有箭头的图例，且想要把箭头数值改为我计算的结果，但更改mangnitudeF后没有得到官网的类似结果。

lanlaniris

看了下感觉没什么问题啊……你单独画一下矢量肩头看下有没有出参考箭头

timevigor

lanlaniris 发表于 2015-12-23 11:06
看了下感觉没什么问题啊……你单独画一下矢量肩头看下有没有出参考箭头

单独画了也没有参考箭头，我还想是不是被bar遮住了，去掉bar也没有

海盗船长

vcres@vcLabelsOn                = True  

lanlaniris

你完整的脚本贴一个上来吧

timevigor

lanlaniris 发表于 2015-12-23 16:11
你完整的脚本贴一个上来吧

;dd1d2
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/cnmap/cnmap.ncl"

;*****************************************************************************
begin
;****************************************************************************
;  read in data
;******************************************************************************
  f1        = addfile("/home/almond/uq25febave.nc","r")
  f2        = addfile("/home/almond/vq25febave.nc","r")
  d1        = f1->qu
  d2        = f2->qv
  ;mm        = sqrt(d1^2+d2^2)                            全风场
  ;exit
;***********************************************
;  creat plot
;***********************************************
  wks = gsn_open_wks("png" ,"q25feb")
  gsn_define_colormap(wks,"BlRe")
;================================================
  res                                = True            
  res@gsnMaximize                    = True
  res@gsnDraw                        = False
  res@gsnFrame                       = False
  ;res@cnFillOn                       = True
  ;res@cnFillPalette                  ="BlWhRe"
;>--------------------------------------------<
;  set for the map
;>--------------------------------------------<
  res@mpMinLatF                      = 0.                        
  res@mpMaxLatF                      = 35.
  res@mpMinLonF                      = 80.
  res@mpMaxLonF                      = 150.
  res@cnLevelSelectionMode           ="ManualLevels"
  res@cnMinLevelValF                 =-300
  res@cnMaxLevelValF                 =300
  res@cnLevelSpacingF                =50
;>--------------------------------------------<
;  set for plot (fillcontour)
;>--------------------------------------------<

  res@lbLabelBarOn                   = True
  res@cnFillOn                       = True  
  
                                          

  res@pmTickMarkDisplayMode          = "Always"                       
  res@gsnAddCyclic                   = False

  res@cnLineLabelPlacementMode       = "Computed"
  res@cnLineLabelDensityF            = 2.0
  res@cnLabelMasking                 = True
  res@cnLineLabelBackgroundColor     = "White"
  ;;res@cnLineLabelFormat             = "@^sg"
  
  res@cnLinesOn                      = True           
  res@cnLineLabelsOn                 = True
  res@cnLineLabelInterval            = 1
  res@cnLineLabelFontHeightF         = 0.01

  plot = gsn_csm_contour_map_ce(wks,d1,res)

;>============================================================<
vcres                             = True
vcres@gsnAddCyclic                = False
vcres@gsnDraw                     = False                                          ; don't draw yet
vcres@gsnFrame                    = False                                          ; don't advance frame yet
vcres@pmTickMarkDisplayMode       = "Always"
vcres@vcLabelsOn                = True                                                   
;vcres@lbLabelBarOn                = True
vcres@vcRefMagnitudeF             = 10.0                                            ; make vectors larger
vcres@vcRefLengthF                = 0.0018                                          ; ref vec length
vcres@vcGlyphStyle                = "CurlyVector"                                   ; turn on curly vectors
vcres@vcMinDistanceF              = 0.017                                           ; thin out vectors
vcres@vcRefAnnoOn                 = False
vcres@vcMonoFillArrowFillColor    = True
vcres@vcLineArrowThicknessF       = 2.0
vcres@vcRefAnnoOn                 =False
vcres@vcVectorDrawOrder           = "PostDraw"
;vcres@tiMainString                ="vaper_flux : kg/(m*s)"
;vcres@gsnLeftString               = ti
vcres@gsnRightString              ="uwf"
plot1=gsn_csm_vector(wks,d1,d2,vcres)
;-------------------------------------------------------------<
;                      add China map
;>------------------------------------------------------------<
cnres                               = True
cnres@china                         = True        ;draw china map or not
cnres@river                         = False       ;draw changjiang&huanghe or not
cnres@province                      = True        ;draw province boundary or not
cnres@nanhai                        = False       ;draw nanhai or not
cnres@diqu                          = False       ;draw diqujie or not

chinamap                            = add_china_map(wks,plot1,cnres)
  overlay(plot,plot1)
  draw (plot)
  frame(wks)
end

timevigor

海盗船长 发表于 2015-12-23 14:07
vcres@vcLabelsOn                = True

试了并没有

centaurus

可能是vcRefAnnoOn=True才会有吧，你设置的False

timevigor

lanlaniris 发表于 2015-12-23 16:11
你完整的脚本贴一个上来吧

vcRefAnnoOn=True,确实是这个问题，多谢

timevigor

lanlaniris 发表于 2015-12-23 16:11
你完整的脚本贴一个上来吧

vcRefAnnoOn=True,确实是这个问题，多谢