moinput: Model Operations Input Data Library

#' fullREMIND
#' 
#' Function that produces the complete regional data set required for the
#' REMIND model.
#' 
#' @param rev data revision which should be used as input (positive numeric).
#' \code{\link{setConfig}} (e.g. for setting the mainfolder if not already set
#' properly).
#' @author Lavinia Baumstark
#' @seealso
#' \code{\link{readSource}},\code{\link{getCalculations}},\code{\link{calcOutput}}
#' @examples
#' 
#' \dontrun{ 
#' fullREMIND()
#' }
#' 
fullREMIND <- function(rev=0) {
  
  rem_years <- seq(2005,2150,5)
  rem_years_hist <- seq(1990,2150,5)
  
  #-------------- macro-economic parameters -----------------------------------------------------------
  calcOutput("Population", years=rem_years_hist,      round=8,  file="f_pop.cs3r")       
  calcOutput("Labour",     years=rem_years,           round=8,  file="f_lab.cs3r")       
  calcOutput("GDPppp",     years=rem_years_hist,      round=8,  file="f_gdp.cs3r") 
  calcOutput("RatioPPP2MER",                          round=8,  file="pm_shPPPMER.cs4r")
  calcOutput("MacroInvestments",                      round=8,  file="p01_boundInvMacro.cs4r")
  calcOutput("SubsStationary", sector = "transport" , round=2,  file="f21_tau_fe_sub_transport.cs4r") 
  calcOutput("TaxesStationary",sector = "transport",  round=2,  file="f21_tau_fe_tax_transport.cs4r")
  calcOutput("SubsStationary", sector = "bit_st",     round=2,  file="f21_tau_fe_sub_bit_st.cs4r") 
  calcOutput("TaxesStationary",sector = "bit_st",     round=2,  file="f21_tau_fe_tax_bit_st.cs4r")
  calcOutput("TaxConvergence",                        round=2,  file="f21_tax_convergence.cs4r")
  calcOutput("TaxLimits", subtype = "maxFeSubsidy",   round=2,  file="f21_max_fe_sub.cs4r")
  calcOutput("TaxLimits", subtype = "maxPeSubsidy",   round=2,  file="f21_max_pe_sub.cs4r")
  calcOutput("TaxLimits", subtype = "propFeSubsidy",  round=2,  file="f21_prop_fe_sub.cs4r")
  calcOutput("SubsStationary", sector = "extraction", round=2,  file="p21_tau_fuEx_sub.cs4r") 
  calcOutput("TaxXport",                              round=2,  file="p21_tau_xpres_tax.cs4r")   # not default, overwritten with 0
  calcOutput("Capital",                               round=6,  file="p29_capitalQuantity.cs4r")
  calcOutput("Capital",   subtype = "CapitalUnit",    round=6,  file="f29_capitalUnitProjections.cs4r")
  calcOutput("FEdemand",  subtype = "FE",             round=8,  file="pm_fe_demand.cs4r")
  calcOutput("FEdemand",  subtype = "ES",             round=6,  file="pm_es_demand.cs4r")
  calcOutput("EnergyEffPaths",                        round=6,  file="p29_efficiency_growth.cs4r")
  calcOutput("Floorspace",                            round=6,  file="p36_floorspace_scen.cs4r")
  calcOutput("WeightNash",                            round=6,  file="p80_eoWeights_fix.cs4r")
  calcOutput("IntertempElastSubst",                   round=6,  file="pm_ies.cs4r")
  calcOutput("TimePref",                              round=6,  file="p23_prtp.cs4r")
  calcOutput("CO2Prices",                             round=2,  file="pm_taxCO2eqHist.cs4r")
  calcOutput("RiskPremium",                           round=6,  file="pm_risk_premium.cs4r")
  calcOutput("NetForeignAsset",                       round=6,  file="pm_nfa_start.cs4r")
  calcOutput("NetForeignAsset",                       round=6,  file="pm_nfa_start.cs4r")
  calcOutput("CalibrationPrices",                     round=6,  file="p29_cesdata_price.cs4r")
  calcOutput("Theil",                                 round=8,  file="f_ineqTheil.cs4r")
  calcOutput("DevelopmentState",                      round=4,  file="f_developmentState.cs3r")
  
  #-------------- energy services parameter -----------------------------------------------------------
  calcOutput("FEdemand", subtype = "EsUeFe_in",       round=8, file="p36_serviceInputs.cs4r")
  calcOutput("FEdemand", subtype = "EsUeFe_out",      round=8, file="p36_serviceOutputs.cs4r")
  
  #-------------- emission parameter ------------------------------------------------------------------
  calcOutput("EconometricEmiParameter",                             round=5, file="p_emineg_econometric.cs3r") 
  calcOutput("EmissionsTe",                                         round=5, file="p_boundEmi.cs4r")
  calcOutput("HistEmissions",subtype="sector",                      round=8, file="p_histEmiSector.cs4r")
  calcOutput("HistEmissions",subtype="MAC",                         round=8, file="p_histEmiMac.cs4r")
  calcOutput("EmiCO2LandUse",                                       round=5, file="p_macPolCO2luc.cs4r")
  calcOutput("MacBaseLandUse",subtype="DirectlyFromMAgPIE",         round=5, file="f_macBaseMagpie.cs4r")
  calcOutput("MacBaseLandUse",subtype="Exogenous",                  round=5, file="f_macBaseExo.cs4r")
  calcOutput("MACCsCO2",                                            round=5, file="p_abatparam_CO2.cs4r", aggregate=FALSE) 
  calcOutput("EmiMac",                                              round=5, file="p_macBase2005.cs4r")    
  calcOutput("EmiMac1990",                                          round=5, file="p_macBase1990.cs4r")
  calcOutput("MACCbaseN2O",                                         round=5, file="p_macBaseVanv.cs4r")   
  calcOutput("MACCsCH4",                                            round=6, file="p_abatparam_CH4.cs4r")
  calcOutput("MACCsN2O",                                            round=6, file="p_abatparam_N2O.cs4r")
  calcOutput("FGas",                                                round=6, file="f_emiFgas.cs4r")
  calcOutput("EmiFossilFuelExtr",                                   round=6, file="p_emiFossilFuelExtr.cs4r")
  calcOutput("Region2MAGICC",                                       round=6, file="p_regi_2_MAGICC_regions.cs3r")
  calcOutput("EmiPollutantExo", subtype="AviationShipping",         round=6, file="f11_emiAPexoGlob.cs4r",aggregate=FALSE)
  calcOutput("EmiPollutantExo", subtype="Waste",                    round=6, file="f11_emiAPexo.cs4r")
  calcOutput("EmiAirPollLandUse",                                   round=6, file="f11_emiAPexoAgricult.cs4r")
  calcOutput('GAINSEmi',       subtype="emissions",                 round=5, file="emi_gains.cs4r")
  calcOutput('GAINSEmi',       subtype="emission_factors",          round=5, file="ef_gains.cs4r")
  calcOutput('GAINSEmi',       subtype="emissions_starting_values", round=5, file="f11_emiAPexsolve.cs4r")
  calcOutput("EmissionFactors",subtype="emission_factors",          round=5, file="f11_emiFacAP.cs4r")
  
  #-------------- energy/technology parameters ---------------------------------------------------------
  calcOutput("PotentialHydro",                        round=3,  file="f_maxProdGradeRegiHydro.cs3r") 
  calcOutput("PotentialWind",                         round=3,  file="f_maxProdGradeRegiWind.cs3r")    
  calcOutput("PotentialGeothermal",                   round=3,  file="f_maxProdGeothermal.cs3r")
  calcOutput("PotentialWeathering",                   round=3,  file="f33_maxProdGradeRegiWeathering.cs3r")
  calcOutput("PotentialWeathering",                   round=3,  file="f33_maxProdGradeRegiWeathering.cs3r")
  calcOutput("CostsTrade",                            round=5,  file="pm_costsPEtradeMp.cs4r")
  calcOutput("CostsTradePeFinancial",                 round=5,  file="pm_costsTradePeFinancial.cs3r")
  calcOutput("ShareCHP",                              round=3,  file="f32_shCHP.cs4r")
  calcOutput("CapacityOffset",                        round=5,  file="p_adj_deltacapoffset.cs4r") 
  calcOutput("CoolingSharesAll",                      round=2,  file="CoolingShares_time.cs4r")
  calcOutput("WaterConsCoef",                         round=3,  file="WaterConsCoef.cs4r",aggregate=FALSE)
  calcOutput("WaterWithCoef",                         round=3,  file="WaterWithCoef.cs4r",aggregate=FALSE)
  calcOutput("IO",   subtype="output",                round=8,  file="f04_IO_output.cs4r")
  calcOutput("IO",   subtype="input",                 round=8,  file="f04_IO_input.cs4r")
  calcOutput("IO",   subtype="trade",                 round=8,  file="f_IO_trade.cs4r")
  calcOutput("ShareIndFE",                            round=3,  file="p37_shIndFE.cs3r")
  calcOutput("Capacity", subtype="capacityByTech",    round=6,  file="p_histCap.cs3r") 
  calcOutput("Capacity", subtype="capacityByPE",      round=6,  file="p_PE_histCap.cs3r") 
  calcOutput("CapacityFactor",                        round=6,  file="f_cf.cs3r")
  calcOutput("StorageFactor",                         round=6,  file="f32_factorStorage.cs4r")
  calcOutput("GridFactor",                            round=6,  file="p32_grid_factor.cs4r")
  calcOutput("FEShares", subtype="ind_coal",          round=5,  file="p_share_ind_fesos.cs4r")
  calcOutput("FEShares", subtype="ind_bio",           round=5,  file="p_share_ind_fesos_bio.cs4r")
  calcOutput("FEShares", subtype="ind_liq",           round=5,  file="p_share_ind_fehos.cs4r")
  calcOutput("BunkersTransportShare",                 round=5,  file="pm_bunker_share_in_nonldv_fe.cs4r")
  calcOutput("Solar",                                 round=5,  file="f_dataRegiSolar.cs3r")
  calcOutput("CapacityEV",                            round=8,  file="pm_boundCapEV.cs4r")
  calcOutput("CapacityNuclear",                       round=5,  file="pm_NuclearConstraint.cs4r")
  calcOutput("CCScapacity",                           round=8,  file="pm_boundCapCCS.cs4r")
  calcOutput("CCSbounds",                             round=8,  file="p_boundCapCCSindicator.cs4r")
  calcOutput("LimitCCS",                              round=8,  file="pm_dataccs.cs3r")
  calcOutput("TranspEff",                             round=8,  file="f35_transp_eff.cs3r")
  calcOutput("VintagesTransport",                     round=8,  file="f35_factorVintages.cs3r")
  calcOutput("BiomassPrices",                                         round=6,  file="f30_bioen_price.cs4r")
  calcOutput("ResFor2ndBioengery",                   years=rem_years, round=5,  file="p30_biolcResidues.cs3r")
  calcOutput("1stBioDem",    subtype ="ethanol_oils",years=rem_years, round=5,  file="p30_bio1stgen.cs3r")
  calcOutput("MAgPIEReport", subtype ="ProductionBiomass",            round=8,  file="p30_biolcProductionLookup.cs4r")
  calcOutput("MAgPIEReport", subtype ="CostTotal",                    round=8,  file="p26_totLUcostLookup.cs4r")
  calcOutput("MAgPIEReport", subtype ="CostMAC",                      round=8,  file="p26_macCostLuLookup.cs4r")
  calcOutput("FossilPolyRent",                        round=8,  file="f31_ffPolyRent.cs4r")
  calcOutput("FossilPolyCumEx",                       round=8,  file="f31_ffPolyCumEx.cs4r")
  calcOutput("FossilExtraction", subtype="FossilExtraction", round=9,  file="f31_ffPolyCoeffs.cs3r")
  calcOutput("FossilExtraction", subtype="UraniumExtraction",round=5,  file="f31_costExPoly.cs3r")
  calcOutput("RLDCCoefficients", subtype = "LoB",            round=6,  file="f32_RLDC_Coeff_LoB.cs3r")
  calcOutput("RLDCCoefficients", subtype = "Peak",           round=6,  file="f32_RLDC_Coeff_Peak.cs3r")
  calcOutput("EarlyRetirementAdjFactor",                               file="p_earlyRetirementAdjFactor.cs3r")
  calcOutput("DiffInvestCosts",  subtype="Invest_Costs",     round=4,  file="p_inco0.cs4r")
  calcOutput("DiffInvestCosts",  subtype="Efficiency",       round=4,  file="pm_eff.cs4r")
  calcOutput("CapacityFactorHist",                           round=4,  file="p_histCapFac.cs4r")
  
  #---------------policy parameters--------------------------------------------------------------------
  calcOutput("EmiTarget", subtype="share_cond" ,          round=4,  file="p45_2005share_target.cs4r") 
  calcOutput("EmiTarget", subtype="multiplier_cond",      round=4,  file="p45_factor_targetyear.cs4r") 
  calcOutput("EmiTarget", subtype="ghg_cond",             round=4,  file="p45_hist_share.cs4r") 
  calcOutput("CapTarget", sources="Rogelj2017", condition="Conditional", round = 4, file = "f40_NDC.cs4r")
  calcOutput("CapTarget", sources="REN21",                round=4,  file="f40_REN21.cs4r")
  calcOutput("CapTarget", sources="NDC+REN21+CHN_NUC",    round=4,  file="f40_NDC+REN21+CHN_NUC.cs4r")
  calcOutput("SharedTarget", subtype="FErenewablesShare", round=3,  file="f40_FE_RenShare.cs4r")
  calcOutput("EffortSharingTarget",                       round=3,  file="p47_EStarget.cs4r")
  calcOutput("EffortSharingRefEmi", subtype="REMIND_GHG", round=6,  file="p47_ES_GHG_referenceEmissions.cs4r")
  calcOutput("EffortSharingRefEmi", subtype="REMIND_CO2", round=6,  file="p47_ES_CO2_referenceEmissions.cs4r")
  
  #-------------- historical data ---------------------------------------------------------------------
  calcOutput("Historical",                            round=5,  file="historical.mif", aggregate="REG+GLO")
  
  #--------------- EDGE Transport ---------------------------------------------------------------------
  sapply(c("value_time", "harmonized_intensities", "price_nonmot",
           "SW", "pref", "UCD_NEC_iso", "logit_exponent"),
         function(stype){
           print(sprintf("Loading %s", stype))
           suppressWarnings(calcOutput("EDGETransport", subtype=stype,
                                       file=paste0(stype, ".cs4r"), aggregate=F))
         })
  
  sapply(c("fe_demand_tech", "fe2es", "esCapCost", "pm_trp_demand"),
         function(stype){
          print(sprintf("Loading %s", stype))
           suppressWarnings(calcOutput("EDGETransport", subtype=stype, round=8,
                                       file=paste0(stype, ".cs4r"), aggregate=T))
         })
  
}