IndexRundeck: P1SoM40Created: Fri May 10 03:30:18 EDT 2024

kppmix

Module:
File: OCNKPP.f
Summary: KPPMIX Main driver subroutine for kpp vertical mixing scheme and interface to greater ocean model
Author :
Version:

Declaration:
Variables:
mdiff : integer, parameter
number of diffusivities for local arrays
Initial Value = 3
ze : real*8, dimension(0:LMO)
GISS vertical layering (m)
zgrid : real*8, dimension(0:km+1)
vertical grid (<= 0) (m)
hwide : real*8, dimension(0:km+1)
layer thicknesses (m)
byhwide : real*8, dimension(0:km+1)
1/layer thicknesses (1/m)
kmtj : integer
number of vertical layers on this row
shsq : real*8, dimension(km)
(local velocity shear)^2 (m/s)^2
dvsq : real*8, dimension(km)
(velocity shear re sfc)^2 (m/s)^2
ustar : real*8
surface friction velocity (m/s)
bo : real*8
surface turbulent buoy. forcing (m^2/s^3)
bosol : real*8
radiative buoyancy forcing (m^2/s^3)
alphadt : real*8, dimension(km)
alpha * DT across interfaces (kg/m^3)
betads : real*8, dimension(km)
beta * DS across interfaces (kg/m^3)
dbloc : real*8, dimension(km)
local delta buoyancy across interfaces (m/s^2)
ritop : real*8, dimension(km)
numerator of bulk Richardson Number (m/s)^2
coriol : real*8
Coriolis parameter (1/s)
ldd : logical
= TRUE for double diffusion
visc : real*8, dimension(0:km+1)
vertical viscosity coefficient (m^2/s)
difs : real*8, dimension(0:km+1)
vertical scalar diffusivity (m^2/s)
dift : real*8, dimension(0:km+1)
vertical temperature diffusivity (m^2/s)
ghats : real*8, dimension(km)
nonlocal transport (s/m^2)
hbl : real*8
boundary layer depth (m)
byhbl : real*8
1/boundary layer depth (1/m)
bfsfc : real*8
surface buoyancy forcing (m^2/s^3)
ws : real*8
momentum velocity scale
wm : real*8
scalar velocity scale
casea : real*8
= 1 in case A; =0 in case B
stable : real*8
1 in stable forcing; 0 in unstable
dkm1 : real*8, dimension(mdiff)
boundary layer difs at kbl-1 level
gat1 : real*8, dimension(mdiff)
shape function at sigma=1
dat1 : real*8, dimension(mdiff)
derivative of shape function at sigma=1
blmc : real*8, dimension(km, mdiff)
boundary layer mixing coefficients
sigma : real*8
normalized depth (d / hbl)
rib : real*8, dimension(2)
bulk Richardson number
kbl : integer
index of first grid level below hbl
kmax : integer
minimum of LSRPD and kmtj, used in swfrac
rigg : real*8
local richardson number
fcon : real*8
function of Rig
fri : real*8
function of Rig
ftop : real*8
function of topography
delta : real*8
fraction hbl lies beteen zgrid neighbors
zehat : real*8
= zeta * ustar**3
ki : integer

mr : integer

ka : integer

ku : integer

kl : integer

iz : integer

izp1 : integer

ju : integer

jup1 : integer

ksave : integer

kn : integer

kt : integer

ratio : real*8

zdiff : real*8

zfrac : real*8

fzfrac : real*8

wam : real*8

wbm : real*8

was : real*8

wbs : real*8

u3 : real*8

bvsq : real*8

delhat : real*8

dvdzup : real*8

dvdzdn : real*8

viscp : real*8

diftp : real*8

visch : real*8

difsh : real*8

f1 : real*8

bywm : real*8

byws : real*8

sig : real*8

a1 : real*8

a2 : real*8

a3 : real*8

gm : real*8

gs : real*8

gt : real*8

dstar : real*8

udiff : real*8

ufrac : real*8

vtsq : real*8

r : real*8

difsp : real*8

difth : real*8

dkmp5 : real*8

Simplex Website Curator: Igor Aleinov — NASA Official: Gavin A. Schmidt

Contact GISS NASA Privacy PolicyAccessibility