windprofile: Calculates wind profile for entire canopy

In ilyamaclean/microclimc: Below or above canopy microclimate modelling with R

Calculates wind profile for entire canopy

Description

calculates wind speed at any given point above or below canopy from wind speed at reference height

Usage

windprofile(
  ui,
  zi,
  zo,
  a = 2,
  PAI,
  hgt,
  psi_m = 0,
  hgtg = 0.05 * hgt,
  zm0 = 0.004
)

Arguments

ui	wind at reference height (m / s)
zi	height of wind measurement (m)
zo	height for which wind is required (m)
a	attenuation coefficient as returned by microctools::attencoef()
PAI	plant area index (m / m)
hgt	height of canopy (m)
psi_m	diabatic correction factor as return by microctools::diabatic_cor()
hgtg	height of ground vegetation layer below canopy (m)
zm0	roughness length (m) of vegetation layer below canopy as returned by microctools::roughlength()

Value

wind speed (m /s) at height zo

Examples

zo<- c(0:200) / 100
uz <- windprofile(2, 2, zo, 2.3, 3, 2)
plot(zo ~ uz, type = "l", xlab = "wind speed", ylab = "height")

ilyamaclean/microclimc documentation built on July 28, 2023, 1:40 a.m.