GrADS填色出现白边问题

Eegle

在画水汽通量散度的时候 出现了白边的问题 在论坛上搜索了一下方法 按照上面做的 第一张图没有白天 后面的几张又出现白边的问题了。是不是循环的时候出错还是？请高手指点 不胜感激。
'reinit'
'enable print E:\zhuang\picture\qdiv.gmf'
'sdfopen E:\zhuang\uwnd.2014.nc'
'sdfopen E:\zhuang\vwnd.2014.nc'
'sdfopen E:\zhuang\rhum.2014.nc'
'sdfopen E:\zhuang\air.2014.nc'
'run E:\zhuang\define_colors.gs'

'set grads off'
'set grid off'
'set poli on'
'set mpdraw on'
'set mproj scaled'
'set mpdset cnworld'
'set xlab on'
'set ylab on'
'set lat 0 60'
'set lon 80 160'
'set xlopts 1 4 0.16'
'set ylopts 1 4 0.16'
'draw map'

h=654
while(h<672)

'set mpdset cnworld'
'set grads off'
'set t 'h''
'set lev 850'
'set gxout shaded'


'define p=lev'
'define e=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhum.3/100'
'define qdiv=hdivg(q*uwnd.1/9.8,q*vwnd.2/9.8)'
'set lat 15 35'
'set lon 105 135'
'set clevs  0 20 30 40 50 60 70 80'
'set ccols  100 81 82 83 84 85 86 87 88 89'

'd qdiv*1e10'
'cbarn'

'set gxout contour'
'set lat 15 35'
'set lon 105 135'
'set cthick 6'
'set ccolor 1'
'set clopts 1 -1 0.16'
'set cmin 0'
'set grads off'
'd qdiv*1e10'

'set gxout barb'
'set digsize 0.06'
'set ccolor 2'
'set cthick 4'  
'set grads off'
'd uwnd.1;vwnd.2'

'q time'
tt=subwrd(result,3)
hour=substr(tt,1,2)
day=substr(tt,4,2)
mon='6'
year='2014'

BT=hour+8      

if(BT>24)
day=day+1
BT=BT-24
endif

'draw title qdiv 850hpa 'year'-'mon'-'day' 'BT':00'
   
'print'
'c'
h=h+1
endwhile
'disable print'
;

river

把你最上面设置的那个经纬度范围放到循环里去

汾西门子

好难哦，就要假期了，想放松的~~

Eegle

汾西门子 发表于 2014-9-28 21:44
好难哦，就要假期了，想放松的~~

哈哈哈哈 放假前搞定 才能过个轻松的假期 不过对于预报员来说 是没有假期的 555555555555

Eegle

river 发表于 2014-9-28 21:05
把你最上面设置的那个经纬度范围放到循环里去

好棒 谢谢 问题解决了  谢谢 真心感谢 每次都这么热心地帮我。

river

Eegle 发表于 2014-9-29 08:47
好棒 谢谢 问题解决了  谢谢 真心感谢 每次都这么热心地帮我。

每次吗？我自己倒是没有太注意