1.3M_ms_high_Z
The test checks the evolution of metal-rich low-mass stars by evolving a 1.3 Msun, metal-rich Z=0.04 model from the pre-main sequence to core hydrogen depletion.
This test case has two parts.
Part 1 (
inlist_zams
) creates the pre-main-sequence model and stops near zams.Part 2 (
inlist_1.3M_ms_high_Z
) evolves the model from zams to when the luminosity reaches log10(L/Lsun) = 0.7.
The final model, click on the image for a larger version, shows
The left plot shows the HR diagram. The yellow curve is a precalculated
HR track loaded from HR_OPAL.dat
, while the green curve is the model calculated.
The right plot shows a profile of the metal mass fraction. The curve rises above the
background metallicity as hydrogen burns to helium and goes below the background metallicity
as core hydrogen depletes.
pgstar commands used for the plots above:
&pgstar
pgstar_interval = 10
! device
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'
! build grid
Grid1_title = '1.3M_ms_high_Z'
Grid1_win_flag = .true.
Grid1_win_width = 15
Grid1_win_aspect_ratio = 0.4 ! aspect_ratio = height/width
Grid1_xleft = 0.05 ! fraction of full window width for margin on left
Grid1_xright = 0.95 ! fraction of full window width for margin on right
Grid1_ybot = 0.08 ! fraction of full window width for margin on bottom
Grid1_ytop = 0.92 ! fraction of full window width for margin on top
Grid1_num_cols = 5 ! divide plotting region into this many equal width cols
Grid1_num_rows = 1 ! divide plotting region into this many equal height rows
Grid1_num_plots = 2 ! <= 10
Grid1_plot_name(1) = 'HR'
Grid1_plot_row(1) = 1 ! number from 1 at top
Grid1_plot_rowspan(1) = 1 ! plot spans this number of rows
Grid1_plot_col(1) = 1 ! number from 1 at left
Grid1_plot_colspan(1) = 2 ! plot spans this number of columns
Grid1_plot_pad_left(1) = 0.025 ! fraction of full window width for padding on left
Grid1_plot_pad_right(1) = 0.05 ! fraction of full window width for padding on right
Grid1_plot_pad_top(1) = 0.05 ! fraction of full window height for padding at top
Grid1_plot_pad_bot(1) = 0.05 ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(1) = 1.0 ! multiply txt_scale for subplot by this
Grid1_plot_name(2) = 'Profile_Panels1'
Grid1_plot_row(2) = 1 ! number from 1 at top
Grid1_plot_rowspan(2) = 1 ! plot spans this number of rows
Grid1_plot_col(2) = 3 ! number from 1 at left
Grid1_plot_colspan(2) = 3 ! plot spans this number of columns
Grid1_plot_pad_left(2) = 0.05 ! fraction of full window width for padding on left
Grid1_plot_pad_right(2) = 0.025 ! fraction of full window width for padding on right
Grid1_plot_pad_top(2) = 0.05 ! fraction of full window height for padding at top
Grid1_plot_pad_bot(2) = 0.05 ! fraction of full window height for padding at bottom
Grid1_txt_scale_factor(2) = 1.0 ! multiply txt_scale for subplot by this
! file output
Grid1_file_flag = .true.
Grid1_file_dir = 'png'
Grid1_file_prefix = 'grid1'
Grid1_file_interval = 1000000 ! output when mod(model_number,Grid1_file_interval)==0
Grid1_file_width = 15 ! negative means use same value as for window
Grid1_file_aspect_ratio = -1 ! negative means use same value as for window
pgstar_model_scale = 1.4
pgstar_model_disp = 1.5
pgstar_model_coord = 0.98
pgstar_model_fjust = 1.0
pgstar_age_scale = 1.4
pgstar_age_disp = 1.5
pgstar_age_coord = -0.04
pgstar_age_fjust = 0.0
! HR diagram
HR_win_flag = .false.
HR_win_aspect_ratio = 1.0
! set static plot bounds
HR_logT_min = 3.55
HR_logT_max = 3.85
HR_logL_min = 0.1
HR_logL_max = 1.0
! show OPAL results
HR_fname = 'HR_OPAL.dat' ! file name for extra HR data
! profile
Profile_Panels1_win_flag = .false.
Profile_Panels1_num_panels = 1
Profile_Panels1_xaxis_name = 'mass'
Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
Profile_Panels1_xmax = -101d0 ! only used if /= -101d0
Profile_Panels1_yaxis_name(1) = 'z_mass_fraction_metals'
Profile_Panels1_ymin(1) = 0.038 ! only used if /= -101d0
Profile_Panels1_ymax(1) = 0.042 ! only used if /= -101d0
Profile_Panels1_other_yaxis_name(:) = ''
/ ! end of pgstar namelist
Last-Updated: 27May2021 (MESA ebecc10) by fxt