Arguments

series of daily minimum temperature (data frame). Must include the following columns (and names): "year", "month", "day" and one or more station id(s), each column one Tn series.

daily_Tx

series of daily maximum temperature (data frame). Must include the following columns (and names): "year", "month", "day" and one or more station id(s), each column one Tx series.

daily_Tm

series of daily mean temperature (data frame). Must include the following columns (and names): "year", "month", "day" and one or more station id(s), each column one Tm series.

daily_P

series of daily precipitation (data frame). Must include the following columns (and names): "year", "month", "day" and one or more station id(s), each column one P series.

first.yr

first year of the period over which indices are calculated

last.yr

last year of the period over which indices are calculated

subs_missing

logical: if TRUE (default), missing values in input series are replaced by corresponding average values for each day in the series.

coeff_HI

Huglin's daylength correction, as long as the number of stations, or a single coeffient; in this case it is recycled on all stations. See details.

quant_Tn_rest

quantile (0..1) for the choice of the position of the representative year in the series, referred to the minimum temperature during "rest" period. Default is 0.5 (median)

quant_Tn_veg

quantile (0..1) for the choice of the position of the representative year in the series, referred to the minimum temperature during the "vegetative" period. Default is 0.5 (median)

quant_Tx_veg

quantile (0..1) for the choice of the position of the representative year in the series, referred to the maximum temperature during the "vegetative" period. Default is 0.5 (median)

indices

set of OIV indices to be listed. Default is all indices (1 to 10). See details for numbered list of indices.

Details

General info about OIV indices can be sought at http://www.oiv.int/oiv/info/enresolution2012?lang=en.
In general, if first.yr or last.yr are NULL (default), the lowest and highest values in series are taken as begin and end of calculation period.
A coherence check is done on input of start / end years.

If any input is missing, corresponding indices are not be calculated. If daily_Tm is missing and both daily_Tn and daily_Tx are passed to the function,
daily_Tm is calculated as the mean of daily_Tn and daily_Tx.

OIV indices are the following:

Tm_veg: 1 - mean temperature during vegetation period. The latter is April - October (N emisphere) or October - April (S emisphere). The case (N or S) is automatically recognised by temperature data.

HI: 4 - heliothermic Huglin index. A degree day summation of the average between mean and maximum temperature above 10 C, corrected by a coefficient of daylength duration.
The coefficient is given by the author in a table according to latitude. If one value is given, it is used for all stations. Default is 1.04 (lat. 44-46). Ref.: Huglin, P., 1978

FSI: 6 - Fregoni's simplified index. Given by the product between the summation of thermal range (Tx - Tn) and the number of days with Tm > 10 C, for 30 days before ripening.
The period before ripening is taken as September (N emisphere) or March (S emisphere). Ref.: Fregoni et Pezzutto, 2000.

BI: 7 - Branas' hydrothermic index. The only index not included in OIV list, used for fungine infection proneness.
It is given by the product of precipitation (mm) by the mean temperature (C) during the period April - October (N emisphere) or October - April (S emisphere). Ref: Eynard and Dal Masso, 1990.

Tn_rest: 8 - minimum temperature during rest period. Useful for assessing winter severity. The rest period is November - March (N emisphere) or May - September (S emisphere).
The case (N or S) is automatically recognised by temperature data.

quant_Tn_rest, quant_Tn_veg, quant_Tx_veg define the statistical rank of the year to be chosen as representative for assessing Tn_rest,
Tn_veg, and Tx_veg, respectfully. 0.5 (default) is the median year, 0 is the minimum (lowest temperature), 1 is the maximum (highest temperature).

The only missing index among those selected by OIV is Riou's Drought Index, which is calculated by function RDI on monthly series.

Value

A table (one line per station) reporting OIV bioclimatic indices. The Branas' index is added.