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