Energy intensity of FCEVs in DEU

[Loisel]

There seem to be quite unexpected efficiency improvements in FCEVs compared to BEVs in DEU (Units: MJ/vkm):

2050,DEU,Compact Car,BEV,1.03322798677885
2055,DEU,Compact Car,BEV,1.03094499484252
2060,DEU,Compact Car,BEV,1.02866200290619
2065,DEU,Compact Car,BEV,1.02637901096987
2070,DEU,Compact Car,BEV,1.02409601903354
2075,DEU,Compact Car,BEV,1.02181302709721
2080,DEU,Compact Car,BEV,1.01953003516088
2085,DEU,Compact Car,BEV,1.01724704322456
2090,DEU,Compact Car,BEV,1.01496405128823
2095,DEU,Compact Car,BEV,1.0126810593519
2100,DEU,Compact Car,BEV,1.01039806741557

vs

2050,DEU,Compact Car,FCEV,1.05547041155214
2055,DEU,Compact Car,FCEV,1.03549958136915
2060,DEU,Compact Car,FCEV,1.01552875118615
2065,DEU,Compact Car,FCEV,0.995557921003148
2070,DEU,Compact Car,FCEV,0.975587090820151
2075,DEU,Compact Car,FCEV,0.955616260637152
2080,DEU,Compact Car,FCEV,0.935645430454153
2085,DEU,Compact Car,FCEV,0.915674600271154
2090,DEU,Compact Car,FCEV,0.895703770088157
2095,DEU,Compact Car,FCEV,0.875732939905158
2100,DEU,Compact Car,FCEV,0.855762109722159

The script to extract these values is:

library(data.table)                                                                                                                                                                                                           
library(edgeTrpLib)                                                                                                                                                                                                           
                                                                                                                                                                                                                              
## use cached input data for speed purpose                                                                                                                                                                                    
setConfig(forcecache=T)                                                                                                                                                                                                       
                                                                                                                                                                                                                              
data_folder <- "/p/tmp/rottoli/ARIADNE_scenarios_EDGE_standalone/EDGE-T"                                                                                                                                                      
                                                                                                                                                                                                                              
inputdata <- loadInputData(data_folder)                                                                                                                                                                                       
                                                                                                                                                                                                                              
region = "DEU"                                                                                                                                                                                                                
loadFactor = inputdata$loadFactor[iso == region]                                                                                                                                                                              
int_dat = inputdata$int_dat[EJ_Mpkm_final > 0 & iso == region & year > 1975]                                                                                                                                                  
                                                                                                                                                                                                                              
dt = loadFactor[int_dat, on=.(year, iso, vehicle_type)]                                                                                                                                                                       
                                                                                                                                                                                                                              
unique(dt[is.na(loadFactor)]$vehicle_type)                                                                                                                                                                                    
                                                                                                                                                                                                                              
dt[, MJ_vkm := EJ_Mpkm_final * loadFactor * 1e6]                                                                                                                                                                              
                                                                                                                                                                                                                              
fwrite(dt[, .(year, iso, vehicle_type, technology, MJ_vkm)], "energy_intensity_edget.csv")

[ahagen-pik]

Energy intensity for BEV and FCEV look like this now in DEU (SSP2, Mix1)

Variable "Energy intensity|...|Compact Car|BEV" "Energy intensity|...|Compact Car|FCEV"
2005     "0.9324612"                            "1.8482122"                            
2010     "0.8695668"                            "1.8571164"                            
2015     "0.7812146"                            "1.8316520"                            
2020     "0.6800309"                            "1.7054838"                            
2025     "0.6370384"                            "1.5306518"                            
2030     "0.6153127"                            "1.3840411"                            
2035     "0.597979"                             "1.274497"                             
2040     "0.5835848"                            "1.2036779"                            
2045     "0.5692542"                            "1.1556837"                            
2050     "0.5577735"                            "1.1237190"                            
2055     "0.5517688"                            "1.1086940"                            
2060     "0.5504266"                            "1.1055069"                            
2070     "0.5504266"                            "1.1055069"                            
2080     "0.5504266"                            "1.1055069"                            
2090     "0.5504266"                            "1.1055069"                            
2100     "0.5504266"                            "1.1055069"                            
2110     "0.5504266"                            "1.1055069"                            
2130     "0.5504266"                            "1.1055069"                            
2150     "0.5504266"                            "1.1055069"                                                   

The script to reproduce these values is here:

library(data.table)

mifmix1 <- fread("path/to/Transport.mif", header = TRUE)
EIdeu <- mifmix1[Variable %like% "Energy intensity" & Variable %like% "Compact Car" & Region ==  "DEU"]
EIdeuBEV_FCEV <- t(EI[Variable %like% "BEV" |Variable %like% "FCEV"][, c(4, 6:25)])
EIdeuBEV_FCEV[1,1] <- str_replace(EIeurBEV_FCEV[1,1], pattern = "Transport\\|Pass\\|Road\\|LDV\\|Four Wheelers\\|Medium", replacement = "...")
EIdeuBEV_FCEV[1,2] <- str_replace(EIeurBEV_FCEV[1,2], pattern = "Transport\\|Pass\\|Road\\|LDV\\|Four Wheelers\\|Medium", replacement = "...")
EIdeuBEV_FCEV