Index

Rundeck: P1SoM40

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

q_eqn_sta

Module: SOCPBL

File: PBL.f

Summary: q_eqn_sta computes the static solutions of the specific humidity, Q, between the surface and the first GCM layer. The boundary condition at the bottom is: kq * dq/dz = cq * usurf * (q - qg). At the top, Q is prescribed. It is called only at the initialization (from within subroutine inits).
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine q_eqn_sta(q, kq, dz, dzh, cq, usurf, qgrnd, qtop, n)

Variables:

u :

z-profle of west-east velocity component

v :

z-profle of south-north velocity component

t :

z-profle of virtual potential temperature

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

z-profle of specific humidity

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

z-profile of moisture diffusivity

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)

cq : real*8, intent(in)

dimensionless moisture flux at surface (dalton number)

usurf : real*8, intent(in)

effective surface velocity

qgrnd : real*8, intent(in)

specific humidity at the ground

qtop : real*8, intent(in)

specific humidity 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

factq : real*8