diff --git a/orszag_tang.athinput b/orszag_tang.athinput
new file mode 100644
index 0000000..275a1e8
--- /dev/null
+++ b/orszag_tang.athinput
@@ -0,0 +1,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
+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
+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       = 200       # Number of cells in each MeshBlock, X1-dir
+nx2       = 200       # 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
+nlim       = -1        # cycle limit
+tlim       = 1.0       # time limit
+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
+
+# 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
+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