Add maxiter arg to find_zero

Author: imreddyTeja

src/Meshes/interval.jl

@@ -296,6 +296,7 @@ function IntervalMesh(
         RootSolvers.SecantMethod(guess₋, guess₊),
         RootSolvers.CompactSolution(),
         RootSolvers.ResidualTolerance(FT_solve(tol)),
+        maxiters = 10_000,
     )
     if h_bottom_sol.converged !== true
         error(
@@ -316,6 +317,7 @@ function IntervalMesh(
         RootSolvers.SecantMethod(guess₋, guess₊),
         RootSolvers.CompactSolution(),
         RootSolvers.ResidualTolerance(FT_solve(tol)),
+        maxiters = 10_000,
     )
     if h_top_sol.converged !== true
         error(
@@ -410,6 +412,7 @@ function IntervalMesh(
         RootSolvers.NewtonsMethodAD(FT_solve(1.0)),
         RootSolvers.CompactSolution(),
         RootSolvers.ResidualTolerance(FT_solve(tol)),
+        maxiters = 10_000,
     )
     if γ_sol.converged !== true
         error(

src/Operators/integrals.jl

@@ -73,6 +73,7 @@ function column_integral_indefinite!(
             RootSolvers.NewtonsMethodAD(ϕ_prev),
             RootSolvers.CompactSolution(),
             RootSolvers.RelativeSolutionTolerance(rtol),
+            maxiters = 10_000,
         )
         ClimaComms.@assert device converged "unable to integrate through \
                                              z = $z with rtol set to $rtol"