From df7eac27c1e547218bd0e768fc6a51c1fba7e082 Mon Sep 17 00:00:00 2001
From: Simone Silvestri <33547697+simone-silvestri@users.noreply.github.com>
Date: Fri, 21 Jul 2023 11:12:34 -0400
Subject: [PATCH] bugfix

---
 src/dynamics/vertical_advection.jl | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/dynamics/vertical_advection.jl b/src/dynamics/vertical_advection.jl
index c74eb8dd2..8b9e7ec38 100644
--- a/src/dynamics/vertical_advection.jl
+++ b/src/dynamics/vertical_advection.jl
@@ -74,8 +74,8 @@ end
 @inline reconstructed_at_face(::UpwindVerticalAdvection{NF, 3}, u, c⁻⁻, c⁻, c⁺, c⁺⁺) where NF = ifelse(u > 0, (2c⁻⁻ + 5c⁻ - c⁺  ) / 6,
                                                                                                               (-c⁻  + 5c⁺ + 2c⁺⁺) / 6)
 
-@inline reconstructed_at_face(::CenteredVerticalAdvection{NF, 2}, u, c⁻⁻, c⁻, c⁺, c⁺⁺) = (c⁻ + c⁺) / 2
-@inline reconstructed_at_face(::CenteredVerticalAdvection{NF, 3}, u, c⁻⁻, c⁻, c⁺, c⁺⁺) = (- c⁻⁻ + 7c⁻ + 7c⁺ - c⁺⁺) / 12
+@inline reconstructed_at_face(::CenteredVerticalAdvection{NF, 2}, u, c⁻⁻, c⁻, c⁺, c⁺⁺) where NF = (c⁻ + c⁺) / 2
+@inline reconstructed_at_face(::CenteredVerticalAdvection{NF, 3}, u, c⁻⁻, c⁻, c⁺, c⁺⁺) where NF = (- c⁻⁻ + 7c⁻ + 7c⁺ - c⁺⁺) / 12
 
 # MULTI THREADED VERSION only writes into layer k
 function _vertical_advection!(  ξ_tend::Grid,           # tendency of quantity ξ at k