Thanks for all responses!

Farhid in which version of ROMS you've implemented it?

Looking in Rutger's version and lmd_vmix.F there is only diffusivity due to the internal waves and not because of surface wind induced waves?

lmd_vmix.F

**Code:**

! Compute interior diffusivity due to wave breaking (Gargett and

! Holloway.

Patrick,

you have put non-breaking surface wind wave effect into kpp in CROCO? Please if you can point me into which subroutine I should look.

My intention is to put it into gls, it should go into shear production term as addition "W"

\begin{equation}

\ni (M^2 + \alpha N^2) + W

\end{equation}

If you have information about wind induced waves from spectral model (or via coupling) it should be straight forward.

W is simple function of H_S is significant wave height, \omega is peak angular frequency, k is wavenumber, z is depth:

\begin{equation}

W = 0.0014 * k (\frac{\omega H_s}{2) e^{kz})^3

\end{equation}

Thanks!

Ivica