Index | Rundeck: P1SoM40 | Created: Fri May 10 03:30:18 EDT 2024 |
kppmix | Module: | File: OCNKPP.f |
Author :
Version:
Declaration:
subroutine kppmix(ldd, ze, zgrid, hwide, kmtj, shsq, dvsq, ustar, bo, bosol, alphadt, betads, dbloc, ritop, coriol, byhwide, visc, difs, dift, ghats, hbl, kbl )
- 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