Fix NaNs check, add more debug info

docs/src/dev/compute_convergence.jl

@@ -169,7 +169,9 @@ function compute_convergence!(
     float_str(x) = @sprintf "%.4f" x
     pow_str(x) = "10^{$(@sprintf "%.1f" log10(x))}"
     function si_str(x)
-        x in (0, Inf, -Inf, NaN) && return string(x)
+        if isnan(x) || x in (0, Inf, -Inf)
+            return string(x)
+        end
         exponent = floor(Int, log10(x))
         mantissa = x / 10.0^exponent
         return "$(float_str(mantissa)) \\times 10^{$exponent}"
@@ -266,6 +268,9 @@ function compute_convergence!(
         callback = scb_cur_avg_sol_and_err,
     )
     if any(isnan, plot2_values)
+        @show default_dt
+        @show count(isnan, plot2_values.u[end])
+        @show length(plot2_values.u[end])
         error("NaN found in plot2_values in problem $(test_name)")
     end
     out_dict[key1][key2]["plot2_values"] = plot2_values