GRIB 懒加载后 diff 报 step 长度不匹配，先检查切片回归

同一份含两个预报时效的 GRIB 文件，直接 open_dataset(...).diff("step") 报 operands cannot be broadcast together with mismatched lengths for dimension 'step': (1, 0)，先 .load() 再差分却正常。若发生在 xarray 从 2025.3.1 更新到 2025.10.1 后，应先检查惰性切片路径。

这条问题涉及 chunks=None 下的 Lazy Indexing 与 cfgrib 后端。完整切片的快捷路径没有规范化后续切片，后端收到省略边界的切片后可能给出错误长度，导致两段相邻样本无法相减。不是所有广播错误都属于这里：先确认加载前后差异与失败维度。

采用包含切片规范化修正的正式版本

xarray 2025.12.0 已修正 LazilyIndexedArray 的切片规范化，后续正式 2026.9.0 同样保留此处理。按项目与 cfgrib 的兼容约束更新后，应保留原来的读取参数和惰性路径验证，不能先把数据全加载再据此判断回归已消失。

首先检查 step 是否确实存在、长度是否符合预期，以及其他维度是否正常。GRIB 解码器可能因变量组合和过滤参数选择不同的数据集；误读文件或少了时效，不能靠版本更新补出缺失样本。

小范围数据可在关闭文件前加载

暂时不能更新时，数据量可控的任务可先选择需要的变量和范围，再加载后差分。下面是针对已含 tp、step 及空间维度文件的结构示例：

import xarray as xr

with xr.open_dataset("forecast.grib2", engine="cfgrib") as ds:
    small = ds["tp"].isel(
        step=slice(0, 2), latitude=slice(0, 20), longitude=slice(0, 20)
    ).load()
result = small.diff("step")

文件名、变量及空间范围需按实际数据替换。所选数据在文件关闭前已加载，之后差分不再依赖原文件的惰性索引。不要对整个多变量大文件照抄 .load()；它可能带来远超原任务的内存占用。若连这一步选择都出现相同后端索引异常，应采用正式修正，不继续拼接切片绕过。

另一种有条件的路径是在读取时使用 chunks="auto"，由可用的 Dask 分块后端处理。它会改变任务图、调度和分块成本，返回对象也不等于已经读完数据；生命周期应覆盖后续真正计算。

比较差分值与坐标，不只比较能否打开

两个时效的一阶差分，差分轴预期剩一个位置；默认 upper 标签保留后一时效的坐标。用可控范围比较加载结果与惰性结果，核对值、dtype、缺失值及其余维度。

若变量表示累计降水，差分只计算相邻时效之差，不会自动处理累计重置或预报周期切换。数据含义仍须按原任务校验。这里的切片修正也不解决高阶 diff 的递归深度问题，两个异常应分别定位。