GrADS显示一页多图时候遇到的奇怪问题

Sky-walker

大家好，我在论坛里学习了站点插值格点显示以及一页多图的方法以后，在自己尝试的时候遇到了一些问题，实在费解，希望大神能帮我解惑，谢谢
遇到的问题是在一页多图的时候，GrADS报错为




我的gs文件是这样的

现在我的显示效果如下


q1区域是显示站点插值图，这部分是可以正常显示的，但是q2,q3,q4区域的图是没法显示


但是q2,q3,q4是可以单独显示成功的，就是和q1一起显示的时候会报这个错误，是不是数据不太匹配？？


Ps:现在我知道问题就是出在插值函数这里，原因可能是因为grid.ctl是我自己制作的格点场，而前面几个test19.ctl应该用的是全球格点场，现在可能问题就是出在这里了，但是这里怎么解决呢，如何让前几个文件用自己的全球格点场，而插值文件用自己制作的格点场呢？？


各位有什么好的建议码，谢谢不吝赐教了！！



那我应该如何设置才能运行成功呢，请大家帮帮我，谢谢

千里草

好像set vpage时需要设置图形的位置，假如你是四幅图，可以这样
'set vpage 1 5.5 0.5 3.5'
意思就是这幅图位置在整个屏幕的位置是：边长从1到5.5，高（宽）从0.5到3.5，以此类推，可以试试看

'set vpage 6 10.5 4 7.5'
'set vpage 3.5 8.5 0.5 3.5'

枫丹白露

从你的描述来看，你应该是用了两个不同维数环境的数据及其描述（.ctl）文件。解决的办法是用'set dfile ...'改变当前维数环境对应的文件，而后从新设置维数环境。比如：
'open data1.ctl'
'open data2.ctl'

'set vpage 'p1
'set lon 'lon1' 'lon2
'set lat 'lat1' lat2
'set t 'time
'set z 'lev
'd date1'

'set vpage 'p2
'set dfile 2'
'set lon '...
'set lat '...
'set z '...
'set t '...
'd data2.2'
如果lon，lat，lev（不是z），time（不是t）的设置是一样的，那么在'set dfile 2'之后不需要重新设置，它们会自动变成第二个数据文件的维数环境。但如果不一样，则必须重新设置。尤其注意lev和time在很多情况下对应的维数环境并不相同。假如第一个ctl文件是从2002年开始逐月，而第二个是从2000年开始逐月，对于'set t 3'，对应第一个文件是2002年3月，通过'set dfile 2'传递到第二个文件，如果只有2000和2001年的数据，那么所有对应2002年3月的数据都是undefied，这里传递的是“2002年3月”这个时间设置，而不是t的值3；lev同理，传递的是某一等高面比如1000bp，而不是这一等高面对应的z的值。

Sky-walker

枫丹白露 发表于 2014-11-28 06:57
从你的描述来看，你应该是用了两个不同维数环境的数据及其描述（.ctl）文件。解决的办法是用'set dfile ... ...

我试一下哈，谢谢！！

hugeal10000

学习中，前辈们厉害！