Index

Rundeck: E6F40

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

k_gcm

Module:

File: ATURB.f

Summary: computes the turbulent stability functions Km (for momentum) and Kh (for heat and moisture), as well as the fluxes (local and non-local). Within the convective PBL, it is according to Holtslag and Boville (1993); within the stable PBL or above the PBL, it is according to Cheng et al. (2002).
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine k_gcm(tvflx, qflx, uflx, vflx, ustar, wstar, dbl, lmonin, ze, lscale, e, qturb, an2, as2, dtdz, dqdz, dudz, dvdz, kh, kq, km, ke, wt, wq, w2, uw, vw, wt_nl, wq_nl, uw_nl, vw_nl, n)

Variables:

tvflx : real*8, intent(in)

virtual potential temperature flux at surface

qflx : real*8, intent(in)

moisture flux at surface

ustar : real*8, intent(in)

friction velocity

wstar : real*8, intent(in)

Deardorff convective velocity scale

dbl : real*8, intent(in)

the height of the PBL (real*8, in meters)

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

turbulent diffusivity for scalars (heat, moisture,...)

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

height (in meters) of layer edges

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

turbulent length scale

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

turbulent kinetic energy

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

sqrt(2*e)

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

g_alpha*dTdz

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

(dUdz)**2+(dVdz)**2

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

dt/dz

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

dq/dz

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

du/dz

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

dv/dz

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

turbulent diffusivity for momentun

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

turbulent diffusivity for e

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

vertical component of 2*e

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

turbulent fluxes

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

turbulent fluxes

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

turbulent fluxes

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

turbulent fluxes

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

non-local part of heat flux wt

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

non-local part of moisture flux wq

wc_nl :

non-local part of tracer flux

n : integer, intent(in)

number of layers

teeny : used from constant

by3 : used from constant

sha : used from constant

kappa : used from socpbl

prt : used from socpbl

ghmin : used from socpbl

ghmax : used from socpbl

d1 : used from socpbl

d2 : used from socpbl

d3 : used from socpbl

d4 : used from socpbl

d5 : used from socpbl

s0 : used from socpbl

s1 : used from socpbl

s2 : used from socpbl

s4 : used from socpbl

s5 : used from socpbl

s6 : used from socpbl

s7 : used from socpbl

s8 : used from socpbl

b1 : used from socpbl

g5 : used from socpbl

k_max : used from socpbl

kmmin : used from socpbl

khmin : used from socpbl

find_phim0 : used from socpbl

find_phih : used from socpbl

uflx : real*8, intent(in)

vflx : real*8, intent(in)

lmonin : real*8, intent(in)

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

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

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

j : integer

loop variable

tmp : real*8

tau : real*8

gm : real*8

gh : real*8

gmmax : real*8

byden : real*8

sm : real*8

sh : real*8

ustar2 : real*8

wstar3 : real*8

zzi : real*8

tau_pt : real*8

w2j : real*8

phih1 : real*8

by_phim1 : real*8

wm1 : real*8

pr1 : real*8

cgh1 : real*8

km_n : real*8

kh_n : real*8

pr : real*8

cgh : real*8

cgq : real*8

zet : real*8

phih : real*8

phim : real*8

by_phim : real*8

wm : real*8

kz : real*8

cgu1 : real*8

cgv1 : real*8