# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
set2 <- createSet(height=30, v.avg=winddata[,6], dir.avg=winddata[,16])
ts <- formatTS(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL, set1, set2)
neubuerg <- clean(neubuerg)
# calculate AEP
neubuerg.wp <- profile(neubuerg, v.set=c(1, 2), dir.set=1)
pw.56 <- readPC(system.file(package="bReeze", "powercurves",
"PowerWind_56_900kW.wtg"))
neubuerg.aep <- aep(neubuerg.wp, pw.56, 71)
# plot AEP
plotAep(neubuerg.aep)
plotAep(neubuerg.aep, show.total=FALSE) # omit total AEP
# change colour and text size
plotAep(neubuerg.aep, col=gray(5:0 / 5), 0.8)
plotAep(neubuerg.aep, col=rainbow(5), cex=1.1)
Run the code above in your browser using DataLab