From 0e2c06754f21bcb5ec650fe0375941b04b8adbcd Mon Sep 17 00:00:00 2001 From: Zhaoyi Shen <11598433+szy21@users.noreply.github.com> Date: Tue, 10 Oct 2023 16:46:39 -0700 Subject: [PATCH] wip --- src/prognostic_equations/advection.jl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/prognostic_equations/advection.jl b/src/prognostic_equations/advection.jl index dcfe9979b5f..cc73904598b 100644 --- a/src/prognostic_equations/advection.jl +++ b/src/prognostic_equations/advection.jl @@ -213,7 +213,6 @@ NVTX.@annotate function explicit_vertical_advection_tendency!(Yₜ, Y, p, t) end end - # FIXME: add vertical advection if p.atmos.turbconv_model isa AdvectiveEDMFX for j in 1:n @. ᶜa_scalar[colidx] = @@ -229,11 +228,12 @@ NVTX.@annotate function explicit_vertical_advection_tendency!(Yₜ, Y, p, t) ) # FIXME: boundary conditions - @. Yₜ.c.sgsʲs.:($$j).h_tot[colidx] -= ᶜinterp( + LBC = Operators.LeftBiasedF2C(; bottom = Operators.SetValue(0)) + @. Yₜ.c.sgsʲs.:($$j).h_tot[colidx] -= LBC( adjoint(ᶠu³ʲs.:($$j)[colidx]) * ᶠgradᵥ(Y.c.sgsʲs.:($$j).h_tot[colidx]), ) - @. Yₜ.c.sgsʲs.:($$j).q_tot[colidx] -= ᶜinterp( + @. Yₜ.c.sgsʲs.:($$j).q_tot[colidx] -= LBC( adjoint(ᶠu³ʲs.:($$j)[colidx]) * ᶠgradᵥ(Y.c.sgsʲs.:($$j).q_tot[colidx]), )