画经度和气压面的风速图时遇到的问题

X1a0wu

风速风量为u，和垂直风速w, 无论怎么改W的速度的大小，都无法在图上反映出垂直速度来
  u=in->u(32,{950:100},{30:40},{110:120})
    v=in->v(32,{950:100},{30:40},{110:120})
    w=in->w(32,{950:100},{30:40},{110:120})
    u1 = u(:,{34.5},:)
    w1 = w(:,{34.5},:)*10
printVarSummary(u)
    q1 =  in->q(32,{950:100},{30:40},{110:120})
    uq1 = q1*u/g
    vq1 = q1*v/g

    lat1 = in->lat({30:40})
    lon1 = in->lon({110:120})

    plot=new(15,graphic)
    scl = 1.e7
    ;---Divergence of moisture flux: uv2dvF => global 'fixed' rectilinear grid

    duvq1  = uv2dv_cfd(uq1, vq1,lat1,lon1,0) *scl   ; (level,lat,lon)
   ;duvq1@long_name = "Divergence of Moisture Flux"
   ; duvq1@units     = "kg/(kg-s)"    ; (1/m)*[(m/s)(g/kg)] => [g/(kg-s)]
    copy_VarCoords(u,duvq1)
    printMinMax(duvq1,0)
    duvq2 = duvq1(:,{34.5},:)
     res = True

  wks = gsn_open_wks("pdf", "w_u_dquv")
    res@gsnDraw           = False        ; don't draw yet
    res@gsnFrame          = False        ; don't advance frame yet
    res@mpFillOn          = False        ; no grey continents
    ; res@cnInfoLabelOn     = False
    ; res@lbLabelBarOn      = False
    res@cnLineLabelsOn    = False
    res@tmXTOn            = False
    res@tmYROn            = False
    res@cnLineThicknessF  = 1.5
    res@cnFillOn   = True
    res@cnLinesOn = False
    res@lbLabelBarOn = False



    gsn_define_colormap(wks,"temp_diff_18lev")
     ;gsn_reverse_colormap(wks)
    ;       set resource for contour plot
    ;>---------------------------------------<
    res@cnSmoothingOn = True ;
    res@cnSmoothingDistanceF   = 0.001
    res@cnLevelSelectionMode   = "ExplicitLevels"
    plot_levels=ispan(-4,4,1)
    res@cnLevels=plot_levels
    res@lbLabelStrings= tostring(plot_levels)
    res@cnInfoLabelOn = False
    res@cnLineLabelsOn          =  False
    res@cnLineLabelFontHeightF = 0.01
    res@cnLineThicknessF = 2.0
    res@cnLineLabelInterval = 1
    res@lbBoxEndCapStyle = "TriangleBothEnds"
    res@gsnScalarContour= True
    ;     set resource for windBarb plot
    ;>---------------------------------------<

    res@gsnAddCyclic = False

    dres= True
        ;dres@gsnDraw           = False        ; don't draw yet
    ;dres@gsnFrame          = False
    dres@vcGlyphStyle               = "CurlyVector"   ;"WindBarb"
    dres@vcRefLengthF               = 0.03
    dres@vcRefAnnoOn               =   True
    dres@vcRefAnnoOrthogonalPosF   = -1.0  
    dres@vcRefMagnitudeF              = 40
    dres@vcRefAnnoString1On= False
    dres@vcRefAnnoString2  ="10m/s"
    dres@vcRefAnnoString2On = True
    ;res@vcWindBarbLineThicknessF  = 1.5
    ;res@vcWindBarbColor = "black"
    dres@vcMinMagnitudeF = 50
    ;res@vcWindBarbTickLengthF = 0.35
    dres@vcMinDistanceF        = 0.05  
    dres@trYReverse        = True

;   dres@tmYLMode              = "Explicit"
;   dres@tmYLValues            = (/90000,80000,50000,30000,10000/)
;   dres@tmYLLabels            = "" + res@tmYLValues/100

    res@trYReverse        = True

    res@tmYLMode              = "Explicit"
    res@tmYLValues            = (/900,800,500,300,100/)
    res@tmYLLabels            = "" + res@tmYLValues

    res@gsnLeftString=name(0)
    plot_uv = gsn_csm_vector(wks,u1,w1,dres)
    plot(0)=gsn_csm_contour(wks,duvq2,res)

    overlay(plot(0),plot_uv)