# load example data
data(winddata)
# create two datasets
set40 <- createSet(height=40, v.avg=winddata[,2], v.std=winddata[,5],
dir.avg=winddata[,14])
set30 <- createSet(height=30, v.avg=winddata[,6], v.std=winddata[,9],
dir.avg=winddata[,16])
# format time stamp
ts <- formatTS(winddata[,1])
# create met mast object
metmast <- createMast(time.stamp=ts, loc=NULL,
desc=NULL, set40=set40, set30=set30)
# plot time series of met mast signals
plotTimeSeries(metmast)
# calculate availability of pairs of wind speed and direction
availability(metmast)
# calculate monthly means of wind speed
monthStats(metmast)
# calculate directional frequency and mean wind speed
frequency(metmast, 1)
# calculate turbulence intensity
turbulence(metmast, 1)
# calculate weibull parameters
wb <- weibull(metmast, 1)
wb
# calculate total wind energy content
energy(wb)
# calculate wind profile
pf <- profile(metmast, v.set=c(1,2), dir.set=1, num.sectors=12)
pf
# import power curve
pc <- readPC(system.file(package="bReeze", "powercurves",
"Enercon_E126_7.5MW.pow"))
pc
# calculate annual energy production
aep <- aep(pf, pc, 135)
aep
# plot AEP
plotAep(aep)
Run the code above in your browser using DataLab