diff --git a/problems/blast.in b/problems/blast.in
new file mode 100644
index 0000000..6047626
--- /dev/null
+++ b/problems/blast.in
@@ -0,0 +1,69 @@
+# problem name and parameters
+#
+problem              = "blast"
+dens                 = 1.00d+00
+ratio                = 1.00d+02
+radius               = 1.00d-01
+csnd                 = 4.08d-01
+gamma                = 1.67d-01
+buni                 = 0.00d+00
+angle                = 0.00d+00
+
+# physics
+#
+equation_system      = "hd"
+equation_of_state    = "adi"
+
+# methods
+#
+time_advance         = "rk2"
+riemann_solver       = "hllc"
+reconstruction       = "tvd"
+limiter              = "mc"
+fix_positivity       = "off"
+clip_extrema         = "off"
+
+# mesh parameters
+#
+rdimx                =  2
+rdimy                =  3
+rdimz                =  1
+xmin                 = -0.5
+xmax                 =  0.5
+ymin                 = -0.75
+ymax                 =  0.75
+zmin                 = -0.5
+zmax                 =  0.5
+
+# refinement control
+#
+ncells               = 8
+nghost               = 2
+minlev               = 1
+maxlev               = 7
+crefmin              = 2.00d-01
+crefmax              = 8.00d-01
+epsref               = 1.00d-02
+refinement_variables = "dens pres"
+
+# boundary conditions
+#
+xlbndry              = "periodic"
+xubndry              = "periodic"
+ylbndry              = "periodic"
+yubndry              = "periodic"
+zlbndry              = "periodic"
+zubndry              = "periodic"
+
+# runtime control parameters
+#
+tmax                 = 1.00d+00
+cfl                  = 4.00d-01
+
+# data output control
+#
+precise_snapshots    = "off"
+snapshot_type        = "p"
+snapshot_interval    = 1.0d-01
+restart_number       = -1
+integrals_interval   = 10
diff --git a/problems/current_sheet.in b/problems/current_sheet.in
new file mode 100644
index 0000000..f5b6093
--- /dev/null
+++ b/problems/current_sheet.in
@@ -0,0 +1,60 @@
+# problem name and parameters
+#
+problem              = "current_sheet"
+resistivity          = 1.00d-05
+alpha                = 1.00d+01
+
+# physics
+#
+equation_system      = "mhd"
+equation_of_state    = "adi"
+
+# methods
+#
+time_advance         = "rk2"
+riemann_solver       = "hlld"
+reconstruction       = "limo3"
+limiter              = "mc"
+fix_positivity       = "off"
+
+# mesh parameters
+#
+rdimx                =  1
+rdimy                =  1
+rdimz                =  1
+xmin                 = -5.0d-01
+xmax                 =  5.0d-01
+ymin                 = -5.0d-01
+ymax                 =  5.0d-01
+zmin                 = -5.0d-01
+zmax                 =  5.0d-01
+
+# refinement control
+#
+ncells               = 8
+nghost               = 2
+minlev               = 1
+maxlev               = 7
+refinement_variables = "dens pres magy"
+
+# boundary conditions
+#
+xlbndry              = "periodic"
+xubndry              = "periodic"
+ylbndry              = "periodic"
+yubndry              = "periodic"
+zlbndry              = "periodic"
+zubndry              = "periodic"
+
+# runtime control parameters
+#
+tmax                 = 1.0d+00
+cfl                  = 4.0d-01
+
+# data output control
+#
+precise_snapshots    = "off"
+snapshot_type        = "p"
+snapshot_interval    = 1.0d-01
+restart_number       = -1
+integrals_interval   = 10
diff --git a/problems/implosion.in b/problems/implosion.in
new file mode 100644
index 0000000..c8fec38
--- /dev/null
+++ b/problems/implosion.in
@@ -0,0 +1,53 @@
+# problem name and parameters
+#
+problem               = "implosion"
+gamma                 =  1.40d+00
+shock_line            =  5.00d-01
+
+# physics
+#
+equation_system       = "hd"
+equation_of_state     = "adi"
+
+# methods
+#
+time_advance          = "rk2"
+riemann_solver        = "hllc"
+reconstruction        = "limo3"
+limiter               = "mc"
+fix_positivity        = "off"
+
+# mesh parameters
+#
+ncells                =  8
+nghosts               =  2
+rdimx                 =  1
+rdimy                 =  1
+xmin                  = -1.00d+00
+xmax                  =  1.00d+00
+ymin                  = -1.00d+00
+ymax                  =  1.00d+00
+
+# refinement control
+#
+maxlev                = 8
+
+# boundary conditions
+#
+#xlbndry               = "reflecting"
+#xubndry               = "reflecting"
+#ylbndry               = "reflecting"
+#yubndry               = "reflecting"
+
+# runtime control parameters
+#
+tmax                  = 1.00d+00
+cfl                   = 4.00d-01
+
+# data output control
+#
+precise_snapshots     = "on"
+snapshot_type         = "p"
+snapshot_interval     = 1.0d-02
+restart_number        = -1
+integrals_interval    = 1
diff --git a/problems/kh.in b/problems/kh.in
new file mode 100644
index 0000000..e31efee
--- /dev/null
+++ b/problems/kh.in
@@ -0,0 +1,62 @@
+# problem name and parameters
+#
+problem              = "kh"
+
+# random number generator parameters
+#
+gentype              = "random"
+
+# physics
+#
+equation_system      = "mhd"
+equation_of_state    = "adi"
+
+# methods
+#
+time_advance         = "rk2"
+riemann_solver       = "hllc"
+reconstruction       = "limo3"
+limiter              = "mc"
+fix_positivity       = "off"
+
+# mesh parameters
+#
+rdimx                =  1
+rdimy                =  1
+rdimz                =  1
+xmin                 = -5.0d-01
+xmax                 =  5.0d-01
+ymin                 = -5.0d-01
+ymax                 =  5.0d-01
+zmin                 = -5.0d-01
+zmax                 =  5.0d-01
+
+# refinement control
+#
+ncells               = 8
+nghost               = 2
+minlev               = 1
+maxlev               = 7
+
+# boundary conditions
+#
+xlbndry              = "periodic"
+xubndry              = "periodic"
+ylbndry              = "periodic"
+yubndry              = "periodic"
+zlbndry              = "periodic"
+zubndry              = "periodic"
+
+# runtime control parameters
+#
+nmax                 = 1000000
+tmax                 = 1.0d+00
+cfl                  = 4.0d-01
+
+# data output control
+#
+precise_snapshots    = "on"
+snapshot_type        = "p"
+snapshot_interval    = 1.0d-01
+restart_number       = -1
+integrals_interval   = 10