R/nr_atm_dep_pasture.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_atm_dep_pasture
#'
#' @param fulldata data to perform actions on
#'
#' @author Jan-Philipp Dietrich
nr_atm_dep_pasture <- function(fulldata) {
    dep_1995<-land(fulldata)[,,"past"]*(setYears(readGDX(fulldata,"f51_n_i_dep_past", format="first_found")[,"y1995",],NULL))
    manure_past<-dimSums(nr_manure(fulldata)[,,c("fuel","grazing_pasture")],dim=c(3.1,3.2))+setNames(nr_manure_recycl_past(fulldata),NULL)
    # assume that atm_dep grows with manure emissions on pasture, as other inputs are remote
    out<-dep_1995*manure_past/setYears(manure_past[,"y1995",],NULL)
    return(out)
}