MESO要素清单：49-153

Contents

MESO要素清单：49-153#

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/ref/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/ref/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

位势高度#

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

温度#

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

东西风#

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

南北风#

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

垂直速度#

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

Note

检索 1000hPa 垂直速度需要增加 stepType="instant" 条件，与输出间隔内最大垂直速度 (编号 19) 区分。

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