athinput.linear_wave1d

<comment>
problem   = 1D Hydro linear wave convergence
reference =
configure = --prob=linear_wave

<job>
problem_id = LinWave   # problem ID: basename of output filenames   

<output1>
file_type  = hst       # History data dump
dt         = 0.01      # time increment between outputs

<output2>
file_type  = vtk       # VTK data dump
variable   = prim      # variables to be output
dt         = 0.01      # time increment between outputs

<time>
cfl_number  = 0.4       # The Courant, Friedrichs, & Lewy (CFL) Number  #> SCALE 0:1.2:0.1
nlim        = -1        # cycle limit
tlim        = 5.0       # time limit                                    #> SCALE 1:100:1
integrator  = vl2       # time integration algorithm
xorder      = 2         # order of spatial reconstruction
ncycle_out  = 1         # interval for stdout summary info
correct_ic  = false     # correct midpoint assumption in initial condition
correct_err = false     # correct midpoint assumption in analytic solution

<mesh>
nx1        = 128       # Number of zones in X1-direction    #> SCALE 16:512:16
x1min      = 0.0       # minimum value of X1                #> ENTRY
x1max      = 1.0       # maximum value of X1                #> ENTRY
ix1_bc     = periodic  # inner-X1 boundary flag
ox1_bc     = periodic  # outer-X1 boundary flag

nx2        = 1         # Number of zones in X2-direction
x2min      = -0.5      # minimum value of X2
x2max      = 0.5       # maximum value of X2

nx3        = 1         # Number of zones in X3-direction
x3min      = -0.5      # minimum value of X3
x3max      = 0.5       # maximum value of X3

num_threads = 1        # maximum number of OMP threads

<hydro>
gamma = 1.666666666666667 # gamma = C_p/C_v
iso_sound_speed = 1.0     # isothermal sound speed

<problem>
compute_error = true  # when 'true' outputs L1 error compared to initial data
wave_flag = 0         # Wave family number (0 - 4 for adiabatic hydro)
amp       = 1.0e-6    # Wave Amplitude
vflow     = 0.0       # background flow velocity
ang_2     = 0.0       # rotation of k around x3 axis
ang_3     = 0.0       # rotation of k around x2 axis