如何使用python从netcdf数据集中找到气候月平均值？

import netCDF4 as nc
import numpy as np

# 打开NetCDF文件
file_path = 'path_to_your_netcdf_file.nc'
dataset = nc.Dataset(file_path)

# 假设时间变量名为'time'，温度变量名为'temperature'
time_var = dataset.variables['time']
temperature_var = dataset.variables['temperature']

# 获取时间和温度数据
times = time_var[:]
temperatures = temperature_var[:]

# 假设时间变量是以小时为单位的，计算每个小时对应的月份
# 这里需要根据你的时间变量的具体格式进行调整
months = (times // 8760) % 12 + 1  # 假设一年有8760小时

# 计算每个月的平均温度
monthly_averages = {}
for month in range(1, 13):
    monthly_temps = temperatures[months == month]
    if len(monthly_temps) > 0:
        monthly_averages[month] = np.mean(monthly_temps)
    else:
        monthly_averages[month] = None

# 打印结果
for month, avg_temp in monthly_averages.items():
    print(f"Month {month}: Average Temperature = {avg_temp}")

# 关闭数据集
dataset.close()