plot_season <- function(x, binwidth=10) {
  require('burnr')
  require('ggplot2')
  stopifnot(is.fhx(x))
  message('plot_season(): This function is experimental and will likely change in the future.')
  x_events <- subset(x, grepl('_fs|_fi', rec_type))
  p <- ggplot2::ggplot(x_events, aes(year, fill = rec_type))
  p <- (p + ggplot2::geom_histogram(binwidth = binwidth)
  brks.major <- NA
  brks.minor <- NA
  yr_range <- diff(range(x_events$year))
  if (yr_range < 100) {
      brks.major = seq(round(min(x_events$year), -1),
                       round(max(x_events$year), -1),
                       10)
      brks.minor = seq(round(min(x_events$year), -1),
                       round(max(x_events$year), -1),
                       5)
  } else if (yr_range >= 100) {
      brks.major = seq(round(min(x_events$year), -2),
                       round(max(x_events$year), -2),
                       100)
      brks.minor = seq(round(min(x_events$year), -2),
                       round(max(x_events$year), -2),
                       50)
  }
  p <- (p + ggplot2::scale_x_continuous(breaks = brks.major, minor_breaks = brks.minor)
          + ggplot2::theme_bw()
          + ggplot2::theme(axis.title.x = ggplot2::element_blank(),
                  legend.title = ggplot2::element_blank(),
                  legend.position = "bottom"))
  p
}