水汽通量散度高度-纬度剖面出错

某某晨

利用hdivg函数计算了水汽通量散度，想要画出高度-纬度剖面图，结果只能出来数字序列，求各位大神指点，以下是gs文件:
'sdfopen   E:\20220817\waper_flux\20220817.nc'
'define qdiv=cdiff(q*u/9.8,q*v/9.8)'*计算水汽通量散度
'set t 1'
'set parea 2 8.5 2 6.5'
'set grads off'
'set grid off'
'set lon 101.5 '
'set lat 35 40'
'set lev 1000 300'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300'
'set gxout contour'
'd qdiv'
'printim E:\20220817\waper_flux\pm1708.png white'
'c'
;

某某晨

用上面的画图语句画比湿可以出来剖面图，是不是hdivg函数的原因呢？如果不用hdivg函数要用什么计算水汽通量散度呢？