# Figure 2 code for
# Inconsistent definitions of GDP: Implications for estimates of decoupling
# created by Gregor Semieniuk, October 2022

# Article: https://doi.org/10.1016/j.ecolecon.2023.108000
# Zenodo repository: hhttps://doi.org/10.5281/zenodo.8381368

library(tidyverse)
library(viridis)

b <- read_csv("Zenodo/2.2_fig2-sourcedata.csv")

pdf("Zenodo/2.3_fig2_plot.pdf",4.5,3.5)
par(mar=c(3.5,3.5,0.3,0.3))
plot(b$yr,b$ey70,type="o",pch=19,frame.plot=F,ann=F,ylim=c(min(b$ey18),100),lwd=3,col=viridis(3)[1])
mtext(side = 1, line = 2.5, "Year")
mtext(side = 2, line = 2.5, "Index, 1930 = 100")
lines(b$yr,b$ey18,type="o",pch=17,col=viridis(3)[2],lwd=3,cex=1.23)
legend("topright",legend = c("US Energy intensity: 1984 GDP","US Energy intensity: 2021 GDP"),bty="n",lty=1,pch=c(19,17),col=viridis(3)[1:2],lwd=3,pt.cex = c(1,1.3))
dev.off()

write_csv(b,"Zenodo/2.4_fig2_datapoints.csv")