AMReX-Astro · zingale · Sep 23, 2023 · Sep 18, 2023 · Sep 18, 2023 · Sep 18, 2023
diff --git a/Docs/source/sdc.rst b/Docs/source/sdc.rst
@@ -47,7 +47,7 @@ The options that affect the nonlinear solve are:
   ``integrator.atol_spec``, ``integrator.rtol_enuc``,
   ``integrator.atol_enuc``.
 
-In all cases, the type of Jacobian (analytic or numerical) is determined by
+  In all cases, the type of Jacobian (analytic or numerical) is determined by
   ``integrator.jacobian``.
 
 

diff --git a/Exec/reacting_tests/bubble_convergence/inputs_2d.128 b/Exec/reacting_tests/bubble_convergence/inputs_2d.128
@@ -29,7 +29,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 
 castro.grav_source_type = 2

diff --git a/Exec/reacting_tests/bubble_convergence/inputs_2d.256 b/Exec/reacting_tests/bubble_convergence/inputs_2d.256
@@ -29,7 +29,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 
 castro.grav_source_type = 2

diff --git a/Exec/reacting_tests/bubble_convergence/inputs_2d.32 b/Exec/reacting_tests/bubble_convergence/inputs_2d.32
@@ -29,7 +29,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 
 castro.grav_source_type = 2

diff --git a/Exec/reacting_tests/bubble_convergence/inputs_2d.512 b/Exec/reacting_tests/bubble_convergence/inputs_2d.512
@@ -29,7 +29,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 
 castro.grav_source_type = 2

diff --git a/Exec/reacting_tests/bubble_convergence/inputs_2d.64 b/Exec/reacting_tests/bubble_convergence/inputs_2d.64
@@ -29,7 +29,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 
 castro.grav_source_type = 2

diff --git a/Exec/reacting_tests/bubble_convergence/job_scripts/cori.bubble_convergence.slurm b/Exec/reacting_tests/bubble_convergence/job_scripts/cori.bubble_convergence.slurm
@@ -21,7 +21,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1"
 
 srun ${CASTRO_EXEC} inputs_2d.64 ${RUNPARAMS} amr.plot_file=bubble_64_sdc4_plt &> 64.out

diff --git a/Exec/reacting_tests/reacting_convergence/convergence_sdc.sh b/Exec/reacting_tests/reacting_convergence/convergence_sdc.sh
@@ -10,7 +10,6 @@ RUNPARAMS="
 castro.sdc_order=2
 castro.time_integration_method=2
 castro.ppm_type=0
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 castro.use_retry=0"
 

diff --git a/Exec/reacting_tests/reacting_convergence/convergence_sdc4.sh b/Exec/reacting_tests/reacting_convergence/convergence_sdc4.sh
@@ -11,7 +11,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1
 castro.use_retry=0"
 

diff --git a/Exec/reacting_tests/reacting_convergence/job_scripts/cori.reacting_convergence.slurm b/Exec/reacting_tests/reacting_convergence/job_scripts/cori.reacting_convergence.slurm
@@ -26,7 +26,6 @@ castro.sdc_order=4
 castro.time_integration_method=2
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=1"
 
 srun ${CASTRO_EXEC} inputs.64 amr.plot_file=react_converge_64_sdc4_plt ${RUNPARAMS}

diff --git a/Exec/science/Detonation/nse_runs/inputs.template.true_sdc b/Exec/science/Detonation/nse_runs/inputs.template.true_sdc
@@ -38,7 +38,6 @@ castro.sdc_quadrature = @@SDC_QUADRATURE@@
 
 castro.limit_fourth_order=1
 castro.use_reconstructed_gamma1=1
-castro.sdc_solve_for_rhoe=1
 castro.sdc_solver=2
 
 

diff --git a/Exec/science/flame/inputs.1d.sdc b/Exec/science/flame/inputs.1d.sdc
@@ -25,7 +25,6 @@ castro.do_react = 1
 castro.time_integration_method = 2
 castro.sdc_order = 4
 
-castro.sdc_solve_for_rhoe = 1
 castro.sdc_solver = 1
 castro.sdc_solver_relax_factor = 1
 

diff --git a/Source/driver/_cpp_parameters b/Source/driver/_cpp_parameters
@@ -267,9 +267,6 @@ sdc_extra                    int           0
 # which SDC nonlinear solver to use?  1 = Newton, 2 = VODE, 3 = VODE for first iter
 sdc_solver                   int           1
 
-# do we solve for (rho e) or (rho E) in the SDC nonlinear solve?
-sdc_solve_for_rhoe           int           1
-
 # for 2-d axisymmetry, do we include the geometry source terms from Bernand-Champmartin?
 use_axisymmetric_geom_source int           1