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
meso3km_grib2_orig_file_path = get_data_file("cma-meso-3km")
meso3km_grib2_orig_file_pathPosixPath('/home/wangdp/project/cedarkit/notebook-project/notebook-devel/repo/data-notebook-project/data-notebook/data/cma-meso-3km/grib2/orig/Z_NAFP_C_BABJ_20260818000000_P_NWPC-GRAPES-3KM-ORIG-02400.grb2')降水¶
总累计降水¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 总累计降水 | rain | APCP | - | (0, 1, 8) | surface | - | 1 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="APCP",
lazy=True
)
assert field.attrs["GRIB_count"] == 1field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="rain",
lazy=True
)
assert field.attrs["GRIB_count"] == 1field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 8},
lazy=True
)
assert field.attrs["GRIB_count"] == 1降雪¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 降雪 | sf | ASNOW | - | (0, 1, 29) | surface | - | 2 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="ASNOW",
lazy=True
)
assert field.attrs["GRIB_count"] == 2field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="sf",
lazy=True
)
assert field.attrs["GRIB_count"] == 2field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 29},
lazy=True
)
assert field.attrs["GRIB_count"] == 2地表温度与辐射通量¶
地表温度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 地表温度 | ts | TMP | t | (0, 0, 0) | surface | - | 3 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TMP",
lazy=True
)
assert field.attrs["GRIB_count"] == 3field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="t",
lazy=True
)
assert field.attrs["GRIB_count"] == 3field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="ts",
lazy=True
)
assert field.attrs["GRIB_count"] == 3field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 0, "parameterNumber": 0},
lazy=True
)
assert field.attrs["GRIB_count"] == 3净长波辐射通量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 净长波辐射通量 | glw | NLWRF | str | (0, 5, 5) | surface | - | 4 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="NLWRF",
lazy=True
)
assert field.attrs["GRIB_count"] == 4field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="str",
lazy=True
)
assert field.attrs["GRIB_count"] == 4field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="glw",
lazy=True
)
assert field.attrs["GRIB_count"] == 4field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 5, "parameterNumber": 5},
lazy=True
)
assert field.attrs["GRIB_count"] == 4净短波辐射通量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 净短波辐射通量 | gsw | NSWRF | ssr | (0, 4, 9) | surface | - | 5 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="NSWRF",
lazy=True
)
assert field.attrs["GRIB_count"] == 5field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="ssr",
lazy=True
)
assert field.attrs["GRIB_count"] == 5field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="gsw",
lazy=True
)
assert field.attrs["GRIB_count"] == 5field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 4, "parameterNumber": 9},
lazy=True
)
assert field.attrs["GRIB_count"] == 5感热通量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 感热通量 | sshf | HFLUX | hflux | (2, 0, 24) | surface | - | 6 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="HFLUX",
lazy=True
)
assert field.attrs["GRIB_count"] == 6field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="hflux",
lazy=True
)
assert field.attrs["GRIB_count"] == 6field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="sshf",
lazy=True
)
assert field.attrs["GRIB_count"] == 6field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 2, "parameterCategory": 0, "parameterNumber": 24},
lazy=True
)
assert field.attrs["GRIB_count"] == 6地表水汽通量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 地表水汽通量 | qfxsfc | - | - | (0, 1, 229) | surface | - | 7 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="qfxsfc",
lazy=True
)
assert field.attrs["GRIB_count"] == 7field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 229},
lazy=True
)
assert field.attrs["GRIB_count"] == 7地形¶
地形高度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 地形高度 | zs | HGT | orog | (0, 3, 5) | surface | - | 8 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="HGT",
lazy=True
)
assert field.attrs["GRIB_count"] == 8field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="orog",
lazy=True
)
assert field.attrs["GRIB_count"] == 8field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="zs",
lazy=True
)
assert field.attrs["GRIB_count"] == 8field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 3, "parameterNumber": 5},
lazy=True
)
assert field.attrs["GRIB_count"] == 8近地面要素¶
2米比湿¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 2米比湿 | q2m | SPFH | 2sh | (0, 1, 0) | heightAboveGround | 2 | 9 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="SPFH",
lazy=True
)
assert field.attrs["GRIB_count"] == 9field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="2sh",
lazy=True
)
assert field.attrs["GRIB_count"] == 9field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="q2m",
lazy=True
)
assert field.attrs["GRIB_count"] == 9field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 0},
lazy=True
)
assert field.attrs["GRIB_count"] == 92米温度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 2米温度 | t2m | TMP | 2t | (0, 0, 0) | heightAboveGround | 2 | 10 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TMP",
level_type="heightAboveGround",
level=2,
lazy=True
)
assert field.attrs["GRIB_count"] == 10field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="2t",
lazy=True
)
assert field.attrs["GRIB_count"] == 10field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="t2m",
lazy=True
)
assert field.attrs["GRIB_count"] == 10field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 0, "parameterNumber": 0},
level_type="heightAboveGround",
level=2,
lazy=True
)
assert field.attrs["GRIB_count"] == 1010米U风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 10米U风 | u10m | UGRD | 10u | (0, 2, 2) | heightAboveGround | 10 | 11 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="UGRD",
lazy=True
)
assert field.attrs["GRIB_count"] == 11field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="10u",
lazy=True
)
assert field.attrs["GRIB_count"] == 11field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="u10m",
lazy=True
)
assert field.attrs["GRIB_count"] == 11field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 2},
lazy=True
)
assert field.attrs["GRIB_count"] == 1110米V风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 10米V风 | v10m | VGRD | 10v | (0, 2, 3) | heightAboveGround | 10 | 12 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="VGRD",
lazy=True
)
assert field.attrs["GRIB_count"] == 12field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="10v",
lazy=True
)
assert field.attrs["GRIB_count"] == 12field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="v10m",
lazy=True
)
assert field.attrs["GRIB_count"] == 12field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 3},
lazy=True
)
assert field.attrs["GRIB_count"] == 12云量¶
总云量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 总云量 | tcc | TCDC | tcc | (0, 6, 1) | entireAtmosphere | - | 13 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TCDC",
lazy=True
)
assert field.attrs["GRIB_count"] == 13field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="tcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 13field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="tcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 13field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 6, "parameterNumber": 1},
lazy=True
)
assert field.attrs["GRIB_count"] == 13低云量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 低云量 | lcc | LCDC | lcc | (0, 6, 3) | entireAtmosphere | - | 14 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="LCDC",
lazy=True
)
assert field.attrs["GRIB_count"] == 14field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="lcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 14field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="lcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 14field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 6, "parameterNumber": 3},
lazy=True
)
assert field.attrs["GRIB_count"] == 14中云量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 中云量 | mcc | MCDC | mcc | (0, 6, 4) | entireAtmosphere | - | 15 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="MCDC",
lazy=True
)
assert field.attrs["GRIB_count"] == 15field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="mcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 15field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="mcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 15field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 6, "parameterNumber": 4},
lazy=True
)
assert field.attrs["GRIB_count"] == 15高云量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 高云量 | hcc | HCDC | hcc | (0, 6, 5) | entireAtmosphere | - | 16 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="HCDC",
lazy=True
)
assert field.attrs["GRIB_count"] == 16field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="hcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 16field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="hcc",
lazy=True
)
assert field.attrs["GRIB_count"] == 16field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 6, "parameterNumber": 5},
lazy=True
)
assert field.attrs["GRIB_count"] == 16整层水量¶
总柱状水汽¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 总柱状水汽 | tvw | TCIWV | - | (0, 1, 64) | atmosphere | - | 17 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TCIWV",
lazy=True
)
assert field.attrs["GRIB_count"] == 17field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="tvw",
lazy=True
)
assert field.attrs["GRIB_count"] == 17field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 64},
lazy=True
)
assert field.attrs["GRIB_count"] == 17总柱状云水¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 总柱状云水 | tcw | TCOLW | - | (0, 1, 69) | atmosphere | - | 18 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TCOLW",
lazy=True
)
assert field.attrs["GRIB_count"] == 18field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="tcw",
lazy=True
)
assert field.attrs["GRIB_count"] == 18field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 69},
lazy=True
)
assert field.attrs["GRIB_count"] == 18总柱状云冰¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 总柱状云冰 | tiw | TCOLI | - | (0, 1, 70) | atmosphere | - | 19 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="TCOLI",
lazy=True
)
assert field.attrs["GRIB_count"] == 19field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="tiw",
lazy=True
)
assert field.attrs["GRIB_count"] == 19field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 70},
lazy=True
)
assert field.attrs["GRIB_count"] == 19边界层与雷达¶
边界层高度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 边界层高度 | pblh | HPBL | blh | (0, 3, 18) | surface | - | 20 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="HPBL",
lazy=True
)
assert field.attrs["GRIB_count"] == 20field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="blh",
lazy=True
)
assert field.attrs["GRIB_count"] == 20field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="pblh",
lazy=True
)
assert field.attrs["GRIB_count"] == 20field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 3, "parameterNumber": 18},
lazy=True
)
assert field.attrs["GRIB_count"] == 20雷达组合反射率¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 雷达组合反射率 | cdbz | CR | - | (0, 16, 224) | surface | - | 21 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="CR",
lazy=True
)
assert field.attrs["GRIB_count"] == 21field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="cdbz",
lazy=True
)
assert field.attrs["GRIB_count"] == 21field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 16, "parameterNumber": 224},
lazy=True
)
assert field.attrs["GRIB_count"] == 21最大统计与强对流诊断¶
最大垂直速度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 最大垂直速度 | wmax#1 | DZDT | wz | (0, 2, 9) | pl | 1000 | 22 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="DZDT",
lazy=True
)
assert field.attrs["GRIB_count"] == 22field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="wz",
lazy=True
)
assert field.attrs["GRIB_count"] == 22field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="wmax#1",
lazy=True
)
assert field.attrs["GRIB_count"] == 22field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 9},
lazy=True
)
assert field.attrs["GRIB_count"] == 22最大10米U风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 最大10米U风 | u10mmax#1 | UGRD | 10u | (0, 2, 2) | heightAboveGround | 10 | 23 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="u10mmax#1",
lazy=True
)
assert field.attrs["GRIB_count"] == 23最大10米V风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 最大10米V风 | v10mmax#1 | VGRD | 10v | (0, 2, 3) | heightAboveGround | 10 | 24 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="v10mmax#1",
lazy=True
)
assert field.attrs["GRIB_count"] == 24最大组合雷达反射率¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 最大组合雷达反射率 | cdbzmax#1 | CR | - | (0, 16, 224) | surface | - | 25 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="cdbzmax#1",
lazy=True
)
assert field.attrs["GRIB_count"] == 25雷达回波顶¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 雷达回波顶 | echotop | RETOP | - | (0, 16, 3) | surface | - | 26 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="RETOP",
lazy=True
)
assert field.attrs["GRIB_count"] == 26field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="echotop",
lazy=True
)
assert field.attrs["GRIB_count"] == 26field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 16, "parameterNumber": 3},
lazy=True
)
assert field.attrs["GRIB_count"] == 262000-5000米最大上升螺旋度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 2000-5000米最大上升螺旋度 | uhmax(2000-5000) | UPHL | - | (0, 7, 15) | heightAboveGroundLayer | 5000, 2000 | 27 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="UPHL",
lazy=True
)
assert field.attrs["GRIB_count"] == 27field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="uhmax(2000-5000)",
lazy=True
)
assert field.attrs["GRIB_count"] == 27field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 7, "parameterNumber": 15},
lazy=True
)
assert field.attrs["GRIB_count"] == 27近地面总体理查森数¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 近地面总体理查森数 | risfc | RI | - | (0, 7, 12) | surface | - | 28 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="RI",
lazy=True
)
assert field.attrs["GRIB_count"] == 28field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="risfc",
lazy=True
)
assert field.attrs["GRIB_count"] == 28field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 7, "parameterNumber": 12},
lazy=True
)
assert field.attrs["GRIB_count"] == 28边界层总体理查森数¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 边界层总体理查森数 | ripbl | RI | - | (0, 7, 12) | mixingLayer | - | 29 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="RI",
level_type="mixingLayer",
level=0,
lazy=True
)
assert field.attrs["GRIB_count"] == 29field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="ripbl",
level_type="mixingLayer",
level=0,
lazy=True
)
assert field.attrs["GRIB_count"] == 29field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 7, "parameterNumber": 12},
level_type="mixingLayer",
level=0,
lazy=True
)
assert field.attrs["GRIB_count"] == 295分钟最大降水¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 5分钟最大降水 | rain5max | - | - | (0, 1, 241) | surface | - | 30 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="rain5max",
lazy=True
)
assert field.attrs["GRIB_count"] == 30field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 1, "parameterNumber": 241},
lazy=True
)
assert field.attrs["GRIB_count"] == 300-600米最大垂直风切变矢量¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 0-600米最大垂直风切变矢量 | shrmax(0-600) | - | - | (0, 2, 237) | heightAboveGroundLayer | 600, 0 | 31 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="shrmax(0-600)",
lazy=True
)
assert field.attrs["GRIB_count"] == 31field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 237},
lazy=True
)
assert field.attrs["GRIB_count"] == 310-600米最大垂直风切变方向¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 0-600米最大垂直风切变方向 | shrmaxdir(0-600) | - | - | (0, 2, 238) | heightAboveGroundLayer | 600, 0 | 32 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="shrmaxdir(0-600)",
lazy=True
)
assert field.attrs["GRIB_count"] == 32field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 238},
lazy=True
)
assert field.attrs["GRIB_count"] == 32最大雷达回波的高度¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 最大雷达回波的高度 | dbzmaxh | - | - | (0, 16, 226) | surface | - | 33 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="dbzmaxh",
lazy=True
)
assert field.attrs["GRIB_count"] == 33field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 16, "parameterNumber": 226},
lazy=True
)
assert field.attrs["GRIB_count"] == 33气压¶
地面气压¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 地面气压 | ps | PRES | sp | (0, 3, 0) | surface | - | 34 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="PRES",
lazy=True
)
assert field.attrs["GRIB_count"] == 34field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="sp",
lazy=True
)
assert field.attrs["GRIB_count"] == 34field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="ps",
lazy=True
)
assert field.attrs["GRIB_count"] == 34field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 3, "parameterNumber": 0},
lazy=True
)
assert field.attrs["GRIB_count"] == 34海平面气压¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 海平面气压 | psl | PRMSL | prmsl | (0, 3, 1) | meanSea | - | 35 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="PRMSL",
lazy=True
)
assert field.attrs["GRIB_count"] == 35field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="prmsl",
lazy=True
)
assert field.attrs["GRIB_count"] == 35field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="psl",
lazy=True
)
assert field.attrs["GRIB_count"] == 35The history saving thread hit an unexpected error (OperationalError('disk I/O error')).History will not be written to the database.
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 3, "parameterNumber": 1},
lazy=True
)
assert field.attrs["GRIB_count"] == 35高度层风¶
东西风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 东西风 | u | UGRD | u | (0, 2, 2) | heightAboveGround | 9层 | 36-44 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="UGRD",
level_type="heightAboveGround",
level=30,
lazy=True
)
assert field.attrs["GRIB_count"] == 36field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="u",
level_type="heightAboveGround",
level=100,
lazy=True
)
assert field.attrs["GRIB_count"] == 39field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="u",
level_type="heightAboveGround",
level=160,
lazy=True
)
assert field.attrs["GRIB_count"] == 42field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 2},
level_type="heightAboveGround",
level=200,
lazy=True
)
assert field.attrs["GRIB_count"] == 44南北风¶
| 要素名 | cemc | wgrib2 | eccodes | GRIB Key | 层次类型 | 层次 | 编号 |
|---|---|---|---|---|---|---|---|
| 南北风 | v | VGRD | v | (0, 2, 3) | heightAboveGround | 9层 | 45-53 |
field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="VGRD",
level_type="heightAboveGround",
level=30,
lazy=True
)
assert field.attrs["GRIB_count"] == 45field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="v",
level_type="heightAboveGround",
level=100,
lazy=True
)
assert field.attrs["GRIB_count"] == 48field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter="v",
level_type="heightAboveGround",
level=160,
lazy=True
)
assert field.attrs["GRIB_count"] == 51field = load_field_from_file(
meso3km_grib2_orig_file_path,
parameter={"discipline": 0, "parameterCategory": 2, "parameterNumber": 3},
level_type="heightAboveGround",
level=200,
lazy=True
)
assert field.attrs["GRIB_count"] == 53