Index

Rundeck: P1SoM40

Created: Fri May 10 03:30:18 EDT 2024

getk

Module: SOCPBL

File: PBL.f

Summary: getk computes the turbulent diffusivities for momentum, Km, for heat, Kh, for moisture, Kq and for kinetic energy, Ke, at the secondary grids, using the GISS second order closure model (Cheng et al., 2000). u,v,t,q,ke are calculated at the primary grid z, while e,lscale,km,kh,gm,gh are calculated at the secondary grid zhat.
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine getk(km, kh, kq, ke, gma, gha, u, v, t, e, lscale, dzh, n)

Variables:

t_real :

z-profiles

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

z-profiles

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

z-profiles

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

z-profiles

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

z-profiles

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

z-profiles

dz :

(j) zhat(j)-zhat(j-1)

dzh : real*8, dimension(n-1), intent(in)

(j) z(j+1)-z(j)

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

turbulent viscosity for u and v equations

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

turbulent conductivity for t and q equations

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

turbulent diffusivity for e equation

gma : real*8, dimension(n-1), intent(out)

normalized velocity gradient, tau**2*as2

gha : real*8, dimension(n-1), intent(out)

normalized temperature gradient, tau**2*an2

tau : real*8

B1*lscale/sqrt(2*e)

as2 : real*8

shear squared, (dudz)**2+(dvdz)**2

an2 : real*8

Brunt-Vaisala frequency, grav/T*dTdz

se :

stability constant for e, adjustable

teeny : used from constant

n : integer, intent(in)

array dimension

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

i : integer

loop variable

dudz : real*8

dvdz : real*8

ell : real*8

den : real*8

qturb : real*8

gh : real*8

gm : real*8

gmmax : real*8

sm : real*8

sh : real*8

sq : real*8

taue : real*8