make surface temperature dependent on surface height

CliMA · Apr 20, 2024 · 98ddd36 · 98ddd36
1 parent 61c6dc0
commit 98ddd36
Showing 1 changed file with 6 additions and 5 deletions.
diff --git a/src/surface_conditions/surface_conditions.jl b/src/surface_conditions/surface_conditions.jl
@@ -335,14 +335,15 @@ function surface_state_to_conditions(
 end
 
 function surface_temperature(::ZonallySymmetricSST, coordinates)
-    FT = eltype(coordinates.lat)
-    T = FT(271) + FT(29) * exp(-coordinates.lat^2 / (2 * 26^2))
+    (; lat, z) = coordinates
+    FT = eltype(lat)
+    T = FT(271) + FT(29) * exp(-coordinates.lat^2 / (2 * 26^2)) - FT(6.5e-3) * z
     return T
 end
 
 function surface_temperature(::ZonallyAsymmetricSST, coordinates)
-    (; lat, long) = coordinates
-    FT = eltype(coordinates.lat)
+    (; lat, long, z) = coordinates
+    FT = eltype(lat)
     #Assume a surface temperature that varies with both longitude and latitude, Neale and Hoskins, 2021
     T =
         (
@@ -355,7 +356,7 @@ function surface_temperature(::ZonallyAsymmetricSST, coordinates)
                 FT(3) * cosd(long + FT(90)) * cospi(FT(0.5) * lat / FT(30))^2 +
                 FT(0)
             ) : FT(0)
-        )
+        ) - FT(6.5e-3) * z
     return T
 end