MESO要素清单：49-153 - 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')

层次¶

21 层，层次列表如下：

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

位势高度¶

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

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

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

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="h",
    level_type="pl",
    level=200,
    lazy=True
)
assert field.attrs["GRIB_count"] == 67

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

温度¶

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

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

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

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

东西风¶

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

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

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

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

南北风¶

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

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

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

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

垂直速度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
垂直速度	w	DZDT	wz	(0, 2, 9)	pl	21层	133-153

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="DZDT",
    level_type="pl",
    level=1000,
    stepType="instant",
    lazy=True
)
assert field.attrs["GRIB_count"] == 133

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="DZDT",
    level_type="pl",
    level=975,
    lazy=True
)
assert field.attrs["GRIB_count"] == 134

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="wz",
    level_type="pl",
    level=850,
    lazy=True
)
assert field.attrs["GRIB_count"] == 138

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="w",
    level_type="pl",
    level=200,
    lazy=True
)
assert field.attrs["GRIB_count"] == 151

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 2,
        "parameterNumber": 9,
    },
    level_type="pl",
    level=100,
    lazy=True
)
assert field.attrs["GRIB_count"] == 153