求助，2ctl左右坐标出图出现错误

追风的阳光

这是gs文件
'reinit'
'clear'
'open D:\zssdata\fnl.ctl'
'open D:\zssdata\pre.ctl'
*'enable print D:\zssdata\svbitt2.gmf'
'set dfile 1'
'set grads off'
'set grid on'
'set lev 1000 100'
'set t 1 13'
'set xlevs 00 06 12 18 00 06 12 18 00 06 12 18 00'
'set ylevs 10 20 30 40 50 60 70 80 90 100'
'define h=100000*hdivg(ugrdprs,vgrdprs)'
'define w=100000*hcurl(ugrdprs,vgrdprs)'
'define svbi=h/w'
'define aa=tloop(aave(svbi,lon=117,lon=120,lat=30,lat=32))'
'set lat 29'
'set lon 117'
'set t 1 13'
'set lev 900 100'
'set zlog on'
'set xyrev off'
'set gxout contour'
'set ccolor 1'
'set cthick 4'
'set xlevs 00 06 12 18 00 06 12 18 00 06 12 18 00'
'set ylevs 975 925 850 700 600 500 400 300 200 100'
'set clopts -1 -1 0.18'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set cint 5'
'd aa'
'set dfile 2'
'set grads off'
'set grid on'
'set ylpos 8.5 l '
'set t 1 4'
'set xlevs 00 06 12 18 00 06 12 18 00 06 12 18 00'
'set ylevs 0 50 100 150 200 250 300 '
'set lon 0'
'set lat 0'
'set t 1 4'
'd a'
*'print'
*'disable print'
'printim d:\zssdata\sbvitaa.gmf'
;
这是pre.ctl
dset D:\zssdata\4zhan.grd
undef -9.99E+08
title a
xdef 1 linear 0.0 1.0
ydef 1 linear 0.0 1.0
tdef 4 linear 00Z11jul2011  1dy
zdef 1 linear 10  1
vars 1
a  0 99  
endvars

2个ctl可分别作图，为何用上述gs文件会出现以下错误：

G-mimzy

在 'd a' 之前加一句 'set lev 10' 试试；因为变量 a 垂直方向只有一层，却按照之前设置的 'set lev 1000 100' 来了

river

a是第二个文件里的，不需要区分吗？应该是d a.2吧，而且display之后也没有添加c命令，如果出图的话不都叠在一起了么

追风的阳光

谢谢图已出来了，但坐标轴没有对齐，我再看看

追风的阳光

river 发表于 2012-12-4 09:33

                               
登录/注册后可看大图

a是第二个文件里的，不需要区分吗？应该是d a.2吧，而且display之后也没有添加c命令，如果出图的话不都叠在 ...

我就是想让俩个图横坐标相同，纵坐标不同，叠加在一起的，文件已用set dfile区分了

river

追风的阳光 发表于 2012-12-4 09:48

                               
登录/注册后可看大图

我就是想让俩个图横坐标相同，纵坐标不同，叠加在一起的，文件已用set dfile区分了

就是说起始时间是一样的吧，那时间长度也应该相等的吧，一个13一个4不一样啊。纵轴的话可以利用set ylpos命令来显示两个纵坐标分别对应你的两套资料

追风的阳光

river 发表于 2012-12-4 09:53

                               
登录/注册后可看大图

就是说起始时间是一样的吧，那时间长度也应该相等的吧，一个13一个4不一样啊。纵轴的话可以利用set ylpos ...

时间间隔一个是6hr一个是1dy，时间轴是对应的

追风的阳光

G-mimzy 发表于 2012-12-4 09:28

                               
登录/注册后可看大图

在 'd a' 之前加一句 'set lev 10' 试试；因为变量 a 垂直方向只有一层，却按照之前设置的 'set lev 1000  ...

这是我添加set lev 10后的出图，但为什么时间轴无法重合呢？

小草

用set parea 控制一下出图的区域，看看能不能重合？

追风的阳光

小草 发表于 2012-12-4 14:02

                               
登录/注册后可看大图

用set parea 控制一下出图的区域，看看能不能重合？

谢谢你，可以了