make_zams

This test case shows an example of creating a 4 Msun, 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 Msun, Z = 0.01 metallicity, pre-main sequence model and evolves the model to the zero age main sequence, terminating when Lnuc / L = 0.9.

../_images/grid_000192.svg

pgstar commands used for the plot above:

&pgstar


  file_device = 'vcps'          ! postscript
  file_extension = 'ps'

  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.