Index

Rundeck: P1SoM40

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

t_eqn

Module: SOCPBL

File: PBL.f

Summary: t_eqn integrates differential eqns for the virtual potential temperature, T, using tridiagonal method over npbl(=8) sublayers between the surface (sublayer 1) and the middle of the first GCM layer (sublayer npbl). The boundary condition at the bottom is: kh * dt/dz = ch * ( usurf*(t1 - tgrnd) +(1+xdelt*qtop)*gusti*tprime ) which includes the effects on the surface flux due to the moist convection wind gustiness and the downdraft temperature perturbation (Redelsperger et al. 2000; Emanuel and Zivkovic 1999), where tprime=tdns-ttop/(1+xdelt*qtop), t1, q1 are the T and Q at the surface, and tdns is the downdraft temperature in K at (i,j), which is calculated in subroutines CONDSE (in CLOUDS2_DRV.f) and PBL (in PBL_DRV.f). At the top, i.e., the middle of the first GCM layer, T is prescribed. ### the following term was removed from BC at the bottom ### + xdelt * t1/(1+xdelt*q1) * kq * dqdz
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine t_eqn(u, v, t0, t, q, z, kh, kq, dz, dzh, ch, usurf, tgrnd, ttop, qtop, dtdt_gcm, dtime, n, dpdxr, dpdyr, dpdxr0, dpdyr0, usurf0, gusti, tprime, tdns, qdns, ddml_eq_1)

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(inout)

z-profle of virtual potential temperature ref. to the surface

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

z-profle of specific humidity

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

z-profle of t at previous time step

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

z-profile of heat conductivity

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

z-profile of moisture diffusivity

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

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)

ch : real*8, intent(in)

dimensionless heat flux at surface (stanton number)

usurf : real*8, intent(in)

effective surface velocity

tgrnd : real*8, intent(in)

virtual potential temperature at the ground

ttop : real*8, intent(in)

virtual potential temperature at the first GCM layer

dtdt_gcm : real*8, intent(in)

ttop tendency from processes other than turbulence

dtime : real*8, intent(in)

time step

n : integer, intent(in)

number of vertical subgrid main layers

sub : real*8, dimension(n)

dia : real*8, dimension(n)

sup : real*8, dimension(n)

rhs : real*8, dimension(n)

qtop : real*8, intent(in)

dpdxr : real*8, intent(in)

dpdyr : real*8, intent(in)

dpdxr0 : real*8, intent(in)

dpdyr0 : real*8, intent(in)

usurf0 : real*8, intent(in)

gusti : real*8, intent(in)

tprime : real*8, intent(in)

tdns : real*8, intent(in)

qdns : real*8, intent(in)

ddml_eq_1 : logical, intent(in)

i : integer

loop variable

facth0 : real*8

facth : real*8

factx : real*8

facty : real*8