! inlist to evolve a 15 solar mass star

! For the sake of future readers of this file (yourself included),
! ONLY include the controls you are actually using.  DO NOT include
! all of the other controls that simply have their default values.

&star_job
  ! see star/defaults/star_job.defaults

  ! begin with a pre-main sequence model
    create_pre_main_sequence_model = .false. !.true. 

  ! save a model at the end of the run
    save_model_when_terminate = .false.
    save_model_filename = '15M_at_CHeB.mod'

  ! display on-screen plots
    pgstar_flag = .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.02

/ ! end of kap namelist


&controls
  ! see star/defaults/controls.defaults

  ! starting specifications
    initial_mass = 15 ! in Msun units
    initial_z = 0.02

  ! when to stop

    ! stop when the star nears ZAMS (Lnuc/L > 0.99)
    Lnuc_div_L_zams_limit = 0.99d0
    ! stop_near_zams = .true.
    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-3
    
    HB_limit = 0.95



  ! wind
  
  ! atmosphere

  ! rotation

  ! element diffusion

  ! mlt
    mixing_length_alpha = alphasmltscale

  ! mixing
    alpha_semiconvection = 0.01d0
    use_Ledoux_criterion = .true.
    
    
  ! timesteps
    time_delta_coeff = timedeltascale
    max_model_number = 2000 

  ! mesh
    mesh_delta_coeff = meshdeltascale
    max_allowed_nz = 16000 ! default 8000

    history_interval = 1 

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


  ! output

/ ! end of controls namelist
