The space differencing in the turbulent kinetic energy equation
(equation (10) in Part I) is evaluated by the forth
order centered scheme for advection terms and the second order
centered scheme for other terms.
In time integration, the forward scheme is adapted for the
are same as those of
(22) and (24).
The finite difference form of the surface flux are as follows.
The bulk Richardson number is calculated as follows.