#' @title reportCroparea
#' @description reports croparea
#'
#' @importFrom magpiesets reporthelper summationhelper
#' @export
#'
#' @param gdx GDX file
#' @param detail if detail=FALSE, the subcategories of groups are not reported (e.g. "soybean" within "oilcrops")
#' @return Croparea as MAgPIE object (million Ha/yr)
#' @author Florian Humpenoeder
#' @examples
#' \dontrun{
#' x <- reportCroparea(gdx)
#' }
#'
reportCroparea <- function(gdx, detail = FALSE) {
x <- NULL
out <- croparea(gdx, level = "regglo", products = "kcr",
product_aggr = FALSE, water_aggr = TRUE)
# Adapting reportingnames for inclusing of fallow land and tree cover in reportLandUse
ReportingnameCroparea <- "Resources|Land Cover|Cropland|Croparea"
out <- reporthelper(x = out, dim = 3.1,
level_zero_name = ReportingnameCroparea, detail = detail)
getNames(out) <- paste(gsub("\\.", "|", getNames(out)), "(million ha)", sep = " ")
out <- summationhelper(out, sep = "+")
out2 <- croparea(gdx, level = "regglo", products = "kcr",
product_aggr = FALSE, water_aggr = FALSE)
total <- dimSums(out2, dim = 3.1)
getNames(total) <- paste(ReportingnameCroparea, "|", reportingnames(getNames(total)),
" (million ha)", sep = "")
out2 <- reporthelper(x = out2, dim = 3.1,
level_zero_name = ReportingnameCroparea, detail = detail)
getNames(out2) <- paste(gsub("\\.", "|", getNames(out2)), "(million ha)", sep = " ")
out2 <- mbind(out2, total)
out2 <- summationhelper(out2, sep = "++")
x <- mbind(out, out2)
return(x)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.