! inlist_co_nova &star_job show_log_description_at_start = .false. !show_net_species_info = .true. !show_net_reactions_info = .true. !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! CHANGE FOR DIFFERENT MDOTS ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! load_saved_model = .true. saved_model_name = 'sewind_cool.mod' save_model_when_terminate = .true. save_model_filename = '1.38Msun_Tc_3e7_mdot_6.0e-7.mod' relax_mass_change = .FALSE. relax_mass_change_min_steps = 20 relax_mass_change_max_yrs_dt = 0.01 relax_mass_change_init_mdot = 5.2d-7 relax_mass_change_final_mdot = 6.0d-7 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! change_net = .true. !new_net_name = 'cno_extras_o18_to_mg26_plus_fe56.net' !new_net_name = 'arnett_to_mg27.net' 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 = .true. set_initial_age = .true. initial_age = 0 ! in years set_initial_model_number = .true. initial_model_number = 0 set_tau_factor = .true. ! set_to_this_tau_factor = 300 ! okay set_to_this_tau_factor = 30 ! relax_tau_factor = .true. ! relax_to_this_tau_factor = 300 ! to get through insane rlo situation ! dlogtau_factor = 0.1 set_initial_dt = .true. years_for_initial_dt = 1d-3 ! pgstar_flag = .true. / ! end of star_job namelist &controls L_div_CpTMdot_limit = -1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! CHANGE FOR DIFFERENT MASSES ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! initial_mass = 1.376 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! initial_z = 0.02d0 mixing_length_alpha = 1.73 mass_change_full_on_dt = 1d-6 ! (seconds) mass_change_full_off_dt = 1d-6 ! (seconds) !rlo_wind_eta = 1 !rlo_wind_min_L = 1d-6 ! only on when L > this limit. (Lsun) !rlo_wind_max_Teff = 1d99 ! only on when Teff < this limit. !rlo_wind_roche_lobe_radius = 100 ! 0.8 ! only on when R > this (Rsun) !rlo_wind_base_mdot = 1d-8 ! base rate of mass loss when R = roche lobe radius (Msun/year) !rlo_wind_scale_height = 1d-2 ! determines exponential growth rate of mass loss (Rsun) !!roche_lobe_xfer_full_on = 0.5d0 ! full accretion when R/RL <= this !!roche_lobe_xfer_full_off = 1.0d0 ! no accretion when R/RL >= this ! NOTE: for super eddington wind, ! we use Ledd averaged over mass to optical depth tau = 100 super_eddington_wind_eta = 1 ! parameter for mass loss driven by super Eddington luminosity !super_eddington_wind_Ledd_factor = 1 ! multiply Ledd by this factor when computing super Eddington wind ! e.g., if this is 2, then only get wind when L > 2*Ledd photosphere_r_upper_limit = 1d2 use_Ledoux_criterion = .false. !alpha_semiconvection = 1d-3 !thermo_haline_coeff = 1d3 logQ_limit = 100 gamma_center_limit = 10000 ! FOR THIS TEST ! star_mass_max_limit = 1.5d0 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! CHANGE FOR DIFFERENT MDOTS ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! mass_change = 6.0d-7 !LOG_directory = 'LOGS_1.20Msun_5.0e-7' !!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !max_timestep = 6.1556d7 ! in seconds. max_timestep <= 0 means no upper limit. min_center_cell_dq = 1d-7 max_center_cell_dq = 1d-6 max_surface_cell_dq = 1d-12 ! 1d-8 photostep = 500 profile_interval = 100 log_cnt = 1 terminal_cnt = 1 write_header_frequency = 10 ! max_num_profile_models = 1500 ! To get details as luminosity is ramping down !report_hydro_solver_progress = .true. ! set true to see info about newton iterations !report_ierr = .true. ! if true, produce terminal output when have some internal error report_why_dt_limits = .false. report_all_dt_limits = .false. write_varcontrol_internals = -1 report_hydro_dt_info = .false. !show_mesh_changes = .true. !mesh_dump_call_number = 79 !trace_evolve = .true. !hydro debugging !hydro_numerical_jacobian = .true. !hydro_check_everything = .true. !hydro_show_rcond = .true. !hydro_inspectB_flag = .true. !hydro_save_numjac_plot_data = .true. !hydro_dump_call_number = 10 min_timestep_limit = 1d-12 !RGB_wind_scheme = '' !AGB_wind_scheme = '' ! 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 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 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_limit = 1 ! limit for magnitude of max change in log10 density at any point delta_lgRho_hard_limit = -1 !relax_hard_limits_after_backup = .false. relax_hard_limits_after_retry = .false. ! limit for magnitude of change in log10(L/Lsun) delta_lgL_limit = 0.05 delta_lgL_hard_limit = 0.15 net_logTcut_lo = 5.3d0 net_logTcut_lim = 5.4d0 accrete_same_as_surface = .false. !accretion_h1 = 0.7000 ! Estimate Nomoto 2007 !accretion_h2 = 0 !accretion_he3 = 0.292910D-04 !accretion_he4 = 0.279971 ! Estimated Nomoto 2007 !accretion_zfracs = 4 ! Lodders 03 accretion_h1 = 0.72901442685D+00 !! Z=0.01 accretion_h2 = 1.4143154258D-05 accretion_he3 = 4.3313428237D-05 accretion_he4 = 0.26092811657D+00 accretion_zfracs = 2 ! GN93; Lodders03=4 !accretion_h1 = 0.000D+00 !! H-ashes !accretion_h2 = 1.4143154258D-05 !accretion_he3 = 4.3313428237D-05 !accretion_he4 = 9.9034D-01 !accretion_zfracs = 2 ! GN93; Lodders03=4 !which_atm_option = 'grey_and_kap' which_atm_option = 'simple_photosphere' !use_artificial_viscosity = .true. !l1_coef = 0.1 ! increase to as much as 1d3 to suppress pulsations !l1_coef = 1000 ! increase to as much as 1d3 to suppress pulsations ! limit for magnitude of max change in log10 temperature at any point ! delta_lgT_limit = 0.05 ! delta_lgT_hard_limit = 1 ! limit for magnitude of max change in log10 temperature at photosphere ! delta_lgTeff_limit = 0.01 ! delta_lgTeff_hard_limit = 0.03 ! change controls during mass ejection by super Eddington wind delta_lgT_limit = 0.5 delta_lgT_hard_limit = 1 delta_lgTeff_limit = 0.10 delta_lgTeff_hard_limit = 0.30 ! FOR DEBUGGING !report_hydro_solver_progress = .true. ! set true to see info about newton iterations !report_ierr = .true. ! if true, produce terminal output when have some internal error !max_years_for_timestep = 3.67628942044319d-05 !report_why_dt_limits = .true. !report_all_dt_limits = .true. !show_mesh_changes = .true. !mesh_dump_call_number = 5189 !okay_to_remesh = .false. !trace_evolve = .true. ! 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 ! hydro debugging !hydro_check_everything = .true. !hydro_inspectB_flag = .true. !hydro_numerical_jacobian = .true. !hydro_save_numjac_plot_data = .true. !small_mtx_decsol = 'lapack' !large_mtx_decsol = 'lapack' !hydro_dump_call_number = 195 varcontrol_target = 2.5d-4 ! 1d-3 !!!!!!!!!!!!!!!!!!!!!!!!!!!!! ! CHANGE FOR DIFFERENT MDOT ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!! !max_age = 300 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 !!! max_model_number = 50001 mesh_delta_coeff = 0.4 !0.25 ! 0.7 ! 1.5 ! 2.0 max_allowed_nz = 20000 ! 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 = .true. ! if false, use aspect ratio > 1 MAIN_win_width = 8 MAIN_win_aspect_ratio = 1.2 ! aspect_ratio = height/width xaxis_by = 'by_log_column_depth' show_main_win_text_info = .true. ! xaxis limits -- to override system default selections main_xmin = 0 !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 = -4.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 = 20 ! 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 = -12 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 = .false. 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 = 'logxm' Profile_xaxis_reversed = .true. Profile_xmin = -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 = 'logL' 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_mass' ! 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 = .true. Abundance_win_width = 7.5 Abundance_win_aspect_ratio = 0.618 ! aspect_ratio = height/width Abundance_xaxis_by = 'by_log_column_depth' ! select xaxis ! power xaxis limits -- to override system default selections Abundance_xmin = 0 ! only used if > -100 Abundance_xmax = -101 ! only used if > -100 ! power yaxis limits -- to override system default selections Abundance_log_mass_frac_min = -6 ! only used if < 0 Abundance_log_mass_frac_max = 0.4 ! 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 ! Abund_Power window -- Abundance and Power together !Abund_Power_win_flag = .true. Abund_Power_win_width = 6 Abund_Power_win_aspect_ratio = 1.2 ! aspect_ratio = height/width Abund_Power_xaxis_by = 'by_log_column_depth' ! select xaxis (same options as for main window) ! not used if == '' Abund_Power_xmin = 0 ! only used if > -100 Abund_Power_xmax = -101 ! only used if > -100 ! file output Abund_Power_file_flag = .false. Abund_Power_file_dir = 'png' Abund_Power_file_prefix = 'abund_power_dmix' Abund_Power_file_cnt = 5 ! output when mod(model_number,Abund_Power_file_cnt)==0 Abund_Power_file_width = -1 ! negative means use same value as for window Abund_Power_file_aspect_ratio = -1 ! negative means use same value as for window ! Abund_Pwr_Dmix window -- Abundance, Power, and Mixing Ds Abund_Pwr_Dmix_win_flag = .false. Abund_Pwr_Dmix_win_width = 6 Abund_Pwr_Dmix_win_aspect_ratio = 1.2 ! aspect_ratio = height/width Abund_Pwr_Dmix_xaxis_by = '' ! select xaxis (same options as for main window) ! not used if == '' Abund_Pwr_Dmix_xmin = -101 ! only used if > -100 Abund_Pwr_Dmix_xmax = -101 ! only used if > -100 ! file output Abund_Pwr_Dmix_file_flag = .false. Abund_Pwr_Dmix_file_dir = 'png' Abund_Pwr_Dmix_file_prefix = 'abund_power_dmix' Abund_Pwr_Dmix_file_cnt = 5 ! output when mod(model_number,Abund_Pwr_Dmix_file_cnt)==0 Abund_Pwr_Dmix_file_width = -1 ! negative means use same value as for window Abund_Pwr_Dmix_file_aspect_ratio = -1 ! negative means use same value as for window ! Abund_Pwr_Dmix_Dyn window -- Abundance, Power, Mixing Ds, and Dynamo Abund_Pwr_Dmix_Dyn_win_flag = .false. Abund_Pwr_Dmix_Dyn_win_width = 6 Abund_Pwr_Dmix_Dyn_win_aspect_ratio = 1.2 ! aspect_ratio = height/width Abund_Pwr_Dmix_Dyn_xaxis_by = '' ! select xaxis (same options as for main window) ! not used if == '' Abund_Pwr_Dmix_Dyn_xmin = -101 ! only used if > -100 Abund_Pwr_Dmix_Dyn_xmax = -101 ! only used if > -100 ! file output Abund_Pwr_Dmix_Dyn_file_flag = .false. Abund_Pwr_Dmix_Dyn_file_dir = 'png' Abund_Pwr_Dmix_Dyn_file_prefix = 'hr_trho' Abund_Pwr_Dmix_Dyn_file_cnt = 5 ! output when mod(model_number,Abund_Pwr_Dmix_Dyn_file_cnt)==0 Abund_Pwr_Dmix_Dyn_file_width = -1 ! negative means use same value as for window Abund_Pwr_Dmix_Dyn_file_aspect_ratio = -1 ! negative means use same value as for window / ! end of pgstar namelist