tke_eps_advect_3d - 3D turbulence advection (Source File: tke_eps_advect_3d.F90)

INTERFACE:

    subroutine tke_eps_advect_3d()

DESCRIPTION:

This routine carries out advection of the prognostic turbulence quantities tke (turbuent kinetic energy, $k$) and eps (lenght scale related turbulence quantity, e.g. dissipation rate of $k$, $\varepsilon$, or turbulent frequency, $\omega=\varepsilon/k$. Here, the TVD advection schemes are used which are also used for the momentum advection. USES:

    use domain, only: imin,imax,jmin,jmax,kmax,az,ax
 #if defined(SPHERICAL) || defined(CURVILINEAR)
    use domain, only: dxv,dyu
 #else
    use domain, only: dx,dy
 #endif
    use m3d, only: turb_adv_split,turb_adv_hor,turb_adv_ver
    use variables_3d, only: tke,eps,dt,uu,vv,ww,hun,hvn,ho,hn
    use advection, only: J7
    use advection_3d, only: do_advection_3d,W_TAG
    use halo_zones, only: update_3d_halo,wait_halo,H_TAG
    use turbulence, only: k_min,eps_min
  $ use omp_lib
    IMPLICIT NONE

REVISION HISTORY:

    Original author(s): Hans Burchard & Karsten Bolding

LOCAL VARIABLES:

    integer                      :: i,j,k
    REALTYPE,dimension(I3DFIELD) :: uuadv,vvadv,wwadv,hoadv,hnadv,huadv,hvadv