Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # Parameters file for core-collapse SN
- # Simulation basics
- basenm = "a140_"
- restart = .false.
- checkpointFileNumber = 0
- plotFileNumber = 0
- run_comment = "a140 1D with radiation"
- #output_directory = "output"
- # IO
- checkpointFileIntervalStep = 0
- checkpointFileIntervalTime = 0.05
- plotFileIntervalStep = 0
- plotFileIntervalTime = 0.
- # wall_clock_checkpoint = 43200.
- wall_clock_time_limit = 21000.
- # Time
- tinitial = 0.0
- tmax = 3.9
- nend = 100000000
- tstep_change_factor = 1.5
- dtinit = 1.E-8
- dtmax = 1.E5
- dtmin = 1.E-20
- # Domain
- geometry = "spherical"
- xmax = 1.9e11
- xmin = 0.0
- xl_boundary_type = "reflect"
- xr_boundary_type = "outflow"
- # Grid/Refinement
- nblockx = 100
- nblocky = 1
- nblockz = 1
- gr_lrefineMaxRedDoByLogR = .true.
- gr_lrefineMaxRedDoByTime = .FALSE.
- gr_lrefineMaxRedLogBase = 10.0
- gr_lrefineMaxRedRadiusFact = 0.15
- gr_lrefineMaxRedTRef = 2.0
- gr_lrefineMaxRedTimeScale = 0.5
- lrefine_max = 3
- lrefine_min = 1
- refine_var_1 = "dens"
- refine_var_2 = "pres"
- refine_var_3 = "none"
- refine_var_4 = "none"
- refine_cutoff_1 = 0.8
- refine_cutoff_2 = 0.8
- refine_cutoff_3 = 0.8
- refine_cutoff_4 = 0.8
- # Simulation
- model_file = "a140.flash"
- nsub = 4
- vel_mult = 1.0
- # Hydro
- useHydro = .TRUE.
- cfl = 0.5
- interpol_order = 1
- updateHydroFluxes = .TRUE.
- eintSwitch = 0.0 # Always use Etot
- convertToConsvdForMeshCalls = .false.
- converttoconsvdinmeshinterp = .false.
- ## SWITCHES SPECIFIC TO THE UNSPLIT HYDRO SOLVER ##
- # I. INTERPOLATION SCHEME:
- order = 3 # Interpolation order (first/second/third/fifth order)
- slopeLimiter = "hybrid" # Slope limiters (minmod, mc, vanLeer, hybrid, limited)
- LimitedSlopeBeta= 1. # Slope parameter for the "limited" slope by Toro
- charLimiting = .true. # Characteristic limiting vs.Primitive limiting
- use_avisc = .true. # use artificial viscosity (originally for PPM)
- cvisc = 0.1 # coefficient for artificial viscosity
- use_flattening = .true. # use flattening (dissipative) (originally for PPM)
- use_steepening = .false. # use contact steepening (originally for PPM)
- use_upwindTVD = .false. # use upwind biased TVD slope for PPM (need nguard=6)
- flux_correct = .true.
- EOSforRiemann = .false.
- transOrder = 1
- use_auxEintEqn = .TRUE.
- hydroComputeDtOption = -1
- # II. RIEMANN SOLVERS:
- RiemannSolver = "hybrid" # Roe, HLL, HLLC, LLF, Marquina
- entropy = .false. # Entropy fix for the Roe solver
- # III. STRONG SHOCK HANDELING SCHEME:
- shockDetect = .true. # Shock Detect for numerical stability
- ## -------------------------------------------------------------##
- # Gravity
- useGravity = .true.
- updateGravity = .TRUE.
- grav_boundary_type = "isolated"
- mpole_3daxisymmetry = .false.
- mpole_dumpMoments = .FALSE.
- mpole_PrintRadialInfo = .false.
- mpole_IgnoreInnerZone = .false.
- mpole_lmax = 16
- mpole_ZoneRadiusFraction_1 = 1.0
- mpole_ZoneExponent_1 = 0.005
- mpole_ZoneScalar_1 = 0.5
- mpole_ZoneType_1 = "logarithmic"
- point_mass = 0.0
- point_mass_rsoft = 0.e0
- use_gravHalfUpdate = .TRUE.
- use_gravConsv = .FALSE.
- use_gravPotUpdate = .FALSE.
- # EOS
- eosMode = "dens_ie"
- eosModeInit = "dens_temp_equi"
- gamma = 1.66666666666667
- # Small numbers
- smallt = 1.0e-20
- smlrho = 1.0e-27
- smltemp = 1.0e-2
- smallp = 1.E-20
- smalle = 1.E-20
- smallu = 1.E-10
- smallx = 1.E-100
- small = 1.E-100
- ### Multigroup diffusion options ###
- rt_useMGD = .true.
- rt_mgdNumGroups = 1
- rt_mgdBounds_1 = 0.0
- rt_mgdBounds_2 = 1.0e14
- op_emitConst = 10.0
- op_transConst = 1.0e-06
- op_absorbConst = 10.0
- # Set boundary conditions:
- rt_mgdXlBoundaryType = "vacuum"
- rt_mgdXrBoundaryType = "vacuum"
- rt_mgdYlBoundaryType = "vacuum"
- rt_mgdYrBoundaryType = "vacuum"
- rt_mgdZlBoundaryType = "vacuum"
- rt_mgdZrBoundaryType = "vacuum"
Advertisement
Add Comment
Please, Sign In to add comment