ncl画垂直剖面图报错

curious

各位大佬，请问用以下画垂直风速和相对湿度的垂直剖面图出现这样的报错是什么原因呢，还有无其他方法去画垂直风速的垂直剖面，官网的方法没咋看懂
代码：
filename = "d:/ncldata/nc/era5"
    f = addfile(filename + ".nc", "r")
   rh=f->r(1,:,:,:)
   v = f->v(1,:,:,:)
   omega=f->w(1,:,:,:)
    level=rh&level
    lev4d=conform_dims(dimsizes(omega),level,2)
    w=omega_to_w(omega,lev4d,rh)
    w=w*100
    wks = gsn_open_wks("png", "d:/ncldata/poumianfnl_20211001_00_00" )
    res = True
    res@gsnDraw = False
    res@gsnFrame = False
    res@gsnLeftString = "relative humidity"
    res@cnFillOn=True
    res@cnLinesOn=False
    res@cnLevelSpacingF = 10
    res@cnFillPalette="GMT_drywet"
    res@cnInfoLabelOn = False
    res@vcMapDirection=False
    res@vcPositionMode="ArrowTail"
    res@vcRefMagnitudeF=2
    res@vcRefLengthF=0.04
    ;绘制坐标标签
    res@tmXBMode="Explicit"
    res@tmXBValues=(/30,35,40,45/)
    res@tmXBLabels=(/"30~S~o~N~N","35~S~o~N~N","40~S~o~N~N","45~S~o~N~N"/)
    plot=gsn_csm_pres_hgt_vector(wks,rh(1,{1000:10},30:45,{113}),v(1,{1000:10},30:45,{113}),w(1,{1000:10},30:45,{113}),res)

                               
登录/注册后可看大图