Esempi in Python per Run.KnownSolution

Esempio n. 1

File: clm_samrai.py Progetto: mahdadt/parflow-360

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

Esempio n. 2

File: octree-simple.py Progetto: mahdadt/parflow-360

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

Esempio n. 3

File: octree-simple.py Progetto: grapp1/kw-intern

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

Esempio n. 4

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

Esempio n. 5

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

Esempio n. 6

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

Esempio n. 7

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

Esempio n. 8

File: water_balance_y.py Progetto: mahdadt/parflow-360

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

Esempio n. 9

File: default_single.py Progetto: grapp1/kw-intern

#-----------------------------------------------------------------------------
# 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()

Esempio n. 10

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

Esempio n. 11

File: var_dz_1D.py Progetto: mahdadt/parflow-360

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

Esempio n. 12

File: small_domain.py Progetto: grapp1/kw-intern

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

Esempio n. 13

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

Esempio n. 14

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

Esempio n. 15

File: default_overland.pfmg_octree.jac.py Progetto: grapp1/kw-intern

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

Esempio n. 16

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

Esempio n. 17

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

Esempio n. 18

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'

Esempio n. 19

File: clm-reuse.py Progetto: mahdadt/parflow-360

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

Esempio n. 20

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

Esempio n. 21

File: default_richards_wells.py Progetto: mahdadt/parflow-360

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

Esempio n. 22

File: richards_box_proctest_vardz.py Progetto: grapp1/kw-intern

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

Esempio n. 23

File: PFCLM_SC.py Progetto: reedmaxwell/ParFlow_example_cases

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

Esempio n. 24

File: pfmg_octree.py Progetto: mahdadt/parflow-360

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

Esempio n. 25

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

Esempio n. 26

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

Esempio n. 27

File: overland_tiltedV_DWE.py Progetto: grapp1/kw-intern

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

Esempio n. 28

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

Esempio n. 29

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

Esempio n. 30

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