GrADS读取若干nc文件转换并合并为二进制文件

速冻。八爪鱼

展望2014 发表于 2014-9-1 08:48
slp就不用循环了吧。还是六楼的方法啊，输出一个变量以后再输出另外一个变量。不是说每个变量都要循环的。

slp在Z层不用循环是肯定的，但是t需要循环，而且其他变量有Z这个维度而slp没有，我在想是不是要用dfile

lqouc

速冻。八爪鱼 发表于 2014-9-1 09:11
slp在Z层不用循环是肯定的，但是t需要循环，而且其他变量有Z这个维度而slp没有，我在想是不是要用dfile

这个用不用dfile都没啥大区别，输出的时候时间是在最外层，其次就是变量，每一个变量有多少层高度互相不影响，你在ctl里面描述清楚就行了。

展望2014

速冻。八爪鱼 发表于 2014-9-1 09:11
slp在Z层不用循环是肯定的，但是t需要循环，而且其他变量有Z这个维度而slp没有，我在想是不是要用dfile

那就把slp套在t循环里面就行了，但是z一定要设定好，也就是说一定要有‘set z ’语句

展望2014

lqouc 发表于 2014-9-1 09:15
这个用不用dfile都没啥大区别，输出的时候时间是在最外层，其次就是变量，每一个变量有多少层高度互相不 ...

正解！

速冻。八爪鱼

lqouc 发表于 2014-9-1 09:15
这个用不用dfile都没啥大区别，输出的时候时间是在最外层，其次就是变量，每一个变量有多少层高度互相不 ...

能帮忙看下gs文件么
"reinit"
"sdfopen /home/air.1948.nc"
"sdfopen /home/slp.1948.nc"
"set gxout fwrite"
"set fwrite /home/combine.1948.grd"
"set x 1 144"
"set y 1 73"


tt=1
while (tt <= 366)
"set t "tt""

zz=1
while (zz <= 17)
"set z "zz""

"d air.1"

zz=zz+1
endwhile

在这里z已经被定义了，然后导致出现错误，

Data Request Warning:  Request is completely outside file limits

slp没有Z维度，这样如何处理之前这个Z的定义？

是不是用dfile 2 然后把维度重新定义一下？但是t循环里面用dfile可以么？

"d slp.2"


tt=tt+1
endwhile

"disable fwrite"
;

速冻。八爪鱼

展望2014 发表于 2014-9-1 09:17
那就把slp套在t循环里面就行了，但是z一定要设定好，也就是说一定要有‘set z ’语句

set z 0?
因为slp这个变量只有lat lon time三个维度，定义Z的话会出现错误:Data Request Warning:  Request is completely outside file limits

lqouc

速冻。八爪鱼 发表于 2014-9-1 09:21
能帮忙看下gs文件么
"reinit"
"sdfopen /home/air.1948.nc"

"d slp.2"前面加一个对level的set，怎么写你看slp的ctiinfo
实在不行用reset清空维度设置

展望2014

速冻。八爪鱼 发表于 2014-9-1 09:22
set z 0?
因为slp这个变量只有lat lon time三个维度，定义Z的话会出现错误:Data Request Warning:  Requ ...

在grads里面任何数据都是4维的。‘set z 1’ 试一下

速冻。八爪鱼

lqouc 发表于 2014-9-1 09:26
"d slp.2"前面加一个对level的set，怎么写你看slp的ctiinfo
实在不行用reset清空维度设置

slp 0 99 sea level pressure(hpa)

速冻。八爪鱼

展望2014 发表于 2014-9-1 09:28
在grads里面任何数据都是4维的。‘set z 1’ 试一下

我在我之前的gs文件红字部分加上了
zz=1
"set z "zz""
"d slp.2"

运行gs出现同样的错误提示。真是不知道怎么办了