将nc文件每月日数据平均为月平均并输出

# -*- coding: utf-8 -*-
"""
Created on Fri Dec  9 11:06:26 2022

@author: ASUS
"""
import netCDF4 as nc
import xarray as xr
from netCDF4 import Dataset
import numpy as np

f=xr.open_dataset('H:\data\CN05data\CN05.1_Tm_1961_2018_daily_025x025.nc') #打开nc文件
tm=f['tm'].loc['1979-01-01':'2018-12-31'] 
lon =tm.lon
lat =tm.lat
tm=tm.loc[tm.time.dt.month.isin([8])]
tm=np.array(tm).reshape(40,31,163,283)
tm=np.nanmean(tm,axis=1)
year = np.arange(1979,2019)
data2_da = xr.DataArray(tm,coords=[year,lat,lon],dims=['time','lat','lon'])
ds = xr.Dataset({'tm':data2_da})
ds.to_netcdf('H:/CN05.1_Tm_1979_2018_monthly_August_025x025.nc')