richards_box_proctest_vardz.Solver.MaxIter = 50000 richards_box_proctest_vardz.Solver.Nonlinear.MaxIter = 100 richards_box_proctest_vardz.Solver.Nonlinear.ResidualTol = 1e-7 richards_box_proctest_vardz.Solver.Nonlinear.EtaChoice = 'EtaConstant' richards_box_proctest_vardz.Solver.Nonlinear.EtaValue = 1e-2 # used to test analytical and FD jacobian combinations richards_box_proctest_vardz.Solver.Nonlinear.UseJacobian = True #pfset Solver.Nonlinear.UseJacobian False richards_box_proctest_vardz.Solver.Nonlinear.DerivativeEpsilon = 1e-14 richards_box_proctest_vardz.Solver.Linear.KrylovDimension = 100 # used to test different linear preconditioners richards_box_proctest_vardz.Solver.Linear.Preconditioner = 'PFMG' ##pfset Solver.Linear.Preconditioner PFMGOctree ##pfset Solver.Linear.Preconditioner MGSemi ##pfset Solver.Linear.Preconditioner SMG richards_box_proctest_vardz.UseClustering = False #----------------------------------------------------------------------------- # Run and Unload the ParFlow output files #----------------------------------------------------------------------------- richards_box_proctest_vardz.run()
#-----------------------------------------------------------------------------
rbpv.Solver = 'Richards'
rbpv.Solver.MaxIter = 50000
rbpv.Solver.Nonlinear.MaxIter = 100
rbpv.Solver.Nonlinear.ResidualTol = 1e-7
rbpv.Solver.Nonlinear.EtaChoice = 'EtaConstant'
rbpv.Solver.Nonlinear.EtaValue = 1e-2
# used to test analytical and FD jacobian combinations
rbpv.Solver.Nonlinear.UseJacobian = True
rbpv.Solver.Nonlinear.DerivativeEpsilon = 1e-14
rbpv.Solver.Linear.KrylovDimension = 100
# used to test different linear preconditioners
rbpv.Solver.Linear.Preconditioner = 'PFMG'
rbpv.UseClustering = False
#-----------------------------------------------------------------------------
# Run and Unload the ParFlow output files
#-----------------------------------------------------------------------------
dir_name = get_absolute_path('test_output/rbpv')
mkdir(dir_name)
rbpv.run(working_directory=dir_name)
# This runs the basic default_richards test case.
# This run, as written in this input file, should take
# 3 nonlinear iterations.
#
# Import the ParFlow TCL package
#
from parflow import Run
default_richards_nocluster = Run("default_richards_nocluster", __file__)
default_richards_nocluster.FileVersion = 4
default_richards_nocluster.UseClustering = False
default_richards_nocluster.Process.Topology.P = 1
default_richards_nocluster.Process.Topology.Q = 1
default_richards_nocluster.Process.Topology.R = 1
#---------------------------------------------------------
# Computational Grid
#---------------------------------------------------------
default_richards_nocluster.ComputationalGrid.Lower.X = -10.0
default_richards_nocluster.ComputationalGrid.Lower.Y = 10.0
default_richards_nocluster.ComputationalGrid.Lower.Z = 1.0
default_richards_nocluster.ComputationalGrid.DX = 8.8888888888888893
default_richards_nocluster.ComputationalGrid.DY = 10.666666666666666
default_richards_nocluster.ComputationalGrid.DZ = 1.0
default_richards_nocluster.ComputationalGrid.NX = 10
default_richards_nocluster.ComputationalGrid.NY = 10
#-----------------------------------------------------------------------------
# This runs the basic default_richards test case.
# This run, as written in this input file, should take
# 3 nonlinear iterations.
#-----------------------------------------------------------------------------
from parflow import Run
from parflow.tools.fs import mkdir, get_absolute_path
drich_n = Run("default_richards_nocluster", __file__)
#---------------------------------------------------------
drich_n.FileVersion = 4
drich_n.UseClustering = False
drich_n.Process.Topology.P = 1
drich_n.Process.Topology.Q = 1
drich_n.Process.Topology.R = 1
#---------------------------------------------------------
# Computational Grid
#---------------------------------------------------------
drich_n.ComputationalGrid.Lower.X = -10.0
drich_n.ComputationalGrid.Lower.Y = 10.0
drich_n.ComputationalGrid.Lower.Z = 1.0
drich_n.ComputationalGrid.DX = 8.8888888888888893
drich_n.ComputationalGrid.DY = 10.666666666666666
