recipe_file <- "recipes/examples/NAO_recipe.yml"

source("modules/Loading/Loading.R")
source("modules/Calibration/Calibration.R")
source("modules/Skill/Skill.R")
source("modules/Saving/Saving.R")
source("tools/prepare_outputs.R")
source("modules/Anomalies/Anomalies.R")
recipe <- prepare_outputs(recipe_file)

#for (smonth in 1:12) {
  data <- load_datasets(recipe)
  gc()
  data <- compute_anomalies(recipe, data)
#  data <- readRDS("../Test_NAOmodule.RDS")
source("modules/Indices/Indices.R")
  nao_s2dv <- Indices(data = data, recipe = recipe)
source("modules/Skill/Skill.R")
  # todo parameter agg to get it from the nao_s2dv?
  skill_metrics <- compute_skill_metrics(recipe = recipe, data = nao_s2dv,
                                         agg = 'region')