[mesa-users] Rotation : relax_initial_omega_div_omega_crit + preserving omega in several inlists

[Héctor MR]

  Dear all,

  I have a question regarding the &star_job command
relax_initial_omega_div_omega_crit. [Disclosure: this is the first time I
am implementing rotation in MESA]

  I am analyzing rotation on an accreting, 1Msun-WD with radius
R=8.25E-3*Rsun, so that the break-up frequency sqrt(G*M/R**3) is ~0.83;
however, when I use new_omega_div_omega_crit = 0.1, I get


  num_steps_to_relax_rotation         100
        relax to omega: wanted-current, current, wanted           1
1.7485786361779884D-05    0.0000000000000000D+00    *1.7485786361779884D-05*
        relax to omega: wanted-current, current, wanted           2
1.7330581050599486D-05    1.7485788969702041D-07    1.7505438940296507D-05
        relax to omega: wanted-current, current, wanted           3
1.7177934014914814D-05    3.5010756599501660D-07    1.7528041580909830D-05
        relax to omega: wanted-current, current, wanted           4
1.7016693789573586D-05    5.2583873499324927D-07    1.7542532524566834D-05
        relax to omega: wanted-current, current, wanted           5
1.6850228116883350D-05    7.0169673253464221D-07    1.7551924849417992D-05
        relax to omega: wanted-current, current, wanted           6
1.6680900334687578D-05    8.7758847614349309D-07    1.7558488810831069D-05
        relax to omega: wanted-current, current, wanted           7
1.6513051889555891D-05    1.0534860275206089D-06    1.7566537917076500D-05


  Rather than, let's say, 0.083, as I was expecting. I have done some
"research" with "grep -r" in $MESA_DIR/star/private:

  /relax.f90: integer function relax_omega_check_model

  /do_one_utils: integer function do_check_limits

  /star_utils.f90: real(dp) function omega_crit(s, k), subroutine
set_surf_avg_rotation_info(s)

  /hydro_vars.f90: subroutine set_rmid_and_Amid        !definition of rmid


  And have not found anything weird. Is there anyone who can tell me what I
am either doing wrong, obviating or misuderstanding? This sums everything
up:

      new_rotation_flag = .true.
      change_initial_rotation_flag = .true.

      !set_initial_omega_div_omega_crit = .true.
      new_omega_div_omega_crit = 0.1

      relax_initial_omega_div_omega_crit = .true.
      num_steps_to_relax_rotation = 100


  Finally, my inlists are split in two pieces. How can I preserve the last
value of omega_div_omega_crit in the first one for the initial angular
velocity in the second one? Via x_ctrl(1)? And if so, how?


  Thank you so much!


  --
  Héctor
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20160616/77a6344e/attachment.html>