Xinxin Yang's avatar
added data_map data
Xinxin Yang committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 
## script: xxx
# nuts: Nomenclature of Territorial Units for Statistics
rm(list = ls())
library(ggtext)
library(rlang)
setwd("D:/public/yang/CAPRIR_Project/caprir")

map_data <- prep_mapdata()
map_data_use <- map_data %>% 
  filter(!grepl("TR.*", CAPRI_NUTS_ID)) %>% # TR: Tuerkiye
  filter(!grepl("UK.*", CAPRI_NUTS_ID))  %>%# UK:
  filter(!grepl("PT20.*", CAPRI_NUTS_ID)) %>% # PT20/PT30
  filter(!grepl("PT30.*", CAPRI_NUTS_ID)) %>%# PT20/PT30
  filter(!grepl("ES70.*", CAPRI_NUTS_ID))  # ES70


file1 <- "res_2_1230cap_after_2014_refdefaultA" # refdefautA
file2 <- "res_2_1230fta_import_bandefaultA" # bandefaultA
setwd("D:/public/yang/CAPRIR_Project/caprir_yang")
gdx.dir <- paste0(getwd(), "/gdx/")

# CERE RAPE SOYA PULS
prod1 = c("CERE","RAPE","SOYA","PULS")

# MAIF OFAR GRAI GRAE
prod2 = c("MAIF", "OFAR","GRAI","GRAE")
 

# products %>% filter(key %in% prod_list)
outdata.dir <- paste0(getwd(), "/outdata.dir")



map_prod_data <- function(benchmark,scenario,prod, saveName){
  
  benchmark <- capri_data(filename = paste0(gdx.dir,benchmark), 
                          selregion = "all", 
                          seldim5 = "", 
                          selcols = prod,
                          selrows = "LEVL" , 
                          simyear = "2030", 
                          scenarioname = "benchmark")
  scenario <- capri_data(filename = paste0(gdx.dir,scenario), 
                         selregion = "all", 
                         seldim5 = "", 
                         selcols = prod,
                         selrows = "LEVL" , 
                         simyear = "2030", 
                         scenarioname = "scenario")
  
  save(benchmark, scenario, file = paste0(outdata.dir,"/",saveName, ".rda"))
  # map_capri(baseline = benchmark, scenario = scenario, comparison = FALSE)
  # map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = FALSE)
  # map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = TRUE)
}

map_prod_data(benchmark = file1, scenario = file2,prod = prod1,saveName = "prod1")
map_prod_data(benchmark = file1, scenario = file2,prod = prod2,saveName = "prod2")

t1 = Sys.time()
rm(benchmark, scenario)
load(file =  paste0(outdata.dir,"/","prod1", ".rda"))
p1 <- map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = TRUE, 11)
p1
ggsave(paste0( outdata.dir,"/cereals.png"), plot = p1,width = 16, height = 9, dpi = 100)
print(Sys.time()-t1)

rm(benchmark, scenario)
load(file =  paste0(outdata.dir,"/","prod2", ".rda"))
map_capri(baseline = benchmark, scenario = scenario, comparison = TRUE, percent_change = TRUE)