Index

Rundeck: P1SoM40

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

t_eqn_sta

Module: SOCPBL

File: PBL.f

Summary: t_eqn_sta computes the static solutions of the virtual potential temperature, T, between the surface and the first GCM layer. The boundary condition at the bottom is: kh * dt/dz = ch * usurf * (t - tg). At the top, T is prescribed. It is called only at the initialization (from within subroutine inits).
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine t_eqn_sta(t, q, kh, kq, dz, dzh, ch, usurf, tgrnd, ttop, n)

Variables:

u :

z-profle of west-east velocity component

v :

z-profle of south-north velocity component

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

z-profle of virtual potential temperature

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

z-profle of specific humidity

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

z-profile of heat conductivity

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

z-profile of specific humidity

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)

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

n : integer, intent(in)

number of vertical main subgrid layers

sub : real*8, dimension(n)

dia : real*8, dimension(n)

sup : real*8, dimension(n)

rhs : real*8, dimension(n)

i : integer

loop variable

facth : real*8