Model not converging

[Francis Oussou]

Hi all,

Please can anyone suggest what to do if the model is not converging?

I increased the iteration and reduced the total residual but it still is not. Kindly find attached the run script and share your view.

Best,

Francis

#-----------------------------------------------------------------------------

# Set solver parameters

#-----------------------------------------------------------------------------

nalohou_run.Solver.TerrainFollowingGrid = True

nalohou_run.Solver = 'Richards'

nalohou_run.Solver.MaxIter = 25000

nalohou_run.OverlandFlowDiffusive = 0

nalohou_run.Solver.Nonlinear.MaxIter = 800

nalohou_run.Solver.Nonlinear.ResidualTol = #1e-10

nalohou_run.Solver.Nonlinear.EtaChoice = 'Walker1'

nalohou_run.Solver.Nonlinear.EtaChoice = 'EtaConstant'

nalohou_run.Solver.Nonlinear.EtaValue = 0.001

nalohou_run.Solver.Nonlinear.UseJacobian = False

nalohou_run.Solver.Nonlinear.DerivativeEpsilon = 1e-16

nalohou_run.Solver.Nonlinear.StepTol = 1e-30

nalohou_run.Solver.Nonlinear.Globalization = 'LineSearch'

nalohou_run.Solver.Linear.KrylovDimension = 50

nalohou_run.Solver.Linear.MaxRestart = 3

nalohou_run.Solver.Linear.Preconditioner = 'PFMG'

nalohou_run.Solver.Linear.Preconditioner.PFMG.MaxIter = 5

nalohou_run.Solver.Linear.Preconditioner.PFMG.Smoother = 'RBGaussSeidelNonSymmetric'

nalohou_run.Solver.Linear.Preconditioner.PFMG.NumPreRelax = 1

nalohou_run.Solver.Linear.Preconditioner.PFMG.NumPostRelax = 1

nalohou_run.Solver.Drop = 1E-20

nalohou_run.Solver.AbsTol = 1E-9

#nalohou_run.Solver.WriteSiloSubsurfData = True

#nalohou_run.Solver.WriteSiloPressure = True

#nalohou_run.Solver.WriteSiloSaturation = True

#nalohou_run.Solver.WriteSiloConcentration = True

#spinup key  ## Read Readme washita

#nalohou_run.OverlandFlowSpinUp = 1

#nalohou_run.OverlandSpinupDampP1 = 1.0

#nalohou_run.OverlandSpinupDampP2 = 0.00001