R/reportPastSoilCarbon.R

In pik-piam/magpie4: MAgPIE outputs R package for MAgPIE version 4.x

#' @title reportPastSoilCarbon
#' @description reports pasture soil carbon
#'
#' @export
#'
#' @param gdx GDX file
#' @return Soil carbon values as magpie objetc
#' @author Marcos Alves
#' @examples
#' \dontrun{
#' x <- reportPastSoilCarbon(gdx)
#' }
#'
reportPastSoilCarbon <- function(gdx) {
  
  x <- NULL
  soilc_cg_target <- NULL
  past_ha_c <- NULL
  # read in data
  try({soilc_cg_target <- readGDX(gdx, "ov_soilc_target", format = "simplest", react = "silent")[, , list("type" = "level")]})
  try({past_ha_c <- readGDX(gdx, "ov_past_area", format = "simplest", react = "silent")[, , list("past_mngt" = "cont_grazing", "type" = "level", "w" = "rainfed")]})
  
  if (!any(c(is.null(soilc_cg_target),is.null(past_ha_c)))) {
  soilc_cg_target_reg <- gdxAggregate(gdx, soilc_cg_target, to = "regglo", weight = past_ha_c, absolute = F)
  soilc_cg_change_reg <- soilc_cg_target_reg - soilc_cg_target_reg[,1,]
  getYears(soilc_cg_change_reg) <- getYears(soilc_cg_target_reg)

  # aggreate and rename
  x <- mbind(x, setNames(soilc_cg_target_reg, "Resources|Soil Carbon Target|Pasture|Continuous grazing|Density (tC per ha)"))
  x <- mbind(x, setNames(soilc_cg_change_reg, "Resources|Soil Carbon Target Change|Pasture|Continuous grazing|Density (tC per ha)"))
  } else { message("Soil carbon was not emulated", appendLF = FALSE) }
  return(x)
}