From 843961eee6fdddfa09d5eec2e025603feec6d271 Mon Sep 17 00:00:00 2001 From: "Navid C. Constantinou" Date: Wed, 7 May 2025 15:25:22 +0300 Subject: [PATCH] Fix mistake in EdsonMomentumStabilityFunctions https://github.com/NOAA-PSL/COARE-algorithm/blob/5b144cf6376a98b42200196d57ae40d791494abe/Matlab/COARE3.6/coare36vn_zrf_et.m#L703 --- .../similarity_theory_turbulent_fluxes.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl b/src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl index c2f4e9f0a..1f9fa4649 100644 --- a/src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl +++ b/src/OceanSeaIceModels/InterfaceComputations/similarity_theory_turbulent_fluxes.jl @@ -277,7 +277,7 @@ These stability functions are obtained by regression to experimental data. The stability parameter for stable atmospheric conditions is defined as ```math dζ = min(ζmax, Aˢζ) -ψₛ = - (Bˢ ζ + Cˢ ( ζ - Dˢ ) ) exp( - dζ) - Cˢ Dˢ +ψₛ = - Bˢ * ζ⁺ - Cˢ * (ζ⁺ - Dˢ) * exp(- dζ) - Cˢ * Dˢ ``` While the stability parameter for unstable atmospheric conditions is calculated @@ -329,7 +329,7 @@ end dζ = min(ζmax, Aˢ * ζ⁺) # Stability parameter for _stable_ atmospheric conditions - ψₛ = - (Bˢ * ζ⁺ + Cˢ * (ζ⁺ - Dˢ)) * exp(- dζ) - Cˢ * Dˢ + ψₛ = - Bˢ * ζ⁺ - Cˢ * (ζ⁺ - Dˢ) * exp(- dζ) - Cˢ * Dˢ # Stability parameter for _unstable_ atmospheric conditions fᵤ₁ = sqrt(sqrt(1 - Aᵘ * ζ⁻))