GFS要素清单：454-471 - 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

gfs_grib2_orig_file_path = get_data_file("cma-gfs")
gfs_grib2_orig_file_path

The history saving thread hit an unexpected error (OperationalError('attempt to write a readonly database')).History will not be written to the database.

PosixPath('/home/wangdp/project/cedarkit/notebook-project/notebook-devel/repo/data-notebook-project/data-notebook/data/cma-gfs/grib2/orig/Z_NAFP_C_BABJ_20260811000000_P_NWPC-GRAPES-GFS-GLB-02400.grib2')

土壤温度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
土壤温度	st(0-10)	TMP	t	(0, 0, 0)	depthBelowLandLayer	0, 0.1	464
土壤温度	st(10-40)	TMP	t	(0, 0, 0)	depthBelowLandLayer	0.1, 0.4	465
土壤温度	st(40-100)	TMP	t	(0, 0, 0)	depthBelowLandLayer	0.4, 1	466
土壤温度	st(100-200)	TMP	t	(0, 0, 0)	depthBelowLandLayer	1, 2	467

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="TMP",
    level_type="depthBelowLandLayer",
    level={
        "first_level": 0,
        "second_level": 0.1,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 464

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="t",
    level_type="depthBelowLandLayer",
    level={
        "first_level": 0.1,
        "second_level": 0.4,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 465

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="st(40-100)",
    lazy=True
)
assert field.attrs["GRIB_count"] == 466

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 0,
        "parameterNumber": 0,
    },
    level_type="depthBelowLandLayer",
    level={
        "first_level": 1,
        "second_level": 2,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 467

土壤湿度¶

要素名	cemc	wgrib2	eccodes	GRIB Key	层次类型	层次	编号
土壤湿度	sw(0-10)	SPFH	q	(0, 1, 0)	depthBelowLandLayer	0, 0.1	468
土壤湿度	sw(10-40)	SPFH	q	(0, 1, 0)	depthBelowLandLayer	0.1, 0.4	469
土壤湿度	sw(40-100)	SPFH	q	(0, 1, 0)	depthBelowLandLayer	0.4, 1	470
土壤湿度	sw(100-200)	SPFH	q	(0, 1, 0)	depthBelowLandLayer	1, 2	471

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="SPFH",
    level_type="depthBelowLandLayer",
    level={
        "first_level": 0,
        "second_level": 0.1,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 468

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="q",
    level_type="depthBelowLandLayer",
    level={
        "first_level": 0.1,
        "second_level": 0.4,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 469

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter="sw(40-100)",
    lazy=True
)
assert field.attrs["GRIB_count"] == 470

field = load_field_from_file(
    gfs_grib2_orig_file_path,
    parameter={
        "discipline": 0,
        "parameterCategory": 1,
        "parameterNumber": 0,
    },
    level_type="depthBelowLandLayer",
    level={
        "first_level": 1,
        "second_level": 2,
    },
    lazy=True
)
assert field.attrs["GRIB_count"] == 471