diff --git a/Project.toml b/Project.toml
index 791aa80534..96c0d0c7ea 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "Oceananigans"
 uuid = "9e8cae18-63c1-5223-a75c-80ca9d6e9a09"
 authors = ["Climate Modeling Alliance and contributors"]
-version = "0.91.3"
+version = "0.91.4"
 
 [deps]
 Adapt = "79e6a3ab-5dfb-504d-930d-738a2a938a0e"
diff --git a/src/TurbulenceClosures/velocity_tracer_gradients.jl b/src/TurbulenceClosures/velocity_tracer_gradients.jl
index 78bedbdb81..be36e6996e 100644
--- a/src/TurbulenceClosures/velocity_tracer_gradients.jl
+++ b/src/TurbulenceClosures/velocity_tracer_gradients.jl
@@ -9,10 +9,10 @@
 
 # Off-diagonal
 @inline ∂x_v(i, j, k, grid, v) = ∂xᶠᶠᶜ(i, j, k, grid, v)
-@inline ∂x_w(i, j, k, grid, w) = ∂xᶠᶜᶜ(i, j, k, grid, w)
+@inline ∂x_w(i, j, k, grid, w) = ∂xᶠᶜᶠ(i, j, k, grid, w)
 
 @inline ∂y_u(i, j, k, grid, u) = ∂yᶠᶠᶜ(i, j, k, grid, u)
-@inline ∂y_w(i, j, k, grid, w) = ∂yᶜᶠᶜ(i, j, k, grid, w)
+@inline ∂y_w(i, j, k, grid, w) = ∂yᶜᶠᶠ(i, j, k, grid, w)
 
 @inline ∂z_u(i, j, k, grid, u) = ∂zᶠᶜᶠ(i, j, k, grid, u)
 @inline ∂z_v(i, j, k, grid, v) = ∂zᶜᶠᶠ(i, j, k, grid, v)