Matlab function to calculate vertical mixing coefficients via the KPP model: KPP.m

Reference for the method (please cite if you use it in published research):
Smyth, W.D., E.D. Skyllingstad, G.B. Crawford and H. Wijesekera, 2002: "Nonlocal fluxes and Stokes drift effects in the K-profile parameterization", Ocean Dynamics 52 (3), 104-115.[PDF]