40-degrees.R

needs(tidyverse, jsonlite, zoo)

url <- 'https://interactive-atlas.ipcc.ch/temporal/serie?ids=19%2C17%2C16&region=AR6&project=EUR-11&scenario=%scenario%&variable=tx40&valueType=VALUE&period=near&temporalFilter=year&baseline=&landMask=false&mountainMask=false'

readModel <- function(scenario) {
  read_json(str_replace(url, '%scenario%', scenario), simplifyVector = T) %>% 
    group_by(year) %>% 
    summarise(
        p10=quantile(value, probs=c(0.1), na.rm=T),
        p25=quantile(value, probs=c(0.25), na.rm=T),
        p50=median(value, na.rm=T),
        p75=quantile(value, probs=c(0.75), na.rm=T),
        p90=quantile(value, probs=c(0.9), na.rm=T)
    ) %>% 
    mutate(
      !!paste0(scenario,'.p10'):=rollapply(p10, 3, mean, align='center',partial=T, fill=NA),
      !!paste0(scenario,'.p25'):=rollapply(p25, 3, mean, align='center',partial=T, fill=NA),
      !!paste0(scenario,'.p50'):=rollapply(p50, 3, mean, align='center',partial=T, fill=NA),
      !!paste0(scenario,'.p75'):=rollapply(p75, 3, mean, align='center',partial=T, fill=NA),
      !!paste0(scenario,'.p90'):=rollapply(p90, 3, mean, align='center',partial=T, fill=NA),
    ) %>% 
    select(-p10, -p25, -p50, -p75, -p90)
       
}

# 
models.hist <- readModel('historical')
models.rcp26 <- readModel('rcp26') %>% filter(year>=2015)
models.rcp45 <- readModel('rcp45') %>% filter(year>=2015)
models.rcp85 <- readModel('rcp85') %>% filter(year>=2015)


models.hist %>%
  full_join(models.rcp26, copy = T, by=c('year'='year')) %>% 
  full_join(models.rcp45, copy = T, by=c('year'='year')) %>% 
  full_join(models.rcp85, copy = T, by=c('year'='year')) %>% write_csv('data.csv')
  

data.csv