make_zams

This test case shows an example of creating a 4 \({\rm M}_\odot\), Z = 0.01 metallicity, pre-main sequence model and evolving it to the zero age main sequence.

This test case has 1 part. Click to see a larger version of a plot.

• Part 1 (inlist_zams) creates a 4 \({\rm M}_\odot\), Z = 0.01 metallicity, pre-main sequence model and evolves the model to the zero age main sequence, terminating when Lnuc / L = 0.9.

pgstar commands used for the plot above:

&pgstar

   file_device = 'vcps'   ! postscript

   pgstar_interval = 10

   Abundance_win_flag = .false.
   Abundance_win_width = 15
   Abundance_win_aspect_ratio = 0.75

   Abundance_title = ''
   Abundance_num_isos_to_show = 6

   Abundance_which_isos_to_show(1) = 'h1'
   Abundance_which_isos_to_show(2) = 'he3'
   Abundance_which_isos_to_show(3) = 'he4'
   Abundance_which_isos_to_show(4) = 'c12'
   Abundance_which_isos_to_show(5) = 'n14'
   Abundance_which_isos_to_show(6) = 'o16'

   num_abundance_line_labels = 4
   Abundance_line_txt_scale_factor = 1.0
   Abundance_legend_max_cnt = 0

   Abundance_xaxis_name = 'mass'
   Abundance_xaxis_reversed = .false.
   Abundance_xmin = 0.0
   Abundance_xmax = 1.0

   Abundance_log_mass_frac_min = -3.0
   Abundance_log_mass_frac_max =  0.3

   HR_win_flag = .false.
   HR_win_width = 12
   HR_win_aspect_ratio = 1.2  ! aspect_ratio = height/width
   HR_title = ''

   HR_logT_min = 3.5
   HR_logT_max = 4.5
   HR_logL_min = 1.0
   HR_logL_max = 3.0

   Text_Summary1_win_flag = .false.
   Text_Summary1_win_width = 10
   Text_Summary1_win_aspect_ratio = 0.15

   Text_Summary1_xleft = 0.01
   Text_Summary1_xright = 0.99
   Text_Summary1_ybot = 0.0
   Text_Summary1_ytop = 1.0
   Text_Summary1_txt_scale = 1.0

   Text_Summary1_num_rows = 1  ! <= 20
   Text_Summary1_num_cols = 3  ! <= 20
   Text_Summary1_name(:,:) = ''

   Text_Summary1_name(1,1) = 'star_mass'
   Text_Summary1_name(1,2) = 'log_LH'
   Text_Summary1_name(1,3) = 'num_zones'

   Grid1_title = 'make_zams'
   Grid1_win_flag = .false.
   Grid1_win_width = 15
   Grid1_win_aspect_ratio = 0.6

   Grid1_plot_name(:) = ''
   Grid1_plot_row(:) = 1            ! number from 1 at top
   Grid1_plot_rowspan(:) = 1        ! plot spans this number of rows
   Grid1_plot_col(:) =  1           ! number from 1 at left
   Grid1_plot_colspan(:) = 1        ! plot spans this number of columns
   Grid1_plot_pad_left(:) = 0.0     ! fraction of full window width for padding on left
   Grid1_plot_pad_right(:) = 0.0    ! fraction of full window width for padding on right
   Grid1_plot_pad_top(:) = 0.0      ! fraction of full window height for padding at top
   Grid1_plot_pad_bot(:) = 0.0      ! fraction of full window height for padding at bottom
   Grid1_txt_scale_factor(:) = 0.7  ! multiply txt_scale for subplot by this

   Grid1_num_cols = 6 ! divide plotting region into this many equal width cols
   Grid1_num_rows = 2 ! divide plotting region into this many equal height rows
   Grid1_num_plots = 10 ! <= 10

   pgstar_show_model_number = .true.
   pgstar_model_scale = 1.0
   pgstar_model_lw = 3
   pgstar_model_disp = 2.0
   pgstar_model_coord = 0.92
   pgstar_model_fjust = 1.0

   pgstar_show_age = .true.
   pgstar_age_scale = 1.0
   pgstar_age_lw = 3
   pgstar_age_disp = 2.0
   pgstar_age_coord = -0.10
   pgstar_age_fjust = 0.0

   Grid1_plot_name(1) = 'Text_Summary1'
   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) = 6        ! plot spans this number of columns

   Grid1_plot_pad_left(1) = -0.06   ! 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.02    ! fraction of full window height for padding at top
   Grid1_plot_pad_bot(1) = 0.39     ! fraction of full window height for padding at bottom
   Grid1_txt_scale_factor(1) = 1.2  ! multiply txt_scale for subplot by this

   Grid1_plot_name(2) = 'Abundance'
   Grid1_plot_row(2) = 1            ! number from 1 at top
   Grid1_plot_rowspan(2) = 2        ! plot spans this number of rows
   Grid1_plot_col(2) =  1           ! 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.10   ! fraction of full window width for padding on right
   Grid1_plot_pad_top(2) = 0.03     ! fraction of full window height for padding at top
   Grid1_plot_pad_bot(2) = 0.03     ! fraction of full window height for padding at bottom
   Grid1_txt_scale_factor(2) = 0.7  ! multiply txt_scale for subplot by this

   Grid1_plot_name(3) = 'HR'
   Grid1_plot_row(3) = 1            ! number from 1 at top
   Grid1_plot_rowspan(3) = 2        ! plot spans this number of rows
   Grid1_plot_col(3) =  5           ! Number from 1 at left
   Grid1_plot_colspan(3) = 3        ! plot spans this number of columns

   Grid1_plot_pad_left(3) = -0.15   ! fraction of full window width for padding on left
   Grid1_plot_pad_right(3) = 0.20   ! fraction of full window width for padding on right
   Grid1_plot_pad_top(3) = 0.03     ! fraction of full window height for padding at top
   Grid1_plot_pad_bot(3) = 0.03     ! fraction of full window height for padding at bottom
   Grid1_txt_scale_factor(3) = 0.7  ! multiply txt_scale for subplot by this

   Grid1_file_flag = .true.
   Grid1_file_dir = 'pgstar_out'
   Grid1_file_prefix = 'grid_'
   Grid1_file_interval = 10000
   Grid1_file_width = -1   ! (inches) negative means use same value as for window

 / ! end of pgstar namelist

Last-Updated: 20Jun2021 (MESA e2acbc2) by fxt.

Last-Run: 22Oct2024 (MESA 9b2017ca) by pmocz on C916PXT6XW in 46 seconds using 8 threads.