This test case exercises radiative levitation and the OP mono opacities in the outer layers of a 0.466 Msun, Z=0.02 metallicity, B-type subdwarf (sdB) model.
This test case has 1 part. Click to see a larger version of a plot.
This test case requires the OP_mono data. See the section
OP mono opacities in
controls.defaults for information on how to obtain this data and set the relevant paths.
Part 1 (
inlist_radiative_levitation) loads a pre-built 0.466 Msun, Z=0.02 metallicity, B-type subdwarf (sdB) model, created from make_sdb with an envelope thickness x_ctrl(1) set to 1e-3. Evolution commences with element diffusion activated (
do_element_diffusion = .true.) and radiative levitation enabled for Fe56 and Ni58 (
diffusion_min_Z_for_radaccel = 26and
diffusion_max_Z_for_radaccel = 28) for temperatures less than 1e7 K (
diffusion_max_T_for_radaccel = 1e7). The evolution terminates after 1 million years. The
run_star_extras.f90then checks that at a mass location of 1e-10 Msun that the mass fraction of Fe56 is greater than 0.02, the mass fraction of Ni58 is greater than 0.007, and that their sum is less than 0.8. If so, a success message is written to the terminal:
found expected effects of radiative levitation
pgstar commands used for the plot above:
&pgstar file_white_on_black_flag = .true. ! white_on_black flags -- true means white foreground color on black background !file_device = 'png' ! png !file_extension = 'png' file_device = 'vcps' ! postscript file_extension = 'ps' pgstar_interval = 10 Profile_Panels3_win_flag = .true. Profile_Panels3_win_width = 10 Profile_Panels3_title = 'radiative_levitation' Profile_Panels3_xaxis_name = 'logxm' Profile_Panels3_xaxis_reversed = .true. Profile_Panels3_xmin = -14 ! only used if /= -101d0 Profile_Panels3_xmax = -1 ! only used if /= -101d0 Profile_Panels3_num_panels = 3 Profile_Panels3_yaxis_name(2) = 'opacity' Profile_Panels3_other_yaxis_name(2) = 'gradr' Profile_Panels3_file_flag = .true. Profile_Panels3_file_dir = 'pgstar_out' Profile_Panels3_file_prefix = 'profile_' Profile_Panels3_file_interval = 10000 Profile_Panels3_file_width = 10 Profile_Panels3_file_aspect_ratio = -1 / ! end of pgstar namelist
Last-Updated: 24Jun2021 (MESA e2acbc2) by fxt.