为何月度水汽通量散度皆为正值？

大老鹰

数据为NCEP 平均月度数值。如果范围扩大，会出现大片区域无法上色现象，提醒“logic error 3或1 in gxshad2 s2follow”。请各位指教
代码如下：

'reinit'
'sdfopen D:\NOAA\ESRL.2.5\air.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\rhum.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\uwnd.mon.ltm.nc'
'sdfopen D:\NOAA\ESRL.2.5\vwnd.mon.ltm.nc'
*输入文件
'set mpdraw off'
*关闭原有底图
'set lon 30 90'
'set lat 0 30'
*设置经纬度
'set lev 850'
t=1
while(t<=12)
*设置层级，时间起始
'set t 't''
'set grads off'
'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.3/9.8,q*vwnd.4/9.8)'
'd qdiv'
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
*水汽通量散度计算
'set gxout shaded'
'run D:\NOAA\200colors.gs'
'set clevs 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1'
*阴影图，设置色带
'set gxout vector'
'd mag(q1,q2)'

'set ccolor 0'
'd uwnd.3*2.5;vwnd.4*2.5'
*矢量图，设置U\V分量，蓝色，U\V分量频率
'set gxout shp'
'set line 1'
'draw shp D:\NOAA\country.shp'
*shp图，白色线形，画shp图 文件路径
'q dim'
a=sublin(result,5)
a=subwrd(a,9)
'draw title 'a
'set rgb 50 255 255 255'
'run cbarn.gs'
*查询维度，导入月份时间作为标题。
'printim D:\NOAA\'%t%'.png x1300 y650 white'
'C'
t=t+1
*输出图片，清屏，循环。
endwhile

夏天夏夏夏

我也出现了这种情况，请问你解决了吗