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
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
|
# NOTE: all the values here must be set unless specifically
# marked as optional
#**********************************
#ALGORITHM CONTROLS (defaults shown; these are both optional)
#**********************************
proj.proj_2 = 1
proj.filter_factor = 0.
ns.init_iter = 2
ns.do_MLsync_proj = 1
#**********************************
#RUN CONTROLS
#**********************************
# MUST SET ONE OR BOTH OF THESE.
max_step = 1 # maximum (integer) number of time steps
stop_time = 2.0 # maximum (real) time (in same units as dt)
amr.n_cell = 16 16 # indices of the domain at level 0
# (this one is optional)
#amr.restart = # set if starting from a restart file, else comment out
# This will default to file "probin" if not set
amr.probin_file = probin_2
#**********************************
#ADAPTIVITY CONTROLS
#**********************************
amr.max_level = 2 # maximum number of levels of refinement
amr.regrid_int = 2 # how often to regrid if max_level > 0
# (this one is optional)
amr.regrid_file = fixed_grids_2 # file name for specifying fixed grids
# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
amr.ref_ratio = 4 2 # integer refinement ratio
#amr.ref_ratio_vect # IntVect refinement ratio
#
# Only call LinOp::ApplyBC() on alternate calls to LinOp::Fsmooth()
# from within LinOp::smooth(). This is not theoretically correct but
# can cut down fairly significantly on communication overhead when
# running in Parallel; especially so for periodic problems.
#
# 0 is false and 1 is true
#
Lp.alternateApplyBC = 0
#**********************************
#VERBOSITY CONTROLS (these are all optional with defaults shown)
#**********************************
ns.v = 1
godunov.v = 0
mac.v = 1
proj.v = 1
diffuse.v = 1
amr.v = 1
mg.v = 1
cg.v = 0
proj.Pcode = 1
#**********************************
#OUTPUT CONTROLS
#**********************************
#THESE FILES WILL BE CREATED ONLY IF YOU SPECIFY THE NAME
#amr.run_log = runlog # name of the run log file
#amr.grid_log = grdlog # name of the grid log file
#amr.data_log = datalog # name of the data log file
#THESE FILES WILL BE CREATED WITH DEFAULT PREFIX IF YOU DONT SPECIFY OTHERWISE
#amr.check_file # prefix of the checkpoint files (default = "chk")
#amr.plot_file # prefix of the plot files (default = "plt")
# (this one is optional - sums will only be done if this is set to positive)
ns.sum_interval = 5 # interval at which to sum the specified quantities
# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
#amr.check_per # (real) period between checkpoint files (in same units as dt)
amr.check_int = 50 # (integer) number of timesteps between checkpoint files
# MUST SET ONLY ONE OF THE TWO BELOW (COMMENT OUT THE OTHER)
#amr.plot_per # (real) period between plot files (in same units as dt)
amr.plot_int = 1 # (integer) number of timesteps between plot files
#**********************************
#TIME STEP CONTROLS
#**********************************
ns.cfl = 0.9 # CFL number used to set dt
# (this is optional)
ns.init_shrink = 0.1 # factor which multiplies the very first time step
#**********************************
#PHYSICS
#**********************************
ns.vel_visc_coef = 0.01
ns.temp_cond_coef = 0.01
ns.scal_diff_coefs = 0.0
ns.do_temp = 1
# (this is optional)
ns.gravity = 9.8 # forcing term is rho * abs(gravity) "down"
#**********************************
#GEOMETRY
#**********************************
geometry.coord_sys = 1 # 0 for x-y (x-y-z), 1 for r-z
geometry.prob_lo = 0. 0. # physical dimensions of the low end of the domain
geometry.prob_hi = 1. 1. # physical dimensions of the high end of the domain
geometry.is_periodic = 0 0
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
# 0 = Interior/Periodic 3 = Symmetry
# 1 = Inflow 4 = SlipWall
# 2 = Outflow 5 = NoSlipWall
ns.lo_bc = 5 1 # boundary conditions on the low end of the domain
ns.hi_bc = 5 2 # boundary conditions on the high end of the domain
#**********************************
#MORE OUTPUT CONTROLS
#**********************************
#
# Select form of FAB output: default is NATIVE
#
# ASCII (this is very slow)
# NATIVE (native binary form on machine -- the default)
# IEEE32 (useful if you want 32bit files when running in double precision)
# 8BIT (eight-bit run-length-encoded)
#
fab.format = NATIVE
#**********************************
#MORE PARALLEL CONTROLS
#**********************************
#
# Initializes DistributionMapping strategy from ParmParse.
#
# ParmParse options are:
#
# DistributionMapping.strategy = ROUNDROBIN
# DistributionMapping.strategy = KNAPSACK
#
# The default strategy is ROUNDROBIN.
#
DistributionMapping.strategy = ROUNDROBIN
DistributionMapping.strategy = KNAPSACK
#**********************************
#GRIDDING CONTROLS (these are optional)
#**********************************
amr.blocking_factor = 1 # the factor by which every patch must be coarsenable
amr.n_error_buf = 2 2 # number of buffer cells around each tagged cell
# necessary for multigrid convergence
mg.usecg = 1
# USE THESE ONLY FOR TESTING
#amr.grid_eff = 0.9 # grid efficiency in the grid creation algorithm
# amr.max_grid_size = 16 # maximum linear dimension of any one grid
#
# StationData.vars -- Names of StateData components to output
# StationData.coord -- BL_SPACEDIM array of Reals
# StationData.coord -- the next one
# StationData.coord -- ditto ...
#
# e.g.
#
#StationData.vars = density xmom ymom zmom
#StationData.coord = 0.1 0.1
#StationData.coord = 0.11 0.21
#StationData.coord = 0.5 0.5
#StationData.coord = 1.23 1.23
#StationData.coord = 0.23 1.53
#StationData.coord = 2.34 2.34
|
