R/nr_snupe.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_snupe
#'
#' @param fulldata data to perform actions on
#' @param level the level
#'
#' @author Jan-Philipp Dietrich
nr_snupe<-function(fulldata,level="reg"){
    budget <- superAggregate(nr_cropland_budget(fulldata),aggr_type="sum",level=level)
    budget_out<- -dimSums(budget[,,c("harvest","ag_residues","bg_residues","fix_crops","seed")],dim=3.1)
    budget_in<-dimSums(budget[,,c("fix_fertilizer","manure","fix_area",      
                                  "ag_residues_rec","bg_residues_rec","atm_deposition","SOM_loss")],dim=3.1)
    out<-budget_out/budget_in
    return(setNames(out,"snupe"))
}