# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], v.std=winddata[,5],
dir.avg=winddata[,14])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set1)
neubuerg <- clean(mast=neubuerg)
# calculate turbulence intensity
neubuerg.ti <- turbulence(mast=neubuerg, turb.set=1, print=FALSE)
# plot turbulence intensity
plotTurbulence(turb=neubuerg.ti)
# change colour, text size etc.
plotTurbulence(turb=neubuerg.ti, cex.axis=0.7, cex.lab=0.9,
circles=c(0.05,0.20,0.05), col="lightgray", col.axis="darkgray",
col.border="gray", col.circle="darkgray", col.cross="darkgray",
col.lab="darkgray", fg=TRUE, lty.circle="dotdash",
lty.cross="longdash", lwd.border=1.2, lwd.circle=1.2,
lwd.cross=1.2, pos.axis=135, sec.space=0.6)
Run the code above in your browser using DataLab