Thursday, January 08, 2009

ncl does not recognize gldas grib files

wgrib gives:

tianyf@tyf ~/gldas]$ wgrib -s GLDAS_NOAH025_M.A200701.001.grb
1:0:d=07010100:NSWRS:sfc:0-3hr ave:NAve=12
2:563720:d=07010100:NLWRS:sfc:0-3hr ave:NAve=12
3:1097066:d=07010100:LHTFL:sfc:0-3hr ave:NAve=12
4:1660786:d=07010100:SHTFL:sfc:0-3hr ave:NAve=12

Undefined parameter table (center 221-221 table 1), using NCEP-opn
5:2224506:d=07010100:GFLUX:sfc:0-3hr ave:NAve=12
6:2727476:d=07010100:LFTX:sfc:anl:NAve=12
7:3291196:d=07010100:4LFTX:sfc:anl:NAve=12
8:3885292:d=07010100:EVP:sfc:anl:NAve=12
9:4388262:d=07010100:SSRUN:sfc:anl:NAve=12
10:4951982:d=07010100:BGRUN:sfc:anl:NAve=12
11:5546078:d=07010100:SNOM:sfc:0-3hr ave:NAve=12
12:5988298:d=07010100:BVF2:sfc:anl:NAve=12
13:6521644:d=07010100:WEASD:sfc:anl:NAve=12
14:7115740:d=07010100:TSOIL:0-4 cm down:anl:NAve=12
15:7649086:d=07010100:TSOIL:0-3 cm down:anl:NAve=12
16:8182432:d=07010100:TSOIL:0-2 cm down:anl:NAve=12
17:8715778:d=07010100:TSOIL:0-1 cm down:anl:NAve=12
18:9218748:d=07010100:SOILM:0-4 cm down:anl:NAve=12
19:9812844:d=07010100:SOILM:0-3 cm down:anl:NAve=12
20:10406940:d=07010100:SOILM:0-2 cm down:anl:NAve=12
21:11001036:d=07010100:SOILM:0-1 cm down:anl:NAve=12
22:11595132:d=07010100:TCDC:sfc:anl:NAve=12
23:12098102:d=07010100:WIND:sfc:anl:NAve=12
24:12540322:d=07010100:TMP:sfc:anl:NAve=12
25:13043292:d=07010100:SPFH:sfc:anl:NAve=12
26:13303258:d=07010100:PRES:sfc:anl:NAve=12
27:13897354:d=07010100:DSWRF:sfc:0-3hr ave:NAve=12
28:14491450:d=07010100:DLWRF:sfc:0-3hr ave:NAve=12

However, when convert it to .nc format using ncl,

[tianyf@tyf ~/gldas]$ ncl_convert2nc GLDAS_NOAH025_M.A200701.001.grb -v weasd
Processing file: GLDAS_NOAH025_M.A200701.001.grb...
Copyright (C) 1995-2007 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 5.0.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
warning:NclGRIB: Unknown grib parameter number detected (155, center 221, table version 1 grib record 5), using default variable name (VAR_155)
warning:NclGRIB: Unknown grib parameter number detected (131, center 221, table version 1 grib record 6), using default variable name (VAR_131)
warning:NclGRIB: Unknown grib parameter number detected (132, center 221, table version 1 grib record 7), using default variable name (VAR_132)
warning:NclGRIB: Unknown grib parameter number detected (235, center 221, table version 1 grib record 9), using default variable name (VAR_235)
warning:NclGRIB: Unknown grib parameter number detected (234, center 221, table version 1 grib record 10), using default variable name (VAR_234)
warning:NclGRIB: Unknown grib parameter number detected (138, center 221, table version 1 grib record 12), using default variable name (VAR_138)
warning:NclGRIB: Unknown grib parameter number detected (204, center 221, table version 1 grib record 27), using default variable name (VAR_204)
warning:NclGRIB: Unknown grib parameter number detected (205, center 221, table version 1 grib record 28), using default variable name (VAR_205)
warning:getfilevartypes: variable (weasd) does not exist in file (GLDAS_NOAH025_M.A200701.001)
fatal:The result of the conditional expression yields a missing value. NCL can not determine branch, see ismissing function
fatal:Execute: Error occurred at or near line 351 in file /tmp/tmp8895.ncl


ncl does handle other grib file correctly,

[tianyf@tyf ~/ds607]$ ncl_convert2nc 19981201.NOAH_d.grb
Processing file: 19981201.NOAH_d.grb...
Copyright (C) 1995-2007 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 5.0.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.

Thus, I can only use wgrib to extract variables from grib files.

wgrib -s GLDAS_NOAH025_M.A200701.001.grb | grep WEASD | wgrib -i -s -text GLDAS_NOAH025_M.A200701.001.grb -o o

No comments: