rscript_get_data_sif_pk_lores.R

library(dplyr)
library(tidyr)
library(purrr)
library(lubridate)
library(magrittr)
library(tidync)

##------------------------------------------------------------------------
## Extract point data and construct separate nested time series data frame
## for each longitde slice
##------------------------------------------------------------------------
dir <- "~/data/gome_2_sif_downscaled/data_halfdeg/"
fileprefix <- "GOME_PK_dcSIF_05deg_8day_"
nclist <- paste0(dir, list.files(dir, pattern = paste0(fileprefix, ".*.nc"), recursive = TRUE))
outdir <- "~/data/gome_2_sif_downscaled/data_tidy_halfdeg/"
varnam <- "SIF"
lonnam <- "lon"
latnam <- "lat"
timenam <- "time"
timedimnam <- "time"

## create files for each longitude slice, containing full time series wrapped for each gridcell (latitude)
rbeni::nclist_to_df(nclist, outdir, fileprefix, varnam, lonnam, latnam, timenam, timedimnam, ncores = 32, single_basedate = TRUE) # 80 cores