/seatac-txx

Code for extreme temperature analysis at Seattle-Tacoma airport

Primary LanguageJupyter NotebookMIT LicenseMIT

This repository contains the details of the code, data processing steps, and software environment associated with the following paper:

Risbey JS, Irving DB, Squire DT, Matear RJ, Monselesan DP, Pook MJ, Ramesh N, Richardson D & Tozer CR (submitted). A large ensemble illustration of how record-shattering heat records can endure. Environmental Research: Climate.

Code, data processing steps, and software environment

The Makefile defines the rules/code/processing steps used to generate most of the resutls in the paper:

  • Figure 1: make plot-historgram CONFIG=workflow_config.mk
  • Figure 2: make plot-return-periods CONFIG=workflow_config.mk
  • Figure 3: make plot-hot-day CONFIG=workflow_config.mk (see find_hottest_model_day.ipynb for information about the hottest day)
  • Figure 4: make plot-z500-rmse CONFIG=workflow_config.mk
  • Figure 5: make plot-water-z500 CONFIG=workflow_config.mk
  • Figure 6: make plot-sample-size-dist CONFIG=workflow_config.mk
  • Figure 7: make plot-by-year CONFIG=workflow_config.mk

The Python scripts called by the Makefile rely on the unseen package: https://github.com/AusClimateService/unseen

The README file for the unseen package describes how to install that package and its dependencies.

Provenance

Each output image file has the command history embedded in the image metadata. It can be viewed by installing exiftool (e.g. $ conda install exiftool) and then running the following at the command line:

$ exiftool path/to/image.png

Data

The observations were daily maximum temperatures at Seattle Tacoma International Airport from the GHCNv2 station dataset, downloaded from the KNMI climate explorer.