calc2ndBioDem <- function(datasource) {
if (datasource == "REMIND") {
x <- readSource("REMIND","extensive")
x <- x[,,"Primary Energy Production|Biomass|Energy Crops (EJ/yr)"]*10^3
x <- collapseNames(x)
first_remind_year <- sort(getYears(x))[1]
x <- time_interpolate(x,seq(1995,2150,5),extrapolation_type = "constant")
# set values in initial years that are not existing in REMIND data to zero
x[,getYears(x)<first_remind_year,]<-0
description <- "2nd generation bioenergy demand for different scenarios taken from R2M41 coupled runs"
} else if (datasource == "REMMAG") {
x <- readSource("REMMAG","biodem")
#harmonize historic period
x[,c(1995,2000,2005,2010),] <- collapseNames(x[,c(1995,2000,2005,2010),"SSP2-Ref-SPA0"])
description <- "2nd generation bioenergy demand for different scenarios taken from R17M3 coupled runs"
} else if (datasource == "SSPResults") {
x<-readSource("SSPResults")
x<- collapseNames(x[,,"Primary Energy|Biomass|Energy Crops (EJ/yr)"])*10^3
description <- "2nd generation bioenergy demand for different scenarios taken from IIASA SSP database"
} else if (datasource == "SSP_and_REM") {
ssp <- calcOutput("2ndBioDem",datasource="SSPResults",aggregate = FALSE)
rem <- calcOutput("2ndBioDem",datasource="REMIND",aggregate = FALSE)
ssp <- time_interpolate(ssp,getYears(rem),extrapolation_type = "constant")
x <- mbind(ssp,rem)
# sort scenarios alphabetically
x <- x[,,sort(getNames(x))]
description <- "2nd generation bioenergy demand for different scenarios taken from R2M41 coupled runs and from IIASA SSP database"
} else {
stop("Unknown datasource",datasource)
}
return(list(x=x, weight=NULL,
description=description,
unit="PJ per year",
note="bioenergy is demanded in the country which is expected to produce the bioenergy (demand after trade)"))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.