各位帮忙看看GS哪里出了问题，小弟实在是找不到

sound

'reinit'
********************************************* data ***
'open e:\ncep\0131\fnl_20150108_12_00.grib1.ctl'
'set lon 101 101'
'set lat 10 30'
********** lev ************
'set lev 1000 200'
'set grid off'
'set timelab off'
'set vpage 0 11 0 8'
'set parea 1 10.5 0 8'
'set grads off'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'

** water vapour **
'define es=6.112*exp(17.2693822*(tmpprs-273.16)/(tmpprs-29.65))'
'define qs=0.622*es/(lev-1.378*es)'
'define q=rhprs*qs/100*1000'  

***  water vapor flux divergence ***
'define qu=ugrdprs*q'
'define qv=vgrdprs*q'
'define qd=hdivg(qu,qv)*1e5/9.8'
*** draw contour  ***
'set gxout contour'
'set clopts 1 4 0.2'
'set csmooth on'
'd qd'
'printim E:\ncep\0131\tp\div-shuiqi08-20.png white'




以上为编写的GS，以及运行错误情况。谢谢各位。

sound

希望给位能给个答案，非常感谢。

深度狙杀

hdivg求水平散度。你设置了'set lon 101 101'，水平方向仅有一维  怎么差分

狂奔的和尚

赞一个！！！！！！！！！！！

sound

sound 发表于 2015-3-10 11:09
希望给位能给个答案，非常感谢。

已解决。找到问题所在了，非常感谢。