MESO要素清单:324-358#
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 |
层次类型 |
层次 |
编号 |
|---|---|---|---|---|---|---|---|
散度 |
RELD |
d |
div |
(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
温度平流#
要素名 |
wgrib2 |
eccodes |
cemc |
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
露点温度#
要素名 |
wgrib2 |
eccodes |
cemc |
GRIB Key |
层次类型 |
层次 |
编号 |
|---|---|---|---|---|---|---|---|
露点温度 |
DPT |
dpt |
td |
(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
温度露点差#
要素名 |
wgrib2 |
eccodes |
cemc |
GRIB Key |
层次类型 |
层次 |
编号 |
|---|---|---|---|---|---|---|---|
温度露点差 |
DEPR |
- |
ttd |
(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
假相当位温#
要素名 |
wgrib2 |
eccodes |
cemc |
GRIB Key |
层次类型 |
层次 |
编号 |
|---|---|---|---|---|---|---|---|
假相当位温 |
EPOT |
papt |
thetaSe |
(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