Here the inlist used to calculate the models discussed in the paper.
&star_job
! MESA VERSION: 5118
! create_pre_main_sequence_model = .true.
! show_log_description_at_start = .false.
!load_saved_model = .true.
!saved_model_name = 'zams.mod'
!save_model_when_terminate = .true.
!save_model_filename = 'zams.mod'
set_rate_n14pg = 'CF88'
set_rate_3a = 'CF88'
set_rate_1212 = 'CF88'
set_rate_c12ag = 'CF88'
kappa_file_prefix = 'gs98'
change_lnPgas_flag = .false.
new_lnPgas_flag = .false.
!---------------------------------------- INIT ROTATION
set_near_zams_surface_rotation_v_steps = 2 !
new_surface_rotation_v = 1 ! 50 ! 150 ! 300 !(km/sec)
!---------------------------------------- PGSTAR OUTPUT
pgstar_flag = .true.
!----------------------------------------- START OF MS
set_initial_age = .true.
initial_age = 0 ! (yrs)
set_initial_model_number = .true.
initial_model_number = 0
/ ! end of star_job namelist
&controls
!---------------------------------------- GLOBAL PARAMS
initial_mass = 1.5 ! 2.0 ! 3.0
initial_z = 0.02
use_Type2_opacities = .true.
!---------------------------------------- RESOLUTION
mesh_delta_coeff = 0.3
varcontrol_target = 1d-4
max_allowed_nz = 10000
!---------------------------------------- MIXING
mixing_length_alpha = 1.6
use_Henyey_MLT = .true.
use_Ledoux_criterion = .true.
alpha_semiconvection = 1d-3
semiconvection_option = 'Langer_85'
thermo_haline_coeff = 2.0
overshoot_f_above_burn_h = 0.2
overshoot_f0_above_burn_h = 0.
overshoot_step_fraction = 1. ! Step function Overshoot
!---------------------------------------- ELEMENT DIFFUSION
do_element_diffusion = .true.
!---------------------------------------- STOP CRITERION
! Lnuc_div_L_upper_limit = 1.00d0 ! defines "zams"
!---------------------------------------- WIND
RGB_wind_scheme = 'Reimers'
AGB_wind_scheme = 'Blocker'
RGB_to_AGB_wind_switch = 1d-4
Reimers_wind_eta = 1d0
Blocker_wind_eta = 5d0
!---------------------------------------- ROTATION
D_DSI_factor = 1.0
D_SH_factor = 1.0
D_SSI_factor = 1.0
D_ES_factor = 1.0
D_GSF_factor = 1.0
D_ST_factor = 1.0
am_nu_ST_factor = 1.0 ! 100 !(TSx100)
am_D_mix_factor = 0.0333
am_nu_factor = 1.0
am_gradmu_factor = 0.05
! set_uniform_am_nu = .true.
! uniform_am_nu = 6d4 ! cm^2/s
!---------------------------------------- OUTPUT
photostep = 50
profile_interval = 10
history_interval = 10
terminal_cnt = 5
write_header_frequency = 5
max_num_profile_models = 1000
!---------------------------------------- BRUNT
calculate_Brunt_N2 = .true.
!---------------------------------------- FGONG
write_pulse_info_with_profile = .true. ! if true, write pulse info file when write profile
pulse_info_format = 'FGONG' ! e.g., 'FGONG', 'OSC', 'GYRE'
add_atmosphere_to_pulse_info = .false.
add_center_point_to_pulse_info = .true. ! add point for r=0
keep_surface_point_for_pulse_info = .false. ! keep k=1
format_for_FGONG_data = '(1P5E16.9,x)' ! '(1P5E18.9,x)'
!---------------------------------------- OTHER
max_tries = 50 ! max number newton iterations before give up
max_tries1 = 500 ! max tries on 1st model
max_tries_for_retry = 50
max_tries_after_backup = 50
max_tries_after_backup2 = 50
/ ! end of controls namelist