!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! inlist_OB
!
! This inlist runs a very high mass, solar-metallicity, star with no termination age
! This is a modified inlist file based on MESA's high_mass test suite
! It DOES NOT require a pre-generated MESA model
! You must change the directory strings to point to
!     - your MESA2HYDRO data directory
!     - your MESA directory
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

&star_job

      mesa_dir = '/home/mjoyce/MESA/mesa-r10398'
      history_columns_file='/home/mjoyce/MESA2HYDRO/data/MESA_controls/history_columns_testsuite.list'
      profile_columns_file='/home/mjoyce/MESA2HYDRO/data/MESA_controls/profile_columns_testsuite.list'

      save_model_when_terminate = .false.
      !save_model_filename = 'final.mod'

      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'profile_OB_solarmet.data'

      create_pre_main_sequence_model = .true.
      show_log_description_at_start = .false.

      pgstar_flag = .false.
      save_pgstar_files_when_terminate = .false.


/ ! end of star_job namelist

&controls
      star_history_name = 'history_OB_solarmet.data'

      use_gold_tolerances = .true.
      use_eosDT2 = .true.

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 0
         max_number_retries = 0 

      max_model_number = 1000  

      initial_mass = 90
      initial_z = 0.02  !1d-5
      Zbase = 0.02      !1d-5
      ! max_age = 5d5 !500,000 if we want this to be on the main sequence
      ! max_age=1d6

      mixing_length_alpha = 2

      conv_core_gap_dq_limit = 0.1

      xa_central_lower_limit_species(1) = 'h1'
      xa_central_lower_limit(1) = 0.5
           
      !min_timestep_limit = 1d0 ! (seconds)

      overshoot_f_above_nonburn_core = 0.014
      overshoot_f_above_nonburn_shell = 0.014
      overshoot_f_below_nonburn_shell = 0.128
      overshoot_f_above_burn_h_core = 0.014
      overshoot_f_above_burn_h_shell = 0.014
      overshoot_f_below_burn_h_shell = 0.014
      overshoot_f_above_burn_he_core = 0.014
      overshoot_f_above_burn_he_shell = 0.014
      overshoot_f_below_burn_he_shell = 0.014            

      overshoot_f0_above_nonburn_core = 0.004
      overshoot_f0_above_nonburn_shell = 0.004
      overshoot_f0_below_nonburn_shell = 0.100
      overshoot_f0_above_burn_h_core = 0.004
      overshoot_f0_above_burn_h_shell = 0.004
      overshoot_f0_below_burn_h_shell = 0.004
      overshoot_f0_above_burn_he_core = 0.004
      overshoot_f0_above_burn_he_shell = 0.004
      overshoot_f0_below_burn_he_shell = 0.004            

      overshoot_scheme(1) = 'exponential'
      overshoot_zone_type(1) = 'nonburn'
      overshoot_zone_loc(1) = 'shell'
      overshoot_bdy_loc(1) = 'bottom'
      overshoot_f(1) = 0.128
      overshoot_f0(1) = 0.100

      overshoot_scheme(2) = 'exponential'
      overshoot_zone_type(2) = 'any'
      overshoot_zone_loc(2) = 'any'
      overshoot_bdy_loc(2) = 'any'
      overshoot_f(2) = 0.014
      overshoot_f0(2) = 0.004

      overshoot_new = .true.


      photo_interval = 50
      profile_interval = 50
      history_interval = 1
      terminal_interval = 10
      write_header_frequency = 100

/ ! end of controls namelist

