将多个NetCDF合并到一个文件中，每个文件的平均值都为平均值

Question

我有3个NetCDF文件，其中包含1981年至2010年之间每月平均每日温度范围（dtr）的数据。一个文件的使用期限为10年，因此每个像素具有10 x 12的值（每月1个值x 10年）。可以下载文件here（需要登录）。

我想将所有这些文件合并为一个文件，其中包含3个文件所覆盖的整个期间的月均dtr。最终的想法是每个像素12个值（30年平均1个每月值）。

我以前没有使用NetCDF格式的经验，也不太了解其结构。我尝试了NetCDF包中的ncdump函数，但没有从中学到很多……

不幸的是，执行以下操作只能计算30年的平均dtr。

"GIS/Clim/BioClim/dtr/" %>%
  list.files(pattern = "\\.nc$", full.names = TRUE) %>%
  raster::stack() %>% 
  raster::overlay(fun = mean)

这是从一个文件中的单个像素提取数据时得到的结果：

library(tidyverse)
library(sf)

p <- tibble(lat = 45.76, lon = -107.12) %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326)

raster::stack("GIS/Clim/dtr/cru_dtr-2001_2010.nc") %>% 
  raster::extract(p, df = TRUE)

# X2001.01.16 X2001.02.15 X2001.03.16 X2001.04.16 X2001.05.16 X2001.06.16
# 1        11.1        11.3        10.8    9.900001    9.800000    9.400001
# X2001.07.16 X2001.08.16 X2001.09.16 X2001.10.16 X2001.11.16 X2001.12.16
# 1   10.100000         8.6   10.200000    9.500000         9.5         9.3
# X2002.01.16 X2002.02.15 X2002.03.16 X2002.04.16 X2002.05.16 X2002.06.16
# 1        10.6        10.7        10.3         9.0    9.600000         9.1
# 
# ...                               ...                               ...
# 
# X2010.07.16 X2010.08.16 X2010.09.16 X2010.10.16 X2010.11.16 X2010.12.16
# 1    9.000000         9.0    8.900001    8.800000         8.3         8.3

如何合并这三个NetCDF文件并计算月平均dtr来创建新的NetCDF？

请注意，Python解决方案也适用于我。