use vertical strain for lilly correction

Author: haakon-e

src/parameterized_tendencies/les_sgs_models/smagorinsky_lilly.jl

@@ -43,7 +43,7 @@ function set_smagorinsky_lilly_precomputed_quantities!(Y, p)
     ᶜθ_v = @. lazy(TD.virtual_pottemp(thermo_params, ᶜts))
     ᶜ∇ᵥθ = @. ᶜtemp_scalar_2 = Geometry.WVector(ᶜgradᵥ(ᶠinterp(ᶜθ_v))).components.data.:1
     ᶜN² = @. ᶜtemp_scalar = grav / ᶜθ_v * ᶜ∇ᵥθ
-    ᶜS_norm = @. ᶜtemp_scalar_2 = √(2 * norm_sqr(ᶜS))
+    ᶜS_norm = projected_strain_rate_norm(ᶜS, Geometry.WAxis())
 
     ᶜRi = @. ᶜtemp_scalar = ᶜN² / (ᶜS_norm^2 + eps(FT))  # Ri = N² / |S|²
     ᶜfb = @. ᶜtemp_scalar = ifelse(ᶜRi ≤ 0, 1, max(0, 1 - ᶜRi / Pr_t)^(1 / 4))