# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2], dir.avg=winddata[,14])
set30 <- createSet(height=30, v.avg=winddata[,6], dir.avg=winddata[,16])
set20 <- createSet(height=20, v.avg=winddata[,10])
ts <- formatTS(time.stamp=winddata[,1])
neubuerg <- createMast(time.stamp=ts, set40, set30, set20)
neubuerg <- clean(mast=neubuerg)
# create profile
neubuerg.wp <- profile(mast=neubuerg, v.set=c(1,2), dir.set=1,
print=FALSE)
# plot all profiles
plotProfile(profile=neubuerg.wp)
# plot only one sector
row.names(neubuerg.wp$profile) # available sectors
plotProfile(profile=neubuerg.wp, sector=3) # ENE by sector number
plotProfile(profile=neubuerg.wp, sector="ene") # ENE by sector code
plotProfile(profile=neubuerg.wp, sector="all") # general profile
# omit 'measured' points
plotProfile(profile=neubuerg.wp, measured=FALSE)
# customize plot
plotProfile(profile=neubuerg.wp, bty="l", bty.leg="o", cex.axis=0.8,
cex.lab=0.9, cex.leg=0.7, col=rainbow(13), col.axis=gray(0.2),
col.box=gray(0.2), col.lab=gray(0.2), col.leg=gray(0.2),
col.ticks=gray(0.2), las=0, lty=c(rep(3,12),1),
lwd=c(rep(1.2,12), 1.7), mar=c(3,3,0.5,0.5), mgp=c(2,0.7,0),
pos.leg="right", xlab="velocity [m/s]", ylab="height [m]",
xlim=c(0,11), ylim=c(0,150), x.intersp=1, y.intersp=1.2)
Run the code above in your browser using DataLab