 
&star_job
 
   show_log_description_at_start = .false.
 
   load_saved_model = .true.
   saved_model_name = 'rse_m8.0.mod'
 
   save_model_when_terminate = .true.
   save_model_filename = 'final.mod'
 
   kappa_file_prefix = 'a09'
   kappa_lowT_prefix = 'lowT_fa05_a09p'
   kappa_CO_prefix = 'a09_co'
   initial_zfracs = 6
 
   num_special_rate_factors = 2
   reaction_for_special_factor(1) = 'r_he4_he4_he4_to_c12'
   reaction_for_special_factor(2) =  'r_c12_ag_o16'
   special_rate_factor(1) = 1.0
   special_rate_factor(2) = 1.0
 
   ! set rate preferences to nacre
        set_rates_preference =.true.
        new_rates_preference = 1 !NACRE RATES
 
   pgstar_flag = .true.
 
/ !end of star_job namelist
 
 
&controls
 
!opacity_factor = 2
!min_logT_for_opacity_factor_off = 0.001d0
!min_logT_for_opacity_factor_on = 0.002d0
!max_logT_for_opacity_factor_on = 5.4d0
!max_logT_for_opacity_factor_off = 5.41d0
 
 
 
   ! approximate Cepheid instability strip edges
 
   ! blue edge for Z=0.02 from Bono et al. 2000 Table 5 page 315
   x_ctrl(2) = 3.6845 ! logT1
   x_ctrl(3) = 4.5 ! logL1
   x_ctrl(4) = 3.8126 ! logT2
   x_ctrl(5) = 2.4 ! logL2
 
   ! red edge for Z=0.02 from Bono et al. 2000 Table 5 page 315
   x_ctrl(6) = 3.5670 ! logT3
   x_ctrl(7) = 4.5 ! logL3
   x_ctrl(8) = 3.8274 ! logT4
   x_ctrl(9) = 2.4 ! logL4
 
   initial_mass = 8.0 ! in Msun units
   ! initial_Y = 0.256d0
   ! initial_Z = 0.008d0
   initial_Y = 0.28d0
   initial_Z = 0.02d0
 
   use_gold_tolerances = .true.
   use_eosDT2 = .true.
   use_eosELM = .false.
 
   use_dedt_form_of_energy_eqn = .true.
   min_cell_energy_fraction_for_dedt_form = 1d-9
   use_eps_correction_for_KE_plus_PE_in_dLdm_eqn = .true.
 
   use_Type2_opacities = .true.
   ! Zbase = 0.008d0
   Zbase = 0.02d0
 
   mixing_length_alpha = 1.8d0
   MLT_option = 'Henyey'
 
   smooth_convective_bdy = .true.
   use_Ledoux_criterion = .true.
   alpha_semiconvection = 0.1
 
   overshoot_f_above_burn_h_core =  0.012
   overshoot_f0_above_burn_h_core = 0.002
 
   ! overshoot_f_above_burn_he_core =  0.012
   ! overshoot_f0_above_burn_he_core = 0.002
 
   overshoot_f_below_nonburn_shell =  0.022
   overshoot_f0_below_nonburn_shell = 0.002
 
   predictive_mix(1) = .true.
   predictive_zone_type(1) = 'any'
   predictive_zone_loc(1) = 'core'
   predictive_bdy_loc(1) = 'any'
 
   predictive_mix(2) = .true.
   predictive_zone_type(2) = 'any'
   predictive_zone_loc(2) = 'surf'
   predictive_bdy_loc(2) = 'any'
 
 
! for core He-burning, to prevent splitting of the core convection zone and/or core breathing pulses
 
   predictive_superad_thresh(1) = 0.01
   predictive_avoid_reversal(1) = 'he4'
 
   varcontrol_target = 5.0d-5
 
   delta_HR_limit = 0.002
 
! stop when the center mass fraction of he4 drops below this limit
    xa_central_lower_limit_species(1) = 'he4'
    xa_central_lower_limit(1) = 1d-4
 
 
