diff --git a/cicecore/cicedynB/general/ice_init.F90 b/cicecore/cicedynB/general/ice_init.F90
index e0ebdfbed..167e4056d 100644
--- a/cicecore/cicedynB/general/ice_init.F90
+++ b/cicecore/cicedynB/general/ice_init.F90
@@ -419,7 +419,7 @@ subroutine input_data
       deltaminEVP = 1e-11_dbl_kind ! minimum delta for viscosities (EVP, Hunke 2001)
       deltaminVP  = 2e-9_dbl_kind  ! minimum delta for viscosities (VP, Hibler 1979)
       capping_method  = 'max'  ! method for capping of viscosities (max=Hibler 1979,sum=Kreyscher2000)
-      maxits_nonlin = 4        ! max nb of iteration for nonlinear solver
+      maxits_nonlin = 10       ! max nb of iteration for nonlinear solver
       precond = 'pgmres'       ! preconditioner for fgmres: 'ident' (identity), 'diag' (diagonal),
                                ! 'pgmres' (Jacobi-preconditioned GMRES)
       dim_fgmres = 50          ! size of fgmres Krylov subspace
@@ -431,7 +431,7 @@ subroutine input_data
       monitor_pgmres = .false. ! print pgmres residual norm
       ortho_type = 'mgs'       ! orthogonalization procedure 'cgs' or 'mgs'
       reltol_nonlin = 1e-8_dbl_kind ! nonlinear stopping criterion: reltol_nonlin*res(k=0)
-      reltol_fgmres = 1e-2_dbl_kind ! fgmres stopping criterion: reltol_fgmres*res(k)
+      reltol_fgmres = 1e-1_dbl_kind ! fgmres stopping criterion: reltol_fgmres*res(k)
       reltol_pgmres = 1e-6_dbl_kind ! pgmres stopping criterion: reltol_pgmres*res(k)
       algo_nonlin = 'picard'        ! nonlinear algorithm: 'picard' (Picard iteration), 'anderson' (Anderson acceleration)
       fpfunc_andacc = 1        ! fixed point function for Anderson acceleration:
diff --git a/configuration/scripts/ice_in b/configuration/scripts/ice_in
index e5cbd6fcc..8262f34ec 100644
--- a/configuration/scripts/ice_in
+++ b/configuration/scripts/ice_in
@@ -167,7 +167,7 @@
     kridge          = 1
     ktransport      = 1
     ssh_stress      = 'geostrophic'
-    maxits_nonlin   = 4
+    maxits_nonlin   = 10
     precond         = 'pgmres'
     dim_fgmres       = 50
     dim_pgmres       = 5
@@ -178,7 +178,7 @@
     monitor_pgmres  = .false.
     ortho_type      = 'mgs'
     reltol_nonlin   = 1e-8
-    reltol_fgmres   = 1e-2
+    reltol_fgmres   = 1e-1
     reltol_pgmres   = 1e-6
     algo_nonlin     = 'picard'
     use_mean_vrel   = .true.
diff --git a/doc/source/user_guide/ug_case_settings.rst b/doc/source/user_guide/ug_case_settings.rst
index 04a275815..b2dd54c33 100644
--- a/doc/source/user_guide/ug_case_settings.rst
+++ b/doc/source/user_guide/ug_case_settings.rst
@@ -469,7 +469,7 @@ dynamics_nml
    "``k1``", "real", "1st free parameter for landfast parameterization", "7.5"
    "``k2``", "real", "2nd free parameter (N/m\ :math:`^3`) for landfast parameterization", "15.0"
    "``maxits_fgmres``", "integer", "maximum number of restarts for FGMRES solver", "1"
-   "``maxits_nonlin``", "integer", "maximum number of nonlinear iterations for VP solver", "4"
+   "``maxits_nonlin``", "integer", "maximum number of nonlinear iterations for VP solver", "10"
    "``maxits_pgmres``", "integer", "maximum number of restarts for PGMRES preconditioner", "1"
    "``monitor_fgmres``", "logical", "write velocity norm at each FGMRES iteration", "``.false.``"
    "``monitor_nonlin``", "logical", "write velocity norm at each nonlinear iteration", "``.false.``"
@@ -482,7 +482,7 @@ dynamics_nml
    "", "``ident``", "Don't use a preconditioner for the FGMRES solver", ""
    "", "``pgmres``", "Use GMRES as preconditioner for FGMRES solver", ""
    "``Pstar``", "real", "constant in Hibler strength formula (N/m\ :math:`^2`)", "2.75e4"
-   "``reltol_fgmres``", "real", "relative tolerance for FGMRES solver", "1e-2"
+   "``reltol_fgmres``", "real", "relative tolerance for FGMRES solver", "1e-1"
    "``reltol_nonlin``", "real", "relative tolerance for nonlinear solver", "1e-8"
    "``reltol_pgmres``", "real", "relative tolerance for PGMRES preconditioner", "1e-6"
    "``revised_evp``", "logical", "use revised EVP formulation", "``.false.``"