新手求助！在画水汽通量散度场时等值线只显示一部分，这是为何？该如何处理？

周一的07

一、给出现有成果图。（450hpa水汽通量散度场）
二、给出我的gs程序。
'reinit'
'open E:\Master_thesis\GrADS\YFY_practice\llj\ctl\sandu\450hpa.ctl'
'set z 1'
'set t 1 12'
'set x 1 241'
'set y 1 161'
'define hdivg450=hdivg(q450*u450/9.8,q450*v450/9.8)*1e5'
'set t 1'
'set parea 1.5 10.0 1.5 8.0'
'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 81 98'
'set lat 27 32'
'set grads off'
'set grid off'
'set xlpos -10'
'set ylpos -10'
'set font 5'
'set gxout contour'
'set ccolor 1'
'set clopts 1 5 0.1'
'd hdivg450'
'run axis.gs -type b  -position o -label on -start 81 -end 98 -suffix `3. -interval 2  -size 0.2 -lsize 0.15 -lthick 0.3 '
'run axis.gs -type L  -position o -start 27 -end 32 -suffix `3. -interval 1 -lsize 0.15 -lthick 0.3'
'set line 2'
'draw shp E:\Master_thesis\GrADS\thesis\arcgis\yarlungfenqv'
'printim E:\Master_thesis\GrADS\thesis\png\sandu\450_1mo.png white'
；

三、设置的范围都对，水汽通量图也能正常出图，排查了几遍，但是还是不知道这个是什么问题引起的，请各位大神指点！感恩！

river

是不是原始资料的范围没这么大啊

灵心流萤

计算的范围一定要比显示的范围大才行，你看下你的数据经纬度范围是多大

周一的07

river 发表于 2021-9-27 22:23
是不是原始资料的范围没这么大啊

您好，我把我的过程和您说一下：
一、我先利用U、V和shum三个数据生成可识别的dat数据，gs程序如下：
'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'
'set gxout fwrite'
'set fwrite E:\Master_thesis\GrADS\thesis\dat\sandu\450hpa.dat'
'set lon 60 120'
'set lat 5 45'
'set z 1 27'
'set t 1 480'
'set lev 450'
'define u450=ave(u.1,t+0,t=480,12)'
'define v450=ave(v.2,t+0,t=480,12)'
'define q450=ave(q.3*1000,t+0,t=480,12)'
'set t 1 12'
'd u450'
'd v450'
'd q450'
'disable fwrite'
return
；
（注：我的原始资料范围确认为经度60~120，纬度5~45，如图所示）
二、为dat数据配备ctl
DSET  E:\Master_thesis\GrADS\YFY_practice\llj\dat\sandu\450hpa.dat
UNDEF  -9999.0
XDEF   241  LINEAR  60   0.25
YDEF   161  LINEAR  5   0.25
ZDEF   1    LINEAR  1    450
TDEF   12  linear 00z0jan1979 1mo
VARS  3
u450 0  99 u450 water all
v450 0  99 v450 water all
q450 0  99 q450 water all
ENDVARS
三、运行我提问时附上的gs程序
麻烦老师再帮我看看，是不是我哪里出错了？感谢！

周一的07

灵心流萤 发表于 2021-9-28 08:44
计算的范围一定要比显示的范围大才行，你看下你的数据经纬度范围是多大

我之前考虑过会不会是这个问题，但是我查了一下，我的原始资料范围确认为经度60~120，纬度5~45，如图所示，所以不知道是哪里出错了。

周一的07

我又重新梳理了一遍，不知道为何，图又能都显示出来了，但还是谢谢各位前辈。

周一的07

周一的07 发表于 2021-9-28 15:03
我又重新梳理了一遍，不知道为何，图又能都显示出来了，但还是谢谢各位前辈。

知道问题了，在写ctl的时候，我把数据的路径输错了。

myj0219

感谢分享！！

但愿

{:eb502:}{:eb502:}