【求助】整层水汽通量程序总是无法运行？

瑾珂

单层的水汽通量是可以画出来的，但运行整层的水汽通量gs总是出错?各位老师能不能帮忙看看是哪里出问题了？
'reinit'
'open J:\0813-032\DJBY\20y_winter\2017_1\fnl_201701.ctl'


i=1
*while(i<=21)
'set t 'i''
'set gxout shaded'
'set mpdset cnworld'
'set map 1 1 9'
'set csmooth on'
'set grads off'
'set grid off'


'set lev 1000 1'
'define tp=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define uw=UGRDprs'
'define vw=VGRDprs'

*------------计算比湿  单位：kg/kg--------------*
'define es=(6.11*exp((17.2693882*(tp-273.16))/(tp-35.86)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define qq=rh*qs/100'
*-------计算水汽通量(单层)  单位：(g*s)/kg=g/(s*hpa*cm)------*
'define qu=uw*qq*1000/9.8'
'define qv=vw*qq*1000/9.8'

'define pres=PRESsfc/100'
'define uqx=vint(pres,qu,300)'
'define vqx=vint(pres,qv,300)'

*------------计算水汽通量散度------*

'define aqdiv=hdivg(uqx,vqx)*100000'

'set grads off'
'set grid off'
'set parea 1 9.5 0.5 8'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set xlint 5'
'set ylint 2'
'set gxout shaded'
'set cint 2'
'd mag(uqx,vqx)'
'run J:\0813-032\2022222snow\gs\output.gs'
'set gxout barb'
'set ccolor 2'
'd uqx;vqx'
'set gxout contour'
'set ccolor 9'
'set cthick 6'
'd aqdiv'
*'cbarm'
*'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbarm 1 1 1'
'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbar_interp1'


'q dim'
it=sublin(result,5)
it=subwrd(it,6)
'draw title 'it'-jifen-waterflux'
'printim  J:\0813-032\DJBY\ca_saf\201701\waterflu\'it'-jifen-waflux.wmf white'
*c
*i=i+1
*endwhile
'disable gxprint'
;

苏拉苏拉

错误提示 'aqu' 这个变量有问题。但是我扫了两眼你的程序，里面并没有aqu这个变量啊

瑾珂

苏拉苏拉 发表于 2022-11-28 15:28
错误提示 'aqu' 这个变量有问题。但是我扫了两眼你的程序，里面并没有aqu这个变量啊

贴错图了，后面我又改了一下，能运行了，就是这个量级有点大

river

瑾珂 发表于 2022-11-28 15:42
贴错图了，后面我又改了一下，能运行了，就是这个量级有点大

量级你不是计算的时候扩大了10的5次方了吗