INTERFACE:
subroutine tke_eps_advect_3d()DESCRIPTION:
This routine carries out advection of the prognostic turbulence quantities tke (turbuent kinetic energy, ) and eps (lenght scale related turbulence quantity, e.g. dissipation rate of , , or turbulent frequency, . 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 NONEREVISION HISTORY:
Original author(s): Hans Burchard & Karsten BoldingLOCAL VARIABLES:
integer :: i,j,k REALTYPE,dimension(I3DFIELD) :: uuadv,vvadv,wwadv,hoadv,hnadv,huadv,hvadv