data-raw/make_extdata/make_dataset_list_for_wc_future.R

In EvolEcolGroup/pastclim: Manipulate Time Series of Palaeoclimate Reconstructions

# start this script from this directory
wc_template <- read.csv("../data_files/worldclim_future_template.csv")
gcm <- c(
  "ACCESS-CM2", "BCC-CSM2-MR", "CMCC-ESM2", "EC-Earth3-Veg", "FIO-ESM-2-0",
  "GFDL-ESM4", "GISS-E2-1-G", "HadGEM3-GC31-LL", "INM-CM5-0", "IPSL-CM6A-LR",
  "MIROC6", "MPI-ESM1-2-HR", "MRI-ESM2-0", "UKESM1-0-LL"
)
scenarios <- c("ssp126", "ssp245", "ssp370", "ssp585")
dates_df <- data.frame(
  orig = c("2021-2040", "2041-2060", "2061-2080", "2081-2100"),
  time_bp = (c(2030, 2050, 2070, 2090) - 1950)
)
resolutions <- c(10, 5) # add 2.5 and 0.5
library(dplyr)
library(stringr)
wc_future <- wc_template[0, ]
for (i_gcm in gcm) {
  for (i_scenario in scenarios) {
    for (i_res in resolutions) {
      wc_this <- wc_template %>%
        mutate(
          dataset = str_replace(dataset, "XXX", i_gcm),
          file_name = str_replace(file_name, "XXX", i_gcm)
        ) %>%
        mutate(
          dataset = str_replace(dataset, "YYY", i_scenario),
          file_name = str_replace(file_name, "YYY", i_scenario)
        ) %>%
        mutate(
          dataset = str_replace(dataset, "ZZZ", as.character(i_res)),
          file_name = str_replace(file_name, "ZZZ", as.character(i_res))
        )
      wc_future <- wc_future %>% bind_rows(wc_this)
    }
  }
}
write.csv(wc_future, "../data_files/worldclim_future_list.csv",
  row.names = FALSE, na = ""
)