restore old defaults for H and He in new arrays for abundance-based timestep controls and fix test inlists

Author: warrickball

star/defaults/controls.defaults

@@ -9485,11 +9485,13 @@
       ! Limit on magnitude of decrease in any cell hydrogen abundance during a single timestep.
       ! dH here is ``abs(xa(j,k) - xa_old(j,k))`` for any cell k and all species j, eg 'h1', 'he4', etc.
       ! special 'species' X (any hydrogen), Y (any helium) and Z (any metals) are
-      ! allowed here. Eg. 'Z' will trigger the timestep control if any metal isotope
+      ! allowed here. E.g. 'Z' will trigger the timestep control if any metal isotope
       ! abundance individually satisfies the conditions below.
       ! Considers all cells except where have convective mixing.
 
-    dX_limit_species(:) = ''
+    dX_limit_species(1) = 'h1'
+    dX_limit_species(2) = 'he4'
+    dX_limit_species(3:) = ''
 
       ! dX_limit_min_X
       ! ~~~~~~~~~~~~~~
@@ -9545,7 +9547,9 @@
 
       ! ::
 
-    dX_div_X_limit_min_X(:) = 1d99
+    dX_div_X_limit_min_X(1) = 1d-3
+    dX_div_X_limit_min_X(2) = 1d-3
+    dX_div_X_limit_min_X(3:) = 1d99
 
 
       ! dX_div_X_limit
@@ -9556,7 +9560,9 @@
 
       ! ::
 
-    dX_div_X_limit(:) = 1d99
+    dX_div_X_limit(1) = 0.9d0
+    dX_div_X_limit(2) = 0.9d0
+    dX_div_X_limit(3:) = 1d99
 
 
       ! dX_div_X_hard_limit

star/test_suite/12M_pre_ms_to_core_collapse/inlist_common

@@ -130,7 +130,7 @@
       delta_lgT_max_limit = 0.003     
       delta_lgT_max_hard_limit = 0.01
 
-      dHe_div_He_limit = -1
+      dX_div_X_limit(2) = -1 ! for he4
 
       dX_nuc_drop_limit = 2d-2
       dX_nuc_drop_limit_at_high_T = -1  ! default = -1 = same as dX_nuc_drop_limit

star/test_suite/20M_pre_ms_to_core_collapse/inlist_common

@@ -133,7 +133,7 @@
       delta_lgT_max_limit = 0.002     
       delta_lgT_max_hard_limit = 0.01
 
-      dHe_div_He_limit = -1
+      dX_div_X_limit(2) = -1 ! for he4
 
       dX_nuc_drop_limit = 2d-2
       dX_nuc_drop_limit_at_high_T = -1  ! default = -1 = same as dX_nuc_drop_limit

star/test_suite/c13_pocket/inlist_c13_pocket

@@ -40,11 +40,6 @@
       ! limit max_model_number as part of test_suite
       max_model_number = 2500
 
-! relax temporal resolution to spend less time in interpulse
-
-      dH_div_H_limit_min_H = 1d-3
-      dH_div_H_limit = 0.9
-
 ! stopping condition
       ! one thermal pulse cycle
       ! see run_star_extras

star/test_suite/c13_pocket/inlist_to_lgL_3.79

@@ -42,9 +42,4 @@
       ! limit max_model_number as part of test_suite
       max_model_number = 2500
 
-! relax temporal resolution to spend less time in interpulse
-
-      dH_div_H_limit_min_H = 1d-3
-      dH_div_H_limit = 0.9
-
 / ! end of controls namelist

star/test_suite/pisn/inlist_common

@@ -129,7 +129,7 @@
       delta_lgT_max_limit = 0.002     
       delta_lgT_max_hard_limit = 0.01
 
-      dHe_div_He_limit = -1
+      dX_div_X_limit(2) = -1 ! for he4
 
       dX_nuc_drop_limit = 5d-1 ! 2d-2
       dX_nuc_drop_limit_at_high_T = -1

star/test_suite/radiative_levitation/inlist_radiative_levitation

@@ -77,9 +77,9 @@
       history_interval = 10
       terminal_interval = 10
       write_header_frequency = 10      
-      
-      dH_div_H_limit_min_H = 1d99 ! disable this
-      dHe_div_He_limit_min_He = 1d99 ! disable this
+
+      dX_div_X_limit_min_X(1) = 1d99 ! disable this for H
+      dX_div_X_limit_min_X(2) = 1d99 ! disable this for He
 
       prune_bad_cz_min_Hp_height = 0.1 ! lower limit on radial extent of cz
 

star/test_suite/split_burn_big_net/inlist_common

@@ -118,7 +118,7 @@
       delta_lgT_max_limit = 0.003     
       delta_lgT_max_hard_limit = 0.01
 
-      dHe_div_He_limit = -1
+      dX_div_X_limit(2) = -1 ! for he4
 
       dX_nuc_drop_limit = 2d-2
       dX_nuc_drop_limit_at_high_T = -1  ! default = -1 = same as dX_nuc_drop_limit

star/test_suite/twin_studies/inlist_common

@@ -123,7 +123,7 @@
       delta_lgT_max_limit = 0.003     
       delta_lgT_max_hard_limit = 0.01
 
-      dHe_div_He_limit = -1
+      dX_div_X_limit(2) = -1 ! for he4
 
       dX_nuc_drop_limit = 2d-2
       dX_nuc_drop_limit_at_high_T = -1  ! default = -1 = same as dX_nuc_drop_limit

star/test_suite/wd_stable_h_burn/inlist_wd_stable_h_burn

@@ -96,7 +96,7 @@
 		xa_function_weight(2) = 20
 		xa_function_param(2) = 1d-2
 
-      !dH_div_H_limit = 0.8d0   ! 10.7 101 3 2137
+      ! dX_div_X_limit(1) = 0.8d0   ! 10.7 101 3 2137
 
       delta_lgL_H_limit = 0.1 ! limit for magnitude of change in lgL_H
       delta_lgL_H_hard_limit = 0.5