-
Notifications
You must be signed in to change notification settings - Fork 4
/
Copy pathinputs.2d.D15-MESA-SF0.8
91 lines (74 loc) · 2.57 KB
/
inputs.2d.D15-MESA-SF0.8
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
# ------------------ INPUTS TO MAIN PROGRAM -------------------
max_step = 0
stop_time = 0.01
# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 0 0
geometry.coord_sys = 1 # R-Z coordinates
geometry.prob_lo = 0.00E+00 -4.00E+09 # Restrict to quadrant
geometry.prob_hi = 4.00E+09 4.00E+09
amr.n_cell = 64 128
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior 3 = Symmetry
# 1 = Inflow 4 = SlipWall
# 2 = Outflow 5 = NoSlipWall
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
castro.lo_bc = 3 2
castro.hi_bc = 2 2
# WHICH PHYSICS
castro.do_hydro = 1
castro.do_grav = 1
castro.do_reflux = 1
castro.do_react = 0
castro.do_sponge = 0
# GRAVITY
gravity.gravity_type = MonopoleGrav
castro.drdxfac = 2
castro.grav_source_type = 4
castro.point_mass = 0.0 #-2.878606188e+32
castro.point_mass_fix_solution = 1
# REACTIONS
castro.react_T_max = 6.0e9
castro.react_T_min = 1.0e8
castro.dtnuc_X = 1.0e-4
castro.dtnuc_e = 1.0e-2
# HYDRO
castro.ppm_type = 1
castro.riemann_solver = 0
castro.allow_negative_energy = 0
castro.density_reset_method = 1
castro.transverse_reset_density = 1
castro.transverse_reset_rhoe = 0
castro.transverse_use_eos = 0
# TIME STEP CONTROL
castro.dt_cutoff = 5.e-20 # level 0 timestep below which we halt
castro.cfl = 0.5 # cfl number for hyperbolic system
castro.init_shrink = 0.1 # scale back initial timestep
castro.change_max = 1.1 # maximum increase in dt over successive steps
castro.use_retry = 1
# DIAGNOSTICS & VERBOSITY
castro.sum_interval = 1 # timesteps between computing mass
castro.v = 1 # verbosity in Castro.cpp
gravity.v = 1 # verbosity in Gravity.cpp
amr.v = 1 # verbosity in Amr.cpp
castro.print_energy_diagnostics = 1
castro.show_center_of_mass = 1
# REFINEMENT / REGRIDDING
amr.max_level = 4 # maximum level number allowed
amr.ref_ratio = 2 2 2 2 # refinement ratio
amr.regrid_int = 2 # how often to regrid
amr.blocking_factor = 4 # block factor in grid generation
amr.max_grid_size = 128
amr.regrid_on_restart = 1
amr.n_error_buf = 2 2 2 2
amr.grid_eff = 0.95
# CHECKPOINT FILES
amr.check_file = chk_chimera_2d_ # root name of checkpoint file
amr.check_int = 1000 # number of timesteps between checkpoints
#amr.check_per = 0.005
# PLOTFILES
amr.plot_file = plt_chimera_2d_
amr.plot_int = 1
#amr.plot_per = 0.01
amr.derive_plot_vars=ALL
# PROBIN FILENAME
amr.probin_file = probin.2d.D15-MESA-SF0.8