&star_job
  ! see star/defaults/star_job.defaults

  ! begin with a pre-main sequence model
    create_pre_main_sequence_model = .true.
  ! reducing the number of preMS relaxation steps from 300 to 100 to save on some computation time
    pre_ms_relax_num_steps = 100

  ! increasing the first time step after a convective core has formed from 0.00001 years to 1 year to speed up the computation
    set_initial_dt = .true.
    years_for_initial_dt = 1d0

  ! display on-screen plots
    pgstar_flag = .true.

  ! wait for confirmation of the user before closing the pgstar plots
    pause_before_terminate = .true.

/ ! end of star_job namelist


&eos
  ! eos options
  ! see eos/defaults/eos.defaults

/ ! end of eos namelist


&kap
  ! kap options
  ! see kap/defaults/kap.defaults
  use_Type2_opacities = .true.
  Zbase = 0.014d0

/ ! end of kap namelist


&controls
  ! see star/defaults/controls.defaults

  ! starting specifications
    initial_mass = 5d0 ! in Msun units

  ! Set initial fractions according to Galactic enrichtment law
  ! Primordial helium abundance (Aver et al. 2021)
  ! Galactic enrichment ratio. Calibrated by requiring that for solar (Xini, Yini, Zini) = (0.710, 0.276, 0.014)
    initial_z = 0.014d0 ! initial metal mass fraction

  ! when to stop
    ! stop when the star nears ZAMS (Lnuc/L > 0.99)
    Lnuc_div_L_zams_limit = 0.99d0
    stop_near_zams = .false.

    ! stop when the center mass fraction of h1 drops below this limit
    xa_central_lower_limit_species(1) = 'h1'
    xa_central_lower_limit(1) = 1d-6

  ! wind

  ! atmosphere

  ! rotation

  ! element diffusion

  ! mlt

  ! mixing
    overshoot_scheme(1) = 'step'
    overshoot_zone_type(1) = 'burn_H'
    overshoot_zone_loc(1) = 'core'
    overshoot_bdy_loc(1) = 'top'
    overshoot_f(1) = 0.30d0
    overshoot_f0(1) = 0.005d0

  ! timesteps

  ! mesh

  ! solver
  ! options for energy conservation (see MESA V, Section 3)
    energy_eqn_option = 'dedt'
    use_gold_tolerances = .true.

/ ! end of controls namelist