! output controls
 
   num_trace_history_values = 2
   trace_history_value_name(1) = 'rel_E_err'
   trace_history_value_name(2) = 'log_rel_run_E_err'
 
   photo_interval = 50
   profile_interval = 50
   history_interval = 1
   terminal_interval = 10
 
   write_profiles_flag = .false.
   profile_header_include_sys_details  = .false. ! useful for py_mesa_reader
 
! use custom opacity factor
 
   use_other_opacity_factor = .false.
 
/ ! end of controls namelist
 
&pgstar
 
   file_device = 'ps'
   file_extension = 'ps'
 
   ! see convective vs radiative regions
   ! Summary_Burn_win_flag = .true.
 
 
   ! HR window -- lg_L vs. lg_Teff
 
   HR_win_flag = .true.
 
   HR_win_width = 6
   HR_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
 
   ! axis limits -- to override system default selections
   ! HR_logT_min = 3.5!3.6 ! -101 ! only used if > -100
   ! HR_logT_max = 3.95 ! only used if > -100
   ! HR_logL_min = 2.7!2.99 ! -101 ! only used if > -100
   ! HR_logL_max = 3.6!3.4 ! -101 ! only used if > -100
 
   show_HR_classical_instability_strip = .false.
   show_HR_target_box = .false.
   HR_target_n_sigma = -3
   HR_target_logL = 0.67d0
   HR_target_logL_sigma = 0.05d0
   HR_target_logT = 3.785d0
   HR_target_logT_sigma = 0.00461d0
 
   ! show mass in the HR diagram
   ! implemented in run_star_extras
   HR_use_decorator = .true.
 
   ! file output
   HR_file_flag = .false.
   HR_file_dir = 'pgstar_out'
   HR_file_prefix = 'hr'
   HR_file_interval = 5 ! output when mod(model_number,HR_file_interval)==0
   HR_file_width = -1 ! negative means use same value as for window
   HR_file_aspect_ratio = -1 ! negative means use same value as for window
 
 
 
 
   Profile_Panels1_win_flag = .true.
 
   Profile_Panels1_win_width = 6
   Profile_Panels1_win_aspect_ratio = 0.75 ! aspect_ratio = height/width
 
   Profile_Panels1_xleft = 0.15
   Profile_Panels1_xright = 0.85
   Profile_Panels1_ybot = 0.15
   Profile_Panels1_ytop = 0.85
   Profile_Panels1_txt_scale = 1.0
   Profile_Panels1_title = 'Profile Panels'
 
   Profile_Panels1_xaxis_name = 'logT'
   Profile_Panels1_xaxis_reversed = .true.
   Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
   Profile_Panels1_xmax = 6.3d0 ! only used if /= -101d0
   Profile_Panels1_xmargin = 0d0
   Profile_Panels1_show_mix_regions_on_xaxis = .true.
 
   Profile_Panels1_yaxis_name(:) = ''
   Profile_Panels1_yaxis_reversed(:) = .false.
   Profile_Panels1_yaxis_log(:) = .false. ! show log10 of abs value
   Profile_Panels1_ymin(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_ymax(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_ycenter(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_ymargin(:) = 0.1
   Profile_Panels1_dymin(:) = -1
 
   Profile_Panels1_other_yaxis_name(:) = ''
   Profile_Panels1_other_yaxis_reversed(:) = .false.
   Profile_Panels1_other_yaxis_log(:) = .false. ! show log10 of abs value
   Profile_Panels1_other_ymin(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_other_ymax(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_other_ycenter(:) = -101d0 ! only used if /= -101d0
   Profile_Panels1_other_ymargin(:) = 0.1
   Profile_Panels1_other_dymin(:) = -1
 
   Profile_Panels1_show_grid = .false.
 
   ! setup default plot
   Profile_Panels1_num_panels = 2
   Profile_Panels1_yaxis_name(1) = 'opacity'
   Profile_Panels1_other_yaxis_name(1) = 'extra_opacity_factor'
   Profile_Panels1_yaxis_name(2) = 'logRho'
   Profile_Panels1_other_yaxis_name(2) = 'logP'
 
/ ! end of pgstar namelist
 
 
 
