diff --git a/.github/workflows/main.yml b/.github/workflows/main.yml
index 7592ed42..54c87997 100644
--- a/.github/workflows/main.yml
+++ b/.github/workflows/main.yml
@@ -62,6 +62,14 @@ jobs:
         cd "$GITHUB_WORKSPACE"/tests/c2001_case0_FD_uniform
         mpirun -np 4 ../../bin/rayleigh.dbg
 
+        # J2011 steady hydro test
+        cd "$GITHUB_WORKSPACE"/tests/j2011_steady_hydro_minimal
+        mpirun -np 4 ../../bin/rayleigh.dbg
+
+        # J2011 steady mhd test
+        cd "$GITHUB_WORKSPACE"/tests/j2011_steady_mhd_minimal
+        mpirun -np 4 ../../bin/rayleigh.dbg
+
         # Generic input test
         cd "$GITHUB_WORKSPACE"
         sh ./tests/generic_input/run_test.sh
diff --git a/tests/j2011_steady_hydro_minimal/Benchmark_Reports/00000020 b/tests/j2011_steady_hydro_minimal/Benchmark_Reports/00000020
new file mode 100644
index 00000000..6d7c7013
--- /dev/null
+++ b/tests/j2011_steady_hydro_minimal/Benchmark_Reports/00000020
@@ -0,0 +1,22 @@
+ ///////////////////////////////////////////////////////////////////////////
+               RAYLEIGH ACCURACY BENCHMARK SUMMARY               
+  
+   Benchmark:  Jones et al. 2001  (Hydrodynamic Case)
+  
+   Radial Resolution      N_R =          128
+   Angular Resolution N_theta =          192
+  
+   Averaging Interval (seconds) :       0.000000
+  
+   Beginning Iteration :           20
+   Ending Iteration    :           20
+   Number of Samples   :            1
+ ----------------------------------------------------------------------------
+   Observable      |    Measured    | Suggested   | % Difference |  Std. Dev.
+ ----------------------------------------------------------------------------
+   Kinetic Energy  :  5.384025E+24  5.570280E+35 -1.000000E+02  0.000000E+00
+   Zonal KE        :  5.518963E+04  6.380990E+34 -1.000000E+02  0.000000E+00
+   Meridional KE   :  4.598309E+07  1.498250E+32 -1.000000E+02  0.000000E+00
+   Entropy         :  8.006096E+05  7.945200E+05  7.664490E-01  0.000000E+00
+   Vphi            : -7.725785E-03  6.902700E+02 -1.000011E+02  0.000000E+00
+   Drift Frequency :  2.510466+273  3.105120E-06  8.084924+280  0.000000E+00
diff --git a/tests/j2011_steady_hydro_minimal/main_input b/tests/j2011_steady_hydro_minimal/main_input
new file mode 100644
index 00000000..85e41239
--- /dev/null
+++ b/tests/j2011_steady_hydro_minimal/main_input
@@ -0,0 +1,34 @@
+&problemsize_namelist
+ n_r = 128
+ n_theta = 192
+ nprow = 2
+ npcol = 2
+/
+&numerical_controls_namelist
+/
+&physical_controls_namelist
+ benchmark_mode = 3
+ benchmark_integration_interval = 20
+ benchmark_report_interval = 20
+/
+&temporal_controls_namelist
+ max_iterations = 20
+ checkpoint_interval = 100000
+ quicksave_interval = 10000
+ num_quicksaves = 2
+/
+&io_controls_namelist
+/
+&output_namelist
+/
+
+&Boundary_Conditions_Namelist
+/
+&Initial_Conditions_Namelist
+/
+&Test_Namelist
+/
+&Reference_Namelist
+/
+&Transport_Namelist
+/
diff --git a/tests/j2011_steady_mhd_minimal/Benchmark_Reports/00000020 b/tests/j2011_steady_mhd_minimal/Benchmark_Reports/00000020
new file mode 100644
index 00000000..481991b7
--- /dev/null
+++ b/tests/j2011_steady_mhd_minimal/Benchmark_Reports/00000020
@@ -0,0 +1,26 @@
+ ///////////////////////////////////////////////////////////////////////////
+               RAYLEIGH ACCURACY BENCHMARK SUMMARY               
+  
+   Benchmark:  Jones et al. 2001  (Steady Dynamo Case)
+  
+   Radial Resolution      N_R =          128
+   Angular Resolution N_theta =          192
+  
+   Averaging Interval (seconds) :       0.000000
+  
+   Beginning Iteration :           20
+   Ending Iteration    :           20
+   Number of Samples   :            1
+ ----------------------------------------------------------------------------
+   Observable      |    Measured    | Suggested   | % Difference |  Std. Dev.
+ ----------------------------------------------------------------------------
+   Kinetic Energy  :  2.305238E+29  8.036230E+36 -1.000000E+02  0.000000E+00
+   Zonal KE        :  7.929201E+26  1.153180E+36 -1.000000E+02  0.000000E+00
+   Meridional KE   :  2.754535E+27  1.015870E+33 -9.999973E+01  0.000000E+00
+   Magnetic Energy :  7.182914E+30  6.133330E+36 -9.999988E+01  0.000000E+00
+   Zonal ME        :  5.878429E+28  4.620460E+36 -1.000000E+02  0.000000E+00
+   Meridional ME   :  1.934906E+28  3.249270E+35 -9.999999E+01  0.000000E+00
+   Entropy         :  5.897223E+05  6.089300E+05 -3.154335E+00  0.000000E+00
+   Vphi            :  7.775163E-02 -2.942200E+03 -1.000026E+02  0.000000E+00
+   Btheta          :  3.824210E-01  2.729200E+02 -9.985988E+01  0.000000E+00
+   Drift Frequency :  1.385626E+37  4.307600E-06  3.216701E+44  0.000000E+00
diff --git a/tests/j2011_steady_mhd_minimal/main_input b/tests/j2011_steady_mhd_minimal/main_input
new file mode 100644
index 00000000..0ffd61ac
--- /dev/null
+++ b/tests/j2011_steady_mhd_minimal/main_input
@@ -0,0 +1,34 @@
+&problemsize_namelist
+ n_r = 128
+ n_theta = 192
+ nprow = 2
+ npcol = 2
+/
+&numerical_controls_namelist
+/
+&physical_controls_namelist
+ benchmark_mode = 4
+ benchmark_integration_interval = 20
+ benchmark_report_interval = 20
+/
+&temporal_controls_namelist
+ max_iterations = 20
+ checkpoint_interval = 100000
+ quicksave_interval  = 25000
+ num_quicksaves = 2
+/
+&io_controls_namelist
+/
+&output_namelist
+/
+
+&Boundary_Conditions_Namelist
+/
+&Initial_Conditions_Namelist
+/
+&Test_Namelist
+/
+&Reference_Namelist
+/
+&Transport_Namelist
+/