&star_job

	read_extra_star_job_inlist(1) = .true.
	extra_star_job_inlist_name(1) = 'inlist_common'

	create_pre_main_sequence_model = .true.
	pre_ms_T_c = 9d5

	! change_mass_years_for_dt = 1d-2
	! set_initial_dt = .true.
	! years_for_initial_dt = 1d-1

	save_model_when_terminate = .true.
	save_model_filename = 'TAMS.mod'

	required_termination_code_string = 'xa_central_lower_limit'

	history_columns_file = ''
	profile_columns_file = ''

    ! Comment to have 1 master history file; Uncomment these if we want to look at stages seperately
    ! set_initial_age = .true.
    ! initial_age = 0
    ! set_initial_model_number = .true.
    ! initial_model_number = 0

    set_initial_cumulative_energy_error = .true.
	new_cumulative_energy_error = 0d0

	pgstar_flag = .false.

/ !end of star_job namelist

&eos

	read_extra_eos_inlist(1) = .true.
	extra_eos_inlist_name(1) = 'inlist_common'

/ !end of eos namelist

&kap

	read_extra_kap_inlist(1) = .true.
	extra_kap_inlist_name(1) = 'inlist_common'

/ !end of kap namelist

&controls

	x_integer_ctrl(1) = 1 ! part number

	read_extra_controls_inlist(1) = .true.
	extra_controls_inlist_name(1) = 'inlist_common'

	log_directory = 'LOGS_MS'
	photo_directory = 'photos_MS'

	x_logical_ctrl(1) = .false.

	max_model_number = 10000

	num_trace_history_values = 3
	trace_history_value_name(1) = 'log_L_div_Ledd'
	trace_history_value_name(2) = 'mass_conv_core'
	trace_history_value_name(3) = 'center_h1'

	xa_central_lower_limit_species(1) = 'h1'
	xa_central_lower_limit(1) = 1d-6

	energy_eqn_option = 'dedt'
	limit_for_rel_error_in_energy_conservation = 1d99
	hard_limit_for_rel_error_in_energy_conservation = 1d99

	convergence_ignore_equL_residuals = .true.

	delta_XH_cntr_limit = 0.0005d0
	delta_lg_XH_cntr_limit = 0.01d0
	mesh_delta_coeff = 1.0

    ! mass_change = -1d-90  ! Rate of accretion (Msun/year). Negative for mass loss. This only applies when the wind_scheme = ''
	
	overshoot_scheme(1) = 'step' !jermyn et al 2022
	overshoot_zone_type(1) = 'any'
	overshoot_zone_loc(1) = 'core'
	overshoot_bdy_loc(1) = 'any'
	overshoot_f(1) = 0.129
	overshoot_f0(1) = 0.0129
	overshoot_D0(1) = 0d0
	overshoot_Delta0(1) = 1d0
	overshoot_scheme(2) = 'exponential'
	overshoot_zone_type(2) = 'nonburn'
	overshoot_zone_loc(2) = 'shell'
	overshoot_bdy_loc(2) = 'any'
	overshoot_f(2) = 0.0174
	overshoot_f0(2) = 0.00174
	overshoot_D0(2) = 0d0
	overshoot_Delta0(2) = 1d0

	do_conv_premix = .true.
	conv_premix_avoid_increase = .true.
	conv_premix_time_factor = 0
	conv_premix_fix_pgas = .true.
	conv_premix_dump_snapshots = .false.

	use_Ledoux_criterion = .true.

	alpha_semiconvection = 0
	semiconvection_option = 'Langer_85 mixing; gradT = gradr'

	num_cells_for_smooth_gradL_composition_term = 20
	threshold_for_smooth_gradL_composition_term = 0.02

	use_superad_reduction = .true.
	superad_reduction_Gamma_limit = 0.5
	superad_reduction_Gamma_limit_scale = 5
	superad_reduction_Gamma_inv_scale = 5
	superad_reduction_diff_grads_limit = 1d-2
	superad_reduction_limit = -1d0

	relax_max_number_retries = 10000

	photo_interval = 500
	profile_interval = 10000
	history_interval = 10
	terminal_interval = 10
	write_header_frequency = 10

	timestep_factor_for_retries = 0.8
	min_timestep_factor = 0.9
	max_timestep_factor = 1.05d0

	retry_hold = 5

	redo_limit = -1

	relax_hard_limits_after_retry = .false.

	scale_max_correction = 0.01d0

	corr_coeff_limit = 1d-1
	ignore_min_corr_coeff_for_scale_max_correction = .true.
	ignore_species_in_max_correction = .true.
	max_resid_jump_limit = 1d99

	restore_mesh_on_retry = .true.
	num_steps_to_hold_mesh_after_retry = 5
	mesh_delta_coeff_for_highT = 1.0d0
	mesh_Pgas_div_P_exponent = 0.5d0

	max_dq = 0.005
	min_dq_for_xa = 1d-5

	max_allowed_nz = 50000
	max_logT_for_k_below_const_q = 100
	max_q_for_k_below_const_q = 0.995
	min_q_for_k_below_const_q = 0.995
	max_logT_for_k_const_mass = 100
	max_q_for_k_const_mass = 0.99
	min_q_for_k_const_mass = 0.99

/ !end of controls namelist

