Merge pull request #761 from pvjeetze/f_MitiConsv

Start script for mitigation and conservation paper

Author: pvjeetze

CHANGELOG.md

@@ -10,6 +10,7 @@ The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/).
 - **scripts** peatland rewetting now automatically considered in `extra/runSEALSallocation.R`
 
 ### added
+- **scripts** added start script for land-based mitigation and biodiversity conservation paper
 - **start_scripts** added `lock_timeout` as option to `start_run` function
 
 ### removed

config/projects/scenario_config_miti_consv.csv

@@ -0,0 +1,2 @@
+;MitiConsv;MitiConsv-NoBuff
+gms$s56_buffer_aff;0.2;0

scripts/start/projects/paper_MitiConsv.R

@@ -0,0 +1,147 @@
+# |  (C) 2008-2024 Potsdam Institute for Climate Impact Research (PIK)
+# |  authors, and contributors see CITATION.cff file. This file is part
+# |  of MAgPIE and licensed under AGPL-3.0-or-later. Under Section 7 of
+# |  AGPL-3.0, you are granted additional permissions described in the
+# |  MAgPIE License Exception, version 1.0 (see LICENSE file).
+# |  Contact: magpie@pik-potsdam.de
+
+# -------------------------------------------------------------
+# description: Land-based mitigation and habitat conservation
+# -------------------------------------------------------------
+
+rev <- "rev9"
+
+cres <- "c200"
+
+######################################
+#### Script to start a MAgPIE run ####
+######################################
+
+library(gms)
+library(lucode2)
+library(magclass)
+library(gdx2)
+
+# Load start functions
+source("scripts/start_functions.R")
+
+# ====================
+# Calibration
+# ====================
+
+source("config/default.cfg")
+
+cfg$title <- "calib_run"
+cfg$output <- c("rds_report", "validation_short")
+cfg$force_replace <- TRUE
+
+cfg$recalibrate <- TRUE
+cfg$recalibrate_landconversion_cost <- TRUE
+
+cfg$best_calib <- TRUE
+cfg$best_calib_landconversion_cost <- FALSE
+
+cfg$gms$s14_use_yield_calib <- 1
+
+# cc is new default
+cfg <- setScenario(cfg, c("nocc_hist", "NPI", "ForestryExo"))
+cfg <- setScenario(cfg, c("MitiConsv"), scenario_config = "config/projects/scenario_config_miti_consv.csv")
+
+# sticky
+cfg$gms$factor_costs <- "sticky_feb18"
+
+# SNV habitat defintion
+cfg$gms$land_snv <- "secdforest, other"
+
+# marginal land scenario
+cfg$gms$c29_marginal_land <- "q33_marginal"
+
+start_run(cfg = cfg)
+calib_tgz <- magpie4::submitCalibration(paste(rev, "MitiConsv", sep = "_"))
+
+# ====================
+# Scenario runs
+# ====================
+
+prefix <- paste(rev, "MitiConsv", cres, sep = "_")
+
+scenarios <- c(
+  "SSP2-REF", "SSP2-PB650-PriceAR", "SSP2-PB650-PriceProt",
+  "SSP2-PB1000-PriceAR", "SSP2-PB1000-PriceProt",
+  "SSP2-PB650-PriceAR-30by30", "SSP2-PB650-PriceProt-30by30",
+  "SSP2-PB650-PriceAR-BH", "SSP2-PB650-PriceProt-BH",
+  "SSP2-PB1000-PriceAR-30by30", "SSP2-PB1000-PriceProt-30by30",
+  "SSP2-PB1000-PriceAR-BH", "SSP2-PB1000-PriceProt-BH"
+)
+
+for (scen in scenarios) {
+  scen <- unlist(strsplit(scen, "-"))
+  ssp <- scen[grepl("SSP", scen)]
+
+  if (length(ssp) == 0) {
+    ssp <- "SSP2"
+  }
+
+  source("config/default.cfg")
+
+  cfg$qos <- "short_highMem"
+
+  cfg$results_folder <- "output/:title::date:"
+
+  cfg$output <- c(
+    "output_check", "extra/disaggregation", "rds_report", "extra/runSEALSallocation"
+  )
+
+  cfg$input["calibration"] <- calib_tgz
+
+  # Climate change switched off for these runs
+  cfg <- setScenario(cfg, c(ssp, "nocc_hist", "NPI", "ForestryExo"))
+  cfg <- setScenario(cfg, c("MitiConsv"), scenario_config = "config/projects/scenario_config_miti_consv.csv")
+
+  # Set path to coupled output
+  pathToCoupledOutput <- "/p/projects/magpie/users/vjeetze/magpie/projects/MitiConsv/C_MitiConsv_Dec24/remind/output/C_rev5_MitiConsv_SSP2-NPi2025-rem-12/REMIND_generic_C_rev5_MitiConsv_SSP2-NPi2025-rem-12.mif"
+
+  # sticky
+  cfg$gms$factor_costs <- "sticky_feb18"
+
+  # SNV habitat defintion
+  cfg$gms$land_snv <- "secdforest, other"
+
+  if ("PriceAR" %in% scen) {
+    cfg$gms$c56_emis_policy <- "redd+natveg_nosoil"
+  }
+
+  if ("PriceProt" %in% scen) {
+    cfg$gms$c56_emis_policy <- "redd+natveg_nosoil"
+    cfg$gms$s56_c_price_induced_AR <- 0
+  }
+
+  if ("PB650" %in% scen) {
+    cfg <- setScenario(cfg, "NDC")
+    # Update path to coupled output
+    pathToCoupledOutput <- "/p/projects/magpie/users/vjeetze/magpie/projects/MitiConsv/C_MitiConsv_Dec24/remind/output/C_rev5_MitiConsv_SSP2-PkBudg650-rem-12/REMIND_generic_C_rev5_MitiConsv_SSP2-PkBudg650-rem-12.mif"
+  }
+
+  if ("PB1000" %in% scen) {
+    cfg <- setScenario(cfg, "NDC")
+    # Update path to coupled output
+    pathToCoupledOutput <- "/p/projects/magpie/users/vjeetze/magpie/projects/MitiConsv/C_MitiConsv_Dec24/remind/output/C_rev5_MitiConsv_SSP2-PkBudg1000-rem-12/REMIND_generic_C_rev5_MitiConsv_SSP2-PkBudg1000-rem-12.mif"
+  }
+
+  if ("30by30" %in% scen) {
+    cfg$gms$c22_protect_scenario <- "30by30"
+  }
+
+  if ("BH" %in% scen) {
+    cfg$gms$c22_protect_scenario <- "BH"
+  }
+
+  # Settings taken from coupled runs
+  cfg$gms$c56_pollutant_prices <- "coupling"
+  cfg$gms$c60_2ndgen_biodem <- "coupling"
+  cfg$path_to_report_ghgprices <- pathToCoupledOutput
+  cfg$path_to_report_bioenergy <- pathToCoupledOutput
+
+  cfg$title <- paste0(prefix, "_", paste(scen, collapse = "-"))
+  start_run(cfg = cfg, codeCheck = FALSE)
+}