No-slip Dirichlet boundary condition

Functions/Subroutines
subroutine	abl_userbc (ix, iy, iz, iside, eg)
	No slip and no penetration condition. More...
real function	abl_pen_k (y_coord, z0)
	Compute the K term in penalty forcing. More...

Detailed Description

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

Note

Be careful in specifying the penalty region, should be away from the wall

abl_userbc()

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

No slip and no penetration condition.

Todo:

Why does u_star -> NaN? Only for statistics: compute u_star assuming log-law