clm_samrai.Mannings.GeomNames = 'domain' clm_samrai.Mannings.Geom.domain.Value = 5.52e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- clm_samrai.PhaseSources.water.Type = 'Constant' clm_samrai.PhaseSources.water.GeomNames = 'domain' clm_samrai.PhaseSources.water.Geom.domain.Value = 0.0 # #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- # clm_samrai.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- # clm_samrai.Solver = 'Richards' clm_samrai.Solver.MaxIter = 500 # clm_samrai.Solver.Nonlinear.MaxIter = 15 clm_samrai.Solver.Nonlinear.ResidualTol = 1e-9 clm_samrai.Solver.Nonlinear.EtaChoice = 'EtaConstant' clm_samrai.Solver.Nonlinear.EtaValue = 0.01 clm_samrai.Solver.Nonlinear.UseJacobian = False clm_samrai.Solver.Nonlinear.DerivativeEpsilon = 1e-12 clm_samrai.Solver.Nonlinear.StepTol = 1e-20
octree.Geom.domain.ICPressure.RefGeom = 'domain' octree.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- octree.PhaseSources.water.Type = 'Constant' octree.PhaseSources.water.GeomNames = 'background' octree.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- octree.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- octree.Solver = 'Richards' octree.Solver.MaxIter = 5 octree.Solver.Nonlinear.MaxIter = 10 octree.Solver.Nonlinear.ResidualTol = 1e-9 octree.Solver.Nonlinear.EtaChoice = 'EtaConstant' octree.Solver.Nonlinear.EtaValue = 1e-5 octree.Solver.Nonlinear.UseJacobian = True octree.Solver.Nonlinear.DerivativeEpsilon = 1e-2
octree_simple.Geom.domain.ICPressure.RefGeom = 'domain' octree_simple.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- octree_simple.PhaseSources.water.Type = 'Constant' octree_simple.PhaseSources.water.GeomNames = 'background' octree_simple.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- octree_simple.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- octree_simple.Solver = 'Richards' octree_simple.Solver.MaxIter = 5 octree_simple.Solver.Nonlinear.MaxIter = 10 octree_simple.Solver.Nonlinear.ResidualTol = 1e-9 octree_simple.Solver.Nonlinear.EtaChoice = 'EtaConstant' octree_simple.Solver.Nonlinear.EtaValue = 1e-5 octree_simple.Solver.Nonlinear.UseJacobian = True octree_simple.Solver.Nonlinear.DerivativeEpsilon = 1e-2 octree_simple.Solver.Linear.KrylovDimension = 10
LWvdz.Mannings.GeomNames = 'domain' LWvdz.Mannings.Geom.domain.Value = 0.00005 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- LWvdz.PhaseSources.water.Type = 'Constant' LWvdz.PhaseSources.water.GeomNames = 'domain' LWvdz.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- LWvdz.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- LWvdz.Solver = 'Richards' LWvdz.Solver.MaxIter = 2500 LWvdz.Solver.TerrainFollowingGrid = True LWvdz.Solver.Nonlinear.MaxIter = 80 LWvdz.Solver.Nonlinear.ResidualTol = 1e-5 LWvdz.Solver.Nonlinear.EtaValue = 0.001 LWvdz.Solver.PrintSubsurf = False
var_dz_1D.Geom.domain.ICPressure.RefGeom = 'domain' var_dz_1D.Geom.domain.ICPressure.RefPatch = 'top' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- var_dz_1D.PhaseSources.water.Type = 'Constant' var_dz_1D.PhaseSources.water.GeomNames = 'domain' var_dz_1D.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- var_dz_1D.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- var_dz_1D.Solver = 'Richards' var_dz_1D.Solver.MaxIter = 2500 var_dz_1D.Solver.Nonlinear.MaxIter = 200 var_dz_1D.Solver.Nonlinear.ResidualTol = 1e-9 var_dz_1D.Solver.Nonlinear.EtaChoice = 'Walker1' var_dz_1D.Solver.Nonlinear.EtaValue = 1e-5 var_dz_1D.Solver.Nonlinear.UseJacobian = True #pfset Solver.Nonlinear.UseJacobian False var_dz_1D.Solver.Nonlinear.DerivativeEpsilon = 1e-10
BC_test.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- BC_test.PhaseSources.water.Type = 'Constant' BC_test.PhaseSources.water.GeomNames = 'domain' BC_test.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- BC_test.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- BC_test.Solver = 'Richards' BC_test.Solver.MaxIter = 2500 BC_test.Solver.Nonlinear.MaxIter = 300 BC_test.Solver.Nonlinear.ResidualTol = 1e-4 BC_test.Solver.Nonlinear.EtaChoice = 'Walker1' BC_test.Solver.Nonlinear.EtaChoice = 'EtaConstant' BC_test.Solver.Nonlinear.EtaValue = 0.001 BC_test.Solver.Nonlinear.UseJacobian = False
dover.Mannings.Geom.right.Value = 5.e-6 dover.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- dover.PhaseSources.water.Type = 'Constant' dover.PhaseSources.water.GeomNames = 'domain' dover.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- dover.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- dover.Solver = 'Richards' dover.Solver.MaxIter = 2500 dover.Solver.Nonlinear.MaxIter = 20 dover.Solver.Nonlinear.ResidualTol = 1e-9 dover.Solver.Nonlinear.EtaChoice = 'EtaConstant' dover.Solver.Nonlinear.EtaValue = 0.01 dover.Solver.Nonlinear.UseJacobian = True dover.Solver.Nonlinear.DerivativeEpsilon = 1e-8 dover.Solver.Nonlinear.StepTol = 1e-20
wby.Mannings.Geom.right.Value = 5.e-6 wby.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- wby.PhaseSources.water.Type = 'Constant' wby.PhaseSources.water.GeomNames = 'domain' wby.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- wby.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- wby.Solver = 'Richards' wby.Solver.MaxIter = 2500 wby.Solver.AbsTol = 1E-12 wby.Solver.Nonlinear.MaxIter = 300 wby.Solver.Nonlinear.ResidualTol = 1e-12 wby.Solver.Nonlinear.EtaChoice = 'Walker1' wby.Solver.Nonlinear.EtaChoice = 'EtaConstant' wby.Solver.Nonlinear.EtaValue = 0.001 wby.Solver.Nonlinear.UseJacobian = False
#----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- dsingle.PhaseSources.water.Type = 'Constant' dsingle.PhaseSources.water.GeomNames = 'background' dsingle.PhaseSources.water.Geom.background.Value = 0.0 dsingle.PhaseConcen.water.tce.Type = 'Constant' dsingle.PhaseConcen.water.tce.GeomNames = 'concen_region' dsingle.PhaseConcen.water.tce.Geom.concen_region.Value = 0.8 dsingle.Solver.WriteSiloSubsurfData = True dsingle.Solver.WriteSiloPressure = True dsingle.Solver.WriteSiloSaturation = True dsingle.Solver.WriteSiloConcentration = True #----------------------------------------------------------------------------- # The Solver Impes MaxIter default value changed so to get previous # results we need to set it back to what it was #----------------------------------------------------------------------------- dsingle.Solver.MaxIter = 5 dsingle.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Run and Unload the ParFlow output files #----------------------------------------------------------------------------- dsingle.run()
rich_fbx.Geom.domain.ICPressure.RefGeom = 'domain' rich_fbx.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- rich_fbx.PhaseSources.water.Type = 'Constant' rich_fbx.PhaseSources.water.GeomNames = 'domain' rich_fbx.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- rich_fbx.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- rich_fbx.Solver = 'Richards' rich_fbx.Solver.MaxIter = 50000 rich_fbx.Solver.Nonlinear.MaxIter = 100 rich_fbx.Solver.Nonlinear.ResidualTol = 1e-6 rich_fbx.Solver.Nonlinear.EtaChoice = 'EtaConstant' rich_fbx.Solver.Nonlinear.EtaValue = 1e-2 rich_fbx.Solver.Nonlinear.UseJacobian = True rich_fbx.Solver.Nonlinear.DerivativeEpsilon = 1e-12
vardz.Geom.domain.ICPressure.RefGeom = 'domain' vardz.Geom.domain.ICPressure.RefPatch = 'top' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- vardz.PhaseSources.water.Type = 'Constant' vardz.PhaseSources.water.GeomNames = 'domain' vardz.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- vardz.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- vardz.Solver = 'Richards' vardz.Solver.MaxIter = 2500 vardz.Solver.Nonlinear.MaxIter = 200 vardz.Solver.Nonlinear.ResidualTol = 1e-9 vardz.Solver.Nonlinear.EtaChoice = 'Walker1' vardz.Solver.Nonlinear.EtaValue = 1e-5 vardz.Solver.Nonlinear.UseJacobian = True vardz.Solver.Nonlinear.DerivativeEpsilon = 1e-10
small_domain.Geom.infiltration.ICPressure.RefPatch = 'infiltration' small_domain.Geom.infiltration.ICPressure.RefGeom = 'domain' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- small_domain.PhaseSources.water.Type = 'Constant' small_domain.PhaseSources.water.GeomNames = 'background' small_domain.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- small_domain.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- small_domain.Solver = 'Richards' small_domain.Solver.MaxIter = 10000 small_domain.Solver.Nonlinear.MaxIter = 15 small_domain.Solver.Nonlinear.ResidualTol = 1e-9 small_domain.Solver.Nonlinear.StepTol = 1e-9 small_domain.Solver.Nonlinear.EtaValue = 1e-5 small_domain.Solver.Nonlinear.UseJacobian = True small_domain.Solver.Nonlinear.DerivativeEpsilon = 1e-7 small_domain.Solver.Linear.KrylovDimension = 25
overland_flux.Mannings.Geom.right.Value = 5.e-6 overland_flux.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- overland_flux.PhaseSources.water.Type = 'Constant' overland_flux.PhaseSources.water.GeomNames = 'domain' overland_flux.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- overland_flux.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- overland_flux.Solver = 'Richards' overland_flux.Solver.MaxIter = 2500 overland_flux.OverlandFlowDiffusive = 0 overland_flux.Solver.Nonlinear.MaxIter = 20 overland_flux.Solver.Nonlinear.ResidualTol = 1e-9 overland_flux.Solver.Nonlinear.EtaChoice = 'EtaConstant' overland_flux.Solver.Nonlinear.EtaValue = 0.01
dover_1.Mannings.Geom.right.Value = 5.e-6 dover_1.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- dover_1.PhaseSources.water.Type = 'Constant' dover_1.PhaseSources.water.GeomNames = 'domain' dover_1.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- dover_1.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- dover_1.Solver = 'Richards' dover_1.Solver.MaxIter = 2500 dover_1.Solver.Nonlinear.MaxIter = 20 dover_1.Solver.Nonlinear.ResidualTol = 1e-9 dover_1.Solver.Nonlinear.EtaChoice = 'EtaConstant' dover_1.Solver.Nonlinear.EtaValue = 0.01 dover_1.Solver.Nonlinear.UseJacobian = False dover_1.Solver.Nonlinear.DerivativeEpsilon = 1e-8
default_overland_pfmg_octree_jac.Mannings.Geom.right.Value = 5.e-6 default_overland_pfmg_octree_jac.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- default_overland_pfmg_octree_jac.PhaseSources.water.Type = 'Constant' default_overland_pfmg_octree_jac.PhaseSources.water.GeomNames = 'domain' default_overland_pfmg_octree_jac.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- default_overland_pfmg_octree_jac.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- default_overland_pfmg_octree_jac.Solver = 'Richards' default_overland_pfmg_octree_jac.Solver.MaxIter = 2500 default_overland_pfmg_octree_jac.Solver.Nonlinear.MaxIter = 20 default_overland_pfmg_octree_jac.Solver.Nonlinear.ResidualTol = 1e-9 default_overland_pfmg_octree_jac.Solver.Nonlinear.EtaChoice = 'EtaConstant' default_overland_pfmg_octree_jac.Solver.Nonlinear.EtaValue = 0.01 default_overland_pfmg_octree_jac.Solver.Nonlinear.UseJacobian = True default_overland_pfmg_octree_jac.Solver.Nonlinear.DerivativeEpsilon = 1e-8
LW_Test.Patch.z_upper.BCPressure.rec.Value = 0.0000 #---------------------------------------------------------------- # Outputs # --------------------------------------------------------------- LW_Test.Solver.PrintSubsurfData = False LW_Test.Solver.PrintPressure = True LW_Test.Solver.PrintSaturation = True LW_Test.Solver.PrintMask = True #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- LW_Test.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- LW_Test.Solver.TerrainFollowingGrid = True LW_Test.Solver.Nonlinear.VariableDz = False LW_Test.Solver.MaxIter = 25000 LW_Test.Solver.Drop = 1E-20 LW_Test.Solver.AbsTol = 1E-8 LW_Test.Solver.MaxConvergenceFailures = 8 LW_Test.Solver.Nonlinear.MaxIter = 80 LW_Test.Solver.Nonlinear.ResidualTol = 1e-6
default_richards_with_silopmpio.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- default_richards_with_silopmpio.PhaseSources.water.Type = 'Constant' default_richards_with_silopmpio.PhaseSources.water.GeomNames = 'background' default_richards_with_silopmpio.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- default_richards_with_silopmpio.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- default_richards_with_silopmpio.Solver = 'Richards' default_richards_with_silopmpio.Solver.MaxIter = 5 default_richards_with_silopmpio.Solver.Nonlinear.MaxIter = 10 default_richards_with_silopmpio.Solver.Nonlinear.ResidualTol = 1e-9 default_richards_with_silopmpio.Solver.Nonlinear.EtaChoice = 'EtaConstant' default_richards_with_silopmpio.Solver.Nonlinear.EtaValue = 1e-5 default_richards_with_silopmpio.Solver.Nonlinear.UseJacobian = True default_richards_with_silopmpio.Solver.Nonlinear.DerivativeEpsilon = 1e-2
sabino.Geom.domain.ICPressure.Value = -10 #---------------------------------------------------------------- # Outputs # ------------------------------------------------------------ sabino.Solver.PrintSubsurf = True sabino.Solver.PrintPressure = True sabino.Solver.PrintSaturation = True sabino.Solver.PrintMask = True #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- sabino.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- sabino.Solver = 'Richards' sabino.Solver.TerrainFollowingGrid = True #----------------------------------------------------------------------------- # Setting up vertical layering #----------------------------------------------------------------------------- sabino.Solver.Nonlinear.VariableDz = True sabino.dzScale.GeomNames = 'domain' sabino.dzScale.Type = 'nzList'
clm.Mannings.GeomNames = 'domain' clm.Mannings.Geom.domain.Value = 1e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- clm.PhaseSources.water.Type = 'Constant' clm.PhaseSources.water.GeomNames = 'domain' clm.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- clm.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- clm.Solver = 'Richards' clm.Solver.MaxIter = 90000 clm.Solver.Nonlinear.MaxIter = 100 clm.Solver.Nonlinear.ResidualTol = 1e-5 clm.Solver.Nonlinear.EtaChoice = 'Walker1' clm.Solver.Nonlinear.EtaValue = 0.01 clm.Solver.Nonlinear.UseJacobian = True clm.Solver.Nonlinear.DerivativeEpsilon = 1e-12 clm.Solver.Nonlinear.StepTol = 1e-30
overland.Mannings.GeomNames = 'domain' overland.Mannings.Geom.domain.Value = 3.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- overland.PhaseSources.water.Type = 'Constant' overland.PhaseSources.water.GeomNames = 'domain' overland.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- overland.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- overland.Solver = 'Richards' overland.Solver.MaxIter = 2500 overland.Solver.Nonlinear.MaxIter = 50 overland.Solver.Nonlinear.ResidualTol = 1e-9 overland.Solver.Nonlinear.EtaChoice = 'EtaConstant' overland.Solver.Nonlinear.EtaValue = 0.01 overland.Solver.Nonlinear.UseJacobian = False overland.Solver.Nonlinear.DerivativeEpsilon = 1e-15
drich.Geom.domain.ICPressure.RefGeom = 'domain' drich.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- drich.PhaseSources.water.Type = 'Constant' drich.PhaseSources.water.GeomNames = 'background' drich.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- drich.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- drich.Solver = 'Richards' drich.Solver.MaxIter = 5 drich.Solver.Nonlinear.MaxIter = 10 drich.Solver.Nonlinear.ResidualTol = 1e-9 drich.Solver.Nonlinear.EtaChoice = 'EtaConstant' drich.Solver.Nonlinear.EtaValue = 1e-5 drich.Solver.Nonlinear.UseJacobian = True drich.Solver.Nonlinear.DerivativeEpsilon = 1e-2
richards_box_proctest_vardz.Geom.domain.ICPressure.RefGeom = 'domain' richards_box_proctest_vardz.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- richards_box_proctest_vardz.PhaseSources.water.Type = 'Constant' richards_box_proctest_vardz.PhaseSources.water.GeomNames = 'domain' richards_box_proctest_vardz.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- richards_box_proctest_vardz.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- richards_box_proctest_vardz.Solver = 'Richards' 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
PFCLM_SC.Mannings.GeomNames = 'domain' PFCLM_SC.Mannings.Geom.domain.Value = 2.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- PFCLM_SC.PhaseSources.water.Type = 'Constant' PFCLM_SC.PhaseSources.water.GeomNames = 'domain' PFCLM_SC.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- PFCLM_SC.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- PFCLM_SC.Solver = 'Richards' PFCLM_SC.Solver.MaxIter = 9000 PFCLM_SC.Solver.Nonlinear.MaxIter = 100 PFCLM_SC.Solver.Nonlinear.ResidualTol = 1e-5 PFCLM_SC.Solver.Nonlinear.EtaChoice = 'Walker1' PFCLM_SC.Solver.Nonlinear.EtaValue = 0.01 PFCLM_SC.Solver.Nonlinear.UseJacobian = False PFCLM_SC.Solver.Nonlinear.DerivativeEpsilon = 1e-12 PFCLM_SC.Solver.Nonlinear.StepTol = 1e-30
pfmg.Geom.domain.ICPressure.RefGeom = 'domain' pfmg.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- pfmg.PhaseSources.water.Type = 'Constant' pfmg.PhaseSources.water.GeomNames = 'background' pfmg.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- pfmg.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- pfmg.Solver = 'Richards' pfmg.Solver.MaxIter = 5 pfmg.Solver.Nonlinear.MaxIter = 10 pfmg.Solver.Nonlinear.ResidualTol = 1e-9 pfmg.Solver.Nonlinear.EtaChoice = 'EtaConstant' pfmg.Solver.Nonlinear.EtaValue = 1e-5 pfmg.Solver.Nonlinear.UseJacobian = True pfmg.Solver.Nonlinear.DerivativeEpsilon = 1e-2
water_balance_y.Mannings.Geom.right.Value = 5.e-6 water_balance_y.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- water_balance_y.PhaseSources.water.Type = 'Constant' water_balance_y.PhaseSources.water.GeomNames = 'domain' water_balance_y.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- water_balance_y.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- water_balance_y.Solver = 'Richards' water_balance_y.Solver.MaxIter = 2500 water_balance_y.Solver.AbsTol = 1E-12 water_balance_y.Solver.Nonlinear.MaxIter = 300 water_balance_y.Solver.Nonlinear.ResidualTol = 1e-12 water_balance_y.Solver.Nonlinear.EtaChoice = 'Walker1' water_balance_y.Solver.Nonlinear.EtaChoice = 'EtaConstant' water_balance_y.Solver.Nonlinear.EtaValue = 0.001 water_balance_y.Solver.Nonlinear.UseJacobian = False
rbp.Geom.domain.ICPressure.RefGeom = 'domain' rbp.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- rbp.PhaseSources.water.Type = 'Constant' rbp.PhaseSources.water.GeomNames = 'domain' rbp.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- rbp.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- rbp.Solver = 'Richards' rbp.Solver.MaxIter = 50000 rbp.Solver.Nonlinear.MaxIter = 100 rbp.Solver.Nonlinear.ResidualTol = 1e-6 rbp.Solver.Nonlinear.EtaChoice = 'EtaConstant' rbp.Solver.Nonlinear.EtaValue = 1e-2 rbp.Solver.Nonlinear.UseJacobian = True rbp.Solver.Nonlinear.DerivativeEpsilon = 1e-12
overland_tiltedV_DWE.Mannings.GeomNames = 'domain' overland_tiltedV_DWE.Mannings.Geom.domain.Value = 3.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- overland_tiltedV_DWE.PhaseSources.water.Type = 'Constant' overland_tiltedV_DWE.PhaseSources.water.GeomNames = 'domain' overland_tiltedV_DWE.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- overland_tiltedV_DWE.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- overland_tiltedV_DWE.Solver = 'Richards' overland_tiltedV_DWE.Solver.MaxIter = 2500 overland_tiltedV_DWE.Solver.Nonlinear.MaxIter = 100 overland_tiltedV_DWE.Solver.Nonlinear.ResidualTol = 1e-9 overland_tiltedV_DWE.Solver.Nonlinear.EtaChoice = 'EtaConstant' overland_tiltedV_DWE.Solver.Nonlinear.EtaValue = 0.01 overland_tiltedV_DWE.Solver.Nonlinear.UseJacobian = False overland_tiltedV_DWE.Solver.Nonlinear.DerivativeEpsilon = 1e-15 overland_tiltedV_DWE.Solver.Nonlinear.StepTol = 1e-20
wbx.Mannings.Geom.right.Value = 5.e-6 wbx.Mannings.Geom.channel.Value = 1.e-6 #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- wbx.PhaseSources.water.Type = 'Constant' wbx.PhaseSources.water.GeomNames = 'domain' wbx.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- wbx.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- wbx.Solver = 'Richards' wbx.Solver.MaxIter = 100 wbx.Solver.AbsTol = 1E-10 wbx.Solver.Nonlinear.MaxIter = 20 wbx.Solver.Nonlinear.ResidualTol = 1e-9 wbx.Solver.Nonlinear.EtaChoice = 'Walker1' wbx.Solver.Nonlinear.EtaChoice = 'EtaConstant' wbx.Solver.Nonlinear.EtaValue = 0.01 wbx.Solver.Nonlinear.UseJacobian = False
richards_FBx.Geom.domain.ICPressure.RefGeom = 'domain' richards_FBx.Geom.domain.ICPressure.RefPatch = 'bottom' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- richards_FBx.PhaseSources.water.Type = 'Constant' richards_FBx.PhaseSources.water.GeomNames = 'domain' richards_FBx.PhaseSources.water.Geom.domain.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- richards_FBx.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- richards_FBx.Solver = 'Richards' richards_FBx.Solver.MaxIter = 50000 richards_FBx.Solver.Nonlinear.MaxIter = 100 richards_FBx.Solver.Nonlinear.ResidualTol = 1e-6 richards_FBx.Solver.Nonlinear.EtaChoice = 'EtaConstant' richards_FBx.Solver.Nonlinear.EtaValue = 1e-2 richards_FBx.Solver.Nonlinear.UseJacobian = True richards_FBx.Solver.Nonlinear.DerivativeEpsilon = 1e-12
crater.Geom.infiltration.ICPressure.RefPatch = 'infiltration' crater.Geom.infiltration.ICPressure.RefGeom = 'domain' #----------------------------------------------------------------------------- # Phase sources: #----------------------------------------------------------------------------- crater.PhaseSources.water.Type = 'Constant' crater.PhaseSources.water.GeomNames = 'background' crater.PhaseSources.water.Geom.background.Value = 0.0 #----------------------------------------------------------------------------- # Exact solution specification for error calculations #----------------------------------------------------------------------------- crater.KnownSolution = 'NoKnownSolution' #----------------------------------------------------------------------------- # Set solver parameters #----------------------------------------------------------------------------- crater.Solver = 'Richards' crater.Solver.MaxIter = 10000 crater.Solver.Nonlinear.MaxIter = 15 crater.Solver.Nonlinear.ResidualTol = 1e-9 crater.Solver.Nonlinear.StepTol = 1e-9 crater.Solver.Nonlinear.EtaValue = 1e-5 crater.Solver.Nonlinear.UseJacobian = True crater.Solver.Nonlinear.DerivativeEpsilon = 1e-7