学生下载了gldas v1 0.25x 0.25度的 3-hourly 数据，据介绍里面是有rainfall数据的。但用grads或者matlab read_grib 读取时，发现没有此项，如图：

（click for enlarged view）

这很奇怪，因为相关说明里分明是有rainfall一项的，如在http://hydro1.sci.gsfc.nasa.gov/dods/GLDAS_NOAH025SUBP_3H.info介绍的，

而且对比两者，里面的参数名也不一致。

仔细研究gldas数据对应的元数据，发现每个文件每项里都有对应的最大值最小值，试着读取并对比，发现第7项4LFTX即是rainf 一项，最大值一致。看样子是read_grib等代码是专门for WMO的，跟gldas的描述不一致。所以在读取gldas等数据时，还需要注意了。

在ftp://hydro1.sci.gsfc.nasa.gov/data/s4pa/GLDAS/GLDAS_NOAH025_3H_E1.002/doc/README.GLDAS2.pdf找到parameter table的相关信息，如下图

注意其中的pds ids，rainfall rate是132，对应于图1的第7个记录，#=132。

大概是因为read_grib的参数表是基于ncep的某种格式，读取gldas需要配置自定义的参数表。