More BC fixes

Author: charleskawczynski

test/Operators/finitedifference/opt_examples.jl

@@ -182,11 +182,12 @@ function alloc_test_operators_in_loops(cfield, ffield)
     end
 end
 function alloc_test_nested_expressions_1(cfield, ffield)
+    FT = Spaces.undertype(axes(cfield))
     (; fx, fy, fz, fϕ, fψ) = ffield
     (; cx, cy, cz, cϕ, cψ) = cfield
     ∇c = Operators.DivergenceF2C()
     wvec = Geometry.WVector
-    LB = Operators.LeftBiasedC2F(; bottom = Operators.SetValue(1))
+    LB = Operators.LeftBiasedC2F(; bottom = Operators.SetValue(FT(1)))
     @. cz = cx * cy * ∇c(wvec(LB(cy))) * ∇c(wvec(LB(cx))) * cϕ * cψ # Compile first
     p = @allocated begin
         @. cz = cx * cy * ∇c(wvec(LB(cy))) * ∇c(wvec(LB(cx))) * cϕ * cψ
@@ -208,12 +209,13 @@ function alloc_test_nested_expressions_2(cfield, ffield)
 end
 
 function alloc_test_nested_expressions_3(cfield, ffield)
+    FT = Spaces.undertype(axes(cfield))
     (; fx, fy, fz, fϕ, fψ) = ffield
     (; cx, cy, cz, cϕ, cψ) = cfield
     Ic = Operators.InterpolateF2C()
     ∇c = Operators.DivergenceF2C()
     wvec = Geometry.WVector
-    LB = Operators.LeftBiasedC2F(; bottom = Operators.SetValue(1))
+    LB = Operators.LeftBiasedC2F(; bottom = Operators.SetValue(FT(1)))
     #! format: off
     @. cz = cx * cy * ∇c(wvec(LB(Ic(fy) * cx))) * ∇c(wvec(LB(Ic(fy) * cx))) * cϕ * cψ # Compile first
     p = @allocated begin