LWvdz.Solver.Nonlinear.Globalization = 'LineSearch' LWvdz.Solver.Linear.KrylovDimension = 80 LWvdz.Solver.Linear.MaxRestarts = 2 LWvdz.Solver.Linear.Preconditioner = 'MGSemi' LWvdz.Solver.Linear.Preconditioner = 'PFMG' LWvdz.Solver.Linear.Preconditioner.PCMatrixType = 'FullJacobian' ##--------------------------------------------------------- # Initial conditions: water pressure #--------------------------------------------------------- # set water table to be at the bottom of the domain, the top layer is initially dry LWvdz.ICPressure.Type = 'HydroStaticPatch' LWvdz.ICPressure.GeomNames = 'domain' LWvdz.Geom.domain.ICPressure.Value = 0.0 LWvdz.Geom.domain.ICPressure.RefGeom = 'domain' LWvdz.Geom.domain.ICPressure.RefPatch = 'z_upper' #spinup key LWvdz.OverlandFlowSpinUp = 1 LWvdz.OverlandSpinupDampP1 = 1.0 LWvdz.OverlandSpinupDampP2 = 0.00001 #----------------------------------------------------------------------------- # Run and do tests #----------------------------------------------------------------------------- LWvdz.run()
LW_var_dz_spinup.Solver.Linear.Preconditioner = 'MGSemi' LW_var_dz_spinup.Solver.Linear.Preconditioner = 'PFMG' LW_var_dz_spinup.Solver.Linear.Preconditioner.PCMatrixType = 'FullJacobian' LW_var_dz_spinup.Solver.WriteSiloPressure = True LW_var_dz_spinup.Solver.WriteSiloSaturation = True ##--------------------------------------------------------- # Initial conditions: water pressure #--------------------------------------------------------- # set water table to be at the bottom of the domain, the top layer is initially dry LW_var_dz_spinup.ICPressure.Type = 'HydroStaticPatch' LW_var_dz_spinup.ICPressure.GeomNames = 'domain' LW_var_dz_spinup.Geom.domain.ICPressure.Value = 0.0 LW_var_dz_spinup.Geom.domain.ICPressure.RefGeom = 'domain' LW_var_dz_spinup.Geom.domain.ICPressure.RefPatch = 'z_upper' #spinup key LW_var_dz_spinup.OverlandFlowSpinUp = 1 LW_var_dz_spinup.OverlandSpinupDampP1 = 1.0 LW_var_dz_spinup.OverlandSpinupDampP2 = 0.00001 #----------------------------------------------------------------------------- # Run and do tests #----------------------------------------------------------------------------- LW_var_dz_spinup.run()