#' 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)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.