Difference between revisions of "Energy balance in Stuttgart"

library(OpasnetUtils)
library(xtable)

dependencies <- data.frame(Name = c(
	"Frac_CHP_renew",
	"Ind_process_heat",
	"Cons_Home_electricity",
	"Cons_Ind_electricity",
	"Cons_Commerce_electricity",
	"Cons_Municip_electricity",
	"Cons_Home_heat",
	"Cons_Ind_heat",
	"Cons_Commerce_heat",
	"Cons_Municip_heat"
), Key = "k5NNrR2Q0VCJharl")

formula <- function(dependencies, ...) {

	ComputeDependencies(dependencies, ...)

	## Define the energy transformation matrix ETM.

	# Get the data from Opasnet Base.

	ETM <- tidy(op_baseGetData("opasnet_base", "Op_en5469"), objname = "ETM", direction = "wide")

	# Add user-defined data to it.

	additions <- data.frame(
		Equation = "Fraction CHP renewable",
		Col = "CHP peat",
		ETMResult = Frac_CHP_renew,
		ETMDescription = "Renewable amount divided by peat. BAU: 78.1/1435.7 = 0.0544 Actual number given by user."
	)
	ETM <- rbind(ETM, additions)

	# Reshape the table to make a matrix, and replace NA with 0.

	ETM <- reshape(ETM, timevar = "Col", idvar = "Equation", v.names = "ETMResult", direction = "wide")
	colnames(ETM) <- gsub("ETMResult.", "", colnames(ETM))

	print(xtable(ETM), type = 'html') # Show the matrix with explanations.

	ETM <- as.matrix(ETM[, !colnames(ETM) %in% c("Equation", "ETMDescription")])

	for(i in 1:ncol(ETM)) {
		ETM[is.na(ETM[, i]), i] <- 0
	}

	# Define the outputs for equations (for most equations, it is a zero sum game).

	out <- c(rep(0, 18), # Define the consumptions (and other known variables).
		Ind_process_heat,
		Cons_Home_electricity,
		Cons_Ind_electricity,
		Cons_Commerce_electricity,
		Cons_Municip_electricity,
		Cons_Home_heat,
		Cons_Ind_heat,
		Cons_Commerce_heat,
		Cons_Municip_heat
	)

	out <- solve(ETM, out) # Solve the matrix.
	out <- data.frame(Energy.class = colnames(ETM), Result = out)

	return(out)
}

energy.balance.Kuopio <- new("ovariable",
	name = "energy.balance.Kuopio",
	output = formula(dependencies)
)

print(xtable(energy.balance.Kuopio@output), type = 'html')
#temp <- EvalOutput(energy.balance.Kuopio)
#print(xtable(temp@output), type = 'html')

objects.put(energy.balance.Kuopio)
#	Ind_process_heat,
#	Cons_Home_electricity,
#	Cons_Ind_electricity,
#	Cons_Commerce_electricity,
#	Cons_Municip_electricity,
#	Cons_Home_heat,
#	Cons_Ind_heat,
#	Cons_Commerce_heat,
#	Cons_Municip_heat
#)

cat("Object energy.balance.Kuopio saved. Write down the key of the model run page for further use.\n")