IDFcurve2 {LPM}

Intensity duration frequency curve for maximum annual rainfall series of different duration

Description

Estimate IDF curve fitting a [mm/h], m ,n, b[h] parameters of maximum annual rainfall series

Usage

IDFcurve2(rain, tc, stvalue1 = 1, stvalue2 = 1, t, Tr = 200, MP = F, Trplot = F)

Arguments

rain

Observed Maximum annual rainfall series [mm] of increasing duration

tc

Time of concentration of Basin [h] , maybe h(t) and i(t) duration must be calculated

stvalue1, stvalue2

Starting values of estimation algorithm. Deault stvalue1=1, stvalue2=1

t

observed rainfall series duration [h] example t=c(1,3,6,12,24) for durations 1,3,6,12,24 hours

Tr

Return period [y]. Default Tr=200

MP

logical: TRUE for 3 parameters formula i=a/(b+t)^m , FALSE for 2 parameters formula i=a*t^(n-1), Default MP=False

Trplot

logical: TRUE for plotting Tr values of a(Tr) parameter. Default Trplot=False

Details

Estimate parameters of Intensity Duration Frequency curves for maximum annual rainfall series of different duration

Value

par

List of estimated parameters: a(Tr), m, b, h(t) [mm], i(t) [mm/h], Offset of least squares optimizer

I

I(t) curve scattered point matrix [mm/h]

Curve

IDF curve scattered point matrix [mm/h]

Note

a(Tr) is defined by Gumbel distribution.

Author(s)

Corrado Tallerini

See Also

IDFcurve

Examples

##     data(milano)
##     ris=IDFcurve2(milano, 1, stvalue1 = 1, stvalue2 = 1,
##     t=c(0.25,0.5,0.75,1,1.25,1.5,2,2.5,3,4,6), Tr = 200, MP=F)
## --> 2 parameters IDF curve estimation of annual maximum rainfall
##     series recorded in Palazzo Marino - Milan (Italy)
##     ris=IDFcurve2(milano, 1, stvalue1 = 1, stvalue2 = 1,
##     t=c(0.25,0.5,0.75,1,1.25,1.5,2,2.5,3,4,6), Tr = 200, MP=T)
## --> 3 parameters IDF curve estimation of annual maximum rainfall
##     series recorded in Palazzo Marino - Milan (Italy)
## --> It's obvious the best performance of the 3 parameters formula