R/reportharvested_area_timber.R

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

#' @title reportharvested_area_timber
#' @description reports MAgPIE harvested area for timber.
#' 
#' @export
#' 
#' @param gdx GDX file
#' @return Area harvested for timber production
#' @author Abhijeet Mishra
#' @examples
#' 
#'   \dontrun{
#'     x <- reportharvested_area_timber(gdx)
#'   }
#' 

reportharvested_area_timber<-function(gdx){
  a <- NULL
  
  timber <- FALSE
  if (as.numeric(readGDX(gdx, "s32_hvarea")) > 0 & as.numeric(readGDX(gdx, "s35_hvarea")) > 0) timber <- TRUE
  
  if(timber){
    a <- harvested_area_timber(gdx,level = "regglo")
    a <- mbind(a,setNames(dimSums(a,dim=3),"Total"))
    getNames(a) <- paste0("Resources|Timber operations|Harvested area for timber|",getNames(a))
    getNames(a) <- paste0(getNames(a)," (mha per yr)")
  } else {message("Disabled (no timber) ", appendLF = FALSE)}
  
  return(a)
}