R/calcValidEmisLucGloGasser.R

In pik-piam/mrvalidation: madrat data preparation for validation purposes

#' @title calcValidEmisLucGloGasser
#' 
#' @description Returns historical LUC emissions
#' 
#' @param subtype Available subtypes are:
#'  \itemize{
#'  \item historical:
#'  \itemize{
#'  \item Gasser_2020
#'  \item LUH2_GCB_2019
#'  \item FRA_2015
#'  }
#'  }
#' @return list of magpie object with data and weight
#' @author Abhijeet Mishra
calcValidEmisLucGloGasser <- function(subtype="Gasser_2020") {
  
  if (subtype %in% c("Gasser_2020","LUH2_GCB_2019","FRA_2015")) {
    out <- collapseNames(readSource("Gasser",subtype = "global",convert=FALSE)[,,subtype]) * 1000 * 44/12 ## Conversion from Pg to Mt and C to CO2
    gross_luc        <- dimSums(out[,,c("Deforestation for cropland","Other deforestation","Other natural land appropriation","Conversions among anthrop biomes")],dim=3)
    regrowth_luc     <- dimSums(out[,,c("Reforestation and afforestation","Other natural land (re)establishment")],dim=3)
    wood_harvest_luc <- dimSums(out[,,c("Wood harvest")],dim=3)
    getNames(gross_luc)        <- "Emissions|CO2|Land|Land-use Change|+|Gross LUC (Mt CO2/yr)"
    getNames(regrowth_luc)     <- "Emissions|CO2|Land|Land-use Change|+|Regrowth (Mt CO2/yr)"
    getNames(wood_harvest_luc) <- "Emissions|CO2|Land|Land-use Change|+|Wood products (Mt CO2/yr)"
    out <- mbind(gross_luc,regrowth_luc,wood_harvest_luc)
    out <- add_dimension(out, dim=3.1, add="scenario", nm="historical")
    out <- add_dimension(out, dim=3.2, add="model",nm=subtype)
    out <- setYears(out,"y2009")
    out <- time_interpolate(dataset = out,interpolated_year = paste0("y",2010:2018),extrapolation_type = "constant",integrate_interpolated_years = TRUE)
    names(dimnames(out))[3] <- "scenario.model.variable"
  } else stop("Invalid subtype. See function description for valid subtypes.") 
  
  return(list(x=out,
              weight=NULL,
              unit="Mt CO2/yr",
              min=-10000,
              description="Historical land-use change CO2 emissions from different sources in Mt CO2/yr")
  )
}