R/nr_som_loss.R

In pik-piam/nitrogen: Nitrogen Library

#' nr_som_loss
#'
#' @param fulldata data to perform actions on
#' @param nutrient the nutrient
#' @param usable ???
#'
#' @author Jan-Philipp Dietrich
nr_som_loss<-function(fulldata, nutrient="nr", usable=FALSE) {
    if (usable==FALSE) {
        out<-superAggregate(readGDX(fulldata, "p51_nr_som", format="first_found"),aggr_type="sum")
    } else {
        out<-superAggregate(readGDX(fulldata, "p51_nr_som_usable", format="first_found"),aggr_type="sum")  
    }
    out<-unwrap(out)
    out[,2:(dim(out)[2]),] <- out[,1:(dim(out)[2]-1),]
    out<-as.magpie(out)
    out[,1,] <- superAggregate(readGDX(fulldata, "f51_nr_som_1995", format="first_found"),aggr_type="sum")
    out[,2,] <- superAggregate(readGDX(fulldata, "f51_nr_som_1995", format="first_found"),aggr_type="sum")
    if (nutrient=="dm"){
        out=out*15
    } else if (nutrient=="nr") {
        # do nothing
    } else {stop("currently only values for nr and dm exist")}
    return(out)
}