 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| • |
KPP
Procedure (part 2)
|
|
|
|
– |
Calculate
surface boundary layer k profiles for T, S, and momentum
|
|
|
|
• |
Vertical
diffusivity for T is parameterized as
|
|
|
|
– |
where is the nonlocal transport term
|
|
|
|
• |
Diffusivity
is parameterized as
|
|
|
|
– |
where
w is a turbulent velocity scale that is a function of the
|
|
|
stability
of the forcing, G is a 3rd order polynomial shape
|
|
|
function,
and sigma is a scale depth varying from 0 to 1 over the
|
|
depth
range H
|
|
|
|
• |
Choose
coefficients of G to match the interior and boundary layer K
|
|
|
profiles,
producing a final K profile with a continuous first vertical
|
|
|
derivative
|
|
|
|
– |
Solve
diffusion equation semi-implicitly with two temporal iterations
|
|
|
|
– |
Diagnose
mixed layer thickness along with T, S, u, and v
|
|