GFS要素清单：46-63 - 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')

UGRD/VGRD¶

等高面u风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
等高面u风	u	UGRD	u	(0, 2, 2)	heightAboveGround	30, 50, 70	46-48
等高面u风	u	UGRD	100u	(0, 2, 2)	heightAboveGround	100	49
等高面u风	u	UGRD	u	(0, 2, 2)	heightAboveGround	120, 140, 160, 180	50-53
等高面u风	u	UGRD	200u	(0, 2, 2)	heightAboveGround	200	54

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="UGRD",
    level_type="heightAboveGround",
    level=30,
    lazy=True
)
assert field.attrs["GRIB_count"] == 46

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="u",
    level_type="heightAboveGround",
    level=50,
    lazy=True
)
assert field.attrs["GRIB_count"] == 47

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="100u",
    lazy=True
)
assert field.attrs["GRIB_count"] == 49

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": 2,
        "parameterNumber": 2,
    },
    level_type="heightAboveGround",
    level=120,
    lazy=True
)
assert field.attrs["GRIB_count"] == 50

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="200u",
    lazy=True
)
assert field.attrs["GRIB_count"] == 54

等高面v风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
等高面v风	v	VGRD	v	(0, 2, 3)	heightAboveGround	30, 50, 70	55-57
等高面v风	v	VGRD	100v	(0, 2, 3)	heightAboveGround	100	58
等高面v风	v	VGRD	v	(0, 2, 3)	heightAboveGround	120, 140, 160, 180	59-62
等高面v风	v	VGRD	200v	(0, 2, 3)	heightAboveGround	200	63

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="VGRD",
    level_type="heightAboveGround",
    level=30,
    lazy=True
)
assert field.attrs["GRIB_count"] == 55

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="v",
    level_type="heightAboveGround",
    level=50,
    lazy=True
)
assert field.attrs["GRIB_count"] == 56

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="100v",
    lazy=True
)
assert field.attrs["GRIB_count"] == 58

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

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="200v",
    lazy=True
)
assert field.attrs["GRIB_count"] == 63