wip

CliMA · Oct 10, 2023 · 0e2c067 · 0e2c067
1 parent 503265d
commit 0e2c067
Showing 1 changed file with 3 additions and 3 deletions.
diff --git 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]),
                 )