MESO要素清单：280-323

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

5 层，层次列表如下：

• 1000

• 925

• 850

• 700

• 500

3 层，层次列表如下：

• 850

• 700

• 500

云量

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
云量	CDCC	ccl	cc	(0, 6, 22)	pl	21层	280-300

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

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

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

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

相对湿度

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
相对湿度	RH	r	rh	(0, 1, 1)	pl	21层	301-321

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

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

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

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

2米相对湿度

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
2米相对湿度	RH	2r	rh2m	(0, 1, 1)	heightAboveGround	2	322

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="RH",
    level_type="heightAboveGround",
    level=2,
    lazy=True
)
assert field.attrs["GRIB_count"] == 322

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="2r",
    lazy=True
)
assert field.attrs["GRIB_count"] == 322

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="rh2m",
    lazy=True
)
assert field.attrs["GRIB_count"] == 322

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

2米露点温度

要素名	wgrib2	eccodes	cemc	GRIB Key	层次类型	层次	编号
2米露点温度	DPT	2d	td2m	(0, 0, 6)	heightAboveGround	2	323

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="DPT",
    level_type="heightAboveGround",
    level=2,
    lazy=True
)
assert field.attrs["GRIB_count"] == 323

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="2d",
    lazy=True
)
assert field.attrs["GRIB_count"] == 323

field = load_field_from_file(
    meso_grib2_orig_file_path,
    parameter="td2m",
    lazy=True
)
assert field.attrs["GRIB_count"] == 323

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