GLDAS v1 0.25度3小时降水数据

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

image
(click for enlarged view)

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

image

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

仔细研究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的相关信息,如下图

image

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

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

1 thought on “GLDAS v1 0.25度3小时降水数据

Leave a Reply

Your email address will not be published. Required fields are marked *