######################################################
# CALCULATE COST DIFFERENTIAL FOR TWO SETS OF INPUTS #
##########################################################################################

# Clear workspace and load necessary packages.
pkg <- c("tidyverse", "magrittr", "simputation", "readxl", "zoo", "tidymodels",
         "collapse")
sapply(pkg, library, char = TRUE)

# Suppress column type guessing messages by 'readr' functions.
options(readr.show_col_types = FALSE)

# Disable 'collapse' namespace to avoid code collisions.
set_collapse(mask = NULL)

# Load base parameters.
source("scripts/base_par.R")
source("scripts/aux_fun.R")

par1 <- par2 <- par
rm(par)

par2$ahv13 <- FALSE

(comp <- 
  compare_projections(par1, par2) %>% select(year, cost_diff))

