[Mesa-users] make_co_wd test case in r15140
Kevin Rorke
kevin.rorke at gmail.com
Mon Jan 10 17:58:15 UTC 2022
Hi everyone
I've been trying to run the make_co_wd test case in r15140 with a 5.7
Msun (or thereabouts) starting mass for some time now. Previously the
run failed during the inlist_remove_env step because of memory
constraints - it seems to require about 16GB of RAM. But even with that
amount of memory, the run eventually fails for any starting mass in that
region.
The only value I've changed is the initial_mass value in inlist_common.
And the environment variables are set to create a new zams model each
time. Can anyone suggest what I might be doing wrong?
Cheers
Kevin
Exceeded max_iters iterations in eval_T_tau_uniform
retry: error in PT_eqns_surf 1887
dt
5.6302130558697700D-21
min_timestep_limit
9.9999999999999995D-21
stopping because of problems dt < min_timestep_limit
terminated evolution: cannot find acceptable model
termination code: min_timestep_limit
DATE: 2022-01-10
TIME: 19:24:41
******************************************************************
failed to create co_wd.mod when running inlist_remove_env_header
******************************************************************
-------------- next part --------------
&star_job
show_log_description_at_start = .false.
change_net = .true.
new_net_name = 'co_burn_extras.net'
set_initial_age = .false.
set_initial_model_number = .false.
pgstar_flag = .true.
/ ! end of star_job namelist
&eos
/ ! end of eos namelist
&kap
Zbase = 0.02d0
kap_file_prefix = 'gs98'
use_Type2_opacities = .true.
/ ! end of kap namelist
&controls
initial_mass = 5.7d0
!initial_mass = 3.1d0 ! this default produces a ~0.6 Msun WD.
! initial_mass in the range 1.0-6.5 Msun will produce CO WDs
! with masses 0.5-1.0 Msun.
! Initial mass around 7 or larger will likely ignite carbon in
! the core, and the test will fail because this would produce an
! ONe WD.
initial_z = 0.02d0
! opacity
! when to stop
max_age = 14d9
! mesh
mesh_delta_coeff = 1.5
! diffusion
!do_element_diffusion = .true.
!do_WD_sedimentation_heating = .true.
!do_diffusion_heating = .false.
!show_diffusion_info = .true. ! terminal output for diffusion
diffusion_v_max = 1d-5
! Bypasses the need to set up diffusion classes.
diffusion_use_full_net = .true.
! Solver controls.
diffusion_use_cgs_solver = .true.
diffusion_use_isolve = .true.
diffusion_rtol_for_isolve = 1d-4
diffusion_atol_for_isolve = 1d-5
diffusion_min_X_hard_limit = -1d-4
diffusion_maxsteps_for_isolve = 1000
diffusion_isolve_solver = 'ros2_solver'
! Timestep controls to prevent steps that are difficult for diffusion.
diffusion_steps_limit = 20
diffusion_steps_hard_limit = 50
diffusion_iters_limit = 50
diffusion_iters_hard_limit = 100
! convection
mixing_length_alpha = 2
use_Ledoux_criterion = .true.
alpha_semiconvection = 0.01
thermohaline_coeff = 2
! timesteps
max_years_for_timestep = 1d9
varcontrol_target = 1d-3
dX_nuc_drop_limit = 1d-2
delta_HR_limit = 0.05
delta_XC_cntr_limit = 0.02d0
delta_XO_cntr_limit = 0.02d0
delta_XC_cntr_hard_limit = 0.05d0
delta_XO_cntr_hard_limit = 0.05d0
! Atm
atm_option = 'T_tau'
atm_T_tau_relation = 'Eddington'
atm_T_tau_opacity = 'iterated'
atm_T_tau_max_iters = 100
! Overshoot
overshoot_scheme(1) = 'exponential'
overshoot_zone_type(1) = 'any'
overshoot_zone_loc(1) = 'any'
overshoot_bdy_loc(1) = 'any'
overshoot_f(1) = 0.014
overshoot_f0(1) = 0.004
! Output
history_interval = 10
terminal_interval = 5
write_header_frequency = 10
/ ! end of controls namelist
&pgstar
Grid2_win_flag = .true.
Grid2_win_width = 10
Grid2_win_aspect_ratio = 0.7 ! aspect_ratio = height/width
!Grid2_plot_name(2) = 'Power'
!Power_xmin = 0.510 ! only used if /= -101d0
!Power_xmax = 0.590 ! only used if /= -101d0
Grid2_file_flag = .true.
Grid2_file_dir = 'png'
Grid2_file_interval = 10 ! output when mod(model_number,Grid2_file_interval)==0
Grid2_file_width = -1 ! (inches) negative means use same value as for window
Grid2_file_aspect_ratio = -1 ! negative means use same value as for window
Grid2_num_plots = 5 ! <= 10
Grid2_plot_name(1) = 'TRho_Profile'
Grid2_plot_row(1) = 1 ! number from 1 at top
Grid2_plot_rowspan(1) = 3 ! 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.0 ! fraction of full window width for padding on left
Grid2_plot_pad_right(1) = 0.0 ! fraction of full window width for padding on right
Grid2_plot_pad_top(1) = 0.0 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(1) = 0.1 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(1) = 0.7 ! multiply txt_scale for subplot by this
Grid2_plot_name(2) = 'HR'
Grid2_plot_row(2) = 4 ! number from 1 at top
Grid2_plot_rowspan(2) = 3 ! plot spans this number of rows
Grid2_plot_col(2) = 1 ! number from 1 at left
Grid2_plot_colspan(2) = 1 ! plot spans this number of columns
Grid2_plot_pad_left(2) = 0.00 ! fraction of full window width for padding on left
Grid2_plot_pad_right(2) = 0.04 ! fraction of full window width for padding on right
Grid2_plot_pad_top(2) = 0.03 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(2) = 0.0 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(2) = 0.65 ! multiply txt_scale for subplot by this
Grid2_plot_name(3) = 'TRho'
Grid2_plot_row(3) = 4 ! number from 1 at top
Grid2_plot_rowspan(3) = 3 ! plot spans this number of rows
Grid2_plot_col(3) = 2 ! number from 1 at left
Grid2_plot_colspan(3) = 1 ! plot spans this number of columns
Grid2_plot_pad_left(3) = 0.04 ! fraction of full window width for padding on left
Grid2_plot_pad_right(3) = 0.00 ! fraction of full window width for padding on right
Grid2_plot_pad_top(3) = 0.03 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(3) = 0.0 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(3) = 0.65 ! multiply txt_scale for subplot by this
Grid2_plot_name(4) = 'Text_Summary1'
Grid2_plot_row(4) = 7 ! number from 1 at top
Grid2_plot_rowspan(4) = 2 ! plot spans this number of rows
Grid2_plot_col(4) = 1 ! number from 1 at left
Grid2_plot_colspan(4) = 4 ! plot spans this number of columns
Grid2_plot_pad_left(4) = -0.08 ! fraction of full window width for padding on left
Grid2_plot_pad_right(4) = 0.0 ! fraction of full window width for padding on right
Grid2_plot_pad_top(4) = 0.08 ! fraction of full window height for padding at top
Grid2_plot_pad_bot(4) = -0.04 ! fraction of full window height for padding at bottom
Grid2_txt_scale_factor(4) = 0.19 ! multiply txt_scale for subplot by this
Grid2_plot_name(5) = 'Profile_Panels3'
Grid2_plot_row(5) = 1 ! number from 1 at top
Grid2_plot_rowspan(5) = 6 ! plot spans this number of rows
Grid2_plot_col(5) = 3 ! number from 1 at left
Grid2_plot_colspan(5) = 2 ! plot spans this number of columns
Grid2_plot_pad_left(5) = 0.15 ! 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
Profile_Panels3_xaxis_name = 'mass'
Profile_Panels3_xaxis_reversed = .false.
Profile_Panels3_xmin = -101d0 ! only used if /= -101d0
Profile_Panels3_xmax = 2.2 ! -101d0 ! only used if /= -101d0
/ ! end of pgstar namelist
More information about the Mesa-users
mailing list