GrADS中close命令对提取中国区域数据的影响

秋含章

'reinit'
'set grid on'
'set grads off'
'sdfopen e:\pr_historical_185001-200512.nc'
'set time jan1961'
'q dims'
a11=sublin(result,5)
b11=subwrd(a11,9)
'set time dec1990'
'q dims'
a12=sublin(result,5)
b12=subwrd(a12,9)
'define p1=ave(pr,t='b11',t='b12')*3600*24'
'close 1'
'sdfopen e:\pr_200601-210112.nc'
'set time jan2071'
'q dims'
c21=sublin(result,5)
d21=subwrd(c21,9)
'set time dec2100'
'q dims'
c22=sublin(result,5)
d22=subwrd(c22,9)
'define p2=ave(pr,t='d21',t='d22')*3600*24'
'define p20=p2-p1'
'close 1'

'set parea 0.5 10.5 1.5 7'
'set lon 70 140'
'set lat 15 55'
'set gxout shaded'
'set mpdset cnworld'
'set display color white'
'set xlint 10'
'set ylint 10'
'set map 1 1 3'

'cnbasemap p20'
'printim d:\gradsout\change-p.gif white'
;
若不加红色“close 1”可以正常绘制中国区域；加了这个命令得到的图如上。由于实际工作中需要处理的文件较多，需要把计算好的关掉。
所以想问问：这个close怎么会对提取中国区域数据产生影响~~~

冬日情愫+

close对数据是不会产生影响的。在这里补充一下，close只能关闭最后打开的那个文件，加入你先后打开的2个文件，只能先关闭第二个文件，在关闭第一个文件，也就是说在两个文件同时都是开着的状态，如果用close  1是是会报错的，只能先close 2。打开多个文件定义和打开一个文件是不同的，其实你没有必要一定要关闭之前的文件，只要在定义变量的时候加上文件号就好，如define p2=ave(pr.2,t=1,t=10)*3600*24

mofangbao

楼主写了这么多有没有问一下自己每一行都是啥意思呢？如果明白了那就有答案了。如果不明白，得去看GrADS官方手册了，都有解释的。

秋含章

mofangbao 发表于 2014-9-24 14:23
楼主写了这么多有没有问一下自己每一行都是啥意思呢？如果明白了那就有答案了。如果不明白，得去看GrADS官 ...

我也用了一阵子GrADS，不过之前画图一般都是打开一个文件，或者虽然有多个文件但是都是对应的元素做运算：所以一般不涉及到关闭文件的问题。
这次要在多个文件里分别取出不同时次的数据做运算。这是刚刚起步，就遇到了个困难。
手册翻了一下（看的是大气所的那本），还是不知道为什么关一下文件会影响图的输出……
这个是删掉CLOSE 1得到的图，也就是我想要得到的图

                               
登录/注册后可看大图

帮忙看看呗~~~~
谢谢啦

冬日情愫+

看了一下楼主的gs文件，我其实奇怪 你加上close 1 难道程序不会提示出错吗，你把两个文件都关闭了，下面还设置经纬度，都没有文件了，设置经纬度还有用吗。比较同意清风的说法。

秋含章

冬日情愫+ 发表于 2014-9-24 17:47
看了一下楼主的gs文件，我其实奇怪 你加上close 1 难道程序不会提示出错吗，你把两个文件都关闭了，下面还 ...

我把经纬度的设置当成单纯的出图设置了……经纬度设置应该是针对文件的操作，是这意思吧？边用边学的，基础确实不扎实~~~
确实报错来的，但我看见图出来了，就没太重视。
谢谢两位啦！

秋含章

冬日情愫+ 发表于 2014-9-24 17:25
close对数据是不会产生影响的。在这里补充一下，close只能关闭最后打开的那个文件，加入你先后打开的2个文 ...

关文件这个事情也是迫不得已啊：因为一些设置，比如set t啥的，都是针对缺省文件的吧（？），如果不关掉，就不让设置其他文件的维度环境了……偏偏我那些维度还都多少有些差别，不能统一设置。