# ------------------ INPUTS TO MAIN PROGRAM ------------------- max_step = 25 # maximum timestep stop_time = 0.04 #max_step = 11 geometry.is_periodic = 0 0 0 geometry.coord_sys = 0 # 0 => cart, 1 => RZ, 2 => Spherical geometry.prob_lo = -4000. 0.0 0.0 geometry.prob_hi = 2000. 187.5 187.5 amr.n_cell = 256 8 8 # REFINEMENT / REGRIDDING amr.max_level = 4 # maximum level number allowed amr.ref_ratio = 2 2 2 2 2 2 # refinement ratio amr.regrid_int = 2 2 2 2 2 2 # how often to regrid amr.blocking_factor = 8 # block factor in grid generation amr.max_grid_size = 64 amr.n_error_buf = 2 2 2 2 2 2 # number of buffer cells in error est amr.n_proper = 1 # default value amr.grid_eff = 0.7 # what constitutes an efficient grid amr.refinement_indicators = tempgrad_rel amr.refine.tempgrad_rel.max_level = 14 amr.refine.tempgrad_rel.relative_gradient = 0.1 amr.refine.tempgrad_rel.field_name = Temp # CHECKPOINT FILES amr.check_file = chk # root name of checkpoint file amr.check_int = 1000 # number of timesteps between checkpoints #amr.restart = chk00011 # PLOTFILES amr.plot_file = plt_ amr.plot_int = 500 # number of timesteps between plot files amr.derive_plot_vars = ALL # PROBLEM PARAMETERS problem.rho0 = 5.45969027775e-13 problem.T0 = 100.0e0 problem.v0 = 588588.075591e0 problem.rho1 = 1.96434778127e-12 problem.T1 = 855.719921848 problem.v1 = 163591.632018e0 problem.idir = 1 # EOS eos.eos_gamma = 1.6666666667 # OPACITY opacity.const_kappa_p = 3.92663697758e-5 opacity.kappa_p_exp_m = 0.0e0 opacity.kappa_p_exp_n = 0.0e0 opacity.const_kappa_r = 0.848902853095e0 opacity.kappa_r_exp_m = 0.0e0 opacity.kappa_r_exp_n = 0.0e0 # 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 = 2 4 4 castro.hi_bc = 2 4 4 # WHICH PHYSICS castro.do_grav = 0 castro.do_hydro = 1 castro.do_radiation = 1 castro.do_reflux = 1 # 1 => do refluxing castro.do_react = 0 # reactions? # hydro cutoff parameters castro.small_dens = 1.e-20 # External source terms castro.add_ext_src=0 # Add external source terms # TIME STEP CONTROL castro.cfl = 0.8 # cfl number for hyperbolic system castro.init_shrink = 0.1 # scale back initial timestep castro.change_max = 1.1 #castro.initial_dt = 0.01 #castro.fixed_dt = 1.e-15 # 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 = 6 # RADIATION LIMITER radiation.limiter = 0 # 0 = no limiter # 2 = correct form of Lev-Pom limiter # 0: f = lambda, 1: f = 1/3, 2: f = 1-2*lambda, 3: f = lambda+(lambda*R)^2 radiation.closure = 1 radiation.fspace_advection_type = 2 radiation.nGroups = 16 radiation.lowestGroupHz = 1.e10 radiation.highestGroupHz = 1.e15 radiation.accelerate = 1 radiation.skipAccelAllowed = 1 radiation.relInTol = 1.e-6 # relative tolerance for inner update loop radiation.absInTol = 0.0 # absolute tolerance for inner update loop radiation.maxInIter = 20 # max iterations for inner update loop radiation.minInIter = 3 # RADIATION TOLERANCES radiation.reltol = 1.e-6 # relative tolerance for implicit update loop radiation.abstol = 0.0 # absolute tolerance for implicit update loop radiation.maxiter = 50 # return after numiter iterations if not converged # 0: both, 1: rhoe, 2: residue of rhoe equation, 3: T radiation.convergence_check_type = 3 # C: consevartive NC: non-conservative # 0: C C ... C, 1: NC NC ... NC C C ... C, 2: NC C NC C ... NC C radiation.matter_update_type = 0 # 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 = 101 102 102 radiation.hi_bc = 101 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 # radiation.hi_bcval = 0 0 0 radiation.lo_bcval0 = 3.27517624962438426E-014 2.82440556862939748E-013 \ 2.42497757527818012E-012 2.06318928406580997E-011 \ 1.72250911764424658E-010 1.38202167847469962E-009 \ 1.01793273622194353E-008 6.19571982820059287E-008 \ 2.42768244637479146E-007 3.56904140550843428E-007 \ 8.25303873091818086E-008 6.56183715461762625E-010 \ 4.87717642913733243E-015 0.0000000000000000 0.0000000000000000 \ 0.0000000000000000 radiation.hi_bcval0 = 2.81241745447867642E-013 2.43427083781028680E-012 \ 2.10589023672203340E-011 1.81989273590520617E-010 \ 1.56933800834926731E-009 1.34726324109544666E-008 \ 1.14601755897558854E-007 9.56358213326894788E-007 \ 7.66592208464726935E-006 5.63476073748010151E-005 \ 3.41365489499386712E-004 1.32317116979096768E-003 \ 1.90152165072088103E-003 4.22478249639437597E-004 \ 3.10661088952744746E-006 1.95926322041350964E-011 # ------------------ 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 = 104 # can be any supported hypre solver flag radsolve.reltol = 1.0e-11 # relative tolerance radsolve.abstol = 0.0 # absolute tolerance (often not necessary) radsolve.maxiter = 200 # linear solver iteration limit radsolve.v = 1 # verbosity hmabec.verbose = 1 # verbosity for HypreMultiABec solvers habec.verbose = 1 # verbosity for HypreABec solvers # # The default strategy is SFC. # DistributionMapping.strategy = ROUNDROBIN DistributionMapping.strategy = KNAPSACK DistributionMapping.strategy = SFC