如何在R中用一维条件子集NetCDF文件

我有一个从1979年到2019年的3维（lon、lat、time）月温度NetCDF文件。我从以下网站下载了该文件：。我想创建一个函数，允许我选择纬度高于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"

---

到目前为止，我只是设法对每年的月份进行子集划分，但现在我需要合并上面解释的条件。非常感谢你

请提供一个带有dput（头部（df，10））的数据帧示例。您的意思是什么？像一小块数据<代码>c（0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831，0.00965203495405831），查看本课程

if (lat >23) { "set all values in lake_temp to NA except in the months December, January, February"}