求教：grads处理fnl数据，图片四周出现空白，求解

繁星一点

   下载的0.5*0.5的fnl 再分析数据 ，利用grads绘图，但是在图片四周出现空白区域，查了好久没查出原因，求大神指教  

gs文件如下：

'reinit'
'open H:\fnl\20140702\201407.ctl'
tt=1
while(tt<2)
'set lat 20 35'
'set lon 105 125'
'set lev 850'
'set t 'tt' '
*****************************************
'set vpage 0.5 11 0.5 8.5 '
'set mpdset cnhimap hunan_q '
'set map 15 1 1'
'draw map'
********************************************
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set xlint 5'
'set ylint 5'
'set gxout shaded '
'set grads off'
'set grid off'
'define_colors'
'set clevs -100 -80 -60 -40 -20 0 20 40 60 80'
'set ccols 26 25 24 23 22 21 41 42 43 44 45'
'define div=hdivg(UGRDprs,VGRDprs)*1e6'
'd div'
'cbarn 0.8 0 5.5 8.2 '
'set gxout contour'
'set clevs -100 -80 -60 -40 -20 0 20 40 60 80'
'define curl=hcurl(UGRDprs,VGRDprs)*1e6'
'd curl'
************************************************
'draw string 0.6 0.6 time='tt''
if(rc='0') break;endif;
'printim H:\20140702\fnl\div\div_850_'tt'_2.gif white'
'print'
'c'
**************************************************
tt=tt+1
endwhile
'disable print'
return


图片如下：可以确定不是我的范围画的太大了，因为我缩小的经纬度范围以后还是出现相同的情况。


请教各位，是不是我的gs文件里面哪一句出了问题？

1092995269

设置两次经纬度，计算的范围大于显示的范围就可以了

繁星一点

1092995269 发表于 2014-11-16 08:03
设置两次经纬度，计算的范围大于显示的范围就可以了

设置了  还是一样出现空白

传说中的谁

你得明白散度的计算方法，某个格点上的散度是需要通过周边几个格点的风向风速来计算得到的。楼上上的已经说的很清楚了，是楼主没有理解，就是计算散度的时候范围略大一些（至少比你想要画图的范围各个方向都增加一个经/纬度），画图的时候再设置成你想要 的范围。

繁星一点

传说中的谁 发表于 2014-11-16 09:06
你得明白散度的计算方法，某个格点上的散度是需要通过周边几个格点的风向风速来计算得到的。楼上上的已经说 ...

我重新画了，计算散度经度范围是100-140E，纬度是10-40N，然后再设置了上述画图范围，还是一样的

传说中的谁

繁星一点 发表于 2014-11-16 22:47
我重新画了，计算散度经度范围是100-140E，纬度是10-40N，然后再设置了上述画图范围，还是一样的

gs呢？

繁星一点

传说中的谁 发表于 2014-11-17 01:18
gs呢？

'reinit'
'open H:\fnl\20140702\201407.ctl'
tt=1
while(tt<2)
'set lat 10 40'
'set lon 100 140'
'set lev 850'
'set t 'tt' '
'define div=hdivg(UGRDprs,VGRDprs)*1e6'
'define curl=hcurl(UGRDprs,VGRDprs)*1e6'
*****************************************
'set vpage 0.5 11 0.5 8.5 '
'set mpdset cnhimap hunan_q '
'set map 15 1 1'
'draw map'
********************************************
'set lat 20 35'
'set lon 105 125'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
'set xlint 5'
'set ylint 5'
'set gxout shaded '
'set grads off'
'set grid off'
'define_colors'
'set clevs -100 -80 -60 -40 -20 0 20 40 60 80'
'set ccols 26 25 24 23 22 21 41 42 43 44 45'
'define div=hdivg(UGRDprs,VGRDprs)*1e6'
'd div'
'cbarn 0.8 0 5.5 8.2 '
'set gxout contour'
'set clevs -100 -80 -60 -40 -20 0 20 40 60 80'
'define curl=hcurl(UGRDprs,VGRDprs)*1e6'
'd curl'
************************************************
'draw string 0.6 0.6 time='tt''
if(rc='0') break;endif;
'printim H:\20140702\fnl\div\div_850_'tt'_2.gif white'
'print'
'c'
**************************************************
tt=tt+1
endwhile
'disable print'
return

======================================================
已经解决了   发现gs文件中脚本编写的问题  标红的那两句重复计算了一次涡度和散度
注释掉以后就好了
非常谢谢！

沸水咖啡

学习一下。。。。。。。。。。。。。。

自由行

受益匪浅，好东东啊

xuehongyang

好东西{:eb502:}{:eb502:}