Fix sign error on stability variable (#579)

* Update similarity_theory_turbulent_fluxes.jl

* run examples

* bump patch release

---------

Co-authored-by: Navid C. Constantinou <[email protected]>

Author: simone-silvestri

Project.toml

@@ -2,7 +2,7 @@ name = "ClimaOcean"
 uuid = "0376089a-ecfe-4b0e-a64f-9c555d74d754"
 license = "MIT"
 authors = ["Climate Modeling Alliance and contributors"]
-version = "0.7.1"
+version = "0.7.2"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"

src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl

@@ -211,7 +211,7 @@ function iterate_interface_fluxes(flux_formulation::SimilarityTheoryFluxes,
 
     # Transfer coefficients at height `h`
     ϰ = flux_formulation.von_karman_constant
-    L★ = ifelse(b★ == 0, Inf, - u★^2 / (ϰ * b★))
+    L★ = ifelse(b★ == 0, Inf, u★^2 / (ϰ * b★))
     form = flux_formulation.similarity_form
 
     χu = ϰ / similarity_profile(form, ψu, Δh, ℓu₀, L★)
@@ -253,7 +253,7 @@ Note that the Monin--Obukhov characteristic length scale is defined
 in terms of ``b_★`` and additionally the Von Karman constant ``ϰ``,
 
 ```math
-L_★ = - u_★² / ϰ b_★ .
+L_★ = u_★² / ϰ b_★ .
 ```
 """
 @inline function buoyancy_scale(θ★, q★, ℂ, 𝒬, g)