Interface to Climatic Research Unit Time-Series Version 3.21 Data

Functions for reading in and manipulating CRU TS3.21: Climatic Research Unit (CRU) Time-Series (TS) Version 3.21 data.

cruts

The gridded Climatic Research Unit (CRU) TS (time-series) 3.21 datasets are month-by-month variations in climate over the period 1901-2012, on high-resolution (0.5 x 0.5 degree) grids, produced by the Climatic Research Unit (CRU) at the University of East Anglia. They are available from here:

http://catalogue.ceda.ac.uk/uuid/ac4ecbd554d0dd52a9b575d9666dc42d

This is an R package for reading in and manipulating CRU TS3.21: Climatic Research Unit (CRU) Time-Series (TS) Version 3.21 data. A test script is below

• First load all required packages

  library(cruts);library(sp);library(raster);library(stringr);library(lubridate);library(ncdf);library(rgdal)

• this dataset contains a SpatialPolygonsDataFram polygon 'ama', I want to crop the CRU TS dataset to this polygon

  load("mydir/polygon_data.RData")

• the unzipped NetCDF file containing the data, I have only tested the package on this file

  fln <- "mydir/cru_ts3.21.1901.2012.tmp.dat.nc"

• this function returns a raster brick with the raw data

  br <- cruts2raster(ncfile=fln,timeRange=c("2004-01-01","2013-12-31"),poly=poly,offset=offset,type="brick")

• this function returns your polygon with the raw data averaged over each region

  temp <- cruts2poly(ncfile=fln,poly=ama,timeRange=c("2004-01-01","2013-12-31"),na.rm=TRUE)

• this function returns your polygon with anomalies either as a raster brick or if a SpatialPolygonsDataFrame is specified, the data averaged over each region

  anom <- getAnomaly(ncfile=fln,poly=ama,timeRange=c("2004-01-01","2013-12-31"),na.rm=TRUE)