R/get_fire_season.R

In ecmwf/caliver: Calibration and Verification of Gridded Model Outputs

#' @title get_fire_season
#'
#' @description Get the season for each date in a vector
#'
#' @param dates vector of dates
#' @param fss Fire Season Start (date in the format Y-m-d)
#' @param fse Fire Season End (date in the format Y-m-d)
#' @param zone this can either: "north", "south" or "tropics"
#' 
#' @return A logical vector, where TRUE corresponds to a date in the fire season
#' and FALSE correspond to a date not in the fire season.
#'
#' @export
#'
#' @examples
#' \dontrun{
#'   # Modify default seasons
#'   seasons <- get_fire_season(dates,
#'                            fss = as.Date("2012-04-01", format = "%Y-%m-%d"),
#'                            fse = as.Date("2012-10-31", format = "%Y-%m-%d"))
#' }
#'

get_fire_season <- function(dates, fss = NULL, fse = NULL, zone = NULL){

  # Convert dates from any year to 2012 dates
  d <- as.Date(strftime(dates, format = "2012-%m-%d"))

  if (!is.null(zone)){

    if (zone == "north"){

      # Fire Season Start
      fss <- as.Date("2012-04-01",  format = "%Y-%m-%d")
      # Fire Season End
      fse <- as.Date("2012-10-31",  format = "%Y-%m-%d")

    }

    if (zone == "south"){

      # Fire Season Start
      fss <- as.Date("2012-10-01", format = "%Y-%m-%d")
      # Fire Season End
      fse <- as.Date("2012-04-30", format = "%Y-%m-%d")

    }

    if (zone == "tropics"){

      # Fire Season Start
      fss <- as.Date("2012-06-01", format = "%Y-%m-%d")
      # Fire Season End
      fse <- as.Date("2012-11-30", format = "%Y-%m-%d")

    }

  }

  # Season
  first_january <- as.Date("2012-01-01", format = "%Y-%m-%d")
  last_december <- as.Date("2012-12-31", format = "%Y-%m-%d")
  if (fss < fse){

    season <- ifelse(d >= fss & d <= fse, TRUE, FALSE)

  } else {

    season <- ifelse(d >= first_january & d <= fse |
                       d >= fss & d <= last_december,
                     TRUE, FALSE)

  }

  return(season)

}