Skip to content

Commit

Permalink
apply autoFormat to reportCapitalStock
Browse files Browse the repository at this point in the history
  • Loading branch information
@fbenke-pik
fbenke-pik committed Oct 7, 2024
1 parent 7511045 commit cc809ab
Showing 1 changed file with 109 additions and 109 deletions.
218 changes: 109 additions & 109 deletions R/reportCapitalStock.R
View file
Original file line number Diff line number Diff line change
Expand Up @@ -15,124 +15,124 @@
#' @author Lavinia Baumstark; Michaja Pehl
#' @seealso \code{\link{convGDX2MIF}}
#' @examples
#'
#' \dontrun{reportCapitalStock(gdx)}
#' \dontrun{
#' reportCapitalStock(gdx)
#' }
#'
#' @export
#' @importFrom gdx readGDX
#' @importFrom magclass getYears mbind setNames
#' @importFrom dplyr tribble
reportCapitalStock <- function(gdx,regionSubsetList=NULL,t=c(seq(2005,2060,5),seq(2070,2110,10),2130,2150)) {

module2realisation <- readGDX(gdx, "module2realisation", react = "silent")
tran_mod = module2realisation[module2realisation$modules == "transport", 2]

pm_conv_cap_2_MioLDV <- 650 # The world has ~715million cars in 2005 (IEA TECO2)

# read sets
teall2rlf <- readGDX(gdx,name=c("te2rlf","teall2rlf"),format="first_found")
teue2rlf <- readGDX(gdx,name=c("teue2rlf", "tees2rlf"),format="first_found")
te <- readGDX(gdx,name=c("te"),format="first_found")
# read variables
vm_cap <- readGDX(gdx,name=c("vm_cap"),field="l",format="first_found")
vm_deltaCap <- readGDX(gdx,name=c("vm_deltaCap"),field="l",format="first_found")
v_investcost <- readGDX(gdx,name=c("vm_costTeCapital","v_costTeCapital","v_investcost"),field="l",format="first_found")
vm_cesIO <- readGDX(gdx, name = 'vm_cesIO', field = 'l')
# read parameters
ppfKap_Ind <- readGDX(gdx, name = 'ppfkap_industry_dyn37', react = 'silent')
steel_process_based <- "steel" %in% readGDX(gdx, "secInd37Prc", react='silent')

# calculate maximal temporal resolution
y <- Reduce(intersect,list(getYears(vm_cap),getYears(v_investcost)))
vm_cap <- vm_cap[,y,]
vm_deltaCap <- vm_deltaCap[,y,]
v_investcost <- v_investcost[,y,]

tmp <- NULL

# ---- report transport capital stocks ----
if (tran_mod == "complex"){
LDV35 <- readGDX(gdx,name=c("LDV35"),format="first_found")
tmp <- mbind(tmp,setNames(dimSums( (vm_cap * v_investcost)[teue2rlf]
,dim=c(3.1,3.2)) * 1000, "Est Capital Stock|ESM|Transp vehic (billion US$2017)"))
tmp <- mbind(tmp,setNames(dimSums( (vm_cap * v_investcost)[teall2rlf][,,LDV35]
,dim=c(3.1,3.2)) * 1000, "Est Capital Stock|ESM|Pet/EV LDV (billion US$2017)"))

tmp <- mbind(tmp,
setNames(
dimSums(mbind(vm_cap * v_investcost),
dim = 3) * 1000,
"Estimated Capital Stock|ESM (billion US$2017)"
)
)

# prepare variables
vm_cap <- vm_cap[teall2rlf]
vm_deltaCap <- vm_deltaCap[teall2rlf]

# build reporting
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,LDV35] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est LDV Stock (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apCarElT"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est EV LDV Stock (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apCarH2T"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est H2 LDV Stock (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apCarPeT"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est ICE LDV Stock (million vehicles)"))

tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,LDV35] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est LDV Sales (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apCarElT"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est EV LDV Sales (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apCarH2T"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est H2 LDV Sales (million vehicles)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apCarPeT"] * pm_conv_cap_2_MioLDV,dim=c(3.1,3.2)),"Est ICE LDV Sales (million vehicles)"))

tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,c("apCarDiT","apcarDiEffT","apcarDiEffH2T")],dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Stock|uedit (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apCarDiT"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Stock|apCarDiT (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apcarDiEffT"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Stock|apcarDiEffT (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_cap[,,"apcarDiEffH2T"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Stock|apcarDiEffH2T (arbitrary unit)"))

tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,c("apCarDiT","apcarDiEffT","apcarDiEffH2T")],dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Sales|uedit (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apCarDiT"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Sales|apCarDiT (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apcarDiEffT"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Sales|apcarDiEffT (arbitrary unit)"))
tmp <- mbind(tmp,setNames(dimSums(vm_deltaCap[,,"apcarDiEffH2T"] ,dim=c(3.1,3.2)), "Services and Products|Transport|non-LDV|Sales|apcarDiEffH2T (arbitrary unit)"))


## add global values
tmp <- mbind(tmp,dimSums(tmp,dim=1))
reportCapitalStock <- function(gdx, regionSubsetList = NULL, t = c(seq(2005, 2060, 5), seq(2070, 2110, 10), 2130, 2150)) {

module2realisation <- readGDX(gdx, "module2realisation", react = "silent")
tran_mod <- module2realisation[module2realisation$modules == "transport", 2]

pm_conv_cap_2_MioLDV <- 650 # The world has ~715million cars in 2005 (IEA TECO2)

# read sets
teall2rlf <- readGDX(gdx, name = c("te2rlf", "teall2rlf"), format = "first_found")
teue2rlf <- readGDX(gdx, name = c("teue2rlf", "tees2rlf"), format = "first_found")
# read variables
vm_cap <- readGDX(gdx, name = c("vm_cap"), field = "l", format = "first_found")
vm_deltaCap <- readGDX(gdx, name = c("vm_deltaCap"), field = "l", format = "first_found")
v_investcost <- readGDX(gdx, name = c("vm_costTeCapital", "v_costTeCapital", "v_investcost"), field = "l", format = "first_found")
vm_cesIO <- readGDX(gdx, name = "vm_cesIO", field = "l")
# read parameters
ppfKap_Ind <- readGDX(gdx, name = "ppfkap_industry_dyn37", react = "silent")
steel_process_based <- "steel" %in% readGDX(gdx, "secInd37Prc", react = "silent")

# calculate maximal temporal resolution
y <- Reduce(intersect, list(getYears(vm_cap), getYears(v_investcost)))
vm_cap <- vm_cap[, y, ]
vm_deltaCap <- vm_deltaCap[, y, ]
v_investcost <- v_investcost[, y, ]

tmp <- NULL

# ---- report transport capital stocks ----
if (tran_mod == "complex") {
LDV35 <- readGDX(gdx, name = c("LDV35"), format = "first_found")
tmp <- mbind(tmp, setNames(dimSums((vm_cap * v_investcost)[teue2rlf],
dim = c(3.1, 3.2)) * 1000, "Est Capital Stock|ESM|Transp vehic (billion US$2017)"))
tmp <- mbind(tmp, setNames(dimSums((vm_cap * v_investcost)[teall2rlf][, , LDV35],
dim = c(3.1, 3.2)) * 1000, "Est Capital Stock|ESM|Pet/EV LDV (billion US$2017)"))

tmp <- mbind(tmp,
setNames(
dimSums(mbind(vm_cap * v_investcost),
dim = 3) * 1000,
"Estimated Capital Stock|ESM (billion US$2017)"
)
)

# prepare variables
vm_cap <- vm_cap[teall2rlf]
vm_deltaCap <- vm_deltaCap[teall2rlf]

# build reporting
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , LDV35] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est LDV Stock (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apCarElT"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est EV LDV Stock (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apCarH2T"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est H2 LDV Stock (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apCarPeT"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est ICE LDV Stock (million vehicles)"))

tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , LDV35] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est LDV Sales (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apCarElT"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est EV LDV Sales (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apCarH2T"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est H2 LDV Sales (million vehicles)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apCarPeT"] * pm_conv_cap_2_MioLDV, dim = c(3.1, 3.2)), "Est ICE LDV Sales (million vehicles)"))

tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , c("apCarDiT", "apcarDiEffT", "apcarDiEffH2T")], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Stock|uedit (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apCarDiT"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Stock|apCarDiT (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apcarDiEffT"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Stock|apcarDiEffT (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_cap[, , "apcarDiEffH2T"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Stock|apcarDiEffH2T (arbitrary unit)"))

tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , c("apCarDiT", "apcarDiEffT", "apcarDiEffH2T")], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Sales|uedit (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apCarDiT"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Sales|apCarDiT (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apcarDiEffT"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Sales|apcarDiEffT (arbitrary unit)"))
tmp <- mbind(tmp, setNames(dimSums(vm_deltaCap[, , "apcarDiEffH2T"], dim = c(3.1, 3.2)), "Services and Products|Transport|non-LDV|Sales|apcarDiEffH2T (arbitrary unit)"))


## add global values
tmp <- mbind(tmp, dimSums(tmp, dim = 1))
}

# ---- report industry energy efficiency capital stocks ----
if (!is.null(ppfKap_Ind) && 0 < length(ppfKap_Ind)) {
mixer <- tribble(
~pf, ~name,
"kap_cement", "Cement",
"kap_chemicals", "Chemicals",
"kap_otherInd", "other")

if (!steel_process_based) {
mixer <- bind_rows(mixer, tribble(
~pf, ~name,
"kap_steel_primary", "Primary Steel",
"kap_steel_secondary", "Secondary Steel"))
}

# ---- report industry energy efficiency capital stocks ----
if (!is.null(ppfKap_Ind) & 0 < length(ppfKap_Ind)) {
mixer <- tribble(
~pf, ~name,
'kap_cement', 'Cement',
'kap_chemicals', 'Chemicals',
'kap_otherInd', 'other')
if (0 != length(setdiff(ppfKap_Ind, mixer$pf))) {
warning(paste("Unknown ppfKap_industry_dyn37 entity.",
"Adjust remind2::reportCapitalStock()"))
}

if (!steel_process_based) {
mixer <- bind_rows(mixer, tribble(
~pf, ~name,
'kap_steel_primary', 'Primary Steel',
'kap_steel_secondary', 'Secondary Steel'))
}

if (0 != length(setdiff(ppfKap_Ind, mixer$pf))) {
warning(paste('Unknown ppfKap_industry_dyn37 entity.',
'Adjust remind2::reportCapitalStock()'))
}

if (0 != length(setdiff(mixer$pf, ppfKap_Ind))) {
warning(paste('Missing ppfKap_industry_dyn37 entity.',
'Adjust remind2::reportCapitalStock()'))
}

eek_Ind <- setNames(vm_cesIO[,y,ppfKap_Ind],
paste0('Capital|Energy Efficiency|Industry|',
mixer[mixer$pf %in% ppfKap_Ind,][['name']],
' (billion US$2017)'))
# add industry EEK and global totals
tmp <- mbind(tmp, mbind(eek_Ind, dimSums(eek_Ind, dim = 1)))
if (0 != length(setdiff(mixer$pf, ppfKap_Ind))) {
warning(paste("Missing ppfKap_industry_dyn37 entity.",
"Adjust remind2::reportCapitalStock()"))
}

# ---- add region aggregates ----
if (!is.null(regionSubsetList))
tmp <- mbind(tmp, calc_regionSubset_sums(tmp, regionSubsetList))
eek_Ind <- setNames(vm_cesIO[, y, ppfKap_Ind],
paste0("Capital|Energy Efficiency|Industry|",
mixer[mixer$pf %in% ppfKap_Ind, ][["name"]],
" (billion US$2017)"))
# add industry EEK and global totals
tmp <- mbind(tmp, mbind(eek_Ind, dimSums(eek_Ind, dim = 1)))
}

# ---- add region aggregates ----
if (!is.null(regionSubsetList))
tmp <- mbind(tmp, calc_regionSubset_sums(tmp, regionSubsetList))

getSets(tmp)[3] <- "variable"
return(tmp)
getSets(tmp)[3] <- "variable"
return(tmp)
}

0 comments on commit cc809ab

Please sign in to comment.