ncl作风矢量的垂直剖面图，不知道问题出在哪里

irides

想画某一垂直剖面的风矢量图。为了合成矢量时 垂直风速能很好地体现出来，在画矢量图时将w扩大了10倍，但画出的图不太对。


ncl脚本：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
begin
a = addfile("cntl.nc","r")
type = "png"
wks = gsn_open_wks(type,"21")
res           = True
time = 37
u = wrf_user_getvar(a,"ua",time)              ; u averaged to mass point
v = wrf_user_getvar(a,"va",time)              ; v averaged to mass point
w = wrf_user_getvar(a,"wa",time)              ; w averaged to mass point
z = wrf_user_getvar(a,"z",time)

;----------------------------插值-------------------------------------------；
startlon = 109.68
startlat = 19.91
endlon   = 110.46
endlat   = 19.24
xy = wrf_user_ll_to_ij(a,(/startlon,endlon/),(/startlat,endlat/),res)
plane  = new(4,float)
plane  = (/ xy(0,0),xy(1,0),xy(0,1),xy(1,1) /)       ;AB经纬度对应网格点
alfa   = atan2((endlat-startlat),(endlon-startlon))    ;计算线段AB与x轴方向的夹角
uv     = u*cos(alfa)+v*sin(alfa)

uvcross = wrf_user_intrp3d(uv,z,"v",plane,alfa,True)      ;xy合成量插值到高度上
wcross  = wrf_user_intrp3d(w,z,"v",plane,alfa,True)       ;w插值到高度上
;-----------------------------画图-------------------------------------------；
vcres                 = True
vcres@vcGlyphStyle    = "LineArrow"
vector = wrf_vector(a,wks,uvcross,(wcross*10),vcres)
cnres                 = True
cnres@cnFillOn        = True
cnres@cnLinesOn       = False
contour = wrf_contour(a,wks,uvcross,cnres)
plot = wrf_overlays(a,wks,(/contour,vector/),True)
end

之前用grads画的图是正确的，二者对比明显是我的ncl脚本的问题。

有木有好心的大神 请教一下哪一步出了问题呢？

xyan88

请问最后怎么解决的啊，可以提供下代码学习吗

freekiller

貌似ncl可以直接出来剖面插值

irides

freekiller 发表于 2017-3-28 22:50
貌似ncl可以直接出来剖面插值

额，不知是一个脚本还是函数呢？

万川一月

w向风为什么要乘以10，不懂

irides

万川一月 发表于 2017-7-15 15:14
w向风为什么要乘以10，不懂

因为与水平方向风速比较的话，垂直方向风速太小了，为了体现垂直方向的运动，会将垂直运动速度扩大

万川一月

噢噢，明白了

halakaka

你好 请问解决了吗？

irides

halakaka 发表于 2017-11-13 17:27
你好 请问解决了吗？

可以参考http://www2.mmm.ucar.edu/wrf/OnL ... f_CrossSection2.htm

halakaka

irides 发表于 2017-11-14 10:28
可以参考http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/CROSS_SECTION/wrf_CrossS ...

嗯嗯~谢啦