-
Notifications
You must be signed in to change notification settings - Fork 99
/
Copy pathinputs_nova_en
110 lines (86 loc) · 3.07 KB
/
inputs_nova_en
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
# ------------------ INPUTS TO MAIN PROGRAM -------------------
max_step = 100000
# PROBLEM SIZE & GEOMETRY
geometry.is_periodic = 1 0
geometry.coord_sys = 0 # 0 => cart, 1 => RZ 2=>spherical
geometry.prob_lo = 0 0
geometry.prob_hi = 8.192e7 4.096e7
amr.n_cell = 256 128
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior 3 = Symmetry
# 1 = Inflow 4 = SlipWall
# 2 = Outflow 5 = NoSlipWall
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
castro.lo_bc = 0 1
castro.hi_bc = 0 2
castro.yl_ext_bc_type = 1
castro.hse_interp_temp = 1
castro.hse_reflect_vels = 1
castro.fill_ambient_bc = 1
castro.ambient_fill_dir = 1
castro.ambient_outflow_vel = 1
# WHICH PHYSICS
castro.do_hydro = 1 #hydrodynamics
castro.do_react = 1 #microphysics
castro.add_ext_src = 0 #external sources
castro.do_grav = 1 #gravity
castro.do_sponge = 0 #using sponge
castro.ppm_type = 1
castro.riemann_solver = 0
castro.use_flattening = 1
castro.small_temp = 1.0e5
castro.small_dens = 1.0e-6
castro.retry_small_density_cutoff = 1.0
gravity.gravity_type = ConstantGrav
gravity.const_grav = -8.84e8
castro.react_rho_min = 1.0e2
castro.react_T_min = 1.0e6
# TIME STEP CONTROL
castro.cfl = 0.5 # cfl number for hyperbolic system
#castro.init_shrink = 0.1 # scale back initial timestep
castro.change_max = 1.1 # max time step growth
# DIAGNOSTICS & VERBOSITY
castro.sum_interval = 1 # timesteps between computing mass
castro.v = 1 # verbosity in Castro.cpp
amr.v = 1 # verbosity in Amr.cpp
# REFINEMENT / REGRIDDING
amr.max_level = 1 # 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 = 32 # block factor in grid generation
amr.max_grid_size = 128
amr.n_error_buf = 2 2 2 2 # number of buffer cells in error est
# CHECKPOINT FILES
amr.check_file = chk # root name of checkpoint file
amr.check_int = 200 # number of timesteps between checkpoints
# PLOTFILES
amr.plot_file = plt # root name of plotfile
amr.plot_int = 500 # number of timesteps between plotfiles
amr.derive_plot_vars = ALL
# DATA LOG
amr.data_log = log
# problem initialization
problem.model_name = "glasner_enriched1600.hse"
problem.apply_vel_field = 0
#problem.velpert_scale = 5.0e5
#problem.velpert_scale_ratio = 7.0
#problem.velpert_amplitude = 7.0e5
#problem.num_vortices = 8
problem.apply_temp_field = 1
problem.width = 1.0e6
problem.amplitude = 0.5
problem.num_vortices = 10
problem.width_ratio = 1.8
# Microphysics
integrator.jacobian = 2
integrator.rtol_spec = 1.e-6
integrator.atol_spec = 1.e-6
integrator.rtol_enuc = 1.e-6
# refinement
amr.refinement_indicators = dens temp
amr.refine.dens.max_level = 1
amr.refine.dens.value_greater = 1.0e2
amr.refine.dens.field_name = density
amr.refine.temp.max_level = 1
amr.refine.temp.value_greater = 1.0e7
amr.refine.temp.field_name = Temp