add new rotating 40Msun to cc test_suite for testing

CODEOWNERS

@@ -6,32 +6,32 @@
 adipls/ @rhdtownsend @warrickball
 astero/ @warrickball
 atm/ @rhdtownsend @warrickball
-autodiff/ 
+autodiff/
 binary/ @orlox
 chem/ @fxt44
-colors/ 
+colors/
 const/ @fxt44
-data/ 
+data/
 docs/ @rhdtownsend
 eos/ @fxt44 @evbauer
 gyre/ @rhdtownsend @warrickball
 hdfio5 @rhdtownsend
-interp_1d/ @rhdtownsend 
-interp_2d/ @rhdtownsend 
+interp_1d/ @rhdtownsend
+interp_2d/ @rhdtownsend
 ionization/ @fxt44
 kap/ @fxt44 @Debraheem
-linters/ 
-turb/ 
+linters/
+turb/
 math/ @rhdtownsend
 mtx/  @mjoyceGR
 net/ @fxt44 @Debraheem
 neu/ @fxt44 @Debraheem
-num/ 
+num/
 rates/ @fxt44 @Debraheem
 #star/
 #star_data/
-#stella/ 
-utils/ @rhdtownsend @warrickball 
+#stella/
+utils/ @rhdtownsend @warrickball
 
 ### Star ###
 
@@ -40,7 +40,7 @@ star/private/phase_separation.f90 @evbauer
 
 ### Test suite cases ###
 
-star/test_suite/1.3M_ms_high_Z @fxt44 
+star/test_suite/1.3M_ms_high_Z @fxt44
 star/test_suite/1.4M_ms_op_mono @fxt44
 star/test_suite/15M_dynamo @fxt44
 star/test_suite/1.5M_with_diffusion @annethoul
@@ -52,20 +52,20 @@ star/test_suite/20M_z2m2_high_rotation
 star/test_suite/5M_cepheid_blue_loop @annethoul @rsmolec
 star/test_suite/7M_prems_to_AGB @annethoul
 star/test_suite/accreted_material_j @orlox
-star/test_suite/adjust_net 
+star/test_suite/adjust_net
 star/test_suite/c13_pocket @fxt44
 star/test_suite/carbon_kh @evbauer
-star/test_suite/cburn_inward 
+star/test_suite/cburn_inward
 star/test_suite/ccsn_IIp @fxt44 @aurimontem
 star/test_suite/check_redo @fxt44
 star/test_suite/check_pulse_atm @warrickball
-star/test_suite/conductive_flame @fxt44 
+star/test_suite/conductive_flame @fxt44
 star/test_suite/conserve_angular_momentum @orlox
 star/test_suite/conv_core_cpm @evbauer
-star/test_suite/custom_colors 
+star/test_suite/custom_colors
 star/test_suite/custom_rates @evbauer  @fxt44
 star/test_suite/high_rot_darkening @orlox
-star/test_suite/dev_ppisn 
+star/test_suite/dev_ppisn
 star/test_suite/diffusion_smoothness @fxt44
 star/test_suite/extended_convective_penetration @fxt44
 star/test_suite/gyre_in_mesa_bcep @rhdtownsend @annethoul
@@ -75,16 +75,16 @@ star/test_suite/gyre_in_mesa_rsg @rhdtownsend @annethoul
 star/test_suite/gyre_in_mesa_spb @rhdtownsend @annethoul
 star/test_suite/gyre_in_mesa_wd @rhdtownsend @annethoul
 star/test_suite/hb_2M @annethoul
-star/test_suite/high_mass @fxt44 
+star/test_suite/high_mass @fxt44
 star/test_suite/high_z @fxt44
 star/test_suite/hot_cool_wind @fxt44
 star/test_suite/hse_riemann @fxt44
 star/test_suite/irradiated_planet @fxt44
 star/test_suite/low_z @aarondotter
-star/test_suite/magnetic_braking 
+star/test_suite/magnetic_braking
 star/test_suite/make_brown_dwarf @fxt44
 star/test_suite/make_co_wd @wmwolf @evbauer @fxt44
-star/test_suite/make_env 
+star/test_suite/make_env
 star/test_suite/make_he_wd @evbauer @fxt44
 star/test_suite/make_metals @fxt44
 star/test_suite/make_o_ne_wd @fxt44 @evbauer
@@ -93,12 +93,12 @@ star/test_suite/make_pre_ccsn_13bvn @Debraheem @aurimontem
 star/test_suite/make_sdb @evbauer
 star/test_suite/make_zams @fxt44
 star/test_suite/make_zams_low_mass @fxt44
-star/test_suite/make_zams_ultra_high_mass @fxt44 
+star/test_suite/make_zams_ultra_high_mass @fxt44
 star/test_suite/ns_c @fxt44
 star/test_suite/ns_h @fxt44
 star/test_suite/ns_he @fxt44
 star/test_suite/other_physics_hooks @fxt44
-star/test_suite/pisn 
+star/test_suite/pisn
 star/test_suite/radiative_levitation @evbauer
 star/test_suite/R_CrB_star @fxt44
 star/test_suite/relax_composition_j_entropy @warrickball
@@ -115,8 +115,8 @@ star/test_suite/rsp_Type_II_Cepheid @rsmolec
 star/test_suite/semiconvection  @annethoul
 star/test_suite/simplex_solar_calibration @mjoyceGR
 star/test_suite/test_case_template @fxt44
-star/test_suite/test_memory 
-star/test_suite/timing 
+star/test_suite/test_memory
+star/test_suite/timing
 star/test_suite/T_tau_gradr @warrickball
 star/test_suite/twin_studies @fxt44
 star/test_suite/wd_acc_small_dm @wmwolf
@@ -130,19 +130,20 @@ 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
 
 
-binary/test_suite/double_bh @orlox 
+binary/test_suite/double_bh @orlox
 binary/test_suite/evolve_both_stars @orlox @mjoyceGR
 binary/test_suite/jdot_gr_check @orlox
 binary/test_suite/jdot_ls_check @orlox
 binary/test_suite/jdot_ml_check @orlox
-binary/test_suite/star_plus_point_mass @orlox 
+binary/test_suite/star_plus_point_mass @orlox
 binary/test_suite/star_plus_point_mass_explicit_mdot @orlox
-binary/test_suite/wind_fed_bhhmxb @orlox 
+binary/test_suite/wind_fed_bhhmxb @orlox
 
 astero/test_suite/astero_adipls @warrickball
 astero/test_suite/astero_gyre @warrickball

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
+/