GFS要素清单：64-223

import numpy as np
import pandas as pd
import xarray as xr

from reki.format.grib.eccodes import load_field_from_file

/home/wangdp/project/cedarkit/notebook-project/notebook-devel/ref/data-notebook-project/data-notebook/.venv/lib/python3.14/site-packages/gribapi/__init__.py:23: UserWarning: ecCodes 2.42.0 or higher is recommended. You are running version 2.34.1
  warnings.warn(

from data_notebook.data import get_data_file

gfs_grib2_orig_file_path = get_data_file("cma-gfs")
gfs_grib2_orig_file_path

PosixPath('/home/wangdp/project/cedarkit/notebook-project/notebook-devel/ref/data-notebook-project/data-notebook/data/cma-gfs/grib2/orig/Z_NAFP_C_BABJ_20260725000000_P_NWPC-GRAPES-GFS-GLB-02400.grib2')

层次

40 层列表，单位百帕 (mb)：

• 1000

• 975

• 950

• 925

• 900

• 850

• 800

• 750

• 700

• 650

• 600

• 550

• 500

• 450

• 400

• 350

• 300

• 275

• 250

• 225

• 200

• 175

• 150

• 125

• 100

• 70

• 50

• 30

• 20

• 10

• 7

• 5

• 3

• 2

• 1.5

• 1

• 0.5

• 0.2

• 0.1

HGT/TMP/UGRD/VGRD

位势高度

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
位势高度	HGT	gh	h	(0, 3, 5)	pl	40层	64-103

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="HGT",
    level_type="pl",
    level=1000,
    lazy=True
)
assert field.attrs["GRIB_count"] == 64

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="gh",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 69

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="h",
    level_type="pl",
    level=1.5,
    lazy=True
)
assert field.attrs["GRIB_count"] == 99

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 3,
        "parameterNumber": 5,
    },
    level_type="pl",
    level=0.1,
    lazy=True
)
assert field.attrs["GRIB_count"] == 103

温度

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
温度	TMP	t	t	(0, 0, 0)	pl	40层	104-143

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="TMP",
    level_type="pl",
    level=1000,
    lazy=True
)
assert field.attrs["GRIB_count"] == 104

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="t",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 109

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="t",
    level_type="pl",
    level=1.5,
    lazy=True
)
assert field.attrs["GRIB_count"] == 139

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 0,
        "parameterNumber": 0,
    },
    level_type="pl",
    level=0.1,
    lazy=True
)
assert field.attrs["GRIB_count"] == 143

东西风

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
东西风	UGRD	u	u	(0, 2, 2)	pl	40层	144-183

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="UGRD",
    level_type="pl",
    level=1000,
    lazy=True
)
assert field.attrs["GRIB_count"] == 144

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="u",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 149

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="u",
    level_type="pl",
    level=1.5,
    lazy=True
)
assert field.attrs["GRIB_count"] == 179

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 2,
        "parameterNumber": 2,
    },
    level_type="pl",
    level=0.1,
    lazy=True
)
assert field.attrs["GRIB_count"] == 183

南北风

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
南北风	VGRD	v	v	(0, 2, 3)	pl	40层	184-223

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="VGRD",
    level_type="pl",
    level=1000,
    lazy=True
)
assert field.attrs["GRIB_count"] == 184

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="v",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 189

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="v",
    level_type="pl",
    level=1.5,
    lazy=True
)
assert field.attrs["GRIB_count"] == 219

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 2,
        "parameterNumber": 3,
    },
    level_type="pl",
    level=0.1,
    lazy=True
)
assert field.attrs["GRIB_count"] == 223