The dynamic kb-equation

INTERFACE:

   subroutine kbeq(nlev,dt,u_taus,u_taub,z0s,z0b,h,NN,SS)

DESCRIPTION:

The transport equation for (half the) buoyancy variance, $k_b= {\langle b'^2 \rangle} /2$, follows from the equation for the buoyancy fluctations (see Sander (1998)). In the case of a Boussinesq-fluid, this equation can be written as

$$\dot{k_b} = {\cal D}_b + P_b - \epsilon_b \; , \quad$$ (158)

where $\dot{k_b}$ denotes the material derivative of $k_b$. $P_b$ is the production of $k_b$ be mean density gradients, and $\epsilon_b$ the rate of molecular destruction. ${\cal D}_b$ represents the sum of the viscous and turbulent transport terms. It is presently evaluated with a simple down gradient model in GOTM.

The production of buoyancy variance by the vertical density gradient is

$$P_b = - {\langle w'b' \rangle} \dfrac{\partial {B}}{\partial {z}} = - {\langle w'b' \rangle} N^2 \quad .$$ (159)

Its computation is discussed in section 4.8.

The rate of molecular destruction, $\epsilon_b$, can be computed from either a transport equation or a algebraic expression, section 4.7.10.

USES:

   use turbulence,   only: Pb,epsb,nuh
   use turbulence,   only: kb,kb_min
   use turbulence,   only: k_ubc, k_lbc, ubc_type, lbc_type
   use util,         only: Dirichlet,Neumann
 
   IMPLICIT NONE

INPUT PARAMETERS:

 
   number of vertical layers
   integer,  intent(in)                :: nlev
 
   time step (s)
   REALTYPE, intent(in)                :: dt
 
   surface and bottom
   friction velocity (m/s)
   REALTYPE, intent(in)                :: u_taus,u_taub
 
   surface and bottom
   roughness length (m)
   REALTYPE, intent(in)                :: z0s,z0b
 
   layer thickness (m)
   REALTYPE, intent(in)                :: h(0:nlev)
 
   square of shear and buoyancy
   frequency (1/s^2)
   REALTYPE, intent(in)                :: NN(0:nlev),SS(0:nlev)

REVISION HISTORY:

   Original author(s): Lars Umlauf
 
   $Log: kbeq.F90,v $
   Revision 1.1  2005-06-27 10:54:33  kbk
   new files needed