Index

Rundeck: P1SoM40

Created: Thu May 9 03:30:17 EDT 2024

e_eqn

Module: SOCPBL

File: PBL.f

Summary: e_eqn integrates differential eqns for the turbulent kinetic energy, e, using tridiagonal method over npbl-1(=7) sublayer edges (i.e., the secondary grids). The boundary condition near the bottom is: e(1)=(1/2)*B1**(2/3)*ustar**2. At the top secondary grid, nearest to the middle of the first GCM layer, e is prescribed.
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine e_eqn(esave, e, u, v, t, km, kh, ke, lscale, dz, dzh, ustar, dtime, n)

Variables:

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

z-profle of west-east velocity component

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

z-profle of south-north velocity component

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

z-profle of virtual potential temperature

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

z-profile of turbulent viscosity

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

z-profile of turbulent conductivity

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

z-profile of turbulent diffusion in eqn for e

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

z-profile of the turbulent dissipation length scale

z :

vertical grids (main, meter)

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

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

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

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

dtime : real*8, intent(in)

time step

ustar : real*8, intent(in)

friction velocity at the surface

n : integer, intent(in)

array dimension

n : integer, intent(in)

array dimension

sub : real*8, dimension(n)

dia : real*8, dimension(n)

sup : real*8, dimension(n)

rhs : real*8, dimension(n)

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

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

j : integer

loop variable

an2 : real*8

dudz : real*8

dvdz : real*8

as2 : real*8

qturb : real*8

ri : real*8

aa : real*8

bb : real*8

cc : real*8

gm : real*8

tmp : real*8