[mesa-users] Fwd: Re: Core resolution: non-convergent behavior increasing mesh_delta_coeff
Mathieu
mathren90 at gmail.com
Mon Aug 24 13:20:01 EDT 2015
Hi everyone,
This is a followup email to describe a solution for improving the
convergence of massive stars structures (computed with mass loss)
differing only in the the spatial resolution (mesh_delta_coeff).
This solution has been found thanks to the help of Frank Timmes. You can
find below my original description of the problem.
To solve this issue, three steps seem to be needed:
* take smaller timesteps, by setting varcontrol_target=1d-4,
dX_nuc_drop_limit = 1d-4, and dX_nuc_drop_limit_at_high_T = 1d-4 (N.B:
for silicon burning you might need to loosen some of these, look at the
inlists of the MESA 3 paper
<http://mesastar.org/results/mesa-iii-instrument-paper-inlists>);
* set the maximum relative mass of cells in the model to max_dq=0.5d-4
in the inlist. This makes the baseline grid have *at least* 20000 zones;
* resolve properly the *outer* portion of the star, to make sure the
amount mass lost at each timestep does not depend on the
mesh_delta_coeff adopted. The way I implemented this is in the routine
routine other_mesh_fcn_data, in my run_star_extras.f (see attachment).
In particular using:
vals(k,2) = 100*weight *
min(s%m(1)/Msun,max(s%m(1)/Msun-0.5, s%m(k)/Msun))
where weight=10, I impose a minimum number of cells in the outermost
0.5Msun of the star, which allows to properly resolve the mass loss
rate, and consequently its effects on the location of the burning
shells, and dimensions of the various cores at oxygen depletion. This
was inspired by the R_function2* entries in the controls namelist.
The attached plots show the evolutionary tracks on the
(rho_c,T_c)-plane, and the abundance profiles at oxygen depletion for
two models computed respectively with mesh_delta_coeff=1.0 and 0.5, with
the attached inlists and run_star_extras.f. This are, so far, my best
proof of the convergence of these models.
I think this shows how important is to check the convergence when
changing the resolution of MESA, and that sometimes problems in very
different regions of the stars might be strictly connected: even if you
just care about the core structure, you might want to be sure to resolve
properly the outermost layers of the star!
If anyone has encountered similar problems, and could suggest better
and/or faster (less cells or larger timesteps) ways to obtain
"converged" solutions, I would be very happy to know. Any other
suggestion is very welcome too.
Best regards,
Mathieu
On 07/24/2015 12:27 PM, Mathieu wrote:
> Hi everyone,
>
> I am trying to run models of massive stars (15 <= M_zams/Msun <= 30)
> until oxygen depletion (X(o16)<0.04 in the center), but I am finding a
> problem with the spatial resolution (mesh gridding) of the inner
> regions, especially during and after neon burning. The problem is
> illustrated in the attached plots.
>
> I compare 30 solar mass models computed starting on the ZAMS with the
> mesa_45.net nuclear reaction network. All models have the same
> physical settings (see attached inlists and run_star_extras.f), except
> for my varying mesh_delta_coeff.*
>
> Increasing the mesh_delta_coeff produces qualitatively and
> quantitatively different tracks on the (central density,central
> temperature)-plane, and significantly different masses of the various
> (He, C/O) cores***(see attached plots). The number of cells at the
> end of each computation shown in the attached plots is, respectively:
>
> mesh_delta_coeff = 1 => nz = 31096
> mesh_delta_coeff = 0.5 => nz = 65763
> mesh_delta_coeff = 0.2 => nz = 142096
>
> I also ran tests without my "customized" mass loss, mesh refinements
> and timestep controls, and different initial masses, both with the
> MESA release 6794 and 7624, and I have found that the discrepancies
> were worse than or comparable to what I am showing here.
> Using my timestep controls (in the routine extras_finish_step) and
> mesh refinements (other_mesh_fcns_data) improves the behavior, but
> only up to neon burning.
>
> I have uploaded here
> <http://www.tapir.caltech.edu/%7Emathren90/materials/MWE.tar.gz>a
> minimal working example with models and photos to restart the runs
> with mesh_delta_coeff=1 and mesh_delta_coeff=0.5, ready to be compiled
> and run (with the release 7624).
>
> This non-convergent behavior when changing the resolution makes a
> detailed study of the core structure in late stages very questionable,
> and I would be grateful if you could help me constrain the mesh
> intelligently, so that the core structure does not depend so
> sensitively on the spatial resolution (provided that enough cells are
> used, of course).
>
> Thanks,
>
> Mathieu
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150824/429b8b50/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: outer_res+max_dq_ab_M.png
Type: image/png
Size: 108211 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150824/429b8b50/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: outer_res+max_dq_Tc_rhoc.png
Type: image/png
Size: 56260 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150824/429b8b50/attachment-0001.png>
-------------- next part --------------
! this is the master inlist that MESA reads when it starts.
! This file tells MESA to go look elsewhere for its configuration
! info. This makes changing between different inlists easier, by
! allowing you to easily change the name of the file that gets read.
&star_job
read_extra_star_job_inlist1 = .true.
extra_star_job_inlist1_name = 'inlist_massive_defaults'
read_extra_star_job_inlist2 = .true.
extra_star_job_inlist2_name = 'inlist_project'
/ ! end of star_job namelist
&controls
read_extra_controls_inlist1 = .true.
extra_controls_inlist1_name = 'inlist_massive_defaults'
read_extra_controls_inlist2 = .true.
extra_controls_inlist2_name = 'inlist_project'
/ ! end of controls namelist
&pgstar
read_extra_pgstar_inlist1 = .false.
extra_pgstar_inlist1_name = 'inlist_massive_defaults'
/ ! end of pgstar namelist
-------------- next part --------------
! inlist_massive_defaults
&star_job
show_log_description_at_start = .false.
!change_initial_net = .true.
!new_net_name = 'approx21_cr60.net'
kappa_file_prefix = 'gs98'
change_v_flag = .true.
new_v_flag = .true.
!new_rotation_flag = .true.
!change_rotation_flag = .false. ! rotation off until near zams
!new_omega_div_omega_crit = 0.5
!near_zams_relax_omega_div_omega_crit = .true.
!num_steps_to_relax_rotation = 50 ! use this many steps to change value
set_initial_age = .true.
initial_age = 0 ! in years
set_initial_model_number = .true.
initial_model_number = 0
!pgstar_flag = .true.
/ ! end of star_job namelist
&controls
!min_T_for_time_averaged_conv_velocity = 3e9
Pextra_factor = -1 ! keep the old form for now
!num_trace_history_values = 1
!trace_history_value_name(1) = 'split_mixing_choice'
!split_mixing_choice = -2
!split_mix_do_burn = .true.
! high center T limit to avoid negative mass fractions
sig_min_factor_for_high_Tcenter = 0.01
! inactive when >= 1d0
! if Tcenter >= Tcenter_min_for_sig_min_factor_full_on,
! then okay to reduce sig by as much as this factor
! as needed to prevent causing negative abundances
Tcenter_min_for_sig_min_factor_full_on = 3.2d9
! if Tcenter >= this, factor = sig_min_factor_for_neg_abundances
! this should be > Tcenter_max_for_sig_min_factor_full_off.
Tcenter_max_for_sig_min_factor_full_off = 2.8d9
! if Tcenter <= this, factor = 1, so has no effect
! this should be < Tcenter_min_for_sig_min_factor_full_on.
! for T > full_off and < full_on, factor changes linearly with Tcenter
velocity_logT_lower_bound = 8
max_dt_yrs_for_velocity_logT_lower_bound = 0.1
mesh_delta_coeff = 1.5
mesh_delta_coeff_for_highT = 2.5
logT_max_for_standard_mesh_delta_coeff = 9.0
logT_min_for_highT_mesh_delta_coeff = 9.5
varcontrol_target = 1d-3
dX_nuc_drop_limit = 1d-2
dX_nuc_drop_limit_at_high_T = 5d-3 ! for center logT > 9.45
delta_Ye_highT_limit = 1d-3
okay_to_reduce_gradT_excess = .true.
gradT_excess_f1 = 1d-4
gradT_excess_f2 = 1d-2
gradT_excess_age_fraction = 0.9d0
gradT_excess_lambda1 = 1.0
gradT_excess_beta1 = 0.35
gradT_excess_lambda2 = 0.5
gradT_excess_beta2 = 0.25
gradT_excess_dlambda = 0.1
gradT_excess_dbeta = 0.1
RGB_wind_scheme = 'Dutch'
AGB_wind_scheme = 'Dutch'
RGB_to_AGB_wind_switch = 1d-4
Dutch_wind_eta = 0.8
include_dmu_dt_in_eps_grav = .true.
use_Type2_opacities = .true.
Zbase = 0.019 !to match Z from Vink et al. 2001
xa_scale = 1d-5
newton_itermin = 2
mixing_length_alpha = 1.5
MLT_option = 'Henyey'
use_Ledoux_criterion = .true.
alpha_semiconvection = 0.01
thermohaline_coeff = 2
thermohaline_option = 'Kippenhahn'
overshoot_f0_above_nonburn_core = 0.0005
overshoot_f0_above_nonburn_shell = 0.0005
overshoot_f0_below_nonburn_shell = 0.0005
overshoot_f0_above_burn_h_core = 0.0005
overshoot_f0_above_burn_h_shell = 0.0005
overshoot_f0_below_burn_h_shell = 0.0005
overshoot_f0_above_burn_he_core = 0.0005
overshoot_f0_above_burn_he_shell = 0.0005
overshoot_f0_below_burn_he_shell = 0.0005
overshoot_f0_above_burn_z_core = 0.0005
overshoot_f0_above_burn_z_shell = 0.0005
overshoot_f0_below_burn_z_shell = 0.0005
overshoot_f_above_nonburn_core = 0.001
overshoot_f_above_nonburn_shell = 0.001
overshoot_f_below_nonburn_shell = 0.001
overshoot_f_above_burn_h_core = 0.001
overshoot_f_above_burn_h_shell = 0.001
overshoot_f_below_burn_h_shell = 0.001
overshoot_f_above_burn_he_core = 0.000
overshoot_f_above_burn_he_shell = 0.000
overshoot_f_below_burn_he_shell = 0.000
overshoot_f_above_burn_z_core = 0.000
overshoot_f_above_burn_z_shell = 0.000
overshoot_f_below_burn_z_shell = 0.000
!min_T_for_acceleration_limited_conv_velocity = 0
screening_mode = 'extended'
tol_residual_norm1 = 1d-5
tol_max_residual1 = 1d-2
iter_for_resid_tol2 = 3
tol_residual_norm2 = 1d99
tol_max_residual2 = 1d99
iter_for_resid_tol3 = 12
min_timestep_limit = 1d-12 ! (seconds)
delta_lgL_He_limit = 0.1 ! limit for magnitude of change in lgL_He
dX_nuc_drop_max_A_limit = 52
dX_nuc_drop_min_X_limit = 1d-4
dX_nuc_drop_hard_limit = 1d99
delta_lgTeff_limit = 0.5
delta_lgL_limit = 0.5
delta_lgRho_cntr_limit = 0.02
!delta_lgE_limit = 0.04
!delta_lgE_hard_limit = 0.06
! mesh_dlog_pp_dlogP_extra = 0.4
! mesh_dlog_cno_dlogP_extra = 0.4
!
! mesh_dlog_burn_n_dlogP_extra = 0.4
! mesh_dlog_3alf_dlogP_extra = 0.4
! mesh_dlog_burn_c_dlogP_extra = 0.4
! xa_function_species(1) = 'h1'
! xa_function_weight(1) = 30
! xa_function_param(1) = 1d-3
! xa_function_species(2) = 'he4'
! xa_function_weight(2) = 30
! xa_function_param(2) = 1d-3
T_mix_limit = 0
/ ! end of controls namelist
&pgstar
!pause = .true.
pgstar_lw = 8
Grid2_win_flag = .true.
Grid2_win_width = 12
Grid2_win_aspect_ratio = 0.5 ! aspect_ratio = height/width
Grid2_num_cols = 7 ! divide plotting region into this many equal width cols
Grid2_num_rows = 8 ! divide plotting region into this many equal height rows
Grid2_num_plots = 5 ! <= 10
Grid2_plot_name(1) = 'TRho'
Grid2_plot_row(1) = 1 ! number from 1 at top
Grid2_plot_rowspan(1) = 6 ! plot spans this number of rows
Grid2_plot_col(1) = 1 ! number from 1 at left
Grid2_plot_colspan(1) = 2 ! plot spans this number of columns
Grid2_plot_pad_left(1) = 0.00 ! fraction of full window width for padding on left
Grid2_plot_pad_right(1) = 0.00 ! fraction of full window width for padding on right
Grid2_plot_pad_top(1) = 0.00 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(1) = 0.0 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(1) = 0.65 ! multiply txt_scale for subplot by this
Grid2_plot_name(2) = 'Text_Summary1'
Grid2_plot_row(2) = 7 ! number from 1 at top
Grid2_plot_rowspan(2) = 2 ! plot spans this number of rows
Grid2_plot_col(2) = 1 ! number from 1 at left
Grid2_plot_colspan(2) = 4 ! plot spans this number of columns
Grid2_plot_pad_left(2) = -0.08 ! fraction of full window width for padding on left
Grid2_plot_pad_right(2) = -0.12 ! fraction of full window width for padding on right
Grid2_plot_pad_top(2) = 0.08 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(2) = -0.04 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(2) = 0.19 ! multiply txt_scale for subplot by this
Grid2_plot_name(3) = 'Abundance'
Grid2_plot_row(3) = 1 ! number from 1 at top
Grid2_plot_rowspan(3) = 3 ! plot spans this number of rows
Grid2_plot_col(3) = 3 ! plot spans this number of columns
Grid2_plot_colspan(3) = 3 ! plot spans this number of columns
Grid2_plot_pad_left(3) = 0.07 ! fraction of full window width for padding on left
Grid2_plot_pad_right(3) = 0.07 ! fraction of full window width for padding on right
Grid2_plot_pad_top(3) = 0.0 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(3) = 0.06 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(3) = 0.65 ! multiply txt_scale for subplot by this
Abundance_xmax = 2.4 ! -101 ! only used if /= -101d0
Power_xmax = 2.4 ! -101 ! only used if /= -101d0
Grid2_plot_name(4) = 'Power'
Grid2_plot_row(4) = 4 ! number from 1 at top
Grid2_plot_rowspan(4) = 3 ! plot spans this number of rows
Grid2_plot_col(4) = 3 ! number from 1 at left
Grid2_plot_colspan(4) = 3 ! plot spans this number of columns
Grid2_plot_pad_left(4) = 0.07 ! fraction of full window width for padding on left
Grid2_plot_pad_right(4) = 0.07 ! fraction of full window width for padding on right
Grid2_plot_pad_top(4) = 0.06 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(4) = 0.0 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(4) = 0.65 ! multiply txt_scale for subplot by this
Grid2_plot_name(5) = 'Profile_Panels1'
Grid2_plot_row(5) = 1 ! number from 1 at top
Grid2_plot_rowspan(5) = 8 ! plot spans this number of rows
Grid2_plot_col(5) = 6 ! number from 1 at left
Grid2_plot_colspan(5) = 2 ! plot spans this number of columns
Grid2_plot_pad_left(5) = 0.05 ! fraction of full window width for padding on left
Grid2_plot_pad_right(5) = 0.03 ! fraction of full window width for padding on right
Grid2_plot_pad_top(5) = 0.0 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(5) = 0.0 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(5) = 0.65 ! multiply txt_scale for subplot by this
TRho_logRho_dlogRho_min = 0.01
TRho_logT_dlogT_min = 0.02
HR_dlogT_min = 0.02
HR_dlogL_min = 0.01
Text_Summary1_name(2,1) = 'star_age'
! file output
!Grid2_file_flag = .true.
Grid2_file_dir = 'png'
Grid2_file_prefix = 'grid_'
Grid2_file_cnt = 10 ! 1 ! output when mod(model_number,Grid2_file_cnt)==0
Grid2_file_width = -1 ! negative means use same value as for window
Grid2_file_aspect_ratio = -1 ! negative means use same value as for window
show_TRho_Profile_eos_regions = .false.
TRho_Profile_xmin = -5.1
TRho_Profile_xmax = 10.2
TRho_Profile_ymin = 5.1
TRho_Profile_ymax = 10.2
Text_Summary1_name(1,3) = 'log_cntr_Rho'
Text_Summary1_name(2,3) = 'log_cntr_T'
Text_Summary1_name(3,3) = 'max_T_mass'
Text_Summary1_name(4,3) = 'max_T_lgRho'
Text_Summary1_name(5,3) = 'max_T_lgT'
!Text_Summary1_name(6,3) = 'center_ye'
!Text_Summary1_name(7,3) = 'iron_core_infall'
!Text_Summary1_name(8,3) = 'non_iron_core_infall'
!Text_Summary1_name(8,1) = 'center o16'
Text_Summary1_name(8,1) = 'species'
!Abundance_win_flag = .true.
Abundance_win_width = 6
Abundance_win_aspect_ratio = 1 ! aspect_ratio = height/width
! Abundance_xaxis_name = 'logR'
! Abundance_xaxis_reversed = .false.
! Abundance_xmin = -4 ! only used if /= -101d0
! Abundance_xmax = -2 ! -101 ! only used if /= -101d0
Abundance_xaxis_name = 'mass'
Abundance_xaxis_reversed = .false.
Abundance_xmin = 0 ! only used if /= -101d0
Abundance_xmax = 2.4 ! -101 ! only used if /= -101d0
Abundance_line_txt_scale_factor = 1.1 ! relative to other text
Abundance_legend_txt_scale_factor = 1.1 ! relative to other text
Abundance_legend_max_cnt = 0
Abundance_log_mass_frac_min = -3 ! only used if < 0
!Abundance_file_flag = .true.
Abundance_file_dir = 'png'
Abundance_file_prefix = 'abund_'
Abundance_file_cnt = 5 ! output when mod(model_number,Abundance_file_cnt)==0
Abundance_file_width = -1 ! (inches) negative means use same value as for window
Abundance_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 = 1 ! aspect_ratio = height/width
Profile_Panels1_title = 'Profile Panels'
Profile_Panels1_num_panels = 3
Profile_Panels1_yaxis_name(1) = 'logT'
Profile_Panels1_other_yaxis_name(1) = 'entropy'
Profile_Panels1_other_ymin(1) = -101d0 ! 0 ! only used if /= -101d0
Profile_Panels1_other_ymax(1) = -101d0 ! only used if /= -101d0
Profile_Panels1_other_ymargin(1) = 0.1
Profile_Panels1_yaxis_name(2) = 'logRho'
Profile_Panels1_other_yaxis_name(2) = 'logR_cm'
Profile_Panels1_other_ymin(2) = -101d0 ! only used if /= -101d0
Profile_Panels1_other_ymax(2) = -101d0 ! only used if /= -101d0
Profile_Panels1_other_ymargin(2) = 0.1
Profile_Panels1_yaxis_name(3) = 'ye'
Profile_Panels1_ymin(3) = -101d0 ! 0.43 ! only used if /= -101d0
Profile_Panels1_ymax(3) = -101d0 ! 0.51 ! only used if /= -101d0
Profile_Panels1_ymargin(3) = 0.1
Profile_Panels1_other_yaxis_name(3) = 'vel_km_per_s'
Profile_Panels1_other_ymin(3) = -101d0 ! -1200 ! only used if /= -101d0
Profile_Panels1_other_ymax(3) = -101d0 ! 300d0 ! only used if /= -101d0
Profile_Panels1_other_ymargin(3) = 0.1 ! 10
! Profile_Panels1_xaxis_name = 'logR'
! Profile_Panels1_xaxis_reversed = .false.
! Profile_Panels1_xmin = -4 ! only used if /= -101d0
! Profile_Panels1_xmax = -2 ! -101d0 ! only used if /= -101d0
Profile_Panels1_xaxis_name = 'mass'
Profile_Panels1_xaxis_reversed = .false.
Profile_Panels1_xmin = 0.0 ! only used if /= -101d0
Profile_Panels1_xmax = 2.4 ! -101d0 ! only used if /= -101d0
Profile_Panels1_xmargin = 0d0
/ ! end of pgstar namelist
-------------- next part --------------
! inlist_project
&star_job
mesa_dir = '' ! empty string defaults to environment variable MESA_DIR
history_columns_file = './history_columns.list' ! if null string, use default.
profile_columns_file = './profile_columns.list' ! if null string, use default.
! show_log_description_at_start = .true.
show_net_species_info = .true.
show_net_reactions_info = .true.
! create_pre_main_sequence_model = .true
change_net = .true.
new_net_name = 'mesa_45.net'
set_rates_preference = .true. ! for use by net + rates modules
new_rates_preference = 2 ! only used if set_rates_preference is true
! 1 = NACRE rates -- this is the default
! 2 = jina reaclib rates
eos_file_prefix = 'mesa'
kappa_file_prefix = 'gs98'
set_initial_age = .true.
initial_age = 0 ! in years
set_initial_model_number = .true.
initial_model_number = 0
!SOLAR CHEMICAL COMPOSITION
change_initial_Z = .true.
new_Z = 0.019d0
change_initial_Y = .true.
new_Y=0.27
save_photo_when_terminate = .true.
pause_before_terminate = .false.
pgstar_flag = .true.
/ !end of star_job
&controls
!enforce homogeneous mesh at high temperatures
use_other_mesh_functions = .true.
!avoid MLT++ oscillations
which_atm_option = 'Eddington_grey'
which_atm_off_table_option = 'grey_and_kap'
Pextra_factor = 1
!report_why_dt_limits = .true.
!report_all_dt_limits = .true.
! starting specifications
initial_mass = 30
!solar chemical composition
initial_Z = 0.02
!initial_Y = 0.27
!Mass Loss
use_other_wind = .true.
RGB_wind_scheme = 'other'
AGB_wind_scheme = 'other'
!hot wind efficiency
x_ctrl(1) = 1.0
!cool wind efficiency
x_ctrl(2) = 1.0
!WR wind efficiency
x_ctrl(3) = 1.0
! when to stop
fe_core_infall_limit = 1d8
!Resolution
mesh_delta_coeff_for_highT = 1.0
logT_max_for_standard_mesh_delta_coeff = 8.0
logT_min_for_highT_mesh_delta_coeff = 8.5
max_allowed_nz=300000000
varcontrol_target = 1d-4
mesh_delta_coeff = 1.0
mesh_max_allowed_ratio=2.5
min_timestep_limit = 1d-12 !in sec
max_dq = 0.5d-4
!nuc sets before Si burning
dX_nuc_drop_limit = 1d-4
dX_nuc_drop_limit_at_high_T = 1d-4
dX_nuc_drop_min_X_limit = 1d-4
dX_nuc_drop_max_A_limit = 52
!Mesh refinement (ready to push the evolution pass O depletion)
mesh_dlog_pp_dlogP_extra = 0.015
mesh_dlog_cno_dlogP_extra = 0.015
mesh_dlog_3alf_dlogP_extra = 0.015
mesh_dlog_burn_c_dlogP_extra = 0.015
mesh_dlog_burn_n_dlogP_extra = 0.015
mesh_dlog_burn_o_dlogP_extra = 0.015
mesh_dlog_burn_ne_dlogP_extra = 0.015
mesh_dlog_burn_na_dlogP_extra = 0.015
mesh_dlog_burn_mg_dlogP_extra = 0.015
mesh_dlog_cc_dlogP_extra = 0.015
mesh_dlog_co_dlogP_extra = 0.015
mesh_dlog_oo_dlogP_extra = 0.015
mesh_dlog_burn_si_dlogP_extra = 0.015
mesh_dlog_burn_s_dlogP_extra = 0.015
mesh_dlog_burn_ar_dlogP_extra = 0.015
mesh_dlog_burn_ca_dlogP_extra = 0.015
mesh_dlog_burn_ti_dlogP_extra = 0.015
mesh_dlog_burn_cr_dlogP_extra = 0.015
mesh_dlog_burn_fe_dlogP_extra = 0.015
mesh_dlog_pnhe4_dlogP_extra = 0.015
mesh_dlog_other_dlogP_extra = 0.015
mesh_dlog_photo_dlogP_extra = 0.015
!first 2 functions set in inlist_massive_defaults
xa_function_species(3) = 'c12'
xa_function_weight(3) = 20
xa_function_param(3) = 1d-2
xa_function_species(4) = 'o16'
xa_function_weight(4) = 20
xa_function_param(4) = 1d-2
xa_function_species(5) = 'ne20'
xa_function_weight(5) = 20
xa_function_param(5) = 1d-2
xa_function_species(6) = 'mg24'
xa_function_weight(6) = 20
xa_function_param(6) = 1d-2
!The following will be useless for the
!wind study, I don't go that far in the evolution
xa_function_species(7) = 'si28'
xa_function_weight(7) = 20
xa_function_param(7) = 1d-2
xa_function_species(8) = 's32'
xa_function_weight(8) = 20
xa_function_param(8) = 1d-2
xa_function_species(9) = 'fe56'
xa_function_weight(9) = 20
xa_function_param(9) = 1d-2
! output to files and terminal
photostep = 10
profile_interval = 10
history_interval = 1
terminal_interval = 10
write_header_frequency = 1
max_num_profile_models = 1000
!GYRE
!write_pulse_info_with_profile = .true.
!pulse_info_format = 'GYRE'
mixing_length_alpha = 2.0
MLT_option = 'Henyey'
alpha_semiconvection = 0.1
!f0 is the fraction of the pressure scale height at which the diffusion coefficient
!is taken to be extended beyond the convective boundary for overshooting
!f is the overshooting parameter for exponential overshooting, see Paxton et al. 2011
overshoot_f_above_nonburn_core = 0.025
overshoot_f0_above_nonburn_core = 0.05
overshoot_f_above_nonburn_shell = 0.025
overshoot_f0_above_nonburn_shell = 0.05
overshoot_f_below_nonburn_shell = 0.025
overshoot_f0_below_nonburn_shell = 0.05
overshoot_f_above_burn_h_core = 0.025
overshoot_f0_above_burn_h_core = 0.05
overshoot_f_above_burn_h_shell = 0.025
overshoot_f0_above_burn_h_shell = 0.05
overshoot_f_below_burn_h_shell = 0.025
overshoot_f0_below_burn_h_shell = 0.05
overshoot_f_above_burn_he_core = 0.025
overshoot_f0_above_burn_he_core = 0.05
overshoot_f_above_burn_he_shell = 0.025
overshoot_f0_above_burn_he_shell = 0.05
overshoot_f_below_burn_he_shell = 0.025
overshoot_f0_below_burn_he_shell = 0.05
overshoot_f_above_burn_z_core = 0.025
overshoot_f0_above_burn_z_core = 0.05
overshoot_f_above_burn_z_shell = 0.025
overshoot_f0_above_burn_z_shell = 0.05
overshoot_f_below_burn_z_shell = 0.025
overshoot_f0_below_burn_z_shell = 0.05
/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: run_star_extras.f
Type: text/x-fortran
Size: 27066 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150824/429b8b50/attachment.bin>
More information about the Mesa-users
mailing list