Try to fix moist gpu example

CliMA · Oct 4, 2023 · 1ec12b9 · 1ec12b9
1 parent 3a80801
commit 1ec12b9
Showing 1 changed file with 2 additions and 3 deletions.
diff --git a/src/prognostic_equations/vertical_diffusion_boundary_layer.jl b/src/prognostic_equations/vertical_diffusion_boundary_layer.jl
@@ -85,15 +85,14 @@ function vertical_diffusion_boundary_layer_tendency!(
     FT = eltype(Y)
     interior_uₕ = Fields.level(Y.c.uₕ, 1)
     interior_coordinates = Fields.level(Fields.coordinate_field(Y.c), 1)
-    z_surface =
-        Fields.level(Fields.coordinate_field(Y.f).z[colidx], Fields.half)
     ᶠp = ᶠρK_E = p.ᶠtemp_scalar
     @. ᶠp[colidx] = ᶠinterp(ᶜp[colidx])
+    ᶜΔz_surface = Fields.Δz_field(interior_coordinates)
     @. ᶠρK_E[colidx] =
         ᶠinterp(Y.c.ρ[colidx]) * eddy_diffusivity_coefficient(
             C_E,
             norm(interior_uₕ[colidx]),
-            interior_coordinates.z[colidx] - z_surface,
+            ᶜΔz_surface[colidx] / 2,
             ᶠp[colidx],
         )