globtemp.R

Function for plotting updated monthly or yearly global temperature anomalies from NASA/GISS GISTEMPv4, Copernicus ERA5, and/or Met Office Hadley Centre HadCRUT4, relative to a desired reference period.

Prerequisites

R packages ggplot2 and lubridate (included in tidyverse)

Supporting functions

• read_era5.R, read_hadcrut4.R, and read_nasa.R: functions for retrieving/reading the most resent data.
• anom2anom.R: function for changing the reference period for timeseries of anomalies.
• distr_col.R: function for distribution of colors blue/red/black for anomalies under/over/equal to desired reference period.
• plot_singular_points.R, plot_singular_bars.R, plot_plural_lines.R: functions for plotting
• top5.R: function for printing the top 5 years with the highest temperatures

User input

Dataset(s) to plot (e.g. c('NASA') or c('NASA','ERA5','HadCRUT')), start of desired reference period (e.g. 1961), end of desired reference period (e.g. 1990), period of the year (e.g. 'Yearly', or 'January', or 'February',...), plotting type ('bars' or 'points') and option (TRUE/FALSE, default = FALSE) to save the plot to a filename (default = NA).

The plotting type is only needed if one single dataset is used as input. Several datasets are always presented as solid lines.

Examples R

globtemp(c('NASA'), 1961, 1990, 'Yearly', 'points')

globtemp(c('NASA'), 1961, 1990, 'Yearly', 'bars')

globtemp(c('NASA','ERA5','HadCRUT'), 1981, 2010, 'September', save_option = TRUE, save_name = "example.png")

Contact

Email: Line Båserud ([email protected])