-
Notifications
You must be signed in to change notification settings - Fork 3
/
orszag_tang.athinput
71 lines (58 loc) · 2.7 KB
/
orszag_tang.athinput
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
# AthenaXXX input file for Orszag-Tang MHD vortex test
<comment>
problem = Orszag-Tang vortex
reference = Orszag,S. & Tang,W., J. Fluid Mech., 90, 129 (1998)
<job>
basename = OrszagTang # problem ID: basename of output filenames
<mesh>
nghost = 3 # Number of ghost cells
nx1 = 200 # Number of zones in X1-direction #> SCALE 50:200:50
x1min = -0.5 # minimum value of X1
x1max = 0.5 # maximum value of X1
ix1_bc = periodic # inner-X1 boundary flag
ox1_bc = periodic # outer-X1 boundary flag
nx2 = 200 # Number of zones in X2-direction #> SCALE 50:200:50
x2min = -0.5 # minimum value of X2
x2max = 0.5 # maximum value of X2
ix2_bc = periodic # inner-X2 boundary flag
ox2_bc = periodic # outer-X2 boundary flag
nx3 = 1 # Number of zones in X3-direction
x3min = -0.5 # minimum value of X3
x3max = 0.5 # maximum value of X3
ix3_bc = periodic # inner-X3 boundary flag
ox3_bc = periodic # outer-X3 boundary flag
<meshblock>
nx1 = 50 # Number of cells in each MeshBlock, X1-dir
nx2 = 50 # Number of cells in each MeshBlock, X2-dir
nx3 = 1 # Number of cells in each MeshBlock, X3-dir
<time>
evolution = dynamic # dynamic/kinematic/static
integrator = rk3 # time integration algorithm
cfl_number = 0.4 # The Courant, Friedrichs, & Lewy (CFL) Number #> SCALE 0:1.2:0.1
nlim = -1 # cycle limit
tlim = 1.0 # time limit #> SCALE 1:100:1
ndiag = 1 # cycles between diagostic output
<mhd>
eos = ideal # EOS type
reconstruct = ppm4 # spatial reconstruction method
rsolver = hlld # Riemann-solver to be used
gamma = 1.666666667 # gamma = C_p/C_v #> SCALE 0.1:2.0:0.05
# Note: OT pgen does not read any parameters from input file
<problem>
pgen_name = orszag_tang # problem generator name
<output1>
file_type = hst # History data dump
data_format = %12.5e # optional data format string
dt = 0.01 # time increment between outputs
<output2>
file_type = tab # Tabular data dump
variable = mhd_w_bcc # variables to be output
data_format = %12.5e # Optional data format string
dt = 0.01 # time increment between outputs
slice_x2 = 0.0 # slice in x2 #> SCALE -0.5:0.5:0.1
slice_x3 = 0.0 # slice in x3
ghost_zones = false # switch to output ghost cells
<output3>
file_type = bin # binary data dump
variable = mhd_w_bcc # variables to be output
dt = 0.01 # time increment between outputs