inst/extdata/saveDataFiles.R

In mengqi-z/plutus: A comprehensive mapping package

# Converting raw data into package data
library(usethis)
library(data.table)
library(XML)

rawDataFolder = "E:/NEXO-UA/Github/plutus/inst/extdata/"

# # exampleDataCSV
# exampleDataCSV = data.table::fread(paste0(rawDataFolder,"exampleData.csv"),header=T)
# exampleDataCSV
# use_data(exampleDataCSV, overwrite=T)
#
#
# # exampleDataR
# exampleDataR = data.frame(sector=c("a","b",'c'),
#                           coefficients=c(1,2,3))
# exampleDataR
# use_data(exampleDataR, overwrite=T)



#-------------------
# Data
#-------------------

# Metis XMl Query Files
xmlFilePath = paste0(rawDataFolder,"ElecQueries.xml")
xmlfile <- XML::xmlTreeParse(xmlFilePath)
xmltop <- XML::xmlRoot(xmlfile)
top <- XML::xmlNode(XML::xmlName(xmltop))
for(i in 1:length(xmltop)){
  top <- XML::addChildren(top, xmltop[[i]])
}
xmlElecQueries <- top
use_data(xmlElecQueries, overwrite=T)

# Capacity factors
data_capfactors <- data.table::fread(file=paste(rawDataFolder,"capacity_factor_by_elec_gen_subsector.csv",sep=""),skip=5,encoding="Latin-1")
data_cap_cost_tech <- data.table::fread(paste(rawDataFolder,"L2233.GlobalTechCapital_elecPassthru.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_cap_cost_cool <- data.table::fread(paste(rawDataFolder,"L2233.GlobalTechCapital_elec_cool.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_cap_cost_int_tech <- data.table::fread(paste(rawDataFolder,"L2233.GlobalIntTechCapital_elec.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_cap_cost_int_cool <- data.table::fread(paste(rawDataFolder,"L2233.GlobalIntTechCapital_elec_cool.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_A23.globaltech_retirement <- data.table::fread(paste(rawDataFolder,"A23.globaltech_retirement.csv",sep=""), skip=1)
data_capac_fac <- data.table::fread(paste(rawDataFolder,"L223.GlobalTechCapFac_elec.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_capac_fac_int <- data.table::fread(paste(rawDataFolder,"L223.GlobalIntTechCapFac_elec.csv", sep=""), skip=1, stringsAsFactors = FALSE)
data_tech_mapping <- data.table::fread(paste(rawDataFolder,"agg_tech_mapping.csv", sep=""), skip=1)
use_data(data_capfactors, overwrite=T)
use_data(data_cap_cost_tech, overwrite=T)
use_data(data_cap_cost_cool, overwrite=T)
use_data(data_cap_cost_int_tech, overwrite=T)
use_data(data_cap_cost_int_cool, overwrite=T)
use_data(data_A23.globaltech_retirement, overwrite=T)
use_data(data_capac_fac, overwrite=T)
use_data(data_capac_fac_int, overwrite=T)
use_data(data_tech_mapping, overwrite=T)

#-------------------
# Example Data
#-------------------

# Example .proj file
# Run this
metis.readgcam(gcamdatabase = "E:/NEXO-UA/GCAM-Workspace/gcam-core_LAC_v02_5Nov2019/output/FinalRuns/IDBNexus_gcam5p3_HadGEM2-ES_rcp8p5",
               paramsSelect = c(
                 "elecNewCapCost",
                 "elecNewCapGW",
                 "elecAnnualRetPrematureCost",
                 "elecAnnualRetPrematureGW",
                 "elecCumCapCost",
                 "elecCumCapGW",
                 "elecCumRetPrematureCost",
                 "elecCumRetPrematureGW"
               ),
               dataProjFile = "E:/NEXO-UA/Github/plutus/inst/extdata/Example_gcam5p3_HadGEM2-ES_rcp8p5.proj",
               saveData = F,
               reReadData = T)
exampleGCAMproj <- rgcam::loadProject('E:/NEXO-UA/Github/plutus/inst/extdata/Example_gcam5p3_HadGEM2-ES_rcp8p5.proj')
use_data(exampleGCAMproj, overwrite = T)