add new rotating 40Msun to cc test_suite for testing

CODEOWNERS

@@ -130,6 +130,7 @@ star/test_suite/wd_stable_h_burn @wmwolf
 
 star/test_suite/12M_pre_ms_to_core_collapse @Debraheem @aurimontem
 star/test_suite/20M_pre_ms_to_core_collapse @Debraheem @aurimontem
+star/test_suite/40M_rotating_pre_ms_to_core_collapse @Debraheem @mathren
 star/test_suite/split_burn_big_net
 star/test_suite/zams_to_cc_80 @Debraheem
 star/test_suite/20M_z2m2_high_rotation

star/test_suite/40M_rotating_pre_ms_to_core_collapse/README.rst

@@ -0,0 +1,44 @@
+.. _40M_rotating_pre_ms_to_core_collapse:
+
+***************************
+40M_rotating_pre_ms_to_core_collapse
+***************************
+
+This test suite evolves a Z=10^-3, 40 |MSun| model from the pre-ms to core collapse.
+For bit for bit convergence, we recomended to run by using the ./run_all script instead of restarting from models,
+see https://github.com/MESAHub/mesa/issues/610.
+
+This test_suite has been tested up to 80 solar masses, up to solar metallicity, with mass loss, and produces reasonable HR-tracks.
+Note that for higher masses at solar metallicity, some combination of Pextra_factor, mass-loss, and/or superadiabatic convection reduction (e.g. mlt++)
+might be necessary to stabilize the surface and avoid numerical issues. See the 80Msun_zams_to_cc test_suite as an example.
+
+For production science we recommend adopting tighter mesh and timestep controls, such as those suggested in the comments of inlist_common.
+
+Physical checks
+===============
+
+None
+
+Inlists
+=======
+
+This test case has seven parts.
+
+* Part 1 (``inlist_make_late_pre_zams``) creates a 40 |Msun|, Z=10^-3 metallicity, pre-main sequence model and evolves it for 100 years.
+
+* Part 2 (``inlist_to_zams``) evolves the model to the zero age main sequence and then relaxes until omega/omega_crit = 0.6.
+
+* Part 3 (``inlist_to_end_core_he_burn``) takes the model to core helium depletion.
+
+* Part 4 (``inlist_remove_envelope``) removes the remianing hydrogen envelope. (optional)
+
+* Part 5 (``inlist_to_end_core_c_burn``) takes the model to core carbon depletion.
+
+* Part 6 (``inlist_to_lgTmax``) evolves the model until the core temperature reaches log T =9.60 (approximately silicon-shell burning)
+
+* Part 7 (``inlist_to_cc``) evolves until core collapse.
+
+
+
+
+Last-Updated: 10Dec2024 by MR

star/test_suite/40M_rotating_pre_ms_to_core_collapse/ck

@@ -0,0 +1,7 @@
+#!/bin/bash
+
+# this provides the definition of check_one
+# check_one
+source "${MESA_DIR}/star/test_suite/test_suite_helpers"
+
+check_one

star/test_suite/40M_rotating_pre_ms_to_core_collapse/clean

@@ -0,0 +1,4 @@
+#!/bin/bash
+
+cd make
+make clean

star/test_suite/40M_rotating_pre_ms_to_core_collapse/gyre.in

@@ -0,0 +1,32 @@
+&model
+/
+
+&mode
+   l = 0
+/
+
+&osc
+   nonadiabatic = .TRUE.
+/
+
+&rot
+/
+
+&num
+   diff_scheme = 'MAGNUS_GL2' 
+/
+
+&scan
+   grid_type = 'LINEAR'
+   freq_min = 0.5
+   freq_max = 5.0
+   freq_min_units = 'ACOUSTIC_DELTA'
+   freq_max_units = 'ACOUSTIC_DELTA'
+   n_freq = 50
+/
+
+&grid
+   w_osc = 10
+   w_exp = 2
+   w_ctr = 10
+/