Davis winds
Monday, September 6th, 2010It appears that the average wind speed in Davis over 20 years is quite stable. There is some variation across years, but each month has a characteristic profile. Right now is September and we can expect afternoon winds that start to taper off at about 5pm. November-January have windy nights, where the overnight wind speed stays higher compared to July-September. I took 20 years worth of hourly weather data from CIMIS (~175k observations) and distilled it down to this plot. Here averages are taken by month, year, and hour. The year and hour makes sense, but month is a bit arbitrary. Why should Aug 31 and Sept 1 be grouped in different months? Sept 1 is more similar to Aug 31 than it is to Sept 28th. Despite the arbitrary nature of month breaks, it is still quite interesting to see how stable the patterns are year to year, hour to hour.
Code
#historic Davis wind speed
#with CIMIS
#Greg Hirson, 2010
library(cimis)
library(lattice)
library(reshape)
library(ggplot2)
#years = 1990:2009
d = mcimisannual("006", 1990:2009, period = "hourly")
closeAllConnections()
#yearly average wind speed by hour and month
d1 = ddply(d, ~ time + months(datetime) + year,
.fun = function(df)mean(df$wind_speed, na.rm=T))
d1$time = d1$time/100
d1$months.datetime. = factor(d1$months.datetime.,
levels = unique(months(d$datetime)))
g1 = ggplot(d1, aes(x = time, y = V1, colour = factor(year)))
g1 + geom_line() + facet_wrap(~months.datetime.) +
scale_x_continuous("Hour") +
scale_y_continuous("Wind Speed (mph)") + labs(colour = "Year")





