关于计算整层水汽通量散度的问题

周一的07

介绍以下我计算整层夏季水汽通量散度的步骤：一、计算出各月的整层水汽输送通量，后计算散度，gs程序如下。（时常为1979-2018年共计480个月）
'reinit'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\U_1979_2018.nc'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\V_1979_2018.nc'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\shum_1979_2018.nc'
'sdfopen E:\Master_thesis\GrADS\thesis\nc_data\sp_1979_2018.nc'
'set gxout fwrite'
'set fwrite E:\Master_thesis\GrADS\thesis\gs\hdivg\nj_hdivg.dat'
n=1
while(n<=480)
'set z 1'
'set t 'n
'set x 1 241'
'set y 1 161'
*====计算水汽通量散度=====================
'define quall=vint(sp.4(z=1)/100,q.3*u.1,100)'
'define qvall=vint(sp.4(z=1)/100,q.3*v.2,100)'
*'d quall'
*'d qvall'
'd hdivg(quall,qvall)'
n=n+1
endwhile
'disable fwrite'
return
；

二、得到各月整层的水汽通量散度后，求出夏季（6、7、8月）的水汽输送通量散度并绘制成图。
我先给出不出图的gs程序：
'reinit'
'open E:\Master_thesis\GrADS\thesis\gs\hdivg\nj_hdivg.ctl'
'set gxout fwrite'
'set fwrite E:\Master_thesis\GrADS\thesis\gs\hdivg\hdivg_summer.dat'
'set z 1'
'set t 1 480'
'set x 1 241'
'set y 1 161'
*======================
'define qujun=ave(hdivg,t=6,t=480,12)'
'define qujul=ave(hdivg,t=7,t=480,12)'
'define quaug=ave(hdivg,t=8,t=480,12)'
'define divsummer=(qujun+qujul+quaug)/3.0'
'set t 1'
'd divsummer'
'disable fwrite'
return

；
以上程序时可以完成计算的，grads没有报错，能成功获得夏季的水汽输送通量散度并存入hdivg_summer.dat；当我尝试直接出图时的gs程序：
'reinit'
'open E:\Master_thesis\GrADS\thesis\gs\hdivg\nj_hdivg.ctl'
'set z 1'
'set t 1 480'
'set x 1 241'
'set y 1 161'
*=========夏季水汽通量散度计算=============
'define qujun=ave(hdivg,t=6,t=480,12)'
'define qujul=ave(hdivg,t=7,t=480,12)'
'define quaug=ave(hdivg,t=8,t=480,12)'
'define divsummer=(qujun+qujul+quaug)/3.0'

*=====出图======
'set t 1'
'set parea 1.5 10 1.5 8.0'
'set mpdset cnworld'
'set map 1 1 3'
'set xlopts 1 3 0.2'
'set ylopts 1 3 0.2'
'set xlint 1'
'set ylint 1'
'set lon 76 103'
'set lat 24 35'
'set grads off'
'set grid off'
'set xlpos -10'
'set ylpos -10'
'set gxout shaded'
'd divsummer'
'cbarn 0.8 0'
'run axis.gs -type b  -position m -label on -start 76 -end 103 -suffix `3. -interval 3  -size 0.2 -lsize 0.15 -lthick 0.3 '
'run axis.gs -type L  -position o -start 24 -end 35 -suffix `3. -interval 2 -lsize 0.15 -lthick 0.3'
'set line 2'
'draw shp E:\Master_thesis\GrADS\thesis\arcgis\yarlung'
'printim E:\Master_thesis\GrADS\thesis\gs\hdivg\hdivg_sum.png white'
；

grads报错情况如下：
Cannot contour grid - all undefined values
Cannot plot color bar: No shading information
Nothing displayed
Nothing displayed


并且我曾成功使用上述出图语句完成出图，但是放在整层的水汽通量散度这里就无法出图；
以上；
我在论坛也借鉴了很多其他的出图程序，但是报错情况都是如此，不知道什么原因，还请老师们指教。

river

你把set t 1 放到 'define divsummer=(qujun+qujul+quaug)/3.0' 前面试试

周一的07

river 发表于 2022-2-14 18:18
你把set t 1 放到 'define divsummer=(qujun+qujul+quaug)/3.0' 前面试试

谢谢您的指导，刚刚试了一下，还是同样的反馈信息，无法出图，不知道是什么原因。

river

我的建议是一步一步来，你出来的资料hdivg_summer.dat有没有配个ctl，看看出的图对不对

周一的07

river 发表于 2022-2-15 21:12
我的建议是一步一步来，你出来的资料hdivg_summer.dat有没有配个ctl，看看出的图对不对

换了个方法弄了，图是出来了，感觉应该是正确的，谢谢river给的思路。

周一的07

换了个方法弄了，图是出来了，感觉应该是正确的，谢谢river给的思路。