# load and prepare data
data(winddata)
set40 <- createSet(height=40, v.avg=winddata[,2], v.max=winddata[,3],
v.min=winddata[,4], v.std=winddata[,5], dir.avg=winddata[,14])
set30 <- createSet(height=30, v.avg=winddata[,6], v.std=winddata[,9],
dir.avg=winddata[,16])
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 time series
plotTimeSeries(mast=neubuerg)
plotTimeSeries(mast=neubuerg, set=1, legend=FALSE)
plotTimeSeries(mast=neubuerg, set="set1", legend=FALSE) # same as above
# change signals
plotTimeSeries(mast=neubuerg, signal=c("v.avg", "v.min", "v.max",
"v.std"))
plotTimeSeries(mast=neubuerg, signal=c("turb.int"))
plotTimeSeries(mast=neubuerg, set=1, signal=c("v.avg", "dir.avg"))
# change time scale
plotTimeSeries(mast=neubuerg, subset=c("2010-01-01 00:10:00", NA))
plotTimeSeries(mast=neubuerg, subset=c("2009-10-11 00:10:00",
"2009-10-11 23:50:00"))
plotTimeSeries(mast=neubuerg, set=1, signal="v.max",
subset=c(NA, "2009-12-27 18:30:00"))
# customize plot
plotTimeSeries(mast=neubuerg, bty="n", bty.leg="o", cex.axis=1.2,
cex.lab=1.4, cex.leg=1.2, col=c("darkblue", "red2", "darkgreen"),
col.axis="darkgray", col.lab="darkgray", col.leg="darkgray",
col.ticks="darkgray", las=0, lty=rep(1.5,3), mar=c(0.5,4,0,0.5),
mgp=c(2,0.5,0), ylab=c("v [m/s]","dir [deg]","ti [-]"), x.intersp=1)
Run the code above in your browser using DataLab