Index

Rundeck: P1SoM40

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

uv_eqn

Module: SOCPBL

File: PBL.f

Summary: uv_eqn integrates differential eqns for mean velocity u and v 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: km * du/dz = cm * usurf * u and km * dv/dz = cm * usurf * v. At the top, i.e., the middle of the first GCM layer, u, v are prescribed.
Author : Ye Cheng/G. Hartke
Version:

---

Declaration:

subroutine uv_eqn(u0, v0, u, v, z, km, dz, dzh, ustar, cm, z0m, utop, vtop, utop_old, vtop_old, dtime, coriol, ug, vg, uocean, vocean, n, dpdxr, dpdyr, dpdxr0, dpdyr0)

Variables:

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

z-profle of west-east velocity component

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

z-profle of south-north velocity component

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

z-profle of u at previous time step

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

z-profle of v at previous time step

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

z-profile of turbulent viscosity

kh :

z-profile of turbulent conductivity

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)

ustar : real*8, intent(in)

friction velocity at the surface

cm : real*8, intent(in)

dimensionless momentum flux at surface (drag coeff.)

ch :

dimensionless heat flux at surface (stanton number)

cq :

dimensionless moisture flux at surface (dalton number)

z0m : real*8, intent(in)

roughness height for momentum (if itype=1 or 2)

utop : real*8, intent(in)

due east component of the wind at the first GCM layer

vtop : real*8, intent(in)

due north component of the wind at the first GCM layer

dtime : real*8, intent(in)

time step

coriol : real*8, intent(in)

the Coriolis parameter

ug : real*8, intent(in)

due east component of the geostrophic wind

vg : real*8, intent(in)

due north component of the geostrophic wind

utop_old : real*8, intent(in)

utop after turb. diffusion during previous timestep

vtop_old : real*8, intent(in)

vtop after turb. diffusion during previous timestep

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)

rhs1 : real*8, dimension(n)

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

uocean : real*8, intent(in)

vocean : real*8, intent(in)

dpdxr : real*8, intent(in)

dpdyr : real*8, intent(in)

dpdxr0 : real*8, intent(in)

dpdyr0 : real*8, intent(in)

iter : integer

loop variable

j : integer

loop variable

i : integer

loop variable

factx : real*8

facty : real*8

dpdx : real*8

dpdy : real*8

usurf : real*8

factor : real*8