#' @title reportNitrogenPollution
#' @description Reports total Nitrogen Pollution as the sum of surplus from cropland, pasture, awms, consumption and non-agricutlural land
#'
#' @importFrom magpiesets reportingnames
#' @export
#'
#' @param gdx GDX file
#' @param dir magpie output directory that contains gridded Nitrogen Data
#' @author Benjamin Leon Bodirsky
#' @seealso
#' \code{\link{NitrogenBudget}}
#'
#' @examples
#' \dontrun{
#' x <- reportNitrogenPollution(gdx)
#' }
#'
reportNitrogenPollution <- function(gdx, dir = ".") {
cropland <- NitrogenBudget(gdx, level = "reg", dir = dir)[,,"surplus"]
pasture <- NitrogenBudgetPasture(gdx, level = "reg", dir = dir)[,,"surplus"]
awms <- dimSums(readGDX(gdx,"ov_manure_confinement")[,,"level"][,,"nr"],dim=3)-dimSums(readGDX(gdx,"ov_manure_recycling")[,,"level"][,,"nr"],dim=3)
consumption <- dimSums(demand(gdx,attributes="nr",product_aggr=TRUE)[,,c("food","other_util","bioenergy","waste")],dim=3)
nonagland <- dimSums(NitrogenBudgetNonagland(gdx, level = "reg", dir = dir)[,,"surplus"],dim=3)
combined <- mbind(
setNames(cropland, "Resources|Nitrogen|Pollution|Surplus|+|Cropland (Mt Nr/yr)"),
setNames(pasture, "Resources|Nitrogen|Pollution|Surplus|+|Pasture (Mt Nr/yr)"),
setNames(awms, "Resources|Nitrogen|Pollution|Surplus|+|Animal Waste Management (Mt Nr/yr)"),
setNames(nonagland, "Resources|Nitrogen|Pollution|Surplus|+|Non-agricultural land (Mt Nr/yr)"),
setNames(consumption, "Resources|Nitrogen|Pollution|Surplus|+|End of life losses (Mt Nr/yr)")
)
combined <- mbind(
setNames(dimSums(combined,dim=3),"Resources|Nitrogen|Pollution|Surplus (Mt Nr/yr)"),
combined
)
landAndAWMS <- dimSums(mbind(cropland, pasture, awms, nonagland), dim = 3)
combined <- mbind(
setNames(landAndAWMS, "Resources|Nitrogen|Nutrient surplus from land and manure management (Mt Nr/yr)"),
combined
)
combined <- mbind(
combined,
setCells(dimSums(combined,dim=1),"GLO")
)
return(combined)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.