flux读取不到数据_科学网—读取GLDAS数据出错问题的解决 - 周锋的博文

读取GLDAS数据出错问题的解决

已有 8053 次阅读

2011-10-18 10:44

||系统分类:科研笔记|

QOCA, mload, wgrib, GLDAS

用wgrib读取GLDAS的grb格式的数据时，会出现如下错误：Undefined parameter table (center 221-221 table 1), using NCEP-opn。

解决方案如下：

wgrib的安装

从ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib/下载最新的安装包，目前最新的是wgrib.tar.v1.8.1.2a，解压后编辑makefile，

CC=gcc

make即可，将生成的可执行程序wgrib拷贝到/usr/local/bin下即可。

终端输入wgrib显示如下：

Portable Grib decoder for NCEP/NCAR Reanalysis etc.

it slices, dices    v1.8.1.0b (5-10) Wesley Ebisuzaki

DWD-tables 2,201-205 (11-28-2005) Helmut P. Frank

spectral: Luis Kornblueh (MPI)

usage: wgrib [grib file] [options]

Inventory/diagnostic-output selections

-s/-v                   short/verbose inventory

-V                      diagnostic output (not inventory)

(none)                  regular inventory

Options

-PDS/-PDS10             print PDS in hex/decimal

-GDS/-GDS10             print GDS in hex/decimal

-verf                   print forecast verification time

-ncep_opn/-ncep_rean    default T62 NCEP grib table

-4yr                    print year using 4 digits

-min                    print minutes

-ncep_ens               ensemble info encoded in ncep format

Decoding GRIB selection

-d [record number|all]  decode record number

-p [byte position]      decode record at byte position

-i                      decode controlled by stdin (inventory list)

(none)                  no decoding

Options

-text/-ieee/-grib/-bin  convert to text/ieee/grib/bin (default)

-nh/-h                  output will have no headers/headers (default)

-dwdgrib                output dwd headers, grib (do not append)

-H                      output will include PDS and GDS (-bin/-ieee only)

-append                 append to output file

-o [file]               output file name, 'dump' is default

Misc

-cmc [file]             use NCEP tables for CMC (dangerous)

即证明正确安装。

用一例子验证下，wgrib GLDAS_NOAH10_M.A201001.001.grb

显示如下：

1:0:d=10010100:NSWRS:kpds5=111:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

2:35306:d=10010100:NLWRS:kpds5=112:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

3:68714:d=10010100:LHTFL:kpds5=121:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

4:104020:d=10010100:SHTFL:kpds5=122:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

Undefined parameter table (center 221-221 table 1), using NCEP-opn

5:139326:d=10010100:GFLUX:kpds5=155:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

6:168938:d=10010100:LFTX:kpds5=131:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

7:202346:d=10010100:4LFTX:kpds5=132:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

8:239550:d=10010100:EVP:kpds5=57:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

9:271060:d=10010100:SSRUN:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

10:304468:d=10010100:BGRUN:kpds5=234:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

11:339774:d=10010100:SNOM:kpds5=99:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

12:371284:d=10010100:BVF2:kpds5=138:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

13:404692:d=10010100:WEASD:kpds5=65:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

14:457078:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12

15:490486:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12

16:523894:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12

17:557302:d=10010100:TSOIL:kpds5=85:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12

18:588812:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12

19:626016:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12

20:663220:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12

21:700424:d=10010100:SOILM:kpds5=86:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12

22:737628:d=10010100:TCDC:kpds5=71:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

23:769138:d=10010100:WIND:kpds5=32:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

24:794956:d=10010100:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

25:826466:d=10010100:SPFH:kpds5=51:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

26:842794:d=10010100:PRES:kpds5=1:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

27:879998:d=10010100:DSWRF:kpds5=204:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

28:917202:d=10010100:DLWRF:kpds5=205:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

出现了这样的错误Undefined parameter table (center 221-221 table 1), using NCEP-opn

无论加 -necp_opn 还是 -necp_rean参数都出现上述错误。

解决方案如下(以我的电脑设置为主，过程一致，你可以参照你的电脑设置，只要你掌握基本的linux知识)：http://disc.sci.gsfc.nasa.gov/hydrology/documentation/hydro_doc.shtml#NOAH下载GRIBTAB for Noah output 表文件 gribtab_GLDAS_NOAH.txt重命名为 gribtab_GLDAS_NOAH 放入root下的data目录，在.bashrc下设置

# GRIBTAB for Noah output

export GRIBTAB=/root/data/gribtab_GLDAS_NOAH

这样，重新运行下 wgrib GLDAS_NOAH10_M.A201001.001.grb

显示：

1:0:d=10010100:SWnet:kpds5=111:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

2:35306:d=10010100:LWnet:kpds5=112:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

3:68714:d=10010100:Qle:kpds5=121:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

4:104020:d=10010100:Qh:kpds5=122:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

5:139326:d=10010100:Qg:kpds5=155:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

6:168938:d=10010100:Snowf:kpds5=131:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

7:202346:d=10010100:Rainf:kpds5=132:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

8:239550:d=10010100:Evap:kpds5=57:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

9:271060:d=10010100:Qs:kpds5=235:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

10:304468:d=10010100:Qsb:kpds5=234:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

11:339774:d=10010100:Qsm:kpds5=99:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

12:371284:d=10010100:AvgSurfT:kpds5=138:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

13:404692:d=10010100:SWE:kpds5=65:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

14:457078:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12

15:490486:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12

16:523894:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12

17:557302:d=10010100:SoilTemp4:kpds5=85:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12

18:588812:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=4:TR=0:P1=0:P2=0:TimeU=1:0-4 cm down:anl:NAve=12

19:626016:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=3:TR=0:P1=0:P2=0:TimeU=1:0-3 cm down:anl:NAve=12

20:663220:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:0-2 cm down:anl:NAve=12

21:700424:d=10010100:SoilMoist4:kpds5=86:kpds6=112:kpds7=1:TR=0:P1=0:P2=0:TimeU=1:0-1 cm down:anl:NAve=12

22:737628:d=10010100:Canopint:kpds5=71:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

23:769138:d=10010100:Wind:kpds5=32:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

24:794956:d=10010100:Tair:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

25:826466:d=10010100:Qair:kpds5=51:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

26:842794:d=10010100:PSurf:kpds5=1:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=12

27:879998:d=10010100:SWdown:kpds5=204:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

28:917202:d=10010100:LWdown:kpds5=205:kpds6=1:kpds7=0:TR=3:P1=0:P2=3:TimeU=1:sfc:0-3hr ave:NAve=12

错误解决，且NOAH提供的28个变量正确，正确解决这些问题，利用qoca下的模块mload将可以利用GLDAS下的数据做各项改正。wgrib的命令如何使用，可参照相关说明文档，在此不做说明！

相关的网址链接：

ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib/

http://www.cpc.ncep.noaa.gov/products/wesley/wgrib.html

http://agdisc.gsfc.nasa.gov/dods/GLDAS_NOAH10_M

http://disc.sci.gsfc.nasa.gov/hydrology/documentation/hydro_doc.shtml#NOAH

http://www.nco.ncep.noaa.gov/pmb/docs/on388/

转载本文请联系原作者获取授权，同时请注明本文来自周锋科学网博客。

链接地址：http://blog.sciencenet.cn/blog-350278-498058.html

上一篇：IGS台站在中国(GMT画图)

下一篇：终于用GMT搞定全球部分IGS台站分布图了

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！