GEPS要素清单：1-85 - 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

geps_grib2_orig_file_path = get_data_file("cma-geps")
geps_grib2_orig_file_path

PosixPath('/home/wangdp/project/cedarkit/notebook-project/notebook-devel/repo/data-notebook-project/data-notebook/data/cma-geps/grib2/orig/Z_NAFP_C_BABJ_20260811000000_P_NWPC-GRAPES-GEPS-GLB-02400-m001.grib2')

等压面：位势高度/温度/风/比湿¶

层次¶

17 层列表（百帕 (hPa)）：

1000
975
950
925
900
850
800
750
700
650
600
550
500
400
300
200
100

位势高度¶

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

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="HGT",
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 1

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="gh",
    level_type="pl",
    level=550,
    lazy=True
)
assert field.attrs["GRIB_count"] == 6

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="h",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 12

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

温度¶

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

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="TMP",
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 18

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="t",
    level_type="pl",
    level=550,
    lazy=True
)
assert field.attrs["GRIB_count"] == 23

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

东西风¶

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

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="UGRD",
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 35

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="u",
    level_type="pl",
    level=550,
    lazy=True
)
assert field.attrs["GRIB_count"] == 40

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

南北风¶

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

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="VGRD",
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 52

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="v",
    level_type="pl",
    level=550,
    lazy=True
)
assert field.attrs["GRIB_count"] == 57

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

比湿¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
比湿	q	SPFH	q	(0, 1, 0)	pl	17层	69-85

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="SPFH",
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 69

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter="q",
    level_type="pl",
    level=550,
    lazy=True
)
assert field.attrs["GRIB_count"] == 74

field = load_field_from_file(
    geps_grib2_orig_file_path,
    parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 0},
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 80