Merge pull request #3119 from CliMA/zs/subsidence

change boundary condition of subsidence
CliMA · Jun 27, 2024 · 3b31125 · 3b31125
2 parents 03fa877 + 47c6042
commit 3b31125
Showing 1 changed file with 3 additions and 2 deletions.
diff --git 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.