get_mask {SeaVal}

Function to create a mask of dry regions from CHIRPS

Description

A gridpoint is masked for a given season (either 'MAM', 'JJAS' or 'OND'), if, on average, less than 10% of the annual total of rainfall occur during the season. This function loads CHIRPS data, and derives this mask as a data table of lon, lat coordinates, only containing the coordinates that shouldn't be masked. You can apply the mask to an existing data table using dt = combine(dt,mask).

Usage

get_mask(
  season,
  clim_years = 1990:2020,
  version = "UCSB",
  resolution = "low",
  us = (resolution == "low")
)

Arguments

season

For which season do you want to calculate the mask? Needs to be either 'MAM', 'JJAS' or 'OND'.

clim_years

Numeric vector of years. Which years should be used to establish the mask?

version, resolution, us

Passed to load_chirps. Which CHIRPS version do you want to use and on what resolution?

Examples

if(interactive()) get_mask('MAM')