This function generates daily climatological maps.

daily_clima(b, dates = NULL, verbose = FALSE)

Arguments

b

RasterBrick/Stack containing the historical observations or a proxy (typically a reanalysis) that is used to derive the climatological information. This needs to contain daily layers for multiple years. b@z should contain dates for comparison (e.g. 2017.01.01).

dates

Dates for which we need to calculate daily climatology. By default, this is a leap year.

verbose

Logical defining whether message should be printed on screen (FALSE by default).

Value

The function returns a RasterBrick (if dates contains one date) or a list of RasterBricks (if dates contains more than one date). Extent, resolution and land-sea mask match those of b. Values are the subset of b related to the given dates.

Details

the climatology is calculated over the period spanning 4 days before and

Examples

if (FALSE) {
  # Generate dummy RasterBrick
  set.seed(0)
  r <- raster(nrows = 2, ncols = 2,
              xmn = 0, xmx = 360, ymn = -90, ymx = 90, vals = 30)
  b <- raster::brick(lapply(1:(365 * 3),
       function(i) raster::setValues(r, runif(n = raster::ncell(r),
                                              min = 0, max = 100))))
  raster::setZ(b) <- seq.Date(from = as.Date("1993-01-01"),
                              to = as.Date("1995-12-31"),
                              by = "day")
  daily_clima(b, as.Date("1996-01-01"))
}