diff --git a/analysis/figures/figure1.pdf b/analysis/figures/figure1.pdf
index fb5fd7355a6afadd80a143d80f403427c4a9c07b..8e1509202fdb71af0e6511d84a6e8ce808782146 100644
Binary files a/analysis/figures/figure1.pdf and b/analysis/figures/figure1.pdf differ
diff --git a/analysis/figures/figure2.pdf b/analysis/figures/figure2.pdf
index 4f0f5c1b3840aa9c06b7e32023ad6187a626484a..e54797b3500c297b0be1ff4b576be0266fd8de9f 100644
Binary files a/analysis/figures/figure2.pdf and b/analysis/figures/figure2.pdf differ
diff --git a/analysis/figures/figure3.pdf b/analysis/figures/figure3.pdf
index 6290bcf36e609833790ef7aa01611aea1afe2996..eb8c64b4683b46acc2febc37c7506fb53e93cbf4 100644
Binary files a/analysis/figures/figure3.pdf and b/analysis/figures/figure3.pdf differ
diff --git a/analysis/figures/figure4.pdf b/analysis/figures/figure4.pdf
index 2d0ab573300965ba5e4b0dc09c444c50087804ba..1103445308200ace19a1c210d10b5caa034b3e50 100644
Binary files a/analysis/figures/figure4.pdf and b/analysis/figures/figure4.pdf differ
diff --git a/analysis/figures/figure5.pdf b/analysis/figures/figure5.pdf
index 14d9701c8881fddbceaf9eb1ae40994f6f2ed3cc..b6e817129b16d7f4a3846f5ecb385e808867d646 100644
Binary files a/analysis/figures/figure5.pdf and b/analysis/figures/figure5.pdf differ
diff --git a/analysis/paper/paper.Rmd b/analysis/paper/paper.Rmd
index 6cb70efd4bac3e16e2a4984a9631fa5b623b54f9..ec9b628624e9ffcc92180dfaca1300a5968c974f 100644
--- a/analysis/paper/paper.Rmd
+++ b/analysis/paper/paper.Rmd
@@ -791,6 +791,35 @@ Our results show that in 2015, higher-income people in higher-income countries h
 
 ```{r}
 
+pdat_int_energy = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
+  ungroup() %>%
+  filter(year==2015) %>%
+  left_join(read_csv(here("analysis/data/derived/sectors_agg_method1_ixi.csv")), 
+            by="sector_agg_id") %>% 
+  mutate(intensity_energy = (total_energy_use_tj)/(total_fd_me)) %>%
+  select(five_sectors, eu_q_rank, 
+         intensity_energy) %>%
+  filter(eu_q_rank == 10) %>%
+  select(-eu_q_rank)
+
+pdat_final_demand = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
+  ungroup() %>%
+  filter(year==2015) %>%
+  left_join(read_csv(here("analysis/data/derived/sectors_agg_method1_ixi.csv")), 
+            by="sector_agg_id") %>%
+  left_join(pdat_int_energy, by="five_sectors") %>%
+  mutate(total_energy_use_tj_new = (total_fd_me)*intensity_energy) %>%
+  mutate(total_energy_use_tj_diff = total_energy_use_tj-total_energy_use_tj_new) %>%
+  select(eu_q_rank,total_energy_use_tj_new) %>%
+  group_by(eu_q_rank) %>%
+  summarise(total_energy_use_tj_new = sum(total_energy_use_tj_new)) %>%
+  mutate(total_blah = total_energy_use_tj_new/33417583,
+         total_blah = total_blah*1000)
+
+ineq_curr_check = pdat_final_demand$total_blah[10]/pdat_final_demand$total_blah[1]
+
+
+
 df_energy_deciles = read_csv(here("analysis/data/raw/energy_figure_data.csv")) 
 
 ineq_curr = df_energy_deciles$pae_energy_use_gj[10]/df_energy_deciles$pae_energy_use_gj[1]