MESO要素清单：280-323 - 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	层次类型	层次	编号
云量	cc	CDCC	ccl	(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

相对湿度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
相对湿度	rh	RH	r	(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米相对湿度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
2米相对湿度	rh2m	RH	2r	(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米露点温度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
2米露点温度	td2m	DPT	2d	(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