1979 年から 2019 年までの 3 つの次元 (経度、緯度、時間) の範囲の月別気温の NetCDF ファイルがあります。 /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"}
これまでのところ、毎年これらの月をサブセット化することができましたが、上記で説明した条件を組み込む必要があります。どうもありがとう！