error langmuir wave coupling
Error in the definition of the unit vector for the calculation of the Stokes drift contribution to Langmuir cells. The norm of the unit vector should also be ^2 in the calculation of zWlc2.
Corrections to be done in zdftke.F90 in the subroutine tke_tke:
Delete the SQRT in z1_norm: z1_norm = 1._wp / MAX( ztauiztaui+ztaujztauj, 1.e-12 ) * tmask(ji,jj,1)