diff --git a/src/utils/abbreviations.jl b/src/utils/abbreviations.jl
index f4aea58534..8b686bc2e9 100644
--- a/src/utils/abbreviations.jl
+++ b/src/utils/abbreviations.jl
@@ -30,10 +30,11 @@ const ᶜadvdivᵥ = Operators.DivergenceF2C(
     top = Operators.SetValue(CT3(0)),
 )
 
-# Subsidence has zero tendency at the top, and has no flux at the bottom.
+# Subsidence has extrapolated tendency at the top, and has no flux at the bottom.
+# TODO: This is not accurate and causes some issues at the domain top.
 const ᶜsubdivᵥ = Operators.DivergenceF2C(
     bottom = Operators.SetValue(CT3(0)),
-    top = Operators.SetDivergence(0),
+    top = Operators.Extrapolate(),
 )
 
 # Precipitation has no flux at the top, but it has free outflow at the bottom.