From 2783afdf317602ef1196db80796577c5d7499b40 Mon Sep 17 00:00:00 2001
From: jaccard <jaccard@pik-potsdam.de>
Date: Wed, 3 Feb 2021 19:37:56 +0100
Subject: [PATCH] edit fig.5

---
 analysis/paper/paper.Rmd | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/analysis/paper/paper.Rmd b/analysis/paper/paper.Rmd
index 78ed6fb..8c0ed4c 100644
--- a/analysis/paper/paper.Rmd
+++ b/analysis/paper/paper.Rmd
@@ -885,8 +885,8 @@ pdat_final_demand = get_sector_summary_by_eu_ntile_direct(eu_q_count) %>%
   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(pae_energy_use_gj = total_energy_use_tj_new/33417583,
-         pae_energy_use_gj = pae_energy_use_gj*1000)
+  mutate(pae_energy_use_tj = total_energy_use_tj_new/33417583,
+         pae_energy_use_gj = pae_energy_use_tj*1000)
 
 df_energy_deciles = pdat_final_demand %>%
   select(eu_q_rank, pae_energy_use_gj)
@@ -958,7 +958,7 @@ c_mean_mer = round((mer[1]+mer[2])/2)
 ## run once to save file
 df_all = NULL
 for (min_energy in seq(from=mer[1], to=mer[2], by=0.25)) {
-  for (mean_energy in seq(from=mea[1], to=mea[2], by=1)) {
+  for (mean_energy in seq(from=mea[1], to=mea[2], by=0.25)) {
     if (min_energy <= mean_energy) {
       df_all = df_all %>%
       bind_rows(df_energy_deciles %>%
@@ -988,7 +988,7 @@ for (min_energy in seq(from=mer[1], to=mer[2], by=0.25)) {
 
 round_by = 10
 
-huh = readRDS(here("analysis/data/raw/scenarios_fine.rds"))
+#huh = readRDS(here("analysis/data/raw/scenarios_fine.rds"))
 
 df_all = df_all %>% #readRDS(here("analysis/data/raw/scenarios_fine.rds"))
   filter(eu_q_rank %in% c(1,10)) %>%
-- 
GitLab