From 3209135bdc344efa37dd1c5c0c51721ae2a3c38b Mon Sep 17 00:00:00 2001 From: Yair Cohen Date: Mon, 6 Jun 2022 11:00:37 -0700 Subject: [PATCH] Update test/relations.jl Co-authored-by: Charles Kawczynski --- test/relations.jl | 12 ++++++------ 1 file changed, 6 insertions(+), 6 deletions(-) diff --git a/test/relations.jl b/test/relations.jl index ed65ee49..7ed2f02a 100644 --- a/test/relations.jl +++ b/test/relations.jl @@ -501,13 +501,13 @@ end # test Clausius Clapeyron relation ts = PhaseEquil_ρeq.(param_set, ρ, e_int, q_tot) - L_v = TD.latent_heat_vapor(param_set, ts) - L_s = TD.latent_heat_sublim(param_set, ts) - λ = TD.liquid_fraction(param_set, ts) - L = λ * L_v + (1 - λ) * L_s - T = TD.air_temperature(param_set, ts) + L_v = TD.latent_heat_vapor.(param_set, ts) + L_s = TD.latent_heat_sublim.(param_set, ts) + λ = TD.liquid_fraction.(param_set, ts) + L = λ .* L_v .+ (1 .- λ) .* L_s + T = TD.air_temperature.(param_set, ts) q_vap_sat = TD.vapor_specific_humidity(param_set, ts) - @test TD.∂q_vap_sat_∂T(param_set, ts) - L^2 * q_vap_sat / (R_v * T^2) + @test all(TD.∂q_vap_sat_∂T(param_set, ts) .≈ L .^ 2 .* q_vap_sat ./ (R_v .* T.^2)) # PhaseEquil (freezing) _T_freeze = FT(ICP.T_freeze(param_set))