File:Winter-NAO-Index.svg

# Build multi languages plots for annual and winter NAO
# based on Hurrell data.

library(tidyverse)
library(glue)

# get data : winter and annual
# add sign column for colors. remove the last row (void)
nao_djfm_hurrell <- "https://climatedataguide.ucar.edu/sites/default/files/nao_station_djfm.txt" %>% 
  read_table(skip = 1, col_types = "id", na = "-999.0", col_names = c("year", "nao")) %>% 
  mutate(sign = if_else(nao < 0, "negative", "positive")) %>% 
  drop_na()

nao_hurrell <- "https://climatedataguide.ucar.edu/sites/default/files/nao_station_annual.txt" %>% 
  read_table(skip = 1, col_types = "id", na = "-999.0", col_names = c("year", "nao")) %>% 
  mutate(sign = if_else(nao < 0, "negative", "positive")) %>% 
  drop_na()


# manage languages
language <- list(
  fr_FR = list(
    winter = list(
      data = nao_djfm_hurrell,
      title = "Indice hivernal de l'oscillation nord-atlantique (ONA)",
      subtitle = "Lisbonne - Stykkishólmur/Reykjavik, décembre à mars",
      caption = "Données : Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nMise à jour régulière. Accédé le",
      x = "Année",
      y = "Différence de pression normalisée (hPa)",
      outDec = ","
    ),
    annual = list(
      data = nao_hurrell,
      title = "Indice annuel de l'oscillation nord-atlantique (ONA)",
      subtitle = "Lisbonne - Stykkishólmur/Reykjavik",
      caption = "Données : Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nMise à jour régulière. Accédé le",
      x = "Année",
      y = "Différence de pression normalisée (hPa)",
      outDec = ","
    )
  ),
  en_US = list(
    winter = list(
      data = nao_djfm_hurrell,            
      title = "North Atlantic Oscillation (NAO) winter index",
      subtitle = "Lisbon - Stykkishólmur/Reykjavík, December to March",
      caption = "NAO Index Data provided by the Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nUpdated regularly. Accessed",
      x = "Year",
      y = "Difference of normalized sea level pressure (hPa)",
      outDec = "."
    ),
    annual = list(
      data = nao_hurrell,
      title = "North Atlantic Oscillation (NAO) annual index",
      subtitle = "Lisbon - Stykkishólmur/Reykjavík",
      caption = "NAO Index Data provided by the Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nUpdated regularly. Accessed",
      x = "Year",
      y = "Difference of normalized sea level pressure (hPa)",
      outDec = "."
    )
  ),
  de_DE = list(
    winter = list(
      data = nao_djfm_hurrell,
      title = "Nordatlantischen Oszillation (NAO) Winter Index",
      subtitle = "Lisbon - Stykkishólmur/Reykjavík, Dezember bis März",
      caption = "Datei : Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nRegelmäßig aktualisiert. Zugänglich am",
      x = "Jahre",
      y = "Differenz der standardisierten Luftdruck (hPa)",
      outDec = ","
    ),
    annual = list(
      data = nao_hurrell,      
      title = "Nordatlantischen Oszillation (NAO) Index",
      subtitle = "Lisbon - Stykkishólmur/Reykjavík",
      caption = "Datei : Climate Analysis Section, NCAR, Boulder, USA, Hurrell (2003)\nRegelmäßig aktualisiert. Zugänglich am",
      x = "Jahre",
      y = "Differenz der standardisierten Luftdruck (hPa)",
      outDec = ","
    )
  )
)

theme_set(theme_bw())
theme_update(plot.caption = element_text(size = 7))
oldDec <- getOption("OutDec")

for (l in names(language)) {
  message(l)
  
  for (t in names(language[[l]])) {
    message(t)
    current <- language[[l]][[t]]
    options(OutDec = current$outDec)
    
    # plot graph
    ggplot(current$data, aes(year, nao)) +
      geom_col(aes(fill = sign)) +
      geom_smooth(span = .15, color = "black", alpha = 0.3) +
      scale_fill_manual(values = c("positive" = "darkorange2",
                                   "negative" = "deepskyblue3")) +
      scale_x_continuous(breaks = seq(1860, max(current$data$year), 20)) +
      guides(fill = FALSE) +
      labs(title = current$title,
           subtitle = current$subtitle,
           caption = glue("{current$caption} {format(Sys.Date(), '%Y-%m-%d')}"),
           x = current$x,
           y = current$y)
    
    ggsave(file = glue("nao_hurrell_{t}_{l}_{Sys.Date()}.svg"), 
           width = 20,
           height = 12.4,
           units = "cm",
           scale = 0.8,
           device = svg)
  }
}

options(OutDec = oldDec)