Index

Rundeck: P1SoM40

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

inits

Module: SOCPBL

File: PBL.f

Summary: inits initializes the winds, virtual potential temperature, and humidity by solving their differential equations for the static solutions, using tridiagonal method over npbl(=8) sublayers between the surface (sublayer 1) and the middle of the first GCM layer (sublayer npbl). (Cheng et a., 2002). It is called by subroutine init_pbl (in PBL_DRV.f), and the latter (init_pbl) is called by subroutine INPUT (in MODELE.f).
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine inits(tgrnd, qgrnd, zgrnd, zgs, ztop, utop, vtop, ttop, qtop, coriol, cm, ch, cq, ustar, lmonin, uocean, vocean, ilong, jlat, itype, dpdxr, dpdyr, dpdxr0, dpdyr0, u, v, t, q, e, ug, vg)

Variables:

n :

number of sub-grid levels for the PBL

tgrnd : real*8, intent(in)

virtual potential temperature of ground,at roughness height

qgrnd : real*8, intent(in)

moisture at the ground, at the roughness height

zgrnd : real*8, intent(in)

zgs : real*8, intent(in)

height of the surface layer (nominally 10 m)

ztop : real*8, intent(in)

height of the first model layer, approx 200 m if lm=9

utop : real*8, intent(in)

x component of wind at the top of the layer

vtop : real*8, intent(in)

y component of wind at the top of the layer

ttop : real*8, intent(in)

virtual potential temperature at the top of the layer

qtop : real*8, intent(in)

moisture at the top of the layer

coriol : real*8, intent(in)

the Coriolis parameter

cm : real*8, intent(out)

dimensionless momentum flux at surface (drag coeff.)

ch : real*8, intent(out)

dimensionless heat flux at surface (stanton number)

cq : real*8, intent(out)

dimensionless moisture flux at surface (dalton number)

bgrid : real*8

log-linear gridding parameter

ustar : real*8, intent(out)

friction speed

ilong : integer, intent(in)

longitude identifier

jlat : integer, intent(in)

latitude identifier

itype : integer, intent(in)

surface type

vg : real*8

passed for scm

ug : real*8

passed for scm

iprint : integer, parameter

longitude for diagnostics
Initial Value = 0

jprint : integer, parameter

latitude for diagnostics
Initial Value = 41

uocean : real*8, intent(in)

vocean : real*8, intent(in)

lmonin : real*8, intent(out)

dpdxr : real*8, intent(in)

dpdyr : real*8, intent(in)

dpdxr0 : real*8, intent(in)

dpdyr0 : real*8, intent(in)

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

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

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

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

gm : real*8, dimension(n-1)

gh : real*8, dimension(n-1)

z : real*8, dimension(n)

dz : real*8, dimension(n)

xi : real*8, dimension(n)

usave : real*8, dimension(n)

vsave : real*8, dimension(n)

tsave : real*8, dimension(n)

qsave : real*8, dimension(n)

zhat : real*8, dimension(n-1)

xihat : real*8, dimension(n-1)

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

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

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

lmonin_dry : real*8

z0m : real*8

z0h : real*8

z0q : real*8

hemi : real*8

psi1 : real*8

psi0 : real*8

psi : real*8

usurf : real*8

tstar : real*8

qstar : real*8

ustar0 : real*8

dtime : real*8

test : real*8

wstar3 : real*8

wstar2h : real*8

usurfq : real*8

usurfh : real*8

ts : real*8

itmax : integer, parameter

Initial Value = 100

iter : integer

loop variable

i : integer

loop variable

w : real*8, parameter

Initial Value = 0.50

tol : real*8, parameter

Initial Value = 1d-3

ierr : integer

dbl : real*8

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

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

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

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

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

dm : real*8