散涡比图形优化

追风的阳光

我想作出散涡比的时间序列图，这是其中一段程序，涡度散度同时大于0的情况，但是图形很不美观
do k=1,nz
do s=1,nt
        do j=1,ny
                do i=1,nx
                if (vor2(i,j,s,k)>=0 .and.div(i,j,s,k)>=0)then
                svbit(i,j,s,k)=div(i,j,s,k)/vor2(i,j,s,k)
                else
                svbit(i,j,s,k)=9.999E+20
                endif
  enddo
  enddo
  enddo
  enddo
这是ctl：
dset d:\zssdata\svbit.grd
undef 9.999E+20
title uwnd
xdef 21 linear 110 1
ydef 11 linear 25 1
zdef 26 levels  
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
tdef 29 linear 00:00Z08jul2011 6hr
vars 1
svbit 26 99 vorticity using hgt
endvars
这是gs：
'reinit'
'clear'
'open D:\zssdata\svbit.ctl'
'enable print D:\zssdata\svbitt2.gmf'
'set grads off'
'set grid on'
'set lev 1000 100'
'set t 9 29'
'set xlevs 00 06 12 18 00 06 12 18 00 06 12 18 00 06 12 18 00 06 12 18 00 '
'define avesvbit=tloop(aave(svbit,lon=117,lon=120,lat=30,lat=32))'
'set xyrev off'
'set gxout contour'
'set ccolor 1'
'set cthick 4'
'set clopts -1 -1 0.18'
'set xlopts 1 5 0.2'
'set cint 10'
'set ylopts 1 5 0.2'
'set parea 1 10 1 7.5'
'set lat 29'
'set lon 117'
'set t 9 29'
'set lev 900 100'
'set zlog on'
'd avesvbit'



最后出图是这样的：

请问这种图形是怎么一回事呢，有什么办法优化它呢？

又是那隻貓

图很不好看0.0

追风的阳光

又是那隻貓 发表于 2013-8-20 05:52

                               
登录/注册后可看大图

图很不好看0.0

恩，怎么调整呢？