如何用NCO/CDO自动批量处理多年NetCDF文件计算月平均？

一、使用CDO（推荐，自动识别时间维度） #

CDO会自动解析文件中的时间戳，无需手动计算时间步长，直接按自然月份分组求平均，操作更省心：

1. 合并所有文件（可选） #

如果你的文件按时间顺序命名（如data_2000_01-03.nc、data_2000_04-06.nc），可以先合并成单个文件统一处理：

cdo mergetime data_*.nc all_data.nc

2. 计算全局月平均 #

用monmean算子直接生成所有月份的平均文件：

cdo monmean all_data.nc monthly_mean.nc

3. 单独处理每个3月文件 #

如果不想合并，可逐个处理每个3月数据文件，再合并结果：

for file in data_*.nc; do
    cdo monmean "$file" "monthly_${file}"
done
# 合并所有单文件生成的月平均结果
cdo mergetime monthly_data_*.nc all_monthly_mean.nc

二、使用NCO #

NCO通过时间属性筛选数据，避开时间步长不一致的问题，两种常用方法：

1. 合并文件后批量计算 #

先合并所有文件：

ncrcat data_*.nc all_data.nc

方法A：循环遍历月份计算

for mon in {1..12}; do
    # 筛选所有年份的第mon个月数据并求平均
    ncra -d time,"month=${mon}" all_data.nc monthly_mean_${mon}.nc
done
# 合并所有月份结果
ncrcat monthly_mean_*.nc all_monthly_mean.nc

方法B：按月份分组批量计算

如果文件的time维度带有month属性（多数气候数据都有），用ncwa直接按月份分组求平均：

ncwa -a time --group=month all_data.nc monthly_mean.nc

2. 处理单个3月文件 #

对每个包含3个月数据的文件，直接筛选对应月份计算：

for file in data_*.nc; do
    # 先通过ncdump -h "$file" | grep month查看文件包含的月份
    # 示例：假设文件含1、2、3月
    ncra -d time,"month=1" "$file" "monthly_1_${file}"
    ncra -d time,"month=2" "$file" "monthly_2_${file}"
    ncra -d time,"month=3" "$file" "monthly_3_${file}"
done
# 合并所有结果
ncrcat monthly_*.nc all_monthly_mean.nc

注意事项 #

• 确保所有NetCDF文件的时间维度属性规范（如units为days since YYYY-MM-DD或标准ISO时间），否则工具无法正确识别时间。
• 若时间格式异常，可先用cdo settime（CDO）或ncatted（NCO）修正时间属性。

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