diff --git a/src/shared/cvmix_kpp.F90 b/src/shared/cvmix_kpp.F90 index 38ee525d..0e132b17 100644 --- a/src/shared/cvmix_kpp.F90 +++ b/src/shared/cvmix_kpp.F90 @@ -2434,6 +2434,7 @@ subroutine cvmix_kpp_compute_turbulent_scales_1d_OBL(sigma_coord, & if (compute_ws) then do kw = 1,n_sigma chi_s = compute_Stokes_chi( xi(kw) , lchi_s=.true. ) + L_StokesL = cvmix_one - xi(kw) !wgl - xi w_s(kw)=compute_phi_inv(zeta(kw),CVmix_kpp_params_in,L_StokesL,lphi_s=.true.) * & vonkar*surf_fric_vel / chi_s end do