curvad=function(nomefile,met1=T,met2=F)
{
out=list()
#Fare attenzion al percorso nelle istruzioni che seguono
#Attenzione: la serie in nomefile deve iniziare dal giorno 1 gennaio e puo' avere lunghezza arbitraria
#ma di almeno 5 anni

riverflow=scan(nomefile)
nanni=trunc(length(riverflow)/365)

if (nanni<5)
{
cat("Attenzione, il numero di anni di osservazione è inferiore a 5. Procedura interrotta", fill=T)
return()
}

duration=seq(1,365*nanni)
duration=duration/nanni

if(met1==T)
{
riverflow1=sort(riverflow,decreasing=T)
out$met1$duration=duration
out$met1$riverflow=riverflow1
}

if(met2==T)
{
duration=seq(1,365)
riverflow2=array(riverflow,dim=c(365,nanni))
for (i in 1:nanni)
riverflow2[,i]=sort(riverflow2[,i],decreasing=T)

avg=rep(0,365)
for (i in 1:365)
avg[i]=mean(riverflow2[i,])

out$met2$duration=duration
out$met2$riverflow=avg
}
return(out)
}