关于 gsn_csm_pres_hgt_vector函数中data属性的问题

追风的阳光

附上脚本，如图，dbz色标有问题，想更改色标，谢谢！
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/wrf/WRFUserARW.ncl"

begin
;
; The WRF ARW input file.  
; This needs to have a ".nc" appended, so just do it.
a=addfile("/home/zssapr/WRF3.5.1/20110710/WRFV3/test/em_real/wrfout_d02_2011-07-05_00:00:00.nc","r")


; We generate plots, but what kind do we prefer?
  type = "x11"
; type = "pdf"
;  type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_CrossSection_vector")
  gsn_define_colormap(wks,"wh-bl-gr-ye-re")
; Set some basic resources
  res = True
;res@MainTitle = "REAL-TIME WRF"
  res@Footer = False
  pltres = True
  ter_res = True
  opts_ter = ter_res
  opts_ter@gsnYRefLine = 0.0
  opts_ter@gsnAboveYRefLineColor = "black"
  opts_ter@gsnDraw = False
  opts_ter@gsnFrame = False
  ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  mdims = getfilevardimsizes(a,"P") ; get some dimension sizes for the file
  nd = dimsizes(mdims)

;---------------------------------------------------------------
  ter  = wrf_user_getvar(a,"ter",0)  
  xlon = wrf_user_getvar(a, "XLONG",0)

   it=161
    tc  = wrf_user_getvar(a,"tc",it)     ; T in C
    rh = wrf_user_getvar(a,"rh",it)      ; relative humidity
    z   = wrf_user_getvar(a, "z",it)     ; grid point height
    dBz =wrf_user_getvar(a,"dbz", it)
    U  = wrf_user_getvar(a,"U",it)     ; U wind
    V = wrf_user_getvar(a,"V",it)      ; V wind
    W :=wrf_user_getvar(a,"W", it)     ;vertical wind
    p = wrf_user_getvar(a, "pressure",it)
    ua =wrf_user_unstagger(U,U@stagger)
    va =wrf_user_unstagger(V,V@stagger)
    wa :=wrf_user_unstagger(W,W@stagger)

    plane = (/ mdims(nd-1)/2, mdims(nd-2)/2 /)    ; pivot point is center of domain (x,y)
    angle=45
     opts = False  
     ua_plane = wrf_user_intrp3d(ua,p,"v",plane,angle,opts)
     va_plane = wrf_user_intrp3d(va,p,"v",plane,angle,opts)
     w_plane =wrf_user_intrp3d(wa,p,"v",plane,angle,opts)
     rh_plane = wrf_user_intrp3d(rh,p,"v",plane,angle,opts)
     dBz_plane = wrf_user_intrp3d(dBz,p,"v",plane,angle,opts)
     tc_plane = wrf_user_intrp3d(tc,p,"v",plane,angle,opts)
     p_plane = wrf_user_intrp3d(p, p,"v",plane, angle,opts)
     ter_plane = wrf_user_intrp2d(ter,plane,angle,opts)
     X_plane = wrf_user_intrp2d(xlon,plane,angle,opts)

     wAve   = avg(w_plane(:,:))           ; used for scaling
     vAve   = avg(va_plane(:,:))
     scale  = fabs(vAve/wAve)
     wscale = w_plane*scale                       ; now scale

     tc_plane!0="level"
     tc_plane&level = p_plane(:, 0)
     copy_VarCoords(tc_plane, va_plane)
     copy_VarCoords(tc_plane, wscale)
      ;printVarSummary(p_plane)
      ;print(p_plane(:,0))
     dBz_plane!0="level"
     dBz_plane&level = p_plane(:, 0)
     copy_VarCoords(dBz_plane, va_plane)
     copy_VarCoords(dBz_plane, wscale)
