•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
–