[mesa-users] astero...
Frank Timmes
fxt44 at mac.com
Tue May 16 10:42:11 EDT 2017
thanks.
see warricks's message:
> you've added the astero stuff to &controls in the usual MESA inlist, rather than to &astero_search_controls.
fxt
> On May 16, 2017, at 7:04 AM, Bruno Lustosa <bruno.moura.pesquisa at gmail.com> wrote:
>
> mesa-r8845
>
>
> ##########################################################################
>
> &star_job
>
> pause_before_terminate = .true.
> create_pre_main_sequence_model = .true.
> pre_ms_relax_num_steps = 50
> pre_ms_T_c = 9e5
>
>
>
> !change_initial_net = .true.
> !new_net_name = 'agb.net'
> !set_abundance = .true.
> set_initial_abundance = .true.
> chem_name = 'li7'
> new_frac = 7.77594697D-09
> !new_frac = 1.06927581475e-08
> !set_abundance_nzlo = -1
> !set_abundance_nzhi = -1
>
>
> change_net=.true.
> eos_file_prefix = 'mesa'
> new_net_name = 'pp_extras.net'
> kap_cache_dir = '/home/bruno/mesa-r8845/data/kap_data/cache'
> kappa_file_prefix = 'OP_gs98'
> dump_missing_metals_into_heaviest = .false.
> set_uniform_initial_composition = .true.
> initial_h1 = 0.6997725177
> !initial y = 0.2797816935
> initial_he3 = 0.0000279782
> initial_he4 = 0.2797537154
> initial_h2 = 0
>
> set_rates_preference = .true.
> new_rates_preference = 1
>
>
> new_rotation_flag = .true.
> change_rotation_flag = .false.
> change_initial_rotation_flag = .true.
>
>
> initial_zfracs = 7
> new_surface_rotation_v = 2 ! solar (km sec^1)
> set_near_zams_surface_rotation_v_steps = 10
>
>
>
> pgstar_flag = .true.
>
>
> / ! end of star_job namelist
>
> &controls
> !max_age = 12000000000
> stop_at_TP = .true.
> ! check for retries and backups as part of test_suite
> ! you can/should delete this for use outside of test_suite
> max_number_backups = 200
> max_number_retries = 500
>
>
> max_model_number = 6000 !958 max
>
>
> initial_mass = 1.08
> initial_z = 0.0204457888
> !initial_z = 0.02
> initial_y = 0.2797816935
> !use_Type2_opacities = .true.
> Zbase = 0.0204457888
>
> am_nu_visc_factor = 0
> am_D_mix_factor = 0.0333333333333333d0
> D_DSI_factor = 0
> D_SH_factor = 1
> D_SSI_factor = 1
> D_ES_factor = 1
> D_GSF_factor = 1
> D_ST_factor = 1
>
> varcontrol_target = 1d-3
> mesh_delta_coeff = 1.5
>
>
> photo_interval = 10
> profile_interval = 50
> history_interval = 1
> terminal_interval = 10
> write_header_frequency = 10
>
> write_pulse_info_with_profile = .true. !.true. ! if true, write pulse info file when write profile
> pulse_info_format = 'FGONG'
>
> log_L_lower_limit = -1
> mixing_length_alpha = 1.9107649497
> MLT_option = 'Cox'
> ! options are:
> ! 'Cox' MLT as developed in Cox & Giuli 1968, Chapter 14.
> ! 'ML1' Bohm-Vitense 1958
> ! 'ML2' Bohm and Cassinelli 1971
> ! 'Mihalas' Mihalas 1978, Kurucz 1979
> ! 'Henyey' Henyey, Vardya, and Bodenheimer 1965
> ! 'Cox' option assumes optically thick material.
> ! the other options are various ways of extending to include optically thin material.
>
> !which_atm_option = 'simple_photosphere'
> smooth_convective_bdy = .true.
> convective_bdy_weight = 1
> D_visc_factor = 0
>
> cool_wind_full_on_T = 9d6
> cool_wind_full_off_T = 1d7
> cool_wind_RGB_scheme = 'Reimers'
> cool_wind_AGB_scheme = 'Blocker'
> RGB_to_AGB_wind_switch = 1d-4
> Reimers_scaling_factor = 0.5d0
> Blocker_scaling_factor = 0.5d0
>
>
> ! 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
> !hydro_show_correction_info = .true.
>
> !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.
>
>
> ! 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
>
> !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
> !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
>
> !trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
> !trace_newton_bcyclic_stephi = 1 ! last model number to trace
>
> !trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
> !trace_newton_bcyclic_iterhi = 2 ! last newton iter to trace
>
> !trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
> !trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
>
> !trace_newton_bcyclic_jlo = 1 ! 1st var to trace
> !trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi
>
> !trace_k = 0
>
> ! misc
> add_atmosphere_to_pulse_info = .true.
> add_center_point_to_pulse_info = .true. ! add point for r=0
> stop_at_TP = .true.
>
>
>
> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
> ! seismic data (frequencies in microHz)
>
> !nu_max = 1010
> ! nu_max is needed when
> ! chi2_seismo_nu_max_fraction > 0 or correction_factor > 0 (see below)
> !nu_max_sigma = -1
> ! nu_max_sigma is needed when chi2_seismo_nu_max_fraction > 0
>
> !delta_nu = 56.28
> !delta_nu_sigma = 1.0
>
> !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
> / ! end of controls namelist
>
>
>
> &pgstar
>
> Power_win_flag = .false.
> Power_win_width = 9
> HR_win_flag = .true.
> HR_win_width = 10
> HR_win_aspect_ratio = 1
> logg_logT_win_flag = .false.
> logg_logT_win_width = 8
> logg_logT_win_aspect_ratio = 1.0
> Abundance_win_flag = .false.
> Abundance_win_width = 6
> Abundance_win_aspect_ratio = 1
> HR_file_flag = .true.
> HR_file_dir = 'png'
> HR_file_prefix = 'HR_'
> HR_title = 'DHR - Bruno'
>
> ! Mode Propagation window
>
> Mode_Prop_win_flag = .false.
>
> Mode_Prop_win_width = 10
> Mode_Prop_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
>
> Mode_Prop_nu_max_obs = 1000
>
> Mode_Prop_xaxis_name = 'r_div_R'
>
> ! power xaxis limits -- to override system default selections
> Mode_Prop_xmin = -101 ! only used if > -100
> Mode_Prop_xmax = 0.3 !-101 ! only used if > -100
>
> ! power yaxis limits -- to override system default selections
> Mode_Prop_ymin = 2 ! -101 ! only used if > -100
> Mode_Prop_ymax = 5 ! -101 ! only used if > -100
>
> ! file output
> Mode_Prop_file_flag = .false.
> Mode_Prop_file_dir = 'pgstar_out'
> Mode_Prop_file_prefix = 'mode_prop'
> Mode_Prop_file_interval = 5 ! output when mod(model_number,Mode_Prop_file_interval)==0
> Mode_Prop_file_width = -1 ! negative means use same value as for window
> Mode_Prop_file_aspect_ratio = -1 ! negative means use same value as for window
>
>
> / ! end of pgstar namelist
> <inlist_2.0>
More information about the Mesa-users
mailing list