fetch3.nhl_transpiration.NHL_functions.solve_Uz#

solve_Uz(z, dz, Cd, a_s, U_top, **kwargs)#

Solves the momentum equation to calculate the vertical wind profile. Applies no-slip boundary condition: wind speed=0 at surface (z = 0). Model for turbulent diffusivity of momentum is from Poggi 2004, eqn 6

Parameters:

  • z (array) – vector of heights [m]

  • dz (float) – Vertical discretization interval [m]

  • Cd (float) – drag coefficient [unitless]

  • a_s (float) – leaf surface area [m2]

  • U_top (float) – Measured wind speed at top of canopy [m s-1]

  • calc_mixing_length (**kwargs to be passed to) –

Returns:

  • Km (array) – turbulent diffusivity of momentum at each height in z [m s-1]

  • U (array) – wind speed at each height in z [m s-1]