Channel flow boundary condition

Variant of Moeng boundary condition which also allows to have a wall function at the top boundary. More...

Functions/Subroutines
subroutine	abl_userbc (ix, iy, iz, iside, eg)
	Stress boundary condition as formulated by Moeng (1984) Also has optional temporal filtering c.f: Yang et al. Physical Review Fluids 2, no. 10 (2017): 104601. https://doi.org/10.1103/PhysRevFluids.2.104601. More...
real function	abl_pen_k (y_coord, z0)
	Compute the K term in penalty forcing. More...

Detailed Description

Variant of Moeng boundary condition which also allows to have a wall function at the top boundary.

Function/Subroutine Documentation

abl_pen_k()

real function abl_pen_k	(	real	y_coord,
		real	z0
	)

Compute the K term in penalty forcing.

Parameters

y_coord	mesh coordinate where the coefficient is evaluated
z0	roughness length

abl_userbc()

subroutine abl_userbc	(	integer	ix,
		integer	iy,
		integer	iz,
		integer	iside,
		integer	eg
	)

Stress boundary condition as formulated by Moeng (1984) Also has optional temporal filtering c.f: Yang et al. Physical Review Fluids 2, no. 10 (2017): 104601. https://doi.org/10.1103/PhysRevFluids.2.104601.

Note

Boundary condition evaluated at higher nodes, may need to be corrected for Ekmann turning

This subroutine MAY NOT be called by every process