Final dataset for the globe and for EU28. Also cifos crop names are added....

Final dataset for the globe and for EU28. Also cifos crop names are added. Issues are the not available disaggregation and the impication on nutrient intake (Rest, tnut, Veggie). Which nutrient content should be taken as a proxy?

Suitability_SPAM/Suitability_SPAM.R

@@ -80,13 +80,36 @@ suit_allvars = suit_agg %>%
 # Mapping -----------------------------------------------------------------
 # CIFOS crop 
 
+# Distinct crops from suitability files 
+dist_suit_crops = suit_allvars %>% 
+  distinct_at(vars(crops), .keep_all = FALSE) %>% 
+  drop_na() %>% write_csv(., "Input_data/distinct_suit_crops.csv")
+# count(dist_suit_crops)
 
-# country names 
+suit_crop = CropList_SPAM_FAO_CIFOS_distinct %>%  
+  dplyr::select(cifos_crops,name_suitability) %>%
+  dplyr::rename(suit_name = name_suitability) %>% 
+  dplyr::right_join(., suit_allvars, by = c("suit_name" = "crops"))  
 
+# country names
 
+# Mapping file
+country_map = read.csv(here::here("Input_data", "Mappings", "dat_map_SPAM_GAUL_CIFOS_FAO.csv")) %>%  
+  as_tibble() %>% 
+  dplyr::select(ADM0_CODE_GAUL, iso3_SPAM, name_cntr_CIFOS2021)
 
-# EU27 - subsetting -------------------------------------------------------
+# joining adm0 code with cifos country names
+suit_country_global = country_map %>%  
+  dplyr::right_join(., suit_crop, by = c("ADM0_CODE_GAUL" = "ADM0_CODE"))
 
+# write_csv(suit_country_global, here::here("Input_data","suit_complete_global.csv"))
 
 
+# EU28 - subsetting -------------------------------------------------------
+EU_iso3 = c("AUT","BEL", "BGR", "HRV", "CZE", "DNK", "EST","FIN", "FRA", "DEU", "GRC", "HUN", "IRL", "ITA","LVA", 
+            "LTU", "LUX", "MLT", "NLD", "POL", "PRT","ROU", "SVK","SVN", "ESP", "SWE", "GBR", "CYP")
+
+# Subset EU28 - ADding yields, renaming, mapping names
+suit_EU28 = suit_country_global %>% filter(., iso3_SPAM  %in% EU_iso3)
+# write_csv(suit_EU28, "Input_data/suitability_data_EU28.csv")