The Model (9)
•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
–