[Mesa-users] Rotating stars with omega/omega_critical

[AVG]

Dear all,

I am trying to simulate a rotating star giving MESA the value of its 
omega/omega_critical, but I failed no matter how I try it.

I have been using the information of the documentation about this 
(https://docs.mesastar.org/en/release-r21.12.1/reference/star_job.html?highlight=omega_crit#set-omega-div-omega-crit), 
but I haven't been able to implement the rotation with it. I have used 
all the possible combinations of the parameters 
'set_omega_div_omega_crit' and 'set_initial_omega_div_omega_crit', with 
and without the 'set_omega_div_omega_crit_step_limit = -1' active, but 
it seems that there is something I am missing.

I attach one of the inlists I have used. I hope anyone can help me with 
this or highlight my failures.

Best regards,

Antonio Varo
-------------- next part --------------

&star_job
  ! see star/defaults/star_job.defaults  
  ! begin with a pre-main sequence model
    create_pre_main_sequence_model = .true.

  ! save a model at the end of the run
    save_model_when_terminate = .false.
    save_model_filename = '15M_at_TAMS.mod'

  ! display on-screen plots
    pgstar_flag = .false.
    
  ! My options
  ! to include rotation in the model
    new_rotation_flag = .true.
    change_rotation_flag = .false.
    change_initial_rotation_flag = .false.

  ! Or in units of Omega_critical (DEFAULT)
    new_omega_div_omega_crit = 0.6
    set_omega_div_omega_crit = .true.
    set_initial_omega_div_omega_crit = .true.
    set_omega_div_omega_crit_step_limit = -1

  ! The following apply when rotation is off and model satisfies the "near zams" test. Each turns on rotation and sets step limit. Only applies if > 0

    set_near_zams_omega_div_omega_crit_steps = -1
    set_near_zams_surface_rotation_v_steps = 1

  ! The following are for settle in the new rotation (DEFAULT)

    num_steps_to_relax_rotation = 100
    relax_omega_max_yrs_dt = 1d4
   
    relax_omega_div_omega_crit = .false.
    relax_initial_omega_div_omega_crit = .false.
    near_zams_relax_omega_div_omega_crit = .false.

    relax_surface_rotation_v = .false.
    relax_initial_surface_rotation_v = .false.
    near_zams_relax_initial_surface_rotation_v = .true.

  ! To include diffusion of angular momentum

    new_D_omega_flag = .false.
    change_D_omega_flag = .false.
    change_initial_D_omega_flag = .false.

    new_am_nu_rot_flag = .false.
    change_am_nu_rot_flag = .false.
    change_initial_am_nu_rot_flag = .false.
    use_D_omega_for_am_nu_rot = .true. 
    
/ ! end of star_job namelist


&eos
  ! eos options
  ! see eos/defaults/eos.defaults

/ ! end of eos namelist


&kap
  ! kap options
  ! see kap/defaults/kap.defaults
  use_Type2_opacities = .true.
  Zbase = 0.02

/ ! end of kap namelist


&controls
  ! see star/defaults/controls.defaults
  ! starting specifications
    initial_mass = 2.9 ! in Msun units
  !  initial_z = 0.02

  ! stop when the center mass fraction of h1 drops below this limit
    xa_central_lower_limit_species(1) = 'h1'
    xa_central_lower_limit(1) = 1d-3
    
  ! stop when total power from reactions burning helium (in Lsun units) is > this.
  !  power_he_burn_upper_limit = 1d3
  max_age = 2d7
  write_pulse_data_with_profile = .true.
  pulse_data_format = 'GYRE'
  
/ ! end of controls namelist