&star_job

      mesa_dir = ''


      save_model_when_terminate = .true.
      !save_model_filename = 'end_h.mod'
      save_model_filename = 'end_oburn.mod'

      !load_saved_model = .true.
      !saved_model_name = 'end_h.mod'

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

      pre_ms_relax_num_steps = 200

      !## reaction network
      change_net = .true.
      new_net_name = 'basic.net' ! basic network; switches to bigger nets when
      ! certain conditions are met in order to reach oxygen depletion with the 
      ! 'auto_extend_net = .true.'
      ! new_net_name = 'o_burn_full.net' ! trunkated reaction network up to ca46
      ! new_net_name = 'approx21_plus_co56.net' ! use smaller network to make test for massive stars only

      !## opacity
      !kappa_file_prefix = 'gs98'
      !initial_zfracs = 3
      kappa_file_prefix = 'a09'
      kappa_lowT_prefix = 'lowT_fa05_a09p'
      initial_zfracs = 6
      kappa_CO_prefix = 'a09_co'

      !## reaction rates
      set_rates_preference = .true. ! for use by net + rates modules; default is false
      new_rates_preference = 2 ! latest jina includes most of the NACRE rates tbc
         ! 1 = NACRE rates -- this is the default
         ! 2 = jina reaclib rates -- to match jina where possible
 
      !jina_reaclib_filename = 'jina_reaclib_results_20171204default'
         !! latest version from JINA Reaclib Database
      !set_rate_c12ag = 'jina reaclib' ! jina_reaclib usually uses Kunz 
         !! one of 'NACRE', 'jina reaclib', 'Kunz', or 'CF88'
         !! note: original CF88 rate is actually multiplied by 1.7 as in Timmes' rates
      !set_rate_n14pg = 'jina reaclib' ! empty string (default) means ignore this control
         !! one of 'NACRE', 'jina reaclib', or 'CF88'
      !set_rate_3a = 'jina reaclib' ! empty string (default) means ignore this control
         !! one of 'NACRE', 'jina reaclib', or 'CF88


      !## modified history.data and profile.data
      history_columns_file = 'history_columns.list'
      profile_columns_file = 'profile_columns.list'

      ! change_v_flag = .true.
      ! new_v_flag = .true.

      !pgstar_flag = .true.

/ ! end of star_job namelist


&controls

      ! steps_before_use_gold_tolerances = 1000
      !use_gold_tolerances = .false.

      extra_terminal_output_file = 'log'

      initial_mass = 100.0
      initial_z = 0.014d0
      !initial_y = 0.986d0
      !initial_he3 = 0d0

      ! atmosphere
      !which_atm_option = 'simple_photosphere'  
      Pextra_factor = 2 ! Pressure at optical depth tau is calculated as `P = tau*g/kap*(1 + Pextra)`

      ! convection
      mixing_length_alpha = 1.67d0
      !MLT_option = 'Henyey'

      !use_simple_es_for_kap= .true.
      use_Type2_opacities = .true.
      Zbase=0.014d0
      okay_to_reduce_gradT_excess = .true.
      !opacity_max = 0.5


      !## controls output
      photo_interval = 100 
      profile_interval = 50 
      history_interval = 1
      terminal_interval = 10
      max_num_profile_models = -1

      !## definition of core boundaries

      ! he core
      he_core_boundary_h1_fraction = 0.01
      ! c core
      c_core_boundary_he4_fraction = 0.01
      ! o core
      o_core_boundary_c12_fraction = 0.01
      ! si core
      si_core_boundary_o16_fraction = 0.01
      ! iron core
      fe_core_boundary_si28_fraction = 0.01

      ! wind mass loss
      ! 'Dutch' scheme combines several Dutch wind schemes: For O-stars the 'Vink' scheme and for
      ! wolf-rayet stars the one from Nugis & Lamers (2000)
      hot_wind_scheme = 'Dutch'
      cool_wind_AGB_scheme = 'Dutch'
      cool_wind_RGB_scheme = 'Dutch'
      Dutch_wind_lowT_scheme = 'de Jager' !'van Loon' gives a too strong mass loss! for Teff<1e4
      Dutch_scaling_factor = 0.49 ! as in GENEC (Maeder & Meynet (2001))

      !max_wind=1.e-99



      !## when to stop

      ! hydrogen depletion
      ! xa_central_lower_limit_species(1) = 'h1'
      ! xa_central_lower_limit(1) = 0.00001

      ! helium depletion
      ! xa_central_lower_limit_species(2) = 'he4'
      ! xa_central_lower_limit(2) = 0.001

      ! oxygen depletion
       xa_central_lower_limit_species(3) = 'o16'
       xa_central_lower_limit(3) = 0.00001
      
      ! silicon depletion
      ! xa_central_lower_limit_species(4) = 'si28'
      ! xa_central_lower_limit(4) = 0.01

      ! collapsing core
      ! fe_core_infall_limit = 1d8 ![cm/sec]

      fe56ec_fake_factor = 1d-4
      !rate_factor_for_fe56ec_to_cr56 does not exist in MESA 8118 anymore



      !## overshoot
      overshoot_f0_above_nonburn_core = 0.004
      overshoot_f_above_nonburn_core = 0.004
      overshoot_f0_above_nonburn_shell = 0.004
      overshoot_f_above_nonburn_shell = 0.004
      overshoot_f0_below_nonburn_shell = 0.004
      overshoot_f_below_nonburn_shell = 0.004     

      overshoot_f0_above_burn_h_core = 0.004
      overshoot_f_above_burn_h_core = 0.004
      overshoot_f0_above_burn_h_shell = 0.004
      overshoot_f_above_burn_h_shell = 0.004
      overshoot_f0_below_burn_h_shell = 0.004
      overshoot_f_below_burn_h_shell = 0.004

      overshoot_f0_above_burn_he_core = 0.004
      overshoot_f_above_burn_he_core = 0.004
      overshoot_f0_above_burn_he_shell = 0.004
      overshoot_f_above_burn_he_shell = 0.004
      overshoot_f0_below_burn_he_shell = 0.004
      overshoot_f_below_burn_he_shell = 0.004

      overshoot_f0_above_burn_z_core = 0.004
      overshoot_f_above_burn_z_core = 0.004
      overshoot_f0_above_burn_z_shell = 0.004
      overshoot_f_above_burn_z_shell = 0.004
      overshoot_f0_below_burn_z_shell = 0.004
      overshoot_f_below_burn_z_shell = 0.004

      !## resolution
      varcontrol_target = 0.3d-3 ! max. change between to steps
      dX_nuc_drop_limit = 0.1 ! limits max change in abundance mass fraction
      min_timestep_limit = 1d-14

/ ! end of controls namelist

&pgstar
         
!pause = .true.
         
Profile_Panels1_xmax = 2.2 ! -101d0 ! only used if /= -101d0
Abundance_xmax = 2.2 ! -101 ! only used if /= -101d0

Profile_Panels1_xaxis_name = 'logR_cm'

!Grid2_win_flag = .true.

! file output
!Grid2_file_flag = .true.
   Grid2_file_dir = 'png'
   Grid2_file_prefix = 'grid_'
   Grid2_file_interval = 5 ! 1 ! output when mod(model_number,Grid2_file_interval)==0
   Grid2_file_width = -1 ! negative means use same value as for window
   Grid2_file_aspect_ratio = -1 ! negative means use same value as for window
      

/ ! end of pgstar namelist
