diff --git a/moment_kinetics/test/fokker_planck_tests.jl b/moment_kinetics/test/fokker_planck_tests.jl
index 8f2dafbd7..cd29e94ca 100644
--- a/moment_kinetics/test/fokker_planck_tests.jl
+++ b/moment_kinetics/test/fokker_planck_tests.jl
@@ -66,8 +66,8 @@ function runtests()
     @testset "Fokker Planck tests" verbose=use_verbose begin
         println("Fokker Planck tests")
         
-        @testset " - test Lagrange-polynomial 2D interpolation" begin
-            println(" - test Lagrange-polynomial 2D interpolation")
+        @testset "Lagrange-polynomial 2D interpolation" begin
+            println("    - test Lagrange-polynomial 2D interpolation")
             ngrid = 9
             nelement_vpa = 16
             nelement_vperp = 8
@@ -140,9 +140,9 @@ function runtests()
             
         end
         
-        @testset " - test weak-form 2D differentiation" begin
+        @testset "weak-form 2D differentiation" begin
         # tests the correct definition of mass and stiffness matrices in 2D
-            println(" - test weak-form 2D differentiation")
+            println("    - test weak-form 2D differentiation")
 
             ngrid = 9
             nelement_vpa = 8
@@ -205,8 +205,8 @@ function runtests()
             finalize_comms!()
         end
         
-        @testset " - test weak-form Rosenbluth potential calculation: elliptic solve" begin
-            println(" - test weak-form Rosenbluth potential calculation: elliptic solve")
+        @testset "weak-form Rosenbluth potential calculation: elliptic solve" begin
+            println("    - test weak-form Rosenbluth potential calculation: elliptic solve")
             ngrid = 9
             nelement_vpa = 8
             nelement_vperp = 4
@@ -353,8 +353,8 @@ function runtests()
             finalize_comms!()                                                                  
         end
         
-        @testset " - test weak-form collision operator calculation" begin
-            println(" - test weak-form collision operator calculation")
+        @testset "weak-form collision operator calculation" begin
+            println("    - test weak-form collision operator calculation")
             ngrid = 9
             nelement_vpa = 8
             nelement_vperp = 4
@@ -510,8 +510,8 @@ function runtests()
             finalize_comms!()                                                                  
         end
         
-        @testset " - test weak-form (slowing-down) collision operator calculation" begin
-            println(" - test weak-form (slowing-down) collision operator calculation")
+        @testset "weak-form (slowing-down) collision operator calculation" begin
+            println("    - test weak-form (slowing-down) collision operator calculation")
             ngrid = 9
             nelement_vpa = 16
             nelement_vperp = 8
@@ -607,8 +607,8 @@ function runtests()
             finalize_comms!()                                                                  
         end
         
-        @testset " - test weak-form Rosenbluth potential calculation: direct integration" begin
-            println(" - test weak-form Rosenbluth potential calculation: direct integration")
+        @testset "weak-form Rosenbluth potential calculation: direct integration" begin
+            println("    - test weak-form Rosenbluth potential calculation: direct integration")
             ngrid = 5 # chosen for a quick test -- direct integration is slow!
             nelement_vpa = 8
             nelement_vperp = 4
diff --git a/moment_kinetics/test/jacobian_matrix_tests.jl b/moment_kinetics/test/jacobian_matrix_tests.jl
index b739bc76a..7e48cc80f 100644
--- a/moment_kinetics/test/jacobian_matrix_tests.jl
+++ b/moment_kinetics/test/jacobian_matrix_tests.jl
@@ -225,7 +225,7 @@ end
 function test_electron_z_advection(test_input; rtol=(2.5e2*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_z_advection"
-    println("    electron_z_advection")
+    println("    - electron_z_advection")
 
     @testset "electron_z_advection" begin
         # Suppress console output while running
@@ -453,7 +453,7 @@ end
 function test_electron_vpa_advection(test_input; rtol=(3.0e2*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_vpa_advection"
-    println("    electron_vpa_advection")
+    println("    - electron_vpa_advection")
 
     @testset "electron_vpa_advection" begin
         # Suppress console output while running
@@ -721,7 +721,7 @@ end
 function test_contribution_from_electron_pdf_term(test_input; rtol=(4.0e2*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_contribution_from_electron_pdf_term"
-    println("    contribution_from_electron_pdf_term")
+    println("    - contribution_from_electron_pdf_term")
 
     @testset "contribution_from_electron_pdf_term" begin
         # Suppress console output while running
@@ -975,7 +975,7 @@ end
 function test_electron_dissipation_term(test_input; rtol=(3.0e0*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_dissipation_term"
-    println("    electron_dissipation_term")
+    println("    - electron_dissipation_term")
 
     @testset "electron_dissipation_term" begin
         # Suppress console output while running
@@ -1207,7 +1207,7 @@ end
 function test_electron_krook_collisions(test_input; rtol=(2.0e1*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_krook_collisions"
-    println("    electron_krook_collisions")
+    println("    - electron_krook_collisions")
 
     @testset "electron_krook_collisions" begin
         # Suppress console output while running
@@ -1445,7 +1445,7 @@ end
 function test_external_electron_source(test_input; rtol=(3.0e1*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_external_electron_source"
-    println("    external_electron_source")
+    println("    - external_electron_source")
 
     @testset "external_electron_source" begin
         # Suppress console output while running
@@ -1718,7 +1718,7 @@ end
 function test_electron_implicit_constraint_forcing(test_input; rtol=(1.5e0*epsilon))
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_implicit_constraint_forcing"
-    println("    electron_implicit_constraint_forcing")
+    println("    - electron_implicit_constraint_forcing")
 
     @testset "electron_implicit_constraint_forcing" begin
         # Suppress console output while running
@@ -1953,7 +1953,7 @@ end
 function test_electron_energy_equation(test_input; rtol=(6.0e2*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_energy_equation"
-    println("    electron_energy_equation")
+    println("    - electron_energy_equation")
 
     @testset "electron_energy_equation" begin
         # Suppress console output while running
@@ -2162,7 +2162,7 @@ end
 function test_ion_dt_forcing_of_electron_ppar(test_input; rtol=(1.5e1*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_ion_dt_forcing_of_electron_ppar"
-    println("    ion_dt_forcing_of_electron_ppar")
+    println("    - ion_dt_forcing_of_electron_ppar")
 
     @testset "ion_dt_forcing_of_electron_ppar" begin
         # Suppress console output while running
@@ -2357,7 +2357,7 @@ end
 function test_electron_kinetic_equation(test_input; rtol=(5.0e2*epsilon)^2)
     test_input = deepcopy(test_input)
     test_input["output"]["run_name"] *= "_electron_kinetic_equation"
-    println("    electron_kinetic_equation")
+    println("    - electron_kinetic_equation")
 
     @testset "electron_kinetic_equation" begin
         # Suppress console output while running