Merge branch 'master' into binaries

problems/current-sheet.in

@@ -0,0 +1,41 @@
+# problem name and parameters
+#
+problem              = "current-sheet"
+resistivity          = 1.00d-05
+
+# physics
+#
+equation_system      = "mhd"
+equation_of_state    = "adi"
+glm_source_terms     = "kepes"
+
+# methods
+#
+time_advance         = "ssprk3(2)4"
+riemann_solver       = "kepes"
+reconstruction       = "gp"
+
+# mesh parameters
+#
+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               = 16
+maxlev               = 8
+refinement_variables = "dens pres magy"
+
+# runtime control parameters
+#
+tmax                 = 1.0d+00
+cfl                  = 5.0d-01
+
+# data output control
+#
+snapshot_interval    = 1.0d-01
+compression_format   = "lz4"

problems/current_sheet.in

@@ -1,61 +0,0 @@
-# problem name and parameters
-#
-problem              = "current_sheet"
-resistivity          = 1.00d-05
-alpha                = 1.00d+01
-
-# physics
-#
-equation_system      = "mhd"
-equation_of_state    = "adi"
-glm_source_terms     = "heglm"
-
-# methods
-#
-time_advance         = "rk2"
-riemann_solver       = "hlld"
-reconstruction       = "limo3"
-limiter              = "mc"
-fix_positivity       = "off"
-
-# mesh parameters
-#
-xblocks              =  1
-yblocks              =  1
-zblocks              =  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
-nghosts              = 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

sources/interpolations.F90

@@ -224,7 +224,7 @@ module interpolations
 
 ! calculate κ = (1 - ν) / ν
 !
-    kappa = min(kappa, (1.0d+00 - cfl) / cfl)
+    kappa = max(kappa, (1.0d+00 - cfl) / cfl)
 
 ! correct central weight
 !