grads绘图不能铺满整个范围

spikeyding

我昨天在坛子里找到了用ncep的资料画水汽通量和水汽通量散度的gs，在这儿先感谢那些乐于分享的大神！
我脚本如下
'reinit'
'open h:/fnl_ding/2013/fnl_20130522_18_00.ctl'
'set mpdset cnworld'
'set lev 1000'
'set t 1'
'set lon 105 115'
'set lat 28 37'
'set gxout shaded'
'set csmooth on'
'define p=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*RHprs/100'
'define qdiv=hdivg(q*UGRDprs/9.8,q*VGRDprs/9.8)'
'd qdiv'
'cbarn'
'printim h:/fnl_ding/picture/div_vf052218_1000.png white x1000 y800'
;

然后画出来的图边上就没有值了

其它变量没问题，但是水汽通量散度画出来就是这样，在此求教各位老师和前辈

spikeyding

对了，忘了说，即便是改了经纬度范围，还是不行

spikeyding

经过师兄的指点，已经知道哪边有问题了，散度的计算用到相邻的格点，所以最边上的格点无法被计算出来。先把区域放大，然后在画图之前缩小就行了。脚本如下
'reinit'
'open h:/fnl_ding/2013/fnl_20130522_18_00.ctl'
'set mpdset cnworld'
'set lev 1000'
'set t 1'
'set lon 100 120'
'set lat 20 45'
'set gxout shaded'
'set csmooth on'
'define p=lev'
'define e=(6.112*exp((17.67*(TMPprs-273.16))/(TMPprs-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*RHprs/100'
'define qdiv=hdivg(q*UGRDprs/9.8,q*VGRDprs/9.8)'
'set lon 105 115'
'set lat 28 37'
'd qdiv'
'cbarn'
'printim h:/fnl_ding/picture/div_vf052218_1000.png white x1000 y800'
;

river

其他变量没问题，那是因为你没画涡度····

mofangbao

其实好多人提到过这个问题了，并且给了方法，可能是你没有注意到，看帖子的同时建议看看看回复，有的回复里面的内容也很精彩！

qq469015280

你把计算的范围扩大，画图的范围缩小，就可以了。

你现在的图没画出来的是边界问题

spikeyding

qq469015280 发表于 2015-4-9 09:59
你把计算的范围扩大，画图的范围缩小，就可以了。

你现在的图没画出来的是边界问题

是的，就是这样

spikeyding

mofangbao 发表于 2015-4-3 16:15
其实好多人提到过这个问题了，并且给了方法，可能是你没有注意到，看帖子的同时建议看看看回复，有的回复里 ...

qlguai

spikeyding 发表于 2015-4-3 13:51
经过师兄的指点，已经知道哪边有问题了，散度的计算用到相邻的格点，所以最边上的格点无法被计算出来。先把 ...

不错，感谢你啦！！！！！！！！！

hugeal10000

插值以后边上有空白！