fMonth.long <- factor(air.long$Month)
meteo <- ggplot(air.long, aes(x = fMonth.long, y = value)) +
    geom_boxplot() + 
    facet_wrap(~variable, nrow = 2)
meteo