# 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(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set1, set2)
neubuerg <- clean(mast=neubuerg)
# calculate AEP
neubuerg.wp <- profile(mast=neubuerg, v.set=c(1, 2), dir.set=1,
print=FALSE)
pw.56 <- readPC(file="PowerWind_56_900kW.wtg")
neubuerg.aep <- aep(profile=neubuerg.wp, pc=pw.56, hub.h=71, print=FALSE)
# plot AEP
plotAep(aep=neubuerg.aep)
# omit total AEP
plotAep(aep=neubuerg.aep, show.total=FALSE)
# change colours
plotAep(aep=neubuerg.aep, col=gray(5:0 / 5))
# change text sizes
plotAep(aep=neubuerg.aep, cex=0.8)
# manual definition of circles
plotAep(aep=neubuerg.aep, circles=c(250, 750, 250))
# plot sectors in foreground
plotAep(aep=neubuerg.aep, fg=TRUE)
# change position of axis labels
plotAep(aep=neubuerg.aep, pos.axis=135)
# no legend
plotAep(aep=neubuerg.aep, width.leg=0)
# freaky
plotAep(aep=neubuerg.aep, border.leg=heat.colors(5), bty.leg="o",
cex.axis=0.5, cex.lab=2, cex.leg=0.5, circles=c(80, 800, 80),
col=rainbow(5), col.axis="green", col.border="orange",
col.circle="purple", col.cross="yellow", col.lab="pink",
col.leg="lightblue", fg=TRUE, lwd.border=2, lwd.circle=3,
lwd.cross=4, lty.circle="12345678", lty.cross="87654321",
sec.space=0.6, title.leg="* WiNd SpEeD *", x.intersp=2, y.intersp=5)
Run the code above in your browser using DataLab