#' @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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.