fp_calculate: Two-dimensional footprint models
In grahamstewart12/tidyflux:

Create a two-dimensional matrix of footprint probabilities.

fp_calculate(timestamp, wd, ustar, l, v_sd, ws, weights = NULL,
  fetch = (z - zd) * 100, grid = 2 * fetch, model = c("KM01", "H00"),
  percent_out = FALSE, z = metadata$tower_height,
  zd = metadata$displacement, zo = metadata$roughness_length,
  raster_out = FALSE, coords, clip_area = NULL)

`timestamp`	A vector of class POSIXct
`wd`	A numeric vector representing wind direction in degrees from North.
`ustar`	A numeric vector representing friction velocity in units ms+1s-1.
`l`	A numeric vector representing Monin-Obukhov length in units m.
`v_sd`	A numeric vector representing the standard deviation of the crosswind component.
`ws`	A numeric vector representing wind speed in units ms+1s-1.
`fetch`	An integer specifying the maximum length in each direction to calculate footprint probabilities. Defaults to 100 * (z - zd).
`grid`	An integer specifying the number of grid cells to be used in calculating footprint probabilities. Defaults to the fetch length (i.e. a 1-m resolution).
`model`	A character string naming the model to be used in footprint calculations. Can be "KM01" for the Kormann & Meixner (2001) model or "H00" for the Hsieh et al. (2000) model.
`z`	The tower height in meters. Defaults to an attempt to retrieve the value from the 'tower_height' variable in the 'metadata' list in the global environment.
`zd`	The zero-plane displacement height in meters. Defaults to an attempt to retrieve the value from the 'displacement' variable in the 'metadata' list in the global environment.
`zo`	The aerodynamic roughness length in meters. Defaults to an attempt to retrieve the value from the 'roughness_length' variable in the 'metadata' list in the global environment.