GLDAS数据（grib格式）用wgrib查看变量时发现部分变量变成var数字

水长东

1.问题简述：我下载了GLDAS的数据，格式为grib。在Linux中用wgrib查看变量时，部分变量名变成了数字。这是为什么？

2.问题详述：展示操作步骤，引出问题。
2.1.我选择了GLDAS数据包，如图1：

图1

变量：不做设置，默认所有变量，如图2：

图2

数据格式：grib
时间：2010.06.01~2010.06.30
空间：不做设置，默认全球
由此获得文件：subset_GLDAS_NOAH025_3H_2.1_20220112_010439.txt

2.2.利用MATLAB进行下载。这里顺便共享给大家一个MATLAB下载数据的程序。
% 气象家园  水长东 %
clear;clc;close all
url=textread('subset_GLDAS_NOAH025_3H_2.1_20220112_010439.txt', '%s');
n=length(url);
for i=1:n
    web(url{i});
end

2.3.把下载到的数据导入到服务器上

2.4.安装wgrib，并在~/.bashrc声明变量，然后source命令使声明的变量立即生效

2.5.输入命令：
wgrib -v GLDAS_NOAH025_3H.A20100601.0000.021.nc4.SUB.grb
获得变量说明，如下图：

图3

每一行文字的含义如下：

图4

我认为在图3中的第15行~第29行的参数名称本来也应当和其他行的参数名称一样，应当是有实际意义的英文简写，为什么变成了var+kpds值？

是wgrib软件安装有问题，还是官网上GLDAS数据本身有问题？

如果本来就应当是参数名称本应当是var+kpds值，README_GLDAS2.pdf上没有相关内容，那么这些var+kpds值与变量的实际含义的对应关系还可以在哪里可以查询？

3.问题的拓展：
用GrADS上获取grb数据的ctl文件（见下面所附内容）后，得到各个参量名称与图2中的各个变量又是怎么对应的？

ctl中为什么有很多变量在官网（图2）上都没有？如以下变量ctl中有而官网中没有：
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
LFTXsfc  0 131,1,0  ** surface Surface lifted index [K]

官网（图2）上有存在的变量在ctl中却没有？如一下变量官网上有而ctl中没有：
SoilMoi0_10cm_inst = Soil moisture content (0-10 cm underground) (kg m-2)
SoilMoi10_40cm_inst = Soil moisture content (10-40 cm underground) (kg m-2)


为什么在2.5节中部分变量被wgrib解析成了var+kpds值，而在ctl中却可以被解析成了有意义的变量？





ctl文件内容：
dset C:/OpenGrADS/datatest/test.grb
index C:/OpenGrADS/datatest/test.grb.idx
undef 9.999E+20
title C:/OpenGrADS/datatest/test.grb
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 255
ydef 600 linear -59.875000 0.25
xdef 1440 linear -179.875000 0.250000
tdef 1 linear 00Z01jun2010 1mo
zdef 1 linear 1 1
vars 36
no4LFTXsfc  0 132,1,0  ** surface Best (4-layer) lifted index [K]
ALBDOsfc  0 84,1,0  ** surface Albedo [%]
BGRUNsfc  0 234,1,0  ** surface Baseflow-groundwater runoff [kg/m^2]
BVF2sfc  0 138,1,0  ** surface Brunt-Vaisala frequency^2 [1/s^2]
CFRZRsfc  0 141,1,0  ** surface Categorical freezing rain [yes=1;no=0]
CICEPsfc  0 142,1,0  ** surface Categorical ice pellets [yes=1;no=0]
CNWATsfc  0 223,1,0  ** surface Plant canopy surface water [kg/m^2]
CRAINsfc  0 140,1,0  ** surface Categorical rain [yes=1;no=0]
DLWRFsfc  0 205,1,0  ** surface Downward long wave flux [W/m^2]
DSWRFsfc  0 204,1,0  ** surface Downward short wave flux [W/m^2]
EVBSsfc  0 199,1,0  ** surface Direct evaporation from bare soil [W/m^2]
EVCWsfc  0 200,1,0  ** surface Canopy water evaporation [W/m^2]
EVPsfc  0 57,1,0  ** surface Evaporation [kg/m^2]
GFLUXsfc  0 155,1,0  ** surface Ground heat flux [W/m^2]
KXsfc  0 133,1,0  ** surface K index [K]
LFTXsfc  0 131,1,0  ** surface Surface lifted index [K]
LHTFLsfc  0 121,1,0  ** surface Latent heat flux [W/m^2]
MCONVsfc  0 135,1,0  ** surface Horizontal moisture divergence [kg/kg/s]
MTHAsfc  0 70,1,0  ** surface Main thermocline anomaly [m]
NLWRSsfc  0 112,1,0  ** surface Net long wave (surface) [W/m^2]
NSWRSsfc  0 111,1,0  ** surface Net short wave (surface) [W/m^2]
PEVPRsfc  0 145,1,0  ** surface Potential evaporation rate [W/m^2]
PRATEsfc  0 59,1,0  ** surface Precipitation rate [kg/m^2/s]
PRESsfc  0 1,1,0  ** surface Pressure [Pa]
PVMWsfc  0 139,1,0  ** surface Potential vorticity (mass-weighted) [1/s/m]
SHTFLsfc  0 122,1,0  ** surface Sensible heat flux [W/m^2]
SNODsfc  0 66,1,0  ** surface Snow depth [m]
SNOMsfc  0 99,1,0  ** surface Snow melt [kg/m^2]
SPFHsfc  0 51,1,0  ** surface Specific humidity [kg/kg]
SSRUNsfc  0 235,1,0  ** surface Storm surface runoff [kg/m^2]
SXsfc  0 134,1,0  ** surface Sweat index [K]
TMPsfc  0 11,1,0  ** surface Temp. [K]
TRANSsfc  0 210,1,0  ** surface Transpiration [W/m^2]
VWSHsfc  0 136,1,0  ** surface Vertical speed shear [1/s]
WEASDsfc  0 65,1,0  ** surface Accum. snow [kg/m^2]
WINDsfc  0 32,1,0  ** surface Wind speed [m/s]
ENDVARS

水长东

已解决，修改或创建gribtab_GLDAS_NOAH.gtb文件（详见https://disc.gsfc.nasa.gov/infor ... ogy%20Documentation），并在~/.bashrc中声明变量：

export GRIBTAB=.../gribtab_GLDAS_NOAH.gtb
(   ...  为路径)

再source ~/.bashrc一下就可以了。

mgyang

楼主好强阿！