如何阻止xarray.Dataset.to_netcdf()强制将time维度置于首位？

xarray 保存 netCDF 文件时，默认会把无限维度（unlimited dimension）放到所有维度的最前面，这是它内部的固定处理逻辑，和原文件的维度顺序没关系。ncdump 是直接读取文件元数据展示维度顺序，所以会看到 time 被移到首位；但 xarray 加载文件时会保留自身 Dataset 对象里的维度顺序，因此读回去后显示正常。

另外你提到原 classic 格式文件的无限维度在末尾，这大概率是旧版本 netCDF 工具创建的——虽然 classic 格式规范建议无限维度作为第一个维度，但有些旧工具没严格执行这个限制，而 xarray 保存时会按更严格的逻辑处理，强制把无限维度前置，这也是顺序变化的原因之一。

不用 ncpdq 或者 transpose，直接给 to_netcdf 传 encoding 参数，强制指定每个变量的维度顺序即可：

# 获取原 Dataset 里的正确维度顺序
dim_order = list(data.dims)

# 给每个变量设置编码，指定维度顺序和原 Dataset 一致
encoding = {var: {"dimensions": dim_order} for var in data.data_vars}

# 保存文件
data.to_netcdf("data_copy.nc", encoding=encoding)

这样保存后的文件用 ncdump -h 查看，维度顺序就和原文件完全一致了。

• format 参数只控制文件存储格式，无法改变 xarray 调整维度顺序的逻辑；
• unlimited_dims='time' 只是明确哪个维度为无限维度，不会改变 xarray 把无限维度前置的默认行为；
• engine='netcdf4' 只是切换底层读写引擎，同样不影响维度顺序的处理。

内容的提问来源于stack exchange，提问作者angeloff