# ------------------ INPUTS TO MAIN PROGRAM ------------------- max_step = 5000000 # maximum timestep stop_time = 6000.0 max_step = 10 geometry.is_periodic = 0 0 0 geometry.coord_sys = 2 # 0 => cart, 1 => RZ, 2 => Spherical geometry.prob_lo = 0.0 0.0 0.0 geometry.prob_hi = 1.0e14 1. 1. amr.n_cell = 102400 8 8 # REFINEMENT / REGRIDDING amr.max_level = 0 # maximum level number allowed amr.ref_ratio = 2 2 2 2 # refinement ratio amr.regrid_int = 2 2 2 2 # how often to regrid amr.blocking_factor = 16 # block factor in grid generation amr.max_grid_size = 64 amr.n_error_buf = 4 4 4 4 # number of buffer cells in error est amr.n_proper = 1 # default value amr.grid_eff = 0.7 # what constitutes an efficient grid # CHECKPOINT FILES amr.check_file = chk # root name of checkpoint file amr.check_int = 100 # number of timesteps between checkpoints # PLOTFILES amr.plot_files_output = 1 # suppress plot files amr.plot_file = plt amr.plot_int = 100 # number of timesteps between plot files amr.derive_plot_vars = ALL # PROBLEM PARAMETERS problem.rbasefac = 0.8e0 problem.rwind0 = 2.0e13 problem.rwind1 = 1.0e14 problem.rhowind1 = 1.e-16 problem.Twind1 = 800.0e0 # filter is used only when rho or time is below. problem.filter_rhomax = 1.e-12 problem.filter_timemax = 4.e0 problem.model_file = "model.input.87a.23" radiation.plot_com_flux = 1 # comoving frame flux radiation.plot_lab_flux = 1 # lab frame flux radiation.plot_lambda = 1 # VERBOSITY amr.v = 1 amr.grid_log = grdlog # name of grid logging file # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< # 0 = Interior 3 = Symmetry # 1 = Inflow 4 = SlipWall # 2 = Outflow 5 = NoSlipWall # >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<< castro.lo_bc = 3 castro.hi_bc = 2 # WHICH PHYSICS castro.do_grav = 1 gravity.gravity_type = MonopoleGrav castro.do_hydro = 1 castro.do_radiation = 1 castro.do_reflux = 1 # 1 => do refluxing castro.do_react = 0 # reactions? #castro.difmag = 0.1 # default: 0.1 # hydro cutoff parameters castro.small_dens = 1.e-18 castro.small_temp = 500.0 castro.small_pres = 1.0e-15 # External source terms castro.add_ext_src=0 # Add external source terms # TIME STEP CONTROL castro.cfl = 0.6 # cfl number for hyperbolic system castro.init_shrink = 0.01 # scale back initial timestep castro.change_max = 1.05 # DIAGNOSTICS & VERBOSITY castro.sum_interval = 1 # timesteps between computing mass castro.v = 1 # ------------------ INPUTS TO RADIATION CLASS ------------------- ##### SolverType ##### # 0: single group diffusion w/o coupling to hydro # 5: SGFLD 6: MGFLD radiation.SolverType = 5 radiation.comoving = 1 radiation.filter_lambda_T = 4 radiation.filter_lambda_S = 0 radiation.filter_prim_int = 1 # filter primitives every ? steps radiation.filter_prim_T = 4 radiation.filter_prim_S = 3 castro.first_order_hydro = 1 # default is 0 (i.e., high-order hydro) # RADIATION TOLERANCES radiation.reltol = 1.e-14 # relative tolerance for implicit update loop radiation.abstol = 0.0 # absolute tolerance for implicit update loop radiation.maxiter = 1000 # return after numiter iterations if not converged # 0: conservative; 1: non-conservative radiation.matter_update_type = 0 # 0: no limiter, 2: app. LP, 12: Bruenn, 22: square root, 32: Minerbo radiation.limiter = 2 radiation.update_limiter = -1000 # 0: f = lambda, 1: f = 1/3, 2: f = 1-2*lambda, 3: f = lambda+(lambda*R)^2 # 4: f = 1/3 + 2/3*(F/cE)^2 # If limiter = 12 or 22, option 3 will not be a monotonic function of R. # So we might want to use option 2 or 4 in these cases. radiation.closure = 3 # RADIATION VERBOSITY radiation.v = 2 # verbosity # We set radiation boundary conditions directly since they do not # correspond neatly to the physical boundary conditions used for the fluid. # The choices are: # 101 = LO_DIRICHLET 102 = LO_NEUMANN # 104 = LO_MARSHAK 105 = LO_SANCHEZ_POMRANING radiation.lo_bc = 102 102 102 radiation.hi_bc = 105 102 102 # For each boundary, we can specify either a constant boundary value # or use a Fortran function FORT_RADBNDRY to specify values that vary # in space and time. # If bcflag is 0 then bcval is used, otherwise FORT_RADBNDRY used: radiation.lo_bcflag = 0 0 0 radiation.hi_bcflag = 0 0 0 # bcval is interpreted differently depending on the boundary condition # 101 = LO_DIRICHLET bcval is Dirichlet value of rad energy density # 102 = LO_NEUMANN bcval is inward flux of rad energy # 104 = LO_MARSHAK bcval is incident flux # 105 = LO_SANCHEZ_POMRANING bcval is incident flux radiation.lo_bcval = 0.0 0.0 0.0 radiation.hi_bcval = 0.0 0.0 0.0 radiation.use_opacity_table_module = 1 # ------------------ INPUTS TO RADIATION SOLVER CLASS ------------------- # solver flag values < 100 use HypreABec, support symmetric matrices only # solver flag values >= 100 use HypreMultiABec, support nonsymmetric matrices # # PFMG does not support 1D. # ParCSR does not work for periodic boundaries. # For MGFLD with accelerate = 2, must use >=100. # # 0 SMG # 1 PFMG (>= 2D only) # 100 AMG using ParCSR ObjectType # 102 GMRES using ParCSR ObjectType # 103 GMRES using SStruct ObjectType # 104 GMRES using AMG as preconditioner # 109 GMRES using Struct SMG/PFMG as preconditioner # 150 AMG using ParCSR ObjectType # 1002 PCG using ParCSR ObjectType # 1003 PCG using SStruct ObjectType radsolve.level_solver_flag = 0 radsolve.reltol = 1.0e-11 # relative tolerance radsolve.abstol = 0.0 # absolute tolerance (often not necessary) radsolve.maxiter = 1000 # linear solver iteration limit radsolve.v = 1 # verbosity hmabec.verbose = 2 # verbosity for HypreMultiABec solvers habec.verbose = 2 # verbosity for HypreABec solvers # # The default strategy is SFC. # DistributionMapping.strategy = ROUNDROBIN DistributionMapping.strategy = KNAPSACK DistributionMapping.strategy = SFC