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

  ! save a model at the end of the run
    save_model_when_terminate = .false.

    load_saved_model = .true.
    load_model_filename = 'middle.mod'

    write_profile_when_terminate = .true.
    filename_for_profile_when_terminate = 'final'

  ! display on-screen plots
    pgstar_flag = .false.
    report_retries = .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
  
  !kap_file_prefix = 'gs98'  
  
  kap_lowT_prefix = 'lowT_fa05_gs98_aFe_p4'
  user_num_kap_lowT_Xs = 10
  user_kap_lowT_Xs = 0.0d0, 0.1d0, 0.2d0, 0.35d0, 0.5d0, 0.7d0, 0.8d0, 0.9d0, 0.95d0, 1.0d0
  user_num_kap_lowT_Zs = 13
  user_kap_lowT_Zs = 0.000d0, 0.0001d0, 0.0003d0, 0.001d0, 0.002d0, 0.004d0, 0.01d0, 0.02d0, 0.03d0, 0.04d0, 0.06d0, 0.08d0, 0.100d0
  user_num_kap_lowT_Xs_for_this_Z = 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 9, 9, 8  

  kap_file_prefix = 'gs98_aFe_p3'
  user_num_kap_Xs = 10
  user_kap_Xs = 0.0d0, 0.1d0, 0.2d0, 0.35d0, 0.5d0, 0.7d0, 0.8d0, 0.9d0, 0.95d0, 1.0d0
  user_num_kap_Zs = 13
  user_kap_Zs = 0.000d0, 0.0001d0, 0.0003d0, 0.001d0, 0.002d0, 0.004d0, 0.01d0, 0.02d0, 0.03d0, 0.04d0, 0.06d0, 0.08d0, 0.100d0
  user_num_kap_Xs_for_this_Z = 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 9, 9, 8
  
  !kap_CO_prefix = 'gs98_aFe_p3_co'
  !user_num_kap_CO_Xs = 10
  !user_kap_CO_Xs = 0.00d0, 0.03d0, 0.10d0, 0.35d0, 0.70d0
  !user_num_kap_CO_Zs = 13
  !user_kap_CO_Zs = 0.000d0, 0.001d0, 0.004d0, 0.010d0, 0.020d0, 0.030d0, 0.050d0, 0.100d0
  !user_num_kap_CO_Xs_for_this_Z = 5, 5, 5, 5, 5, 5, 5, 5

  use_Type2_opacities = .false.
  !Zbase = 0.02
  
  include_electron_conduction = .true.
  
/ ! end of kap namelist


&controls
    ! atm
    atm_option = 'T_tau'
    atm_T_tau_relation = 'Eddington'
    atm_T_tau_opacity = 'varying'
    atm_build_tau_outer = 1d-4
    
    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
  mixing_length_alpha = 1.9195933809594632
  log_directory = 'LOGS_alpha_0.31' 

  do_element_diffusion = .true.

  ! when to stop

  ! stop when the temperature goes below this limit
   log_g_lower_limit = 2.50
   log_L_upper_limit = 1.2
   
  ! solver
     ! options for energy conservation (see MESA V, Section 3)
     energy_eqn_option = 'dedt'
     use_gold_tolerances = .true.

  ! output
    history_interval = 1
    profile_interval = 1
    max_num_profile_models = 10000

  ! timestepping
   time_delta_coeff = 0.5

  ! envelope overshoot
    overshoot_scheme(1) = 'exponential'
    overshoot_zone_type(1) = 'any'
    overshoot_zone_loc(1) = 'shell'
    overshoot_bdy_loc(1) = 'bottom'
  
    overshoot_f0(1) = 0.005
    overshoot_f(1) = 0.01


/ ! end of controls namelist
