GLDAS文件用gards批量读取(已解决)

wenzerg

    求助：近日下载了2009年全年的GLDAS三小时一次的grb格式文件：
文件名如下类似：

单个时次的文件已经用grib2ctl做好了ctl文件和gribmap.exe做好idx文件，可以用grads打开，现在的问题是：
如何用单个ctl打开多时次的文件名无规则的grb文件？多个时次的idx文件也生成不了

我生成的ctl如下：
dset g:/data/GLDAS/3h/GLDAS_NOAH025SUBP_3H.A2009001.0000.001.2012139232511.grb
index g:/data/GLDAS/3h/GLDAS_NOAH025SUBP_3H.A2009001.0000.001.2012139232511.grb.idx

options  byteswapped template
undef 9.999E+20
title g:/data/GLDAS/3h/GLDAS_NOAH025SUBP_3H.A2009001.0000_01.grb
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 0
ydef 600 linear -59.875000 0.25
xdef 1440 linear -179.875000 0.250000
tdef 2 linear 21Z31dec2008 3hr
zdef 1 linear 1 1
vars 28
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
BGRUNsfc  0 234,1,0  ** surface Baseflow-groundwater runoff [kg/m^2]
BVF2sfc  0 138,1,0  ** surface Brunt-Vaisala frequency^2 [1/s^2]

.........
.........
ENDVARS

GLDAS文件命名规则太奇怪了，红色部分的序号是日历日，不是常用的年月日，也就是2009年最后一个文件是：GLDAS_NOAH025SUBP_3H.A2009365.2100.001.2010048205332.grb

若把ctl文件第一行改为：dset g:/data/GLDAS/3h/GLDAS_NOAH025SUBP_3H.A2009%d3.%h200.001.2012139232511.grb
还是打开出错。。。。。。。。
这样的文件如何用grads用一个ctl控制？

##############分割线###################问题解决如下###################
2楼说的对，只要文件名命名有规则并且是按“年-月-日-时-分“这种命名的话，就可以批量读取了。经过探索，困难就出现在如何批量重命名上。
找了好久，找到了款“拖把更名器”的工具，网上有很多，我就不贴下载地址了。
然后自己写了个文本文件，文件名写成如下规则：
GLDAS_2009-01-01_00_00
GLDAS_2009-01-01_03_00
GLDAS_2009-01-01_06_00
GLDAS_2009-01-01_09_00
GLDAS_2009-01-01_12_00
GLDAS_2009-01-01_15_00
GLDAS_2009-01-01_18_00
GLDAS_2009-01-01_21_00
GLDAS_2009-01-02_00_00
GLDAS_2009-01-02_03_00
GLDAS_2009-01-02_06_00
GLDAS_2009-01-02_09_00
GLDAS_2009-01-02_12_00
GLDAS_2009-01-02_15_00
GLDAS_2009-01-02_18_00
GLDAS_2009-01-02_21_00
GLDAS_2009-01-03_00_00
GLDAS_2009-01-03_03_00
GLDAS_2009-01-03_06_00
GLDAS_2009-01-03_09_00
GLDAS_2009-01-03_12_00
GLDAS_2009-01-03_15_00


......
......
GLDAS_2009-12-31_21_00
用拖把更名器中“调用文本文件命名”功能，把上述文件名全部更换。


基本快大功告成，再把ctl文件稍作修改：红色部分为修改为能批量生成idx文件的地方：
dset g:/data/GLDAS/3h/GLDAS_2009-%m2-%d2_%h2_00.grb
index g:/data/GLDAS/3h/GLDAS_2009-3h.idx
options  byteswapped template
undef 9.999E+20
title g:/data/GLDAS/3h/GLDAS_2009-3h.grb
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 0
ydef 600 linear -59.875000 0.25
xdef 1440 linear -179.875000 0.250000
tdef 2920 linear 00Z01jan2009 3hr
zdef 1 linear 1 1
vars 28


最后执行 “！gribmap -e -i G:\data\GLDAS\3h\GLDAS_2009-3h.ctl" 生成idx文件，大功已告成！

贴两幅GLDAS的图留做纪念！是noah模式生成的全球0.25度陆面同化数据。分辨率还是很高的

nuist88513

没用过，帮顶下~

river

它的命名应该是为了比较方便看出资料的更多内容吧，你可以修改成有规则的文件名然后再用一个ctl描述吧

yushaotian

受用了，顶一个！

幻想的春

这个挺好的
意思就是可以批量处理GRADS的数据了？

humes

谢谢楼主，认真学习一下！

shrily521

vars 28，是说明变量时28个吧，求详细解释这28个量的具体名称。或者给个介绍了的网址也行啊。我最近在用这个gldas的数据，但是还不太了解这些数据，求指导。谢谢

shrily521

index g:/data/GLDAS/3h/GLDAS_NOAH025SUBP_3H.A2009001.0000.001.2012139232511.grb.idx
这个是生成.idx数据的意思？这个文件时用来干什么的？

shrily521

还想问下楼主，那个用grads怎么处理？要是已经生成了index和ctl怎么画图呢？

272363310

楼主好样的要是把代码贴出来更好