Index

Rundeck: E6F40

Created: Fri May 17 03:30:27 EDT 2024

e_gcm

Module:

File: ATURB.f

Summary: finds the turbulent kinetic energy (e, in m^2/s^2). Within the PBL, e is determined according to the parameterization of the Large Eddy Simulation (LES) data (Moeng and Sullivan, 1994), above the PBL, e is calculated by the second order closure model of Cheng et al. (2002).
Author : Ye Cheng
Version:

---

Declaration:

subroutine e_gcm(tvflx, wstar, ustar, dbl, lmonin, ze, g_alpha, an2, as2, lscale, e, n)

Variables:

lmonin : real*8, intent(in)

Monin-Obukov length

g_alpha : real*8, dimension(n), intent(in)

grav*alpha

teeny : used from constant

by3 : used from constant

kappa : used from socpbl

emax : used from socpbl

rimax : used from socpbl

b1 : used from socpbl

c1 : used from socpbl

c2 : used from socpbl

c3 : used from socpbl

c4 : used from socpbl

c5 : used from socpbl

gm_at_rimax : used from socpbl

emin : used from socpbl

find_phim0 : used from socpbl

n : integer, intent(in)

array dimension

tvflx : real*8, intent(in)

wstar : real*8, intent(in)

ustar : real*8, intent(in)

dbl : real*8, intent(in)

an2 : real*8, dimension(n), intent(in)

as2 : real*8, dimension(n), intent(in)

lscale : real*8, dimension(n), intent(in)

ze : real*8, dimension(n+1), intent(in)

j : integer

loop variable

e : real*8, dimension(n), intent(out)

ri : real*8

gm : real*8

aa : real*8

bb : real*8

cc : real*8

phim : real*8

tmp : real*8

wstar3 : real*8

ustar3 : real*8

zj : real*8

zet : real*8

eps : real*8

ej : real*8

kz : real*8