grads求助！request beyond file limits

寄居蟹

'reinit'
'open G:\con.ctl'
'open G:\ens.ctl'
'open G:\dp.ctl'
'enable print G:\dp_uv.gmf'

'set lev 850'
'define du=u.2-u'
'define dv=v.2-v'

'set dfile 3'
'set lev 1000'
'define dp=dpave.3'

'set mpdset cnworld'
'set map 15 1 6'
'set grads off'
'set grid off'
'set lat -20 60'
'set lon 30 160'

'set gxout contour'
'd dp'
'set gxout vector'
'd du;dv'
'print'
'disable print'
;
三个ctl文件的格点，时间，和层次都是一样的。但运行出现以下问题.请问问题出在哪？谢谢！

river

set dfile 3之后你设置的是set lev 1000，但是du和dv是850hPa上的，所以就超出范围了

guobch

路过   

寄居蟹

river 发表于 2013-9-6 20:05

                               
登录/注册后可看大图

set dfile 3之后你设置的是set lev 1000，但是du和dv是850hPa上的，所以就超出范围了

那请问怎样把1000pha的dp变量和850hpa的风场叠加呢？

lemon_tree

寄居蟹 发表于 2013-9-9 09:29

                               
登录/注册后可看大图

那请问怎样把1000pha的dp变量和850hpa的风场叠加呢？

建议每画一个变量之前都设置一下层数试试

寄居蟹

lemon_tree 发表于 2013-9-9 16:21

                               
登录/注册后可看大图

建议每画一个变量之前都设置一下层数试试

可是我后面加了，也还是相同的错误

叶の幽灵

用括号 d dp(z=1)
        'd du(z=3);dv(z=3)'
试试

寄居蟹

叶の幽灵 发表于 2013-9-10 10:23

                               
登录/注册后可看大图

用括号 d dp(z=1)
        'd du(z=3);dv(z=3)'
试试

也还是不行啊，出来还是同样的错误

叶の幽灵

寄居蟹 发表于 2013-9-10 11:09

                               
登录/注册后可看大图

也还是不行啊，出来还是同样的错误

这个真不知道了。。。好好对照一下ctl的数据范围吧

寄居蟹

叶の幽灵 发表于 2013-9-10 14:15

                               
登录/注册后可看大图

这个真不知道了。。。好好对照一下ctl的数据范围吧

分开单独每个ctl都是可以画图的，应该ctl范围没错吧？