Make sure we reach error

Author: charleskawczynski

examples/hybrid/driver.jl

@@ -1,21 +1,21 @@
 # Test-specific definitions (may be overwritten in each test case file)
 # TODO: Allow some of these to be enironment variables or command line arguments
 
-import ClimaCore
-ClimaCore.DebugOnly.call_post_op_callback() = true
-function ClimaCore.DebugOnly.post_op_callback(result, args...; kwargs...)
-    has_nans = result isa Number ? isnan(result) : any(isnan, parent(result))
-    has_inf = result isa Number ? isinf(result) : any(isinf, parent(result))
-    if has_nans || has_inf
-        has_nans && println("NaNs found!")
-        has_inf && println("Infs found!")
-        # Let's define the stack trace so that we know where this came from
-        st = stacktrace()
-        ClimaCore.DebugOnly.print_depth_limited_stack_trace(st; maxtypedepth=2)
+# import ClimaCore
+# ClimaCore.DebugOnly.call_post_op_callback() = true
+# function ClimaCore.DebugOnly.post_op_callback(result, args...; kwargs...)
+#     has_nans = result isa Number ? isnan(result) : any(isnan, parent(result))
+#     has_inf = result isa Number ? isinf(result) : any(isinf, parent(result))
+#     if has_nans || has_inf
+#         has_nans && println("NaNs found!")
+#         has_inf && println("Infs found!")
+#         # Let's define the stack trace so that we know where this came from
+#         st = stacktrace()
+#         ClimaCore.DebugOnly.print_depth_limited_stack_trace(st; maxtypedepth=2)
 
-        error("Halting execution.")
-    end
-end
+#         error("Halting execution.")
+#     end
+# end
 
 upwinding_mode = :none
 horizontal_mesh = nothing # must be object of type AbstractMesh

examples/hybrid/sphere/held_suarez_rhoe.jl

@@ -10,8 +10,7 @@ horizontal_mesh = cubed_sphere_mesh(; radius = R, h_elem = 4)
 npoly = 4
 z_max = FT(30e3)
 z_elem = 10
-# t_end = FT(60 * 60 * 24 * 10)
-t_end = FT(60 * 60 * 24)
+t_end = FT(60 * 60 * 24 * 10)
 dt = FT(400)
 dt_save_to_sol = FT(60 * 60 * 24)
 dt_save_to_disk = FT(0) # 0 means don't save to disk
@@ -27,7 +26,7 @@ additional_cache(ᶜlocal_geometry, ᶠlocal_geometry, dt) = merge(
     held_suarez_cache(ᶜlocal_geometry),
 )
 function additional_tendency!(Yₜ, Y, p, t)
-    @show t
+    t > 805400.0 && @show t
     hyperdiffusion_tendency!(Yₜ, Y, p, t)
     sponge && rayleigh_sponge_tendency!(Yₜ, Y, p, t)
     held_suarez_tendency!(Yₜ, Y, p, t)