Simple Precpitation Intensity Index

Description

This function computes the climdex index SDII.

Usage

climdex.sdii(ci)

Arguments

Details

climdex.sdii computes the climdex index SDII, or Simple Precipitation Intensity Index. This is defined as the sum of precipitation in wet days (days with preciptitation over 1mm) during the year divided by the number of wet days in the year.

Value

A vector containing the value of the index for each year.

Note

fclimdex rounds to 1 decimal place, whereas climdex.sdii does not. This results in some small differences.

References

http://etccdi.pacificclimate.org/list_27_indices.shtml

See Also

climdexInput.raw, climdexInput.csv.

Examples

library(PCICt)

## Create a climdexInput object from some data already loaded in and
## ready to go.

## Parse the dates into PCICt.
tmax.dates <- as.PCICt(do.call(paste, ec.1018935.tmax[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
tmin.dates <- as.PCICt(do.call(paste, ec.1018935.tmin[,c("year",
"jday")]), format="%Y %j", cal="gregorian")
prec.dates <- as.PCICt(do.call(paste, ec.1018935.prec[,c("year",
"jday")]), format="%Y %j", cal="gregorian")

## Load the data in.
ci <- climdexInput.raw(ec.1018935.tmax$MAX_TEMP,
ec.1018935.tmin$MIN_TEMP, ec.1018935.prec$ONE_DAY_PRECIPITATION,
tmax.dates, tmin.dates, prec.dates, base.range=c(1971, 2000))

## Create a timeseries of annual SDII values.
sdii <- climdex.sdii(ci)