inst/examples/test.generate.R

In RMRAINGEN: RMRAINGEN (R Multi-site RAINfall GENeretor): a package to generate daily time series of rainfall from monthly mean values

# TODO: Add comment
# 
# Author: ecor
###############################################################################


library(RMRAINGEN)

set.seed(125)
data(trentino)

year_min <- 1961
year_max <- 1990

period <- PRECIPITATION$year>=year_min & PRECIPITATION$year<=year_max
station <- names(PRECIPITATION)[!(names(PRECIPITATION) %in% c("day","month","year"))]
prec_mes <- PRECIPITATION[period,station]

## removing nonworking stations (e.g. time series with NA)
accepted <- array(TRUE,length(names(prec_mes)))
names(accepted) <- names(prec_mes)
for (it in names(prec_mes)) {
	accepted[it]  <- (length(which(!is.na(prec_mes[,it])))==length(prec_mes[,it]))
}

prec_mes <- prec_mes[,accepted]
## the dateset is reduced!!!
prec_mes <- prec_mes[,1:3]

coeff <- CoeffYWeq(data=prec_mes,p=3,tolerance=0.001)

generation <- generate(coeff,n=10,names=names(prec_mes))

origin <- paste(year_min,1,1,sep="-")


coeff_monthly <- CoeffYWeq(data=prec_mes,p=2,tolerance=0.001,sample="monthly",origin=origin)


generation_monthly <- generate(coeff_monthly,year_min=year_min,year_max=year_max,names=names(prec_mes),precipitation.indicator=TRUE)