change boundary condition of subsidence

CliMA · Jun 27, 2024 · 47c6042 · 47c6042
1 parent 03fa877
commit 47c6042
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.