# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2], v.std=winddata[,5])
set30 <- createSet(height=30, v.avg=winddata[,6], v.std=winddata[,9])
set20 <- createSet(height=20, v.avg=winddata[,10], v.std=winddata[,13])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set40, set30, set20)
neubuerg <- clean(mast=neubuerg)
# plot
plotTurbIEC(mast=neubuerg, set=1)
plotTurbIEC(mast=neubuerg, set="set1") # same as above
# customize plot
plotTurbIEC(mast=neubuerg, set=1, bty="l", cex.axis=0.8, cex.lab=0.9,
cex.leg=0.7, col.axis="darkblue", col.box="lightblue", col.lab=
"darkblue", col.leg="darkblue", col.ticks="darkblue", las=0,
leg.text=c("IEC class A", "IEC class B", "IEC class C", "measured"),
mar=c(3,3,0.5,0.5), mgp=c(1.8,0.5,0), pos.leg="top", xlab="v [m/s]",
ylab="ti [-]", xlim=c(0,25), ylim=c(0,0.5), x.intersp=1, y.intersp=1)
# customize bars
plotTurbIEC(mast=neubuerg, set=1, col="gray", border="black", space=0.6)
# customize lines
plotTurbIEC(mast=neubuerg, set=1, line=gray(1:3 / 10), lty=2:4,
lwd=0.5:2.5)
Run the code above in your browser using DataLab