;***********************************************
        ; create plot
        ;***********************************************
        res                 = True                    ; plot mods desired
        res@tiMainString    = "Pressure/Height Vector" ; title
        res@tiYAxisString           = "Pressure (mb)"
        res@cnLineLabelsOn  = False              ; turn off line labels
        res@cnFillOn        = False               ; turn on color fill
        res@cnLineColor = "white"
        res@lbLabelStride   = 2                  ; every other color
        ;res@ContourParameters       = (/ 0., 50., 2. /)
        res@gsnSpreadColors = True               ; use full range of color map
        res@cnMissingValPerimOn     = True
        res@cnMissingValFillColor   = 0
        res@cnMissingValFillPattern = 11
        res@tmYLMode                = "Explicit"

        res@gsnSpreadColors = True               ; use full range of color map

        res@vcRefMagnitudeF = 3.0                ; define vector ref mag
        res@vcRefLengthF    = 0.045              ; define length of vec ref
        res@vcGlyphStyle    = "CurlyVector"      ; turn on curley vectors
        res@vcMinDistanceF  = 0.03               ; thin out vectors
        res@vcMapDirection  = False
     plot_vector  = gsn_csm_pres_hgt_vector(wks,tc_plane(2:,:),va_plane(2:,:),wscale(2:,:),res )

;**************************************************
     ;dBz plot
；    opts_dBz=True
；    opts_dBz@cnFillOn        = True               ; turn on color fill
  ；    opts_dBz@lbLabelStride   = 2                  ; every other color
  ；      opts_dBz@ContourParameters       = (/ 0., 50., 2. /)
  ；      opts_dBz@cnMissingValFillColor   = 0
   ；     opts_dBz@cnMissingValFillPattern = 11
   ；     opts_dBz@gsnDraw = False
   ；     opts_dBz@gsnFrame = False
    ；    contour_dBz = wrf_contour(a,wks,dBz_plane(2:,:),opts_dBz)

    ；    plot = wrf_overlays(a,wks,(/plot_vector,contour_dBz/),True)

end

longlivehj

楼主，到底啥问题啊？
标题所说的函数参数属性问题？overlay的问题？色标的问题？
从你给出的图来看，啥问题也没有！

追风的阳光

longlivehj 发表于 2014-10-14 22:11
楼主，到底啥问题啊？
标题所说的函数参数属性问题？overlay的问题？色标的问题？
从你给出的图来看，啥 ...

我想修改dBz填色的色标，即色标从dBz=0开始。我在gsn_csm_pres_hgt_vector的res中无论加入res@ContourParameters       = (/ 0., 50., 2. /)或是res@cnLevelSelectionMode ="Manual"， res@cnLevelSpacingF =2，res@cnMinLevelValF = 0.，res@cnMaxLevelValF =50.都没有用。于是另写contour_dBz想新赋属性与plot_vector叠加（overlays）发现没有用，contour_dBz不显示，报错

longlivehj

追风的阳光 发表于 2014-10-15 09:14
我想修改dBz填色的色标，即色标从dBz=0开始。我在gsn_csm_pres_hgt_vector的res中无论加入res@ContourPar ...

1. ContourParameters是wrf绘图函数里面的属性，在gsn_csm_pres_hgt_vector中使用当然没有作用。
2. cnLevelSelectionMode的方式仍然没有效果，看看是不是dBz数据本身的原因，比如为constant，或者不在指定等值线范围内等等，最好用printMinMax查看一下。
3. 将gsn_csm_*出来的图和wrf_*出来的图用wrf_overlays叠加，难免会出现问题。尽量只选择一套系统，就是要么全部用gsn_csm_*，用overlay叠加；要么全部用wrf_*，用wrf_overlays叠加。

追风的阳光

longlivehj 发表于 2014-10-15 09:41
1. ContourParameters是wrf绘图函数里面的属性，在gsn_csm_pres_hgt_vector中使用当然没有作用。
2. cnL ...

dbz范围是在等值线范围内的，-30到39dbz。dbz单独用gsn_csm_contour画再overlay叠加，依旧不显示。
contour_dBz = gsn_csm_contour(wks,dBz_plane(2:,:),opts_dBz)
plot_vector  = gsn_csm_pres_hgt_vector(wks,tc_plane(2:,:),va_plane(2:,:),wscale(2:,:),res )   
overlay(contour_dBz,plot_vector)                       ; overlay the U-wind plot on the temperature plot
draw(contour_dBz)出来的图形是这个：

等值线依旧无法设置（opts_dBz@cnLevelSpacingF =2有用， opts_dBz@cnMinLevelValF = 0
        opts_dBz@cnMaxLevelValF =50无用），dbz图像倒了个。。。。
如果是
overlayplot_vector,contour_dBz)                    
draw(plot_vector)
那么图形只有vector，没有dBz图像