Add missing L_StokesL computation

CVMix · Oct 30, 2024 · a8d8125 · a8d8125
1 parent 5e3d0ce
commit a8d8125
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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