[Mesa-users] Reactions
MEGAN ALLEN
M.ALLEN-2018 at hull.ac.uk
Thu Oct 20 10:37:58 UTC 2022
Hi,
This was the inlist I used.
Thanks,
Megan
________________________________
From: Francis Timmes <fxt44 at mac.com>
Sent: 19 October 2022 07:01
To: MEGAN ALLEN <M.ALLEN-2018 at hull.ac.uk>
Cc: Frank Timmes <fxt44 at mac.com>; mesa-users <mesa-users at lists.mesastar.org>
Subject: Re: [Mesa-users] Reactions
hi negan,
its difficult to assess with just a plot.
i'll echo rob's question from his 05sep2022 reply to your previous message
https://linkprotect.cudasvc.com/url?a=https%3a%2f%2flists.mesastar.org%2fpipermail%2fmesa-users%2f2022-September%2f014010.html&c=E,1,ELB0WHnXV9oAXwzBiRIfphYSayY5RdEEoAcM5q6u7iR_JT4P5PRa702pJ2Sw-1kQdlkZxxJhkEfDkaLyBHjkTPMVz40r7gX3noxsn6ZQOh7GrO-FXEqg-Cjp4g,,&typo=1
> Can you show your inlists?
fxt
> On Oct 17, 2022, at 3:43 AM, MEGAN ALLEN via Mesa-users <mesa-users at lists.mesastar.org> wrote:
>
> Hi,
>
> I have asked about this previously but I'm having trouble with a set of reactions. I'm using revision 10108.
> The only difference between the two models in the HR-diagram attached is in the https://linkprotect.cudasvc.com/url?a=https%3a%2f%2fAGB.net&c=E,1,LhBJaRmeVE7-22EgX8Q3pVb30-5h44RYlHkTPFpAji2n4qYHdrCgFjgv5gHioZPANY8yq-mLHILBtHyMhdF95d1ad6cbquju6Xl0rzFB8yajhjQ,&typo=1 model there are the following reactions:
> r_c12_pg_n13
> r_n13_wk_c13
> r_c13_pg_n14
>
> In rc12_to_n14 the 3 reactions above are substituted for that single reaction.
>
> I was trying to replicate an inlist that used revision 3709 and the model using rc12_to_n14 follows the HR diagram from the inlist I tried to update, however using this reaction I can't look for C13 pockets which is why I needed the set of 3 reactions but they produce a different result.
>
> Everything else about the two inlists and reactions are identical and as far as I can see both the set of 3 reactions and rc12_to_n14 should produce the same result so I was hoping to find out why there is such a significant difference.
>
> Thanks,
> Megan
> <HRD_ReactionComparison.png>_______________________________________________
> mesa-users at lists.mesastar.org
> https://linkprotect.cudasvc.com/url?a=https%3a%2f%2flists.mesastar.org%2fmailman%2flistinfo%2fmesa-users&c=E,1,7E0F5ZXS50-Tc2HI9qiUP5CVjtVHVbQHZpWdmNOKgjQsTXT52uA0RrztmWtRUXrFIc4P3kys8GVlu9u6DwUwq0yOCFRgDBsl0gOUmJ4mDHkVL7N_oEBmXzpB6Q0,&typo=1
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20221020/510a5a01/attachment.htm>
-------------- next part --------------
&star_job\
create_pre_main_sequence_model = .true.
save_model_number = 100
save_model_filename = '100.mod'
save_model_when_terminate = .true.
save_model_filename = 'final.mod'
initial_zfracs = 2
kappa_file_prefix = 'gn93'
kappa_CO_prefix = 'gn93_co'
change_net = .true.
change_initial_net = .true.
new_net_name = 'agb.net'
file_for_uniform_xa = 'falk_xa.list.agb'
set_uniform_xa_from_file = .true.
set_rates_preference = .true.
new_rates_preference = 1
set_rate_c12ag = 'Kunz'
set_rate_n14pg = 'jina reaclib'
set_rate_3a = 'jina reaclib'
change_v_flag = .true.
! change_initial_v_flag = .false.
new_v_flag = .true.
! end of star_job namelist
&controls
max_model_number = 10000
initial_mass = 1.65
initial_z = 0.01d0
varcontrol_target = 1d-4
mesh_delta_coeff = 0.5
split_mixing_choice = 0
T_mix_limit = 1d4
photo_interval = 500
profile_interval = 100
max_num_profile_models = 100
history_interval = 1
terminal_interval = 10
write_header_frequency = 10
max_allowed_nz = 5000
mixing_length_alpha = 1.73
xa_function_species(1) = 'h1'
xa_function_weight(1) = 5
xa_function_param(1) = 3d-4
xa_function_species(2) = 'he4'
xa_function_weight(2) = 5
xa_function_param(2) = 1d-2
xa_function_species(3) = 'c13'
xa_function_weight(3) = 5
xa_function_param(3) = 1d-4
xa_function_species(4) = 'n14'
xa_function_weight(4) = 4
xa_function_param(4) = 1d-6
varcontrol_target = 1d-4
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
delta_lgL_He_limit = 0.01d0
!delta_lgL_He_hard_limit = 0.05
lgL_He_burn_min = 2.0
DUP_varcontrol_factor = 1d0
dH_limit_min_H = 1.0d-2
dH_limit = 0.1
dH_div_H_limit_min_H = 1d-5
dH_div_H_limit = 0.5
!log_L_lower_limit = -1
!smooth_convective_bdy = .true.
!convective_bdy_weight = 1
limit_overshoot_Hp_using_size_of_convection_zone = .false.
overshoot_below_noburn_shell_factor = 9
ovr_below_burn_he_shell_factor = 0.5714
mass_for_overshoot_full_on = 1.8
mass_for_overshoot_full_off = 1.1
overshoot_f_above_nonburn_core = 0.014
overshoot_f0_above_nonburn_core = 0.014
overshoot_f_above_nonburn_shell = 0.014
overshoot_f0_above_nonburn_shell = 0.014
overshoot_f_below_nonburn_shell = 0.014
overshoot_f0_below_nonburn_shell = 0.014
overshoot_f_above_burn_h_core = 0.014
overshoot_f0_above_burn_h_core = 0.014
overshoot_f_above_burn_h_shell = 0.014
overshoot_f0_above_burn_h_shell = 0.014
overshoot_f_below_burn_h_shell = 0.014
overshoot_f0_below_burn_h_shell = 0.014
overshoot_f_above_burn_he_core = 0.014
overshoot_f0_above_burn_he_core = 0.014
overshoot_f_above_burn_he_shell = 0.014
overshoot_f0_above_burn_he_shell = 0.014
overshoot_f_below_burn_he_shell = 0.014
overshoot_f0_below_burn_he_shell = 0.014
overshoot_f_above_burn_z_core = 0.014
overshoot_f0_above_burn_z_core = 0.014
overshoot_f_above_burn_z_shell = 0.014
overshoot_f0_above_burn_z_shell = 0.014
overshoot_f_below_burn_z_shell = 0.014
overshoot_f0_below_burn_z_shell = 0.014
!overshoot_D2_above_nonburn = 100000000000d0
!overshoot_D2_below_nonburn = 100000000000d0
!overshoot_D2_above_burn_h = 100000000000d0
!overshoot_D2_below_burn_h = 100000000000d0
!overshoot_D2_above_burn_he = 100000000000d0
!overshoot_D2_below_burn_he = 100000000000d0
!overshoot_D2_above_burn_z = 100000000000d0
!overshoot_D2_below_burn_z = 100000000000d0
!overshoot_f2_above_nonburn = 0.27d0
!overshoot_f2_below_nonburn = 0.27d0
!overshoot_f2_above_burn_h = 0.27d0
!overshoot_f2_below_burn_h = 0.27d0
!overshoot_f2_above_burn_he = 0.27d0
!overshoot_f2_below_burn_he = 0.27d0
!overshoot_f2_above_burn_z = 0.27d0
!overshoot_f2_below_burn_z = 0.27d0
which_atm_option = 'simple_photosphere'
!cool_wind_full_on_T = 9.99d9
!hot_wind_full_on_T = 1d10
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.04d0
! super_eddington_scaling_factor = 1
use_Type2_opacities = .true.
kap_Type2_full_off_X = 0.99d0
kap_Type2_full_on_X = 0.98d0
kap_Type2_full_off_dZ = 0.00000d0
kap_Type2_full_on_dZ = 0.0000d0
Zbase = 0.01 !0.01
}
More information about the Mesa-users
mailing list