MESO要素清单：31-48 - 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

meso_grib2_orig_file_path = get_data_file("cma-meso-1km")
meso_grib2_orig_file_path

PosixPath('/home/wangdp/project/cedarkit/notebook-project/notebook-devel/repo/data-notebook-project/data-notebook/data/cma-meso-1km/grib2/orig/Z_NAFP_C_BABJ_20260725000000_P_CEMC-CMA-MESO-1KM-ORIG-02400.grb2')

等高面u风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
等高面u风	u	UGRD	u	(0, 2, 2)	heightAboveGround	30, 50, 70	31-33
等高面u风	u	UGRD	100u	(0, 2, 2)	heightAboveGround	100	34
等高面u风	u	UGRD	u	(0, 2, 2)	heightAboveGround	120, 140, 160, 180	35-38
等高面u风	u	UGRD	200u	(0, 2, 2)	heightAboveGround	200	39

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

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

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

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

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

等高面v风¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
等高面v风	v	VGRD	v	(0, 2, 3)	heightAboveGround	30, 50, 70	40-42
等高面v风	v	VGRD	100v	(0, 2, 3)	heightAboveGround	100	43
等高面v风	v	VGRD	v	(0, 2, 3)	heightAboveGround	120, 140, 160, 180	44-47
等高面v风	v	VGRD	200v	(0, 2, 3)	heightAboveGround	200	48

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

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

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

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

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