如何处理R中的netCDF数据，以便从月度数据中找到年平均温度？

我正在尝试处理一些气候数据。我有一个关于纬度、经度、时间和地表温度的.nc文件，从01-850到12-1049，每个月都有这些数据。我需要将这些月数据转换为年数据，并找到每个纬度，经度的年平均温度

ncin <- nc_open("tas_Amon_MIROC-ES2L_past1000_r1i1p1f2_gn_085001-104912.nc")

lon <- ncvar_get(ncin,"lon")
lon <- ifelse(lon > 180, -(360 - lon), lon)

lat <- ncvar_get(ncin, "lat")

time <- format(as.Date(time, origin = ("850-1-1"), tz = "UTC"),"%Y-%m" )

tas.array <- ncvar_get(ncin,"tas")
tas.array <- tas.array - 273.15

ncin