为什么我画的地形图是这个样子？假相当位温却没有了。

天边的彩虹

我用的地形资料是top.ctl。利用一个楼主提供的gs,绘制的图形是这样子。求高手解决。
'reinit'
'open g:/gradsdata/20140525/fnl0525.ctl'
'open g:/gradsdata/dixing/top.ctl'
'enable print g:/gradsdata/20140525/thetase.gmf'
'set lat 10 50'
'set lon 90 140'
'set lev 1000 100'
'set xlopts 1 4 0.2'
'set ylopts 1 4 0.2'
'set grid on'
'set grads off'
'set z 1 19'
'set ylevs 1000  925  850   700 600 500 400  300 200'
'set zlog on '

'set dfile 2'
'set t 1'
'set lat 20 40'
'set lon 111.5'
'set lev 1000'
'set xlint 5'
'set ylint 100'
'define l1=smth9(ave(top,lon=111.5,lon=111.5))'
'define l2=top*0.0'
'set vrange 0 8000'
'set ylab off'
'set gxout linefill'
'set lfcols 1 2'
'd l1;l2'

t=1
*时间
while(t<=13)
'set t 't
aaa=t
'set gxout contour'
'set dfile 1'
'set lat 20 40'
'set lon 111.5'
'set z 1 19'
'set ylevs 1000  925  850  700 600 500 400  300 200'
'set zlog on '
'set clab on'
'set clopts -1 -1 0.1'
'set cint 2'
'set clab on'
'define t=tmpprs'
'define rh=rhprs'
'define prs=lev'
'define es=(6.112*exp(17.67*(t-273.15)/(t-29.65)))'
'define q=rh*(0.62197*es/(prs-es))/100.'
'define e=prs*q/(0.62197+q)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(t)-log(e)-4.805)'
'define theta=t*pow((1000/prs),(0.2854*(1.0-0.28*q)))'
'define eqt9=theta*exp(((3376./tlcl)-2.54)*q*(1.0+0.81*q))'
'set xlpos -20'                                                                           
'd eqt9'

river

把这句'set vrange 0 8000'去掉看看

天边的彩虹

我使用了'set vpage 'set parea ，解决了这个不一致的问题。

luojianyu

请问是用1°×1°grib数据画的吗？可以分享top.ctl'吗？作剖面图叠加地形是我一直没解决的问题，求教了！