MESO要素清单：154-279 - 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	层次类型	层次	编号
比湿	q	SPFH	q	(0, 1, 0)	pl	21层	154-174

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

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

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

云水混合比¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
云水混合比	qc	CLMR	clwmr	(0, 1, 22)	pl	21层	175-195

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

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

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

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

雨水混合比¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
雨水混合比	qr	RWMR	rwmr	(0, 1, 24)	pl	21层	196-216

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

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

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

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

冰水混合比¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
冰水混合比	qi	ICMR	icmr	(0, 1, 23)	pl	21层	217-237

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

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

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

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

雪水混合比¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
雪水混合比	qs	SNMR	snmr	(0, 1, 25)	pl	30层	238-258

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

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

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

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

霰¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
霰	graupel	GRLE	grle	(0, 1, 32)	pl	21层	259-279

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

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

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

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