[Mesa-users] Setting Omega/Omega_C and equatorial velocity

[AVG]

Dear all,

I am trying to replicate the results of the Table 1 of [1], specifically 
those of 2.9 solar masses.

In order to do that, I have simulated three stars whose rotation starts 
at ZAMS (with 2.9 solar masses and omega/omega_C = 0.2, 0.4 and 0.6), 
and I have simulated their profiles with GYRE at the age of 20 Gyears. 
As I found many discrepancies, I wonder if my problem could be solved by 
setting a value for equatorial velocity of rotation (so I could set this 
magnitude to the values of Table 1 of [1]).

How could I set an equatorial velocity to be used alongside the 
new_omega_div_omega_crit (taking into account that rotation only starts 
at ZAMS but not before)?


[1] Salmon et al. 2014, A&A, 
https://www.aanda.org/articles/aa/pdf/2014/09/aa23259-13.pdf

Best regards,

Antonio Varo