使用ncra --mro计算季节平均仅返回首条记录值的问题

1. 你当前使用的NCO 4.6.3为2016年发布的旧版本，存在已知的--mro参数配合手动设置的记录维度处理bug，这是你将time设为记录维度后输出文件损坏、无法读取的核心原因。
2. 原始文件中time维度默认不是记录类型（unlimited维度），--mro的子周期平均功能仅对记录维度生效，你执行第一条命令时收到的警告已经明确提示了该限制，因此ncra没有执行3个月平均操作，仅提取了每个子周期的第一个元素，导致结果和抽取每年9月的输出完全一致。

按实现成本从低到高推荐以下三种方案：

• 方案一：升级NCO到4.9.0及以上的稳定版本。升级完成后依次执行以下命令即可得到正确结果：

  # 将time设置为记录维度
  ncks --mk_rec_dmn time original.nc4 record.nc4
  # 按每3个时间步为周期计算平均
  ncra --mro -d time,,,3,3 record.nc4 record_average.nc4
  

• 方案二：使用CDO工具完成计算，无需转换维度，操作更简单：

  # 按连续3个时间步为单位计算平均
  cdo timselmean,3 original.nc4 average.nc4
  

• 方案三：如果无法升级NCO也不想安装新工具，可通过拆分-平均-拼接的手动流程实现：

  # 假设总共有n年数据，将n替换为实际年份数
  for i in {0..n-1}
  do
    # 提取对应年份的三个月数据
    ncks -d time,$((i*3)),$((i*3+2)) original.nc4 year_${i}.nc4
    # 计算单年季度平均
    ncra year_${i}.nc4 year_${i}_avg.nc4
  done
  # 拼接所有年份的平均结果
  ncrcat year_*_avg.nc4 average.nc4
  # 清理临时文件
  rm year_*.nc4
  

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