兰州渭源路什字公交站的广告牌。
疯狂到连字都不会写了!(CRZILY?)
学生下载了gldas v1 0.25x 0.25度的 3-hourly 数据,据介绍里面是有rainfall数据的。但用grads或者matlab read_grib 读取时,发现没有此项,如图:
这很奇怪,因为相关说明里分明是有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需要配置自定义的参数表。