scientimate.surfaceroughness

ustar, z0, d = scientimate.surfaceroughness(z, u, delta=None, dispout='no')

Description

Calculate shear velocity and surface roughness from a given velocity profile using Karimpour et al. (2012) method

Inputs

z

Distance from a surface (elevation, height) in (m)

u

Velocity at z in (m/s)

delta=max(z)

Boundary layer height in (m)

dispout=’no’

Define to display outputs or not (‘yes’: display, ‘no’: not display)

Outputs

z0

Surface roughness in (m)

ustar

Shear velocity (u*) in (m/s)

d

Zero plane displacement distance in (m)

Note: Above values are for a logarithmic velocity profile as:

u=(u*/K)*ln((z-d)/z0)

Examples

import scientimate as sm
import numpy as np

z=np.arange(0.1,1,0.05)
u=2/0.4*np.log((z-0.003)/0.002)
ustar,z0,d=sm.surfaceroughness(z,u,np.max(z),'yes')

References

Karimpour, A., Kaye, N. B., & Baratian-Ghorghi, Z. (2012). Modeling the neutrally stable atmospheric boundary layer for laboratory scale studies of the built environment. Building and Environment, 49, 203-211.