R/nr_fixation_total.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_fixation_total
#'
#' @param fulldata data to perform actions on
#' @param som_loss ???
#'
#' @author Jan-Philipp Dietrich
nr_fixation_total<-function(fulldata,som_loss=TRUE){
    budget <- nr_cropland_budget(fulldata)
    if(som_loss==TRUE){
        out<-dimSums(budget[,,c("fix_fertilizer","fix_crops","fix_area","SOM_loss")],dim=3.1)
    } else {
        out<-dimSums(budget[,,c("fix_fertilizer","fix_crops","fix_area")],dim=3.1)
    }
    return(setNames(out,"fixation"))
}