MESO要素清单：324-358 - 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

5 层，层次列表如下：

1000
925
850
700
500

3 层，层次列表如下：

散度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
散度	div	RELD	d	(0, 2, 13)	pl	5层	324-328

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

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

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="div",
    level_type="pl",
    level=700,
    lazy=True
)
assert field.attrs["GRIB_count"] == 327

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

温度平流¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
温度平流	tadv	-	-	(0, 0, 224)	pl	3层	329-331

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

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="tadv",
    level_type="pl",
    level=700,
    lazy=True
)
assert field.attrs["GRIB_count"] == 330

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

露点温度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
露点温度	td	DPT	dpt	(0, 0, 6)	pl	21层	332-352

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

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

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

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

温度露点差¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
温度露点差	ttd	DEPR	-	(0, 0, 7)	pl	3层	353-355

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

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="ttd",
    level_type="pl",
    level=700,
    lazy=True
)
assert field.attrs["GRIB_count"] == 354

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

假相当位温¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
假相当位温	thetaSe	EPOT	papt	(0, 0, 3)	pl	3层	356-358

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

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

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="thetaSe",
    level_type="pl",
    level=700,
    lazy=True
)
assert field.attrs["GRIB_count"] == 357

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