Update Citation

CITATION.cff

@@ -118,6 +118,12 @@ authors:
     affiliation: "Potsdam Institute for Climate Impact Research"
     email: [email protected]
 
+  - family-names: Steinhauser
+    given-names: Jan
+    orcid: https://orcid.org/0000-0002-5989-6855
+    affiliation: "Potsdam Institute for Climate Impact Research"
+    email: [email protected]
+
   - family-names: Lotze-Campen
     given-names: Hermann
     orcid: https://orcid.org/0000-0002-0003-5508

config/default.cfg

@@ -1390,7 +1390,7 @@ cfg$gms$s56_cprice_red_factor <- 1    # def = 1
 # * Note: For best consistency it is recommended to use trajectories from the most recent
 # *     coupled REMIND-MAgPIE runs. Currently, this is R32M46.
 # * Available options:
-# *     none: no 2nd generation bioenergy demand 
+# *     none: no 2nd generation bioenergy demand
 # * R32M46: Coupled REMIND-MAgPIE runs with REMIND 3.2.1.dev412 and MAgPIE 4.6.11
 # *   NPi: Current policies; limit peak warming to 3.0°C
 # *   NDC: Nationally determined contributions; limit peak warming to 2.5°C
@@ -1489,9 +1489,9 @@ cfg$gms$c56_cprice_aff <- "secdforest_vegc"   # def = secdforest_vegc
 # * should not be higher than s32_planing_horizon
 cfg$gms$s56_c_price_exp_aff <- 50   # def = 50
 
-# * Share of carbon credits for re/afforestation projects set aside in a buffer to 
+# * Share of carbon credits for re/afforestation projects set aside in a buffer to
 # * account for potential failure of permanent carbon storage.
-# * Alternatively, this parameter can be interpreted as a reduction factor 
+# * Alternatively, this parameter can be interpreted as a reduction factor
 # * of the carbon price for non-permanent CDR.
 # * Values > 0 will reduce the incentive for c-price-induced re/afforestation
 cfg$gms$s56_buffer_aff <- 0.5   # def = 0.5

modules/60_bioenergy/1st2ndgen_priced_feb24/equations.gms

@@ -66,10 +66,12 @@ q60_res_2ndgenBE(i2) ..
   =g=
   sum(ct,i60_res_2ndgenBE_dem(ct,i2));
 
-*' Finally, an incentive is provided for the production of 1st generation
-*' bioenergy from oils and ethanol even beyond the exogeneous minimum demand.
-*' The incentive is kept low, but should provide a more realistic
-*' overproduction from couple products.
+*' Finally, an incentive is provided for the production of 1st and 2nd generation
+*' bioenergy beyond the exogeneous minimum demand. 1st generation bioenergy can be incentivized
+*' mass- or energy-based, the former, for comparability, in line with other realizations and constant over time. 
+*' The energy-based incentive can take different forms and is applied to both 1st and 2nd generation. 
+*' Combined with low or fade-out exogenous demands, this is useful to assess bioenergy production potentials, however
+*' the endogenous technological change in $13\_TC may react very strongly and create a positive feedback loop
 
 q60_bioenergy_incentive(i2).. vm_bioenergy_utility(i2)
   =e= sum((ct,k1st60), vm_dem_bioen(i2,k1st60) * (-i60_1stgen_bioenergy_subsidy_tdm(ct)))

scripts/output/extra/disaggregation.R

@@ -38,6 +38,11 @@ land_hr_split_file <- file.path(outputdir, "cell.land_split_0.5.mz")
 land_hr_shr_split_file <- file.path(outputdir, "cell.land_split_0.5_share.mz")
 luh_side_layers <- file.path(outputdir, "luh2_side_layers_0.5.mz")
 bii_hr_out_file <- file.path(outputdir, "cell.bii_0.5.mz")
+peatland_v2_hr_file <- file.path(outputdir, "f58_peatland_area_0.5.mz")
+peatland_on_intact_hr_file <- file.path(outputdir, "f58_peatland_intact_0.5.mz")
+peatland_on_degrad_hr_file <- file.path(outputdir, "f58_peatland_degrad_0.5.mz")
+peatland_hr_out_file <- file.path(outputdir, "cell.peatland_0.5.mz")
+peatland_hr_share_out_file <- file.path(outputdir, "cell.peatland_0.5_share.mz")
 
 cfg <- gms::loadConfig(file.path(outputdir, "config.yml"))
 
@@ -54,6 +59,11 @@ if (length(map_file) > 1) {
 #  Output functions
 # -----------------------------------------
 
+.fixCoords <- function(x) {
+  getSets(x, fulldim = FALSE)[1] <- "x.y.iso"
+  return(x)
+}
+
 .writeDisagg <- function(x, file, comment, message) {
   base::message(message)
   x <- .fixCoords(x)
@@ -84,7 +94,8 @@ if (length(map_file) > 1) {
   map <- readRDS(map_file)
 
   # create full time series
-  land_consv_hr <- new.magpie(map[, "cell"], getYears(land_consv_lr), getNames(wdpa_hr))
+  land_consv_hr <- new.magpie(map[, "cell"], getYears(land_consv_lr), getNames(wdpa_hr),
+                             sets = c("x.y.iso", "year", getSets(wdpa_hr, fulldim = FALSE)[3]))
   land_consv_hr[, getYears(land_consv_hr), ] <- wdpa_hr[, nyears(wdpa_hr), ]
   land_consv_hr[, getYears(wdpa_hr), ] <- wdpa_hr
 
@@ -93,7 +104,8 @@ if (length(map_file) > 1) {
       consv_prio_all <- read.magpie(consv_prio_hr_file)
       consv_prio_hr <- new.magpie(
         cells_and_regions = map[, "cell"],
-        names = getNames(consv_prio_all, dim = 2), fill = 0
+        names = getNames(consv_prio_all, dim = 2), fill = 0,
+        sets = c("x.y.iso", "year", "data")
       )
       iso <- readGDX(gdx, "iso")
       consv_iso <- readGDX(gdx, "policy_countries22")
@@ -301,6 +313,7 @@ land_hr <- interpolateAvlCroplandWeighted(
   snv_pol_shr = snv_pol_shr,
   snv_pol_fader = snv_pol_fader
 )
+land_hr <- .fixCoords(land_hr)
 
 # Write output
 .writeDisagg(land_hr, land_hr_out_file,
@@ -535,4 +548,3 @@ rm(land_hr, peat_hr)
 gc()
 
 message("Finished disaggregation")
-