! 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.

  ! save a model at the end of the run
    save_model_when_terminate = .true.
    save_model_filename = 'output.mod'
    write_profile_when_terminate = .true.
`   filename_for_profile_when_terminate = 'final'



  ! display on-screen plots
    pgstar_flag = .false.

/ ! 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
  use_Type2_opacities = .true.
  Zbase = 0.02

/ ! end of kap namelist


&controls
    ! 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 = 1.0
    initial_z = 0.02,
    initial_y = 0.248

    z_fraction_li = 0.0
    z_fraction_be = 0.0
    z_fraction_b = 0.0
    z_fraction_c = 0.17101576181026934
    z_fraction_n = 0.0500958905046917
    z_fraction_o = 0.46510371751090424
    z_fraction_f = 2.965957036304816e-05
    z_fraction_ne = 0.10432253561618107
    z_fraction_na = 0.002113469299743889
    z_fraction_mg = 0.03973311349463119
    z_fraction_al = 0.0034239791414352938
    z_fraction_si = 0.042848083394414496
    z_fraction_p = 0.00037536542434060276
    z_fraction_s = 0.029472738473157494
    z_fraction_cl = 0.000482029049607548
    z_fraction_ar = 0.004314713997081053
    z_fraction_k = 0.00022162130846910976
    z_fraction_ca = 0.003947871545639328
    z_fraction_sc = 2.859200656782797e-06
    z_fraction_ti = 0.00021552278077822373
    z_fraction_v = 2.1904480814983374e-05
    z_fraction_cr = 0.001045747674447158
    z_fraction_mn = 0.0005798689006202789
    z_fraction_fe = 0.07593487242689285
    z_fraction_co = 0.00021077326444727607
    z_fraction_ni = 0.00448790113041295
    z_fraction_cu = 0.0
    z_fraction_zn = 0.0

  mixing_length_alpha = 2.0,

  log_directory = 'LOGS_alpha_0.0' 

  do_element_diffusion = .true.

  ! when to stop

    ! stop when the star nears ZAMS (Lnuc/L > 0.99)
    !Lnuc_div_L_zams_limit = 0.99d0
    ! stop_near_zams = .true.
    max_age = 4.57e9


  ! 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
  ! 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
