Code for reproducibility of the ATLAS boxplots and scatterplots of model projections#

´´´sh reproducibility/projections/boxplots_TandP.R ‘’’ is the star script of the ATLAS!!

It is flexible and easy to use. There is no special requirements for executing the script (only an R environment with a few packages), as it uses R functions and data stored in this repository.

The R function used is computeFigures (available at ‘’’sh ./datasets/aggregated-regionally/scripts/computeFigures.R’’’) and the data used are the regional means available at ‘’’sh ./datasets-aggregated-regionally/data/ ‘’’.

It allows for the selection of:

  • Seasons

  • Reference period

  • Surface (Land, sea or both)

  • Region/s from the ‘’’sh ./reference-regions/’’’ [Iturbide et al., 2020].

For example, the following parameter configuration,

# select seasons, use c(12,1,2) for winter
scatter.seasons <- list(c(12, 1, 2), 6:8)
# select reference period
ref.period <- 1995:2014
# select the area, i.e. "land", "sea" or "landsea"
area <- "land"
# Select reference regions.  Select the CORDEX domain to be considered
regions <- c("ECA", "EAS"); cordex.domain <- "EAS"

Will result in the following boxplots and scatterplots:

References#

IGutierrezA+20

Maialen Iturbide, José M. Gutiérrez, Lincoln M. Alves, Joaquín Bedia, Ruth Cerezo-Mota, Ezequiel Cimadevilla, Antonio S. Cofiño, Alejandro Di Luca, Sergio Henrique Faria, Irina V. Gorodetskaya, Mathias Hauser, Sixto Herrera, Kevin Hennessy, Helene T. Hewitt, Richard G. Jones, Svitlana Krakovska, Rodrigo Manzanas, Daniel Martínez-Castro, Gemma T. Narisma, Intan S. Nurhati, Izidine Pinto, Sonia I. Seneviratne, Bart van den Hurk, and Carolina S. Vera. An update of IPCC climate reference regions for subcontinental analysis of climate model data: definition and aggregated datasets. Earth System Science Data, 12(4):2959–2970, November 2020. doi:10.5194/essd-12-2959-2020.