GFS要素清单：64-223 - CEMC Data Guide

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/repo/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/repo/data-notebook-project/data-notebook/data/cma-gfs/grib2/orig/Z_NAFP_C_BABJ_20260811000000_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¶

位势高度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
位势高度	h	HGT	gh	(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

The history saving thread hit an unexpected error (OperationalError('attempt to write a readonly database')).History will not be written to the database.

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

温度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
温度	t	TMP	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

东西风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
东西风	u	UGRD	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

南北风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
南北风	v	VGRD	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