IDFcurve2 | R Documentation |
Estimate IDF curve fitting a [mm/h], m ,n, b[h] parameters of maximum annual rainfall series
IDFcurve2(rain, tc, stvalue1 = 1, stvalue2 = 1, t, Tr = 200, MP = F, Trplot = F)
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 |
Estimate parameters of Intensity Duration Frequency curves for maximum annual rainfall series of different duration
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] |
a(Tr) is defined by Gumbel distribution.
Corrado Tallerini
IDFcurve
## 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
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.