
&star_job
      show_log_description_at_start = .false.

      load_saved_model = .true.
      load_model_filename = 'start_RC_noMloss.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'end_core_he_burn_noMloss.mod'
      required_termination_code_string = 'xa_central_lower_limit'

      pgstar_flag = .true.

/ ! end of star_job namelist

&eos


/ ! end of eos namelist

&kap
      Zbase = 0.02d0

      kap_file_prefix = 'gs98'
      use_Type2_opacities = .true.
      kap_lowT_prefix = 'lowT_fa05_gs98'

/ ! end of kap namelist

&controls      

      xa_central_lower_limit_species(1) = 'he4'
      xa_central_lower_limit(1) = 0.01

      energy_eqn_option = 'eps_grav'
      use_gold2_tolerances = .true.

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    ! To stop the happening of the core breathing pulses, by neglecting the
    ! eps_g term. It yields models that are less physically consistent (since the
    ! luminosity caused by core contraction is ignored), but that do not exhibit
    ! core breathing pulses. See the note in the Convective Premixing section.
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

      eps_grav_factor = 0.05d0 ! if 0, the models do not converge

    ! Ignore the energy residuals that are caused by neglecting eps_g
    ! If not, the models do not converge
      max_abs_rel_run_E_err = -10d0
      limit_for_rel_error_in_energy_conservation = 1d-1
      hard_limit_for_rel_error_in_energy_conservation = 1d-1
      
      convergence_ignore_equL_residuals = .true.

      ! limit max_model_number as part of test_suite
      max_model_number = 18000
         
      initial_mass = 1.0
      initial_z = 0.02d0

      use_Ledoux_criterion = .true. 

      varcontrol_target = 1d-3
      delta_lgL_He_limit = 0.01d0

    ! Convective boundary mixing prescriptions:

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    ! 1. OVERSHOOTING
    ! Required parameters define overshooting above the convective He burning core
    ! (gradT = gradr in the overshoot region)
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

    !overshoot_scheme(1) = 'step'
    !overshoot_zone_type(1) = 'burn_He'
    !overshoot_zone_loc(1) = 'core'
    !overshoot_bdy_loc(1) = 'top'
    !overshoot_f(1) = 1.01              ! in unit of Hp
    !overshoot_f0(1) = 0.01             ! in unit of Hp

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    ! 2. PENETRATIVE CONVECTION
    ! This scheme needs the overshooting parameters to stay active!
    ! (gradT = grada in the overshoot region)
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

    x_logical_ctrl(1) = .false.

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    ! 3. CONVECTIVE PREMIXING
    ! Convective premixing ensures that gradr = grada at the inner boundary of
    ! the convective core. It leads to a semi-convective region (gradr = grada) at
    ! the external boundary of the overshoot region in some cases. It may yields to
    ! potential core breathing pulses (CBP). See the note in the Convective Premixing 
    ! section.
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

    do_conv_premix = .false.

    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 
    ! 4. MAXIMAL OVERSHOOT (BONUS EXERCICE)
    ! The scheme is not directly implemented in the MESA,
    ! but the "maximal overshoot" structure is rather reached through the
    ! use of the predictive mixing scheme. 
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 

       predictive_mix(1) = .true.
    predictive_zone_type(1) = 'any'
    predictive_zone_loc(1) = 'core'
    predictive_bdy_loc(1) = 'top'
    predictive_superad_thresh(1) = 5d-3
    predictive_avoid_reversal(1) = 'he4' ! the minimize the occurrence of CBP

  ! output    

      photo_interval = 50
      profile_interval = 10
      history_interval = 5
      terminal_interval = 10
      write_header_frequency = 10

  ! change the LOGS directory
    log_directory = 'max_overshoot/LOGS'

  ! change the history filename
   star_history_name = 'history_1M_maxOV.data'

/ ! end of controls namelist


&pgstar
   
/ ! end of pgstar namelist
