! 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 = .true.
    pre_ms_T_c = 5d5

  ! save a model at the end of the run
    save_model_when_terminate = .true.
    save_model_filename = 'zams.mod'

  ! display on-screen plots
    pgstar_flag = .false.
    report_retries = .true.

/ ! en of star_job namelist

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

/ ! end of eos namelist


&kap
  ! kap options
  ! see kap/defaults/kap.defaults
  kap_file_prefix = 'gs98'
  use_Type2_opacities = .false.
  ! Zbase = 0.02

/ ! end of kap namelist


&controls
    max_number_retries = 10000
    
    ! gyre stuff
    write_pulse_data_with_profile = .true.
    pulse_data_format = 'GYRE'
    add_atmosphere_to_pulse_data = .true.
    add_center_point_to_pulse_data = .true.
    keep_surface_point_for_pulse_data = .true.
    add_double_points_to_pulse_data = .true.
    get_delta_nu_from_scaled_solar = .true.


  ! see star/defaults/controls.defaults

  ! starting specifications
    initial_mass = 0.7779340840137416
    initial_z = 0.00017742603878168832
    initial_y = 0.24556187651072667

    z_fraction_li = 0.0
    z_fraction_be = 0.0
    z_fraction_b = 0.0
    z_fraction_c = 0.099760636286597
    z_fraction_n = 0.029223025171423916
    z_fraction_o = 0.5539529530759232
    z_fraction_f = 1.7301666115942203e-05
    z_fraction_ne = 0.12425137555602701
    z_fraction_na = 0.0012328749109602726
    z_fraction_mg = 0.04732337052269569
    z_fraction_al = 0.0019973500346744623
    z_fraction_si = 0.051033396286330045
    z_fraction_p = 0.00021896632904370327
    z_fraction_s = 0.035102945639339383
    z_fraction_cl = 0.00028118767643664657
    z_fraction_ar = 0.002516952877250199
    z_fraction_k = 0.0001292809651783743
    z_fraction_ca = 0.004702037456881943
    z_fraction_sc = 1.6678911567704386e-06
    z_fraction_ti = 0.0002566943164981959
    z_fraction_v = 1.2777798493536745e-05
    z_fraction_cr = 0.0006100282938470821
    z_fraction_mn = 0.00033826174778478686
    z_fraction_fe = 0.044295982484074826
    z_fraction_co = 0.00012295284803509143
    z_fraction_ni = 0.0026179801652320505
    z_fraction_cu = 0.0
    z_fraction_zn = 0.0

  mixing_length_alpha = 1.9195933809594632

  log_directory = 'LOGS_alpha_to_zams_0.31' 

  do_element_diffusion = .false.

  ! when to stop

  ! at zams
    stop_near_zams = .true.
    Lnuc_div_L_zams_limit = 0.99d0

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

  ! output
    history_interval = 1
    profile_interval = 0

/ ! end of controls namelist
