求指导，关于水汽通量剖面图

咿呀呀

想做一个水汽通量垂直剖面图，但是总运行错误。
代码如下：
'reinit'
'set grid off'
'set timelab off'
'set grads off'
'sdfopen e:\201409\air.2014.nc'
'sdfopen e:\201409\rhum.2014.nc'
'sdfopen e:\201409\shum.2014.nc'
'sdfopen e:\201409\uwnd.2014.nc'
'sdfopen e:\201409\vwnd.2014.nc'
'set t 992 1044'
'set lev 1001 101'
'define p=lev'
'define e=(6.112*exp((17.67*(air.1-273.16))/(air.1-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhum.2/100'
'define qdiv=hdivg(q*uwnd.4/9.8,q*vwnd.5/9.8)'
'set lon 108'
'set lat 34'
'd qdiv'
'enable print e:\201409\1\qdiv.gmf'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set gxout shaded'
'set clab on '
'set cint 0.3'
'set cmax 0'
'set rbcols auto'
'set csmooth on'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'd qdiv'
'draw  title \sep.6-18 qdiv'
'draw map'
'draw xlab \time'
'draw ylab \presure'
'print'
'disable print'
'printim e:\201409\1\qdiv1.png white'
'c'
错误贴图，求指导。

zhouling

楼主还是把错误贴出来比较好。

咿呀呀

错误如下：

panghuaji

你是画一点的剖面图吗？
可能是时间设置问题，'set t 992 1044'是一个时间段，而程序中没有循环语句，建议改为'set t 992'先试试能否出结果。
也可能是数据问题，你设置是画一个点的剖面，根据您的经纬度，如果你的数据格点大于1*1，则你画的那个点可能就没数据，建议先画一条线的剖面看看。

咿呀呀

panghuaji 发表于 2014-10-30 11:08
你是画一点的剖面图吗？
可能是时间设置问题，'set t 992 1044'是一个时间段，而程序中没有循环语句，建议 ...

我试了，还是不出图。我用同样的代码，垂直速度图都出来了。我现在的代码是想作图，横坐标是日期，纵坐标是高度场。

独行侠

水汽通量剖面图你为什么会用到求散度的函数呢，之所以不出图很有可能是温度出了问题，nc文件你的air不需要减273.16

白令海

我也遇到了同样的问题，如果是在平面图上画某一层的水汽通量散度就不会出问题。但是，画剖面图就会问题了。都过了快四年了，不知道楼主的问题解决了没有。

269772539

我考虑是不是系统计算智能在某一层上计算，即为二维数据的计算？