R/nr_harvest.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_harvest
#'
#' @param fulldata data to perform actions on
#' @param nutrient the nutrient
#'
#' @author Jan-Philipp Dietrich

nr_harvest<-function(fulldata, nutrient="nr") {
    nr<-collapseNames(readGDX(fulldata, "im_attributes_harvest", format="first_found")[,,nutrient])
    nr<-production(fulldata, level="reg", crop_aggr=FALSE,water="sum") * collapseNames(nr[,,-which(getNames(nr)=="pasture")])
    return(nr)
}