困扰！怎样用GrADS将多个nc文件合并（为了将其中的相同要素化成一个序列）

晚晴山

             最近用wget批量下载了若干.nc的文件，但是下下来才发现这些文件是一天一个压缩包！现在想要将连续十天的（比如温度）化成一个时间序列，批量读入是没有问题的，问题是GrADS没有数组的概念，怎样将10个不同文件中的某个要素量进行合并咧。我语言表达不好，不知道大家能不能听懂。总之之前我都是用笨方法，就是先用GrADS提取，然后再用fortran合并。不过现在文件实在太多了。。。。。。。。。。
            跪求会的大神点拨一下。

mofangbao

先打开一个准备要写入的文件，可以是grads的二进制文件，也可以是nc格式，具体语法要自己上官网看，然后用循环依次打开每个要合并的文件（注意及时关闭），设置当前写入的维度输出就行了，论坛也有类似的例子，要自己搜索一下。最方便的还是用MeteoInfo合并了

晚晴山

mofangbao 发表于 2015-8-5 21:55
先打开一个准备要写入的文件，可以是grads的二进制文件，也可以是nc格式，具体语法要自己上官网看，然后用 ...

好，谢谢版主这么晚还回复我~

兰溪之水

写ctl批量读取不就好了？

晚晴山

兰溪之水 发表于 2015-8-5 22:08
写ctl批量读取不就好了？

其实是不行的。。。。。。

兰溪之水

晚晴山 发表于 2015-8-5 22:11
其实是不行的。。。。。。

单个文件能读？

晚晴山

兰溪之水 发表于 2015-8-6 08:01
单个文件能读？

主要的问题出在用的MERRA的资料 它的xdef不是线性的 而是像 xdef levels -180 -179.333 -178.667........一共540个。。。。。。。写死个我了要。昨晚尝试了用MI合并，但是文件太多 太大，一个晚上都还没合并完成 ，今天看来要硬着头皮把x维的540个levels一个一个敲进去了！！！

四叶草

本人感觉最简单就是用set gxout fwrite 的命令提取，不必再用fortran编码读写。最近我刚做，我是一次性合并18个文件中的任意变量，只要成功提取一个时次（层次），其他时次（层次）的就是复制、粘贴就能搞定。方法虽笨，但比较好理解。

晚晴山

四叶草 发表于 2015-8-6 10:24
本人感觉最简单就是用set gxout fwrite 的命令提取，不必再用fortran编码读写。最近我刚做，我是一次性合并 ...

能否gs给我看一下，或发我邮箱627967828@qq.com   万谢~

四叶草

我发邮件给你，我这边帖子不能添加附件……