请教一下求平均水汽通量中风场平均怎么求

白鹿青崖

毕业论文的时候求水汽通量遇到了问题，用典型月份代表的。后来想了很久，也和同学讨论过，但是在求多年平均的时候风场出现了问题。没学过GrADS，都是自己看的手册，难免疏漏，望各位不吝赐教，批评指正，不胜感激。
reinit
sdfopen E:\shum.mon.mean.nc
sdfopen E:\uwnd_mon_mean.nc
sdfopen E:\vwnd_mon_mean.nc
enable print E:\SHUIQI850.gmf
set grads off
set grid off
set lev 850
set lon 114 136
set lat 38 56
define q=ave(shum.1,t=7,t=816,12)
求风场平均(uwnd.2,vwnd.3应为和shum。1对应的平均风场)
define a=sqrt(uwnd.2*uwnd.2+vwnd.3*vwnd.3)
define q1=q/(1+q)
set cmin 0.1
set cint 0.05
d a*q1/9.8
set gxout vector
run cbarn.gs
print
disable print

river

什么问题，有提示否？

白鹿青崖

river 发表于 2016-7-23 08:18
什么问题，有提示否？

没有错误提示，就是平均风场那个地方求不出来了。画出的图与实际情况不符，所以想请教一下平均风场究竟应该怎么求。多谢！

白鹿青崖

river 发表于 2016-7-23 08:18
什么问题，有提示否？

又奋斗了一天，好像风场没啥问题了，可是水汽数值上好像不太对
sdfopen E:\uwnd_mon_mean.nc
sdfopen E:\vwnd_mon_mean.nc
enable print E:\SHUIQIjijiepingjun.gmf
set grads off
set grid off
set mpdset cnworld
set lon 30 160
set lat 0 80
set lev 850
define q1=ave(shum.1,t=6,t=816,12)
define q2=ave(shum.1,t=7,t=816,12)
define q3=ave(shum.1,t=8,t=816,12)
define q=q1+q2+q3
define qu=uwnd.2*q/9.8
define qv=vwnd.3*q/9.8
d qu*6;qv*6
print
disable print

白鹿青崖

river 发表于 2016-7-23 08:18
什么问题，有提示否？

对不起，新手不太会用，刚才弄错了，程序是这个……
reinit
sdfopen E:\shum.mon.mean.nc
sdfopen E:\uwnd_mon_mean.nc
sdfopen E:\vwnd_mon_mean.nc
enable print E:\SHUIQIx1000,9.gmf
set grads off
set grid off
set mpdset cnworld
set gxout shaded
set lon 114 136
set lat 38 56
set lev 1000
set cmin 0.1
set cint 0.05
define shum=ave(shum.1,t=9,t=816,12)
define u=ave(uwnd.2,t=9,t=816,12)
define v=ave(vwnd.3,t=9,t=816,12)
define a=sqrt(u*u+v*v)
define q1=shum/(1+shum)
d a*q1/9.8
d u;v
run cbarn.gs
print
disable print