rsp_BLAP

This test case checks the non-linear pulsation evolution of a 0.36 Msun, Teff = 26,000 K, L = 320 Lsun, Z = 0.05 metallicity model - a blue large-amplitude pulsator model originally contributed by Alfred Gautschy.

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

  • Part 1 (inlist_rsp_BEP) creates the initial 0.36 Msun, Teff = 26,000 K, L = 320 Lsun, Z = 0.05 metallicity model, and writes the results of conducting a linear nonadiabatic stability analysis to the LOGS directory (see Section 2.2 of MESA V for details). The evolution with RSP then begins. After 10 periods, the run_star_extras.f90 checks if the energy conservation is less than 1e-5 and if fundamental period is within 1% of the expected 0.05543 day period. If these values are within bounds, then a message is written to the terminal and the run terminates:

rel_run_E_err  -7.5272028903215587E-009
good match for period   5.5242272536049263E-002   5.5430000000000000E-002
../_images/grid_0007318.svg

pgstar commands, in addition to those in inlist_rsp_pgstar_default, 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 = 100

     pgstar_age_scale = 0.8
     pgstar_age_lw = 3
     pgstar_age_disp = 3.9
     pgstar_age_coord = -0.11
     pgstar_age_fjust = 0.0

     pgstar_model_disp = 3.9

     History_Panels2_txt_scale = 0.7
     Profile_Panels2_txt_scale = 0.6
     logL_R_txt_scale = 0.7
     logL_v_txt_scale = 0.7
     logL_Teff_txt_scale = 0.7

      Grid2_win_flag = .true.
      Grid2_win_width = 12
      Grid2_title = 'rsp_blap - blue large-amplitude pulsator'
      Grid2_txt_scale_factor(:) = 1.0

       Grid2_file_flag = .true.
       Grid2_file_dir = 'pgstar_out'
       Grid2_file_prefix = 'grid_'
       Grid2_file_interval = 10000
       Grid2_file_width = -1
       Grid2_file_aspect_ratio = -1

/ ! end of pgstar namelist

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