&star_job mesa_data_dir = '/home/swj/mesa/mesa_code/mesa3709/data' create_pre_main_sequence_model = .true. !save_model_number = !save_model_filename = '' !load_saved_model = .true. !saved_model_name = '' !saved_model_name = '' change_net = .true. ! new_net_name = 'eccc_apecxnet_nsource_URCA.net' ! new_net_name = 'FULL_03.net' new_net_name = 'approx21.net' !set_uniform_xa_from_file = .true. !file_for_uniform_xa = '' set_rates_preference = .true. ! for use by net + rates modules new_rates_preference = 2 ! 1 = NACRE rates -- this is the default ! 2 = jina reaclib rates -- to match jina where possible set_rate_c12ag = 'Kunz' ! empty string means ignore this control ! one of 'NACRE', 'Buchmann', 'Kunz', or 'CF88' ! note: original CF88 rate is actually multiplied by 1.7 as in Timmes' rates set_rate_n14pg = 'Imbriani' ! empty string means ignore this control ! one of 'NACRE', 'Imbriani', or 'CF88' set_rate_3a = 'Fynbo' ! empty string means ignore this control ! one of 'NACRE', 'Fynbo', or 'CF88' use_se_output = .true. log_columns_file = 'log_columns_SJ.list' set_uniform_initial_composition = .true. initial_h1 = 7.0644343473E-01 initial_h2 = 1.3705268516E-05 initial_he3 = 4.5399908474E-05 initial_he4 = 2.7349746009E-01 initial_zfracs = 2 change_v_flag = .true. new_v_flag = .false. / ! end of star_job namelist &controls initial_mass = 12.0 initial_z = 0.02 mixing_length_alpha = 1.6 operator_coupling_choice = 0 ! 0 -- fully coupled log_center_density_limit = 11. min_timestep_limit = 1d-12 ! (seconds) !report_why_dt_limits = .true. !report_all_dt_limits = .true. ! alternative correction strategy - for runs encountering several backups: tol_correction_norm_alt = 1d-2 tol_max_correction_alt = 2d0 max_number_backups = -1 ! ignore if < 0 max_number_retries = -1 ! ignore if < 0 max_backups_in_a_row = 100 ! if do more than this many without a successful step, then terminate the run. ! some logistics max_model_number = 400000 ! negative means no maximum ! lower limits on central mass fractions ! stop when central abundance drops below this limit ! can have up to num_xa_central_limits of these (see star_def for value) !xa_central_lower_limit_species(1) = 'he4' ! isotope name as defined in chem_def !xa_central_lower_limit(1) = 1d-5 !xa_central_lower_limit_species(1) = 'o16' ! isotope name as defined in chem_def !xa_central_lower_limit(1) = 1d-3 photostep = 50 profile_interval = 50 max_num_profile_models = 10000 ! maximum number of saved profiles log_cnt = 1 terminal_cnt = 10 write_header_frequency = 10 ! mesh max_allowed_nz = 10000 ! maximum number of grid points allowed ! resolve the C13 pockets - not for massive stars... !xa_function_species(1) = 'h1' ! name of nuclide as defined in chem_def ! xa_function_weight(1) = 5 ! xa_function_param(1) = 3d-4 !xa_function_species(2) = 'he4' ! name of nuclide as defined in chem_def ! xa_function_weight(2) = 5 ! xa_function_param(2) = 1d-2 !xa_function_species(3) = 'c13' ! name of nuclide as defined in chem_def ! xa_function_weight(3) = 4 ! xa_function_param(3) = 1d-3 !xa_function_species(4) = 'n14' ! name of nuclide as defined in chem_def ! xa_function_weight(4) = 3 ! xa_function_param(4) = 1d-2 mesh_dlog_pp_dlogP_extra = 0.15 mesh_dlog_cno_dlogP_extra = 0.15 mesh_dlog_3alf_dlogP_extra = 0.15 mesh_dlog_burn_c_dlogP_extra = 0.15 mesh_dlog_burn_n_dlogP_extra = 0.15 mesh_dlog_burn_o_dlogP_extra = 0.15 mesh_dlog_burn_ne_dlogP_extra = 0.08 ! more neon burning resolution mesh_dlog_burn_na_dlogP_extra = 0.15 mesh_dlog_burn_mg_dlogP_extra = 0.15 mesh_dlog_cc_dlogP_extra = 0.15 mesh_dlog_co_dlogP_extra = 0.15 mesh_dlog_oo_dlogP_extra = 0.15 mesh_dlog_burn_si_dlogP_extra = 0.15 mesh_dlog_burn_s_dlogP_extra = 0.15 mesh_dlog_burn_ar_dlogP_extra = 0.15 mesh_dlog_burn_ca_dlogP_extra = 0.15 mesh_dlog_burn_ti_dlogP_extra = 0.15 mesh_dlog_burn_cr_dlogP_extra = 0.15 mesh_dlog_burn_fe_dlogP_extra = 0.15 mesh_dlog_pnhe4_dlogP_extra = 0.15 mesh_dlog_other_dlogP_extra = 0.15 mesh_dlog_photo_dlogP_extra = 1 xa_function_species(1) = 'ne20' xa_function_weight(1) = 20 xa_function_param(1) = 1d-2 xa_function_species(2) = 'si30' xa_function_weight(2) = 20 xa_function_param(2) = 1d-2 xa_function_species(3) = 's34' xa_function_weight(3) = 20 xa_function_param(3) = 1d-2 xa_function_species(4) = 'cr52' xa_function_weight(4) = 20 xa_function_param(4) = 1d-2 xa_function_species(5) = 'fe56' xa_function_weight(5) = 20 xa_function_param(5) = 1d-2 xa_function_species(6) = 'si28' xa_function_weight(6) = 20 xa_function_param(6) = 1d-2 xa_function_species(7) = 's32' xa_function_weight(7) = 20 xa_function_param(7) = 1d-2 ! timesteps ! general varcontrol_target = 1d-4 ! this is the target value for relative variation in the structure from one model to the next. ! the default timestep adjustment is to increase or reduce the timestep depending on whether ! the actual variation was smaller or greater than this value. ! DISABLE THESE sum_xa_tolerance = 1d99 min_xa_hard_limit = -1d99 delta_lgRho_cntr_limit = 0.05 delta_lgRho_cntr_hard_limit = 0.1 delta_lgT_cntr_limit = 0.01 delta_lgT_cntr_hard_limit = 0.02 delta_lgTeff_limit = 0.01 delta_lgTeff_hard_limit = 0.02 delta_lgL_limit = 0.1 delta_lgL_hard_limit = 0.2 ! TP - again, NOT for massive stars! !delta_lgL_He_limit = 0.01 ! for TP resolution !delta_lgL_He_hard_limit = 0.05 !lgL_He_burn_min = 2.0 ! ignore changes in lgL_He if value ! is less than this ! DUP DUP_varcontrol_factor = 1.0 ! multiply varcontrol_target by this when TP_state == 2, i.e., during dredge up during thermal pulse ! time step controls for interpulse - NOT for massive stars.... !dH_limit_min_H = 1.e-2 !dH_limit = 0.1 !dH_div_H_limit_min_H = 1d-5 !dH_div_H_limit = 0.5 dX_limit_min_X = 5d-2 ! general burning timestep limiters for massive stars dX_limit = 2d-2 ! general burning timestep limiters for massive stars ! limits on max drop in abundance mass fraction from burning with possible mixing inflow ! this considers both nuclear rate and offsetting effect of mixing inflow !dX_nuc_drop_min_X_limit = 5d-2 ! SJONES !dX_nuc_drop_limit only for X > dX_nuc_drop_min_X_limit !dX_nuc_drop_max_A_limit = 52 ! SJONES !dX_nuc_drop_limit only for species with A <= dX_nuc_drop_max_A_limit !dX_nuc_drop_limit = 1.d-3 ! 1d-2 ! SJONES ! Overshooting mass_for_overshoot_full_on = 1.8 ! Msun units mass_for_overshoot_full_off = 0. ! Msun units overshoot_f_above_nonburn = 0 overshoot_f_below_nonburn = 0 overshoot_f_above_burn_h = 0 overshoot_f_below_burn_h = 0 overshoot_f_above_burn_he = 0 overshoot_f_below_burn_he = 0 overshoot_f_above_burn_z = 0 overshoot_f_below_burn_z = 0 ! C13 pocket overshoot_below_noburn_factor = 0.3 ! He-shell flash convection zone ovr_below_burn_he_factor = 0.143 ! atmosphere option which_atm_option = 'simple_photosphere' ! 'simple_photosphere' ! don't integrate, just estimate for tau=2/3 ! mass loss RGB_wind_scheme = 'Dutch' AGB_wind_scheme = 'Dutch' RGB_to_AGB_wind_switch = 1d-4 !Reimers_wind_eta = 0.5d0 !Blocker_wind_eta = 0.05d0 Dutch_wind_eta = 0.8 ! opacities use_CO_enhanced_opacities = .true. base_Z = 0.02 !base_fC = 0.17123643187d0 ! base value of XC/Z !base_fO = 0.4812175d0 ! base value of XO/Z / ! end of controls namelist