# load and prepare data
data(winddata)
set1 <- createSet(height=40, v.avg=winddata[,2], v.std=winddata[,5],
dir.avg=winddata[,14])
ts <- formatTS(winddata[,1])
neubuerg <- createMast(time.stamp=ts, loc=NULL, desc=NULL, set1)
neubuerg <- clean(neubuerg)
# calculate turbulence intensity
neubuerg.ti <- turbulence(neubuerg, 1)
# plot turbulence intensity
plotTurbulence(neubuerg.ti)
# change colour and text size
plotTurbulence(neubuerg.ti, col="darkgray", cex=1.3)
plotTurbulence(neubuerg.ti, col=3, cex=0.9)
Run the code above in your browser using DataLab