make_zams_ultra_high_mass

This test case shows an example of creating a 250 Msun, Z = 1e-4 metallicity, model close to the main sequence.

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

  • Part 1 (inlist_zams_ultra_high_mass) loads the pre-built 90 Msun, Z = 1e-4 metallicity, model late_pre_zams_90M.mod. This initial model is then relaxed to 250 Msun using relax_initial_mass_scale = .true.. At the 250 Msun point, the model already exceeds the termination criteria of Lnuc/L = 0.9 and there is no evolution in the HR diagram:

../_images/grid_000075.svg

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

 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 = -5.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 = 4.0
 HR_logT_max = 5.3
 HR_logL_min = 6.0
 HR_logL_max = 7.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 = .true.
 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.