我有一个从 1979 年到 2019 年具有 3 个维度（经度、纬度、时间）的月度温度 NetCDF 文件。我从这个网站下载了该文件：https ://cds.climate.copernicus.eu/cdsapp#!/dataset /reanalysis-era5-single-levels-monthly-means?tab=form。我想创建一个函数，允许我选择纬度高于 23 度的所有网格单元，对于这些网格单元，只需使用例如 12 月、1 月、2 月的月份来计算年平均值。任何想法将不胜感激。
到目前为止，我只是设法单独执行任务，但我正在寻找更多类似的东西：
mynetcdf <- 'lake.nc'
ncFile <- nc_open(mynetcdf)
lake_temp <- ncvar_get(ncFile,"p0001")
lon <- ncvar_get(ncFile,"longitude")
lat <- ncvar_get(ncFile,"latitude")
time <- ncvar_get)ncFile,"time"
我正在寻找类似的东西：
if (lat >23) { "set all values in lake_temp to NA except in the months December, January, February"}
到目前为止，我只是设法对每年的这些月份进行了子集化，但现在我需要合并上述条件。太感谢了！