! inlist_swwind &star_job show_log_description_at_start = .false. load_saved_model = .true. saved_model_name = 'sewind.mod' change_net = .true. new_net_name = 'wd_accr.net' set_rates_preference = .true. ! for use by net + rates modules new_rates_preference = 1 ! 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. !change_v_flag = .true. !new_v_flag = .false. set_to_this_tau_factor = 300 ! tau_factor of 1 puts outer cell at photosphere. set_tau_factor = .true. ! change tau_factor without reconverging. !pgstar_flag = .true. / ! end of star_job namelist &controls ! check for retries and backups as part of test_suite ! you can/should delete this for use outside of test_suite max_number_backups = 32800 max_number_retries = 3280000 initial_z = 0.02d0 ! initial metallicity mixing_length_alpha = 1.73 ! NOTE we need to relax correction tolerances and remove checks for residuals tol_correction_norm = 1d-3 tol_max_correction = 1d-1 max_iter_for_resid_tol1 = 999 tol_residual_norm1 = 1d99 tol_max_residual1 = 1d99 gamma_center_limit = 10000 max_model_number = 54001 !max_age = 9.445d6 !!! M(WD_PROG)=3; M(DONOR)=2.5; SHOULD REACH 1.0085 Msun !!! !max_age = 8.021d6 !!! M(WD_PROG)=7; M(DONOR)=6; SHOULD REACH 1.2015 Msun !!! max_age = 8.021d5 !!! M(WD_PROG)=7; M(DONOR)=6; SHOULD REACH 1.2015 Msun !!! correction_xa_limit = 5d-3 L_div_CpTMdot_limit = -1 ! disable for this test min_center_cell_dq = 1d-7 max_center_cell_dq = 1d-6 max_surface_cell_dq = 1d-8 photostep = 500 profile_interval = 100 log_cnt = 1 terminal_cnt = 1 write_header_frequency = 10 min_timestep_limit = 1d-25 accrete_same_as_surface = .false. !accretion_h1 = 0.70644343473D+00 !! Z=0.02 !accretion_h2 = 1.3705268516D-05 !accretion_he3 = 4.5399908474D-05 !accretion_he4 = 0.27349746009D+00 accretion_h1 = 0.72901442685D+00 !! Z=0.01 accretion_h2 = 1.4143154258D-05 accretion_he3 = 4.3313428237D-05 accretion_he4 = 0.26092811657D+00 !accretion_h1 = 0.74932831975D+00 !! Z=0.001 !accretion_h2 = 1.4537251426D-05 !accretion_he3 = 4.1435596024D-05 !accretion_he4 = 0.24961570740D+00 !accretion_h1 = 0.75135970904D+00 !! Z=0.0001 !accretion_h2 = 1.457666114D-05 !accretion_he3 = 4.1247812802D-05 !accretion_he4 = 0.24848446649D+00 accretion_zfracs = 2 ! GN93; Lodders03=4 !xa_function_species(1) = 'h1' ! name of nuclide as defined in chem_def !xa_function_weight(1) = 20 !xa_function_param(1) = 1d-6 !xa_function_species(2) = 'he4' ! name of nuclide as defined in chem_def !xa_function_weight(2) = 20 !xa_function_param(2) = 1d-2 ! resolve the C13 pockets xa_function_species(1) = 'h1' ! name of nuclide as defined in chem_def xa_function_weight(1) = 20 xa_function_param(1) = 1d-6 xa_function_species(2) = 'he4' ! name of nuclide as defined in chem_def xa_function_weight(2) = 20 xa_function_param(2) = 1d-2 xa_function_species(3) = 'c13' ! name of nuclide as defined in chem_def xa_function_weight(3) = 5 xa_function_param(3) = 1d-4 xa_function_species(4) = 'n14' ! name of nuclide as defined in chem_def xa_function_weight(4) = 4 xa_function_param(4) = 1d-6 delta_lgL_H_limit = 0.05 ! limit for magnitude of change in lgL_H delta_lgL_H_hard_limit = 0.5 lgL_H_burn_min = 1.5 ! ignore changes in lgL_H if value is less than this delta_lgL_He_limit = 0.025 ! limit for magnitude of change in lgL_He delta_lgL_He_hard_limit = 1 lgL_He_burn_min = 2.5 ! ignore changes in lgL_He if value is less than this 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 relax_hard_limits_after_retry = .false. mesh_delta_coeff = 0.4 !mesh_delta_coeff = 0.2 max_allowed_nz = 8000 ! maximum number of grid points allowed varcontrol_target = 1d-2 super_eddington_wind_eta = 1. !super_eddington_wind_eta = 1.d-2 super_eddington_wind_Ledd_factor = 1. ! 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 ! Overshooting mass_for_overshoot_full_on = 1.8 ! Msun units mass_for_overshoot_full_off = 0.3 ! Msun units overshoot_f_above_nonburn = 0.014 overshoot_f_below_nonburn = 0.014 overshoot_f_above_burn_h = 0.014 overshoot_f_below_burn_h = 0.014 overshoot_f_above_burn_he = 0.011 overshoot_f_below_burn_he = 0.011 overshoot_f_above_burn_z = 0.014 overshoot_f_below_burn_z = 0.014 !overshoot_D2_above_nonburn = 100000000000 overshoot_D2_below_nonburn = 100000000000 !overshoot_D2_above_burn_h = 100000000000 overshoot_D2_below_burn_h = 100000000000 !overshoot_D2_above_burn_he = 100000 overshoot_D2_below_burn_he = 100000 !overshoot_D2_above_burn_z = -0.0001 overshoot_D2_below_burn_z = -0.0001 !overshoot_f2_above_nonburn = 0.25 overshoot_f2_below_nonburn = 0.25 !overshoot_f2_above_burn_h = 0.25 overshoot_f2_below_burn_h = 0.25 !overshoot_f2_above_burn_he = 0.14 overshoot_f2_below_burn_he = 0.14 !overshoot_f2_above_burn_z = -0.00001 overshoot_f2_below_burn_z = -0.00001 ! C13 pocket overshoot_below_noburn_factor = 1.0 net_logTcut_lo = 0 net_logTcut_lim = 0 !mass_change = 1.0d-7 ! rate of accretion (Msun/year) mass_change = 3.162277660168379d-07 ! rate of accretion (Msun/year) (=d-6.5) !mass_change = 3.407009466358378d-07 ! rate of accretion (Msun/year); tentativo caso 1.25 Msun !mass_change = 3.30d-07 ! rate of accretion (Msun/year); tentativo caso 1.25 Msun !mass_change = 3.651741272548377d-07 ! rate of accretion (Msun/year); tentativo caso 1.25 e 1.31 Msun !mass_change = 4.641588833612782d-07 ! rate of accretion (Msun/year); tentativo caso 1.31 Msun ! atmosphere option which_atm_option = 'simple_photosphere' ! mass loss !RGB_wind_scheme = 'Reimers' !AGB_wind_scheme = 'Blocker' !RGB_to_AGB_wind_switch = 1d-4 !Reimers_wind_eta = 0.5d0 !Blocker_wind_eta = 0.08d0 !!Blocker_wind_eta = 0.01d0 !super_eddington_wind_eta = 1 !supersonic_wind_eta = 1 ! opacities use_CO_enhanced_opacities = .true. base_Z = 0.01 base_fC = 0.17123643187d0 ! base value of XC/Z base_fO = 0.4812175d0 ! base value of XO/Z max_logT_for_special_kap = 5.7 special_kap_use_CO_enhanced = .true. special_kap_Z_base = 0.0001 !opacity_max = 110. !opacity_factor = 0.9 / ! end of controls namelist &pgstar ! top level controls !pause = .true. ! if true, the code waits for user to enter a RETURN on the command line ! main window MAIN_win_flag = .true. show_TRho_Profile_with_main = .true. show_HR_TRho_with_main = .false. ! if false, use aspect ratio > 1 MAIN_win_width = 7 MAIN_win_aspect_ratio = 1.2 ! aspect_ratio = height/width xaxis_by = 'by_logxm' show_main_win_text_info = .true. ! xaxis limits -- to override system default selections main_xmin = -10 !0.4 !main_xmax = 0.5 logxq_cutoff = -12.1 ! min value when using logxm for xaxis !main_logL_max = 7 !main_logL_min = 0 ! control for plot showing abundances log_mass_frac_ymax = 0.7 ! making this > 0 helps readability log_mass_frac_ymin = -7.5 num_abundance_line_labels = 5 log_abund_vary = 0.1 ! don't show species if it would vary by < this in plot ! file output !MAIN_file_flag = .true. MAIN_file_dir = 'png' MAIN_file_prefix = 'main' MAIN_file_cnt = 1 ! output when mod(model_number,main_file_cnt)==0 ! TRho Profile window -- current model in T-Rho plane TRho_Profile_win_flag = .false. TRho_Profile_win_width = 7.5 TRho_Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width show_TRho_Profile_text_info = .false. show_TRho_Profile_legend = .true. show_HR_TRho_with_TRho_Profile = .false. show_TRho_Profile_burn_labels = .true. TRho_Profile_show_logQ_limit = .true. !show_TRho_Profile_cross_hair = .true. show_TRho_Profile_mass_locs = .false. show_TRho_Profile_kap_regions = .false. show_TRho_Profile_eos_regions = .true. show_TRho_Profile_degeneracy_line = .true. show_TRho_Profile_Pgas_Prad_line = .true. show_TRho_Profile_burn_lines = .true. !show_TRho_Profile_annotation1 = .true. !show_TRho_Profile_annotation2 = .true. !show_TRho_Profile_annotation3 = .true. ! axis limits TRho_Profile_xmin = -10 TRho_Profile_xmax = 10 TRho_Profile_ymin = 3.0 TRho_Profile_ymax = 10 ! file output !TRho_Profile_file_flag = .true. TRho_Profile_file_dir = 'png' TRho_Profile_file_prefix = 'trho_profile' TRho_Profile_file_cnt = 2 ! output when mod(model_number,TRho_Profile_file_cnt)==0 ! Profile window !Profile_win_flag = .true. Profile_win_width = 7.5 Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width show_Profile_legend = .false. Profile_legend_coord = 0.65 Profile_legend_fjust = 0.0 Profile_legend_disp1 = -7.5 Profile_legend_del_disp = -1.5 Profile_show_mass_locs = .false. Profile_show_decorated_line = .false. show_HR_TRho_with_Profile = .false. show_Profile_text_info = .false. Profile_text_info_xfac = 0.78 ! controls x location Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text Profile_text_info_yfac = 0.95 ! controls y location of 1st line Profile_text_info_dyfac = -0.04 ! controls line spacing show_Profile_cross_hair = .false. show_Profile_annotation1 = .false. show_Profile_annotation2 = .false. show_Profile_annotation3 = .false. ! axis choices ! to get a list of axis identifiers, edit your star_job controls ! to set show_profile_column_numbers = .true. ! you can use any identifier for x or y Profile_xaxis_name = 'mass' !Profile_xaxis_reversed = .true. Profile_xmin = 0.9999 !-12 ! only used if > -100 !Profile_xmax = 1.2 !-101 ! only used if > -100 Profile_yaxis_name = 'entropy' !'grada' Profile_yaxis_reversed = .false. Profile_ymin = -111 !5.4 ! only used if > -100 Profile_ymax = -111 !8.6 ! only used if > -100 Profile_dymin = 0.025 Profile_other_yaxis_name = 'extra_heat' Profile_other_yaxis_reversed = .false. Profile_other_ymin = -101 ! only used if > -100 Profile_other_ymax = -101 ! only used if > -100 ! file output !Profile_file_flag = .true. Profile_file_dir = 'png' Profile_file_prefix = 'profile' Profile_file_cnt = 5 ! output when mod(model_number,Profile_file_cnt)==0 ! Convection window -- history of convection and more. CONV_win_flag = .false. CONV_win_width = 9 CONV_win_aspect_ratio = 0.62 ! aspect_ratio = height/width ! axis limits CONV_xmax = -1 ! step number. negative means use default. CONV_xmin = -1 ! step number. negative means use default. CONV_max_width = 0 ! only used if > 0. causes xmin to move with xmax. CONV_mmax = 1.00007 !-1 ! (Msun units) negative means use start initial mass CONV_mmin = 1 !0 ! (Msun units) CONV_lgLmax = -101 ! only used if > -100; (L in Lsun units) CONV_lgLmin = -2 ! only used if > -100; (L in Lsun units) CONV_show_burn = .true. CONV_show_mixing = .true. CONV_show_log_radius = .true. CONV_show_luminosities = .true. CONV_show_mass_boundaries = .true. show_CONV_annotation1 = .false. show_CONV_annotation2 = .false. show_CONV_annotation3 = .false. ! file output CONV_file_flag = .false. CONV_file_dir = 'pgstar_out' CONV_file_prefix = 'conv' CONV_file_cnt = 5 ! output when mod(model_number,CONV_file_cnt)==0 CONV_file_width = -1 ! negative means use same value as for window CONV_file_aspect_ratio = -1 ! negative means use same value as for window ! Surface History window Surf_Hist_win_flag = .false. Surf_Hist_win_width = 6.5 Surf_Hist_win_aspect_ratio = 1.2 ! aspect_ratio = height/width ! axis limits Surf_Hist_xmax = -1 ! step number. negative means use default. Surf_Hist_xmin = -1 ! step number. negative means use default. Surf_Hist_max_width = 0 ! only used if > 0. causes xmin to move with xmax. ! abundance limits Surf_Hist_mass_frac_min = 0 Surf_Hist_mass_frac_max = 1.1 ! > 0 helps legibility show_Surf_Hist_annotation1 = .false. show_Surf_Hist_annotation2 = .false. show_Surf_Hist_annotation3 = .false. ! file output Surf_Hist_file_flag = .false. Surf_Hist_file_dir = 'pgstar_out' Surf_Hist_file_prefix = 'conv' Surf_Hist_file_cnt = 5 ! output when mod(model_number,Surf_Hist_file_cnt)==0 Surf_Hist_file_width = -1 ! negative means use same value as for window Surf_Hist_file_aspect_ratio = -1 ! negative means use same value as for window ! Power window Power_win_flag = .false. Power_win_width = 8.5 Power_win_aspect_ratio = 0.618 ! aspect_ratio = height/width Power_xaxis_by = 'by_logxm' ! power xaxis limits -- to override system default selections Power_xmin = -101 ! only used if > -100 Power_xmax = -101 ! only used if > -100 ! power yaxis limits -- to override system default selections Power_ymin = -101 ! only used if > -100 Power_ymax = -101 ! only used if > -100 show_Power_cross_hair = .false. ! file output Power_file_flag = .false. Power_file_dir = 'pgstar_out' Power_file_prefix = 'power' Power_file_cnt = 5 ! output when mod(model_number,Power_file_cnt)==0 Power_file_width = -1 ! negative means use same value as for window Power_file_aspect_ratio = -1 ! negative means use same value as for window ! Abundance window Abundance_win_flag = .false. Abundance_win_width = 8.5 Abundance_win_aspect_ratio = 0.618 ! aspect_ratio = height/width Abundance_xaxis_by = 'by_logxm' ! power xaxis limits -- to override system default selections Abundance_xmin = -101 ! only used if > -100 Abundance_xmax = -101 ! only used if > -100 ! power yaxis limits -- to override system default selections Abundance_log_mass_frac_min = 1 ! only used if < 0 Abundance_log_mass_frac_max = 1 ! only used if < 0 show_Abundance_cross_hair = .false. ! file output Abundance_file_flag = .false. Abundance_file_dir = 'pgstar_out' Abundance_file_prefix = 'abund' Abundance_file_cnt = 5 ! output when mod(model_number,Abundance_file_cnt)==0 Abundance_file_width = -1 ! negative means use same value as for window Abundance_file_aspect_ratio = -1 ! negative means use same value as for window / ! end of pgstar namelist