if (FALSE) {
# load example data
data("winddata", package="bReeze")
# create two datasets
set40 <- set(height=40, v.avg=winddata[,2], v.std=winddata[,5],
dir.avg=winddata[,14])
set30 <- set(height=30, v.avg=winddata[,6], v.std=winddata[,9],
dir.avg=winddata[,16])
# format time stamp
ts <- timestamp(timestamp=winddata[,1])
# create met mast object
metmast <- mast(timestamp=ts, set40=set40, set30=set30)
# plot time series of met mast signals
plot(metmast)
# calculate frequency and mean wind speed per wind direction sector
freq <- frequency(mast=metmast, v.set=1)
# plot frequency
plot(freq)
# calculate availability of pairs of wind speed and direction
availability(mast=metmast)
# calculate monthly means of wind speed
month.stats(mast=metmast)
# calculate turbulence intensity
turbulence(mast=metmast, turb.set=1)
# calculate weibull parameters
wb <- weibull(mast=metmast, v.set=1)
# calculate total wind energy content
energy(wb=wb)
# calculate wind profile
pf <- windprofile(mast=metmast, v.set=c(1,2), dir.set=1)
# import power curve
pc <- pc("Enercon_E126_7.5MW.pow")
# calculate annual energy production
aep <- aep(profile=pf, pc=pc, hub.h=135)
# plot AEP
plot(aep)
}
Run the code above in your browser using DataLab