Multi-annual predictions of hot, dry and hot-dry compound extremes

The zip file contains several folders with the script used for the computation of hot (tx90p) dry (SPI and SPEI) and hot-dry compound (HDSPI and HDSPEI) extremes. There is a different script for the different dataset used, specifically the decadal hindcasts (dcppA) the observations and the reanalysis datasets (obs-rea) and the historical simulations (hist). The files are in R(.R), python(.py) and bash script(.sh) format.

For TX90p:

1. We run the python script tx90p_(nameofdataset).py.
2. We run the script yearsum_seamask_regrid_(nameofdataset)_tx.sh.

For SPI and SPEI:

1. We run the R script Script_Accum_pr-wb_(nameofdataset).R
2. We run the R script Accum_to_SPI-SPEI_(nameofdataset).R
3. To compute the dailycopies of the monthly SPI and SPEI values, which will be used in the detection of the compound extremes, we run the script dailycopies_SPI-SPEI_obs-rea.sh for the observation and reanalysis and dailycopies_mod_SPI-SPEI_(nameofdataset).sh for the models.
4. We run the script yearsum_seamask_regrid_(nameofdataset)_droughts.sh.

For HDSPI and HDSPEI:

1. We run the python script HDSPI-HDSPEI_(nameofdataset).py.
2. We run the script yearsum_remask_regrid_(nameofdataset)_HD.sh.