【大神求助！】画风场的时间高度剖面图

vangie24

用18UTC的v减去一天内v的平均值代表风场日变化，在两广区域，把1.1-12.31（排除2.28的情况）共365天日变化都作59年（1958-2016）平均，再取区域平均，得出日变化的时间层次剖面图，下面是我写的gs文件，画出来的图跟实际情况很不符合，请大家帮忙找出问题！改过很多次都没办法，求大神指点！
*********************原始数据v.ctl*********************
dset ^%y4/v%y4%m2.dr
undef 9.999E+20
title anl_p125_vgrd
options yrev zrev template
ydef 145 linear -90.000000 1.25
xdef 288 linear 0.000000 1.250000
tdef 120000 linear 00Z01jan1958 6hr
zdef 37 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100 70 50 30 20 10 7 5 3 2 1
vars 1
v 37 34,100,0 ** v wind [m/s]
ENDVARS
*********************处理数据生成dat（以1、2月为例）*********************
'reinit'
'open v.ctl'

'set gxout fwrite'
imon.1='jan'
imon.2='feb'
'set lon 105 117.5'
'set lat 21 27.25'
'set z 1 27'
h=1
while(h<=2)
'set fwrite vdaily'imon.h'.dat'
if(h=1);day=31;endif
if(h=2);day=28;endif
id=1
while(id<=day)
m=1
while(m<=27)
'set z 'm''
yr=1958
while(yr<=2016)
'set time 06Z'id''imon.h''yr''
if(yr=1958);'define v1mean=0.0';endif
'define v1=v(t+2)-ave(v,t+0,t+3)'
'define v1mean=v1mean+v1'
yr=yr+1
endwhile
'define vd=v1mean/59.0'
'd vd'
m=m+1
endwhile
id=id+1
endwhile
'disable fwrite'
h=h+1
endwhile
*********************一月份的vdailyjan.dat为例写vdailyjan.ctl*********************
dset vdailyjan.dat
undef 9.999E+20
title anl_p125_vgrd
ydef 6 linear 21.000000 1.25
xdef 11 linear 105.000000 1.250000
tdef 31 linear 01jan1958 1dy
zdef 27 levels
1000 975 950 925 900 875 850 825 800 775 750 700 650 600 550 500 450 400 350 300 250 225 200 175 150 125 100
vars 1
vd 27 34,100,0 ** v wind [m/s]
ENDVARS
*********************画图的gs*********************
'reinit'
'open vdailyjan.ctl'
'set gxout shaded'
'color.gs -2 2 0.2 -kind dblue->mblue->sblue->white->yellow->sred->mred'
'set x 1'
'set y 1'
'set z 1 27'
'set t 1 31'
'define a=aave(vd,lon=105,lon=117.5,lat=21,lat=27.25)'
'd a'
'cbarn.gs'
'printim t1.png x1600 y1236 white'

'c'
********************图*********************

tangjie199365

请问 里面vd是什么意思 ？在wrf里面是u还是v